A method and system for comparing natural gas sources

Abstract

This invention belongs to the field of natural gas exploration technology, and discloses a method and system for natural gas source correlation. The method includes: obtaining natural gas samples from a target gas reservoir in a deep-to-ultra-deep high-evolutionary stratigraphic study area; determining the measured values ​​of Ar component abundance in air standard samples; and... 40 Ar / 36 Ar isotope measurements; determination of Ar component abundance measurements in natural gas samples and 40 Ar / 36 Ar isotope measurements; determination of Ar abundance correction values ​​in natural gas samples; determination of the radiogenic origin of natural gas samples. 40 Ar abundance value; determining the age and dominant source rocks of natural gas in the target gas reservoir. This invention enriches and improves the existing technical indicator system for natural gas genetic identification and source rock correlation, providing effective methods and technical means for the genetic identification and source rock correlation of complex natural gas in deep-to-ultra-deep layers, and providing important technical support for guiding natural gas exploration.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of natural gas exploration, and particularly relates to a natural gas source comparison method and system. BACKGROUND

[0002] With the current natural gas exploration and development expanding from the middle-shallow layer to the deep-ultra deep layer of the old high evolution layer series, the exploration objects are increasingly complex and difficult, and the difficulty and risk of exploration and discovery are increasing.

[0003] Currently, natural gas origin identification and gas source comparison mainly use multiple methods such as natural gas components and isotopes, light hydrocarbon components and isotopes, biomarkers in condensate oil and reservoir asphalt, and non-hydrocarbon gas components and isotopes. With the current natural gas exploration and development expanding from the middle-shallow layer to the deep-ultra deep layer of the old high evolution layer series, the exploration objects are increasingly complex and difficult, and the difficulty and risk of exploration and discovery are increasing, and the discovery of conventional single origin gas reservoirs is becoming less and less, and deep high evolution complex multi-source gas reservoirs are increasing. In particular, deep-ultra deep natural gas is mostly high-mature high-over mature natural gas, the components are mainly methane, the content of ethane and above heavy hydrocarbons is extremely low, the dry coefficient is high, and the carbon and hydrogen isotopes of alkane gas are reversed or inverted, and the useful information that can be used is less and less, the origin is very complex, and the conventional origin identification and gas source comparison method is difficult to identify, which restricts the determination of the origin of deep-ultra deep high evolution complex natural gas, the determination of the main hydrocarbon source rock layer series, and the selection of exploration direction.

[0004] In summary, the current natural gas origin identification and gas source comparison method provides less useful information in the deep-ultra deep old high evolution layer series, and cannot complete the determination of the origin of deep-ultra deep high evolution complex natural gas, the determination of the main hydrocarbon source rock layer series, and the exploration. SUMMARY

[0005] In view of the above problems, the application provides a natural gas source comparison method and system, which adopts the following technical scheme:

[0006] A natural gas source comparison method, comprising the following steps: obtaining a natural gas sample of a target gas reservoir in a deep-ultra deep high evolution layer series research area; determining a measurement value of Ar component abundance in an air standard sample and 40 Ar / 36 Ar isotope measurement value; determining a measurement value of Ar component abundance in the natural gas sample and 40 Ar / 36 Ar isotope measurement value; determining a corrected value of Ar component abundance in the natural gas sample by the air standard Ar component abundance correction coefficient and the measurement value of Ar component abundance in the natural gas sample; and determining a correlation between the Ar component abundance corrected value in the natural gas sample and the Ar component abundance in the air standard sample. 40 Ar / 36The Ar isotope actual value and the Ar component abundance correction value determine the radioactivity origin of the natural gas sample 40 The Ar abundance value; according to the radioactivity origin of the natural gas sample 40 The Ar abundance value and the known geological data of the target gas reservoir determine the age of the hydrocarbon source rock and the main hydrocarbon source rock of the natural gas in the target gas reservoir.

[0007] Further, the measured value of the Ar component abundance in the air standard sample and 40 Ar / 36 The Ar isotope measured value includes the following steps:

[0008] The air standard sample is purified to obtain rare gases in the air standard sample;

[0009] The Ar component abundance in the rare gases in the air standard sample is measured multiple times, and the average value is taken to obtain the measured value of the Ar component abundance in the air standard sample;

[0010] The Ar component abundance in the rare gases in the air standard sample is measured multiple times, and the average value is taken to obtain the measured value of the Ar component abundance in the air standard sample; 40 Ar / 36 The Ar isotope value, and the average value is taken to obtain the 40 Ar / 36 Ar isotope measured value.

[0011] Further, it further includes the following steps:

[0012] According to the average value of the multiple Ar component abundance measured values in the air standard sample and the standard gas Ar component abundance correction value, the air standard sample Ar component abundance correction coefficient is determined.

[0013] Further, it further includes the following steps:

[0014] According to the average value of the multiple 40 Ar / 36 Ar isotope measured value in the air standard sample and the standard gas 40 Ar / 36 Ar isotope actual value, the air standard sample 40 Ar / 36 Ar isotope correction coefficient.

[0015] Further, the measured value of the Ar component abundance in the natural gas sample and 40 Ar / 36 The Ar isotope measured value includes the following steps:

[0016] The natural gas sample is purified to obtain rare gases in the natural gas sample;

[0017] The abundance of Ar component in rare gases in natural gas samples was measured multiple times, and the average value was taken to obtain the measured value of Ar component abundance in natural gas samples.

[0018] Multiple measurements of rare gases in natural gas samples 40 Ar / 36 Ar isotope values ​​were taken, and the average value was used to obtain the Ar isotope values ​​in the natural gas sample. 40 Ar / 36 Ar isotope measurements.

[0019] Furthermore, based on the natural gas samples 40 Ar / 36 Actual Ar isotope values ​​and Ar component abundance correction values ​​were used to determine the radiogenic origin of the natural gas sample. 40 Ar abundance values ​​are as follows:

[0020] C 40 Ar*=C 40 Ar 校 ×( 40 Ar / 36 Ar 实 -295.5) / ( 40 Ar / 36 Ar 实 +1.188)

[0021] In the formula, C 40 Ar* indicates the radiogenic origin of the natural gas sample. 40 Ar abundance value, C 40 Ar 校 This indicates the correction value for the abundance of the Ar component in the natural gas sample. 40 Ar / 36 Ar 实 Indicating natural gas sample 40 Ar / 36 Actual values ​​of Ar isotopes.

[0022] Furthermore, the known geological data of the target gas reservoir include the average porosity of the target gas reservoir, the average potassium content of the possible source rocks of the target gas reservoir, the average potassium content of the known source rocks in the known gas reservoirs in the study area, and the actual stratigraphic age of the known source rocks of the known gas reservoirs.

[0023] Furthermore, based on the radiogenic nature of the natural gas samples... 40 Based on Ar abundance values ​​and known geological data of the target gas reservoir, determine the age of the source rocks and the dominant source rocks of the natural gas in the target gas reservoir, including the following steps:

[0024] According to the obtained porosity average value of the known gas reservoir in the same research area as the target gas reservoir, the average potassium content of the known source rock system in the known gas reservoir, and the radioactivity genesis of multiple natural gas samples of the known gas reservoir 40 The average value of Ar abundance, and the balance compensation coefficient of the research area are determined.

[0025] According to the radioactivity genesis of natural gas in the target gas reservoir 40 The average value of Ar abundance, the average potassium content of the known source rock system in the target gas reservoir, and the balance compensation coefficient of the research area are used to estimate the age of the source rock of natural gas in the target gas reservoir.

[0026] The estimated age of the source rock of natural gas in the target gas reservoir is compared with the actual stratigraphic age of the source rock existing in the research area, and the main source rock of natural gas in the target gas reservoir is determined.

[0027] Further, according to the obtained porosity average value of the known gas reservoir in the same research area as the target gas reservoir, the average potassium content of the known source rock system in the known gas reservoir, and the radioactivity genesis of multiple natural gas samples of the known gas reservoir 40 The average value of Ar abundance, and the balance compensation coefficient of the research area are determined.

[0028]

[0029] In the formula, C 40 Ar* 平均 represents the radioactivity genesis of natural gas in the known gas reservoir 40 The average value of Ar abundance, [K] represents the average potassium content of the known source rock system in the known gas reservoir, represents the porosity average value of the known gas reservoir, M represents the balance compensation coefficient of the research area, and X K represents the abundance of K in nature 40 K, and λ K represents 40 K decays into 40 Ar, and λ represents 40 K decays into 40 Ar and 40 Ca, and t represents the age of the known source rock of the known gas reservoir in the research area.

[0030] The application also provides a natural gas source comparison system, comprising:

[0031] A sampling container is used to obtain a natural gas sample of a target gas reservoir in a deep-ultra deep high evolution layer research area.

[0032] A quadrupole mass spectrometer is used to determine the measurement value of the Ar component abundance in the air standard sample and the natural gas sample.

[0033] Rare gas isotope mass spectrometer for determining the isotopic composition of air standards and natural gas samples 40 Ar 36 Ar isotope measurements;

[0034] A first calculation module is configured to determine the Ar component abundance correction value in the natural gas sample by correcting the Ar component abundance measured value in the natural gas sample by the air standard Ar component abundance correction coefficient;

[0035] A second calculation module is configured to determine the radiogenic 40 Ar 36 Ar isotope actual value and the Ar component abundance correction value, to determine the radiogenic 40 Ar abundance value of the natural gas sample;

[0036] A third calculation module is configured to determine the age of the source rock of the natural gas in the target gas reservoir and the main source rock of the natural gas in the target gas reservoir according to the radiogenic 40 Ar abundance value of the natural gas sample and the known geological data of the target gas reservoir.

[0037] Further, the known geological data of the target gas reservoir includes the average porosity of the target gas reservoir, the average potassium content of the possible source rock system of the target gas reservoir, the average potassium content of the known source rock system in the known gas reservoir in the same research area as the target gas reservoir, and the actual stratigraphic age of the known source rock of the known gas reservoir.

[0038] Further, the third calculation module is specifically configured to:

[0039] determine the balance compensation coefficient of the research area according to the average porosity of the known gas reservoir in the same research area as the target gas reservoir, the average potassium content of the known source rock system in the known gas reservoir, and the average value of the radiogenic 40 Ar abundance of the plurality of natural gas samples of the known gas reservoir;

[0040] estimate the age of the source rock of the natural gas in the target gas reservoir according to the radiogenic 40 Ar abundance value of the natural gas in the target gas reservoir, the average potassium content of the known source rock system in the target gas reservoir, and the balance compensation coefficient of the research area;

[0041] compare the estimated age of the source rock of the natural gas in the target gas reservoir with the actual stratigraphic age of the source rock existing in the research area to determine the main source rock of the natural gas in the target gas reservoir.

[0042] The present application enriches and perfects the existing natural gas origin identification and gas source comparison technical index system, can provide effective method and technical means for deep-ultra deep high evolution complex natural gas origin identification and gas source comparison, provide important technical support for determining deep-ultra deep high evolution complex natural gas origin and source, determining main hydrocarbon source rock layer and its contribution, and guiding deep-ultra deep natural gas exploration.

[0043] Other features and advantages of the present application will be set forth in the following description, and in part will become apparent to those skilled in the art upon examination of the following or can be learned by practice of the present application. The objects and other advantages of the present application can be realized and attained by the structure particularly pointed out in the description and appended claims. BRIEF DESCRIPTION OF DRAWINGS

[0044] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description are some embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor on the basis of these drawings.

[0045] Figure 1 A flowchart of a natural gas source comparison method according to an embodiment of the present application is shown;

[0046] Figure 2 A diagram of Ar abundance values of several target gas reservoirs in a certain basin determined by the method according to an embodiment of the present application is shown. 40 A diagram of Ar abundance values of several target gas reservoirs in a certain basin determined by the method according to an embodiment of the present application is shown. DETAILED DESCRIPTION

[0047] In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the following will combine the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are some embodiments of the present application, but not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the protection scope of the present application.

[0048] It should be noted that the terms "first", "second" and the like in the present application are used to distinguish similar objects, and do not necessarily have to describe a specific order or sequence.

[0049] The present application first introduces the radioactive origin argon component:

[0050] Argon has three natural isotopes: 36 Ar, 38 Ar and 40Ar, the isotopic abundance of each is 0.337%, 0.063%, 99.600% respectively, wherein, 40 Ar is absolutely dominant, 40 Ar is formed by the decay of 40 K, 40 K has two branch decay modes, and Ar is produced by the decay of B, 40 Ca, and Ar is produced by electron capture, 40 Ar. Natural gas radioactivity causes 40 Ar is mainly derived from source rock 40 Decay of K.

[0051] Based on the above theory, the abundance of the radioactivity caused argon component can provide necessary information for determining the cause and source of deep-ultra deep high evolution complex natural gas, and therefore, the embodiment of the application provides a natural gas source comparison method and system, based on the abundance of the radioactivity caused argon component in natural gas, to provide an effective technical means for identifying the cause of deep-ultra deep high evolution complex natural gas and comparing gas sources.

[0052] As shown in Figure 1 , a natural gas source comparison method comprises the following steps:

[0053] S1, obtaining a natural gas sample of a target gas reservoir of a deep-ultra deep high evolution layer series research area through a sampling container.

[0054] For example, the sampling container can be a stainless steel high-pressure cylinder with double valves.

[0055] Step S1 is specifically as follows: a stainless steel high-pressure cylinder with double valves is selected to collect a natural gas sample of a deep-ultra deep high evolution gas reservoir, a mechanical pump is used to vacuum the stainless steel high-pressure cylinder to a set value (for example, the vacuum degree is 10 -1 Pa) before sampling, and the stainless steel high-pressure cylinder is repeatedly flushed with natural gas of the target gas reservoir for 4-6 times for 5-10 minutes during the sampling process, and the sampling is stopped when the gas pressure in the stainless steel high-pressure cylinder is between 3-6 MPa, and the natural gas sample of the target gas reservoir is obtained.

[0056] S2, determining the measurement value of the Ar component abundance and 40 Ar / 36 Ar isotope measurement value of the air standard sample through a quadrupole mass spectrometer and a noble gas isotope mass spectrometer, comprising the following steps:

[0057] S21, connecting the steel cylinder containing the air standard sample with a natural gas noble gas sampling device, purifying the active gas such as N2 and O2 outside the noble gas in the air standard sample by using a natural gas noble gas pretreatment and sampling device, and obtaining the noble gas in the air standard sample.

[0058] S22, send the rare gas in the air standard sample into the four-stage mass spectrometer, measure the Ar component abundance by using the four-stage mass spectrometer, measure 3 times and take the average of the measurement values, and obtain the measurement value of the Ar component abundance in the air standard sample.

[0059] S23, send the rare gas in the air standard sample into the rare gas isotope mass spectrometer, measure the Ar isotope value by using the rare gas isotope mass spectrometer, measure 3 times and take the average of the measurement values, and obtain the Ar isotope measurement value in the air standard sample. 40 Ar / 36 40 Ar / 36

[0060] S3, determine the Ar component abundance correction coefficient of the air standard sample according to the average of the measurement values of the plurality of Ar component abundances in the air standard sample and the Ar component abundance correction value (fixed value 0.934%) of the standard gas, and the specific steps are as follows:

[0061] Repeat steps S21 and S22 6 times, calculate the average of the Ar component abundance measurement values in the air standard sample for 6 times, and determine the Ar component abundance correction coefficient of the air standard sample according to the Ar component abundance correction value of the standard gas and the average of the Ar component abundance measurement values in the air standard sample for 6 times.

[0062] S4, determine the Ar isotope correction coefficient of the air standard sample according to the average of the plurality of Ar isotope measurement values in the air standard sample and the Ar isotope actual value (fixed value 295.5) of the standard gas, and the specific steps are as follows: 40 Ar / 36 40 Ar / 36 40 Ar / 36

[0063] Repeat steps S21 and S23 6 times, calculate the average of the Ar isotope measurement values in the air standard sample for 6 times, and determine the Ar isotope correction coefficient of the air standard sample according to the Ar isotope actual value of the standard gas and the average of the Ar isotope measurement values in the air standard sample for 6 times. 40 Ar / 36 40 Ar / 36 40 Ar / 36 40 Ar / 36

[0064] S5, determine the Ar component abundance measurement value and the Ar isotope measurement value in the natural gas sample by using the four-stage mass spectrometer and the rare gas isotope mass spectrometer, including the following steps: 40 Ar / 36 ​​​​​​​​​​

[0065] S51, connect the sampling container (stainless steel high-pressure cylinder) containing the natural gas sample to the natural gas rare gas sampling device through the pressure reducing valve, and use the natural gas rare gas sampling device to purify the active gases such as CH4 and CO2, N2 in the natural gas sample, so as to obtain the rare gas in the natural gas sample.

[0066] S52, send the rare gas in the natural gas sample into the quadrupole mass spectrometer, and use the quadrupole mass spectrometer to measure the Ar component abundance, repeat the measurement for 3 times and take the average value, so as to obtain the Ar component abundance measurement value in the natural gas sample.

[0067] S53, send the rare gas in the natural gas sample into the rare gas isotope mass spectrometer, and use the rare gas isotope mass spectrometer to measure the Ar isotope value, repeat the measurement for 3 times and take the average value, so as to obtain the Ar isotope measurement value in the natural gas sample. 40 Ar / 36 40 Ar / 36 .

[0068] S6, determine the Ar component abundance correction value in the natural gas sample by the air standard Ar component abundance correction coefficient and the Ar component abundance measurement value in the natural gas sample, and mark it as C 40 Ar 校 .

[0069] S7, determine the Ar isotope actual value in the natural gas sample by the air standard Ar isotope correction coefficient and the Ar isotope measurement value in the natural gas sample, and mark it as 40 Ar / 36 40 Ar / 36 40 Ar / 36 40 Ar 36 实 .

[0070] S8, determine the radioactive origin Ar abundance value of the natural gas sample according to the Ar isotope actual value and the Ar component abundance correction value in the natural gas sample, and mark it as C 40 Ar / 36 40 Ar 40

[0071] C 40 Ar*=C 40 Ar 校 ×( 40 Ar / 36 Ar 实 -295.5) / (​​​​​​​40 Ar / 36 Ar 实 +1.188) (1)

[0072] S9. Based on the radiogenic origin of natural gas samples 40 Based on the Ar abundance value and known geological data of the target gas reservoir, determine the age of the source rocks of the natural gas in the target gas reservoir and the main source rocks.

[0073] Among them, the known geological data of the target gas reservoir includes the average porosity of the known gas reservoir, the average potassium content of the possible source rocks of the target gas reservoir, the average potassium content of the known source rocks in the known gas reservoir in the study area, and the actual stratigraphic age of the known source rocks of the known gas reservoir.

[0074] For example, when practical conditions permit, sampling wells for natural gas from the target gas reservoir may involve source rock samples, and potassium content analysis of the source rocks may be performed or potassium content data of the relevant source rock strata of the target gas reservoir may be collected.

[0075] Step S9 is as follows:

[0076] S91. Based on the average porosity of known gas reservoirs in the same study area as the target gas reservoir. The average potassium content [K] of known source rocks in the gas reservoir, and the radiogenicity of multiple natural gas samples from the known gas reservoir. 40 The average abundance of Ar is denoted as C. 40 Ar* 平均 The equilibrium compensation coefficient M for the study area is determined as follows:

[0077]

[0078] In the formula, X K In nature 40 K abundance was 1.167 × 10⁻⁶. -4 ;λ K express 40 K decays 40 The decay constant of Ar is 5.81 × 10⁻⁶. -11 ;λ represents 40 K decays 40 Ar and 40 The total decay constant of Ca is 5.543 × 10⁻⁶. -10 t represents the age of the source rocks of the known gas reservoirs in the study area.

[0079] S92. Based on the actual stratigraphic ages of the known source rocks in the known gas reservoir, the average potassium content [K] of the known source rock series in the known gas reservoir, and the equilibrium compensation coefficient M obtained for the study area, determine the radiogenic origin of the natural gas in the target gas reservoir. 40Ar abundance value, denoted as C 40 Ar*.

[0080] S93, according to the target gas reservoir natural gas in radioactivity cause 40 Ar abundance value C 40 Ar*, the average value of the potassium content of the known source rock system in the target gas reservoir [K], and the balance compensation coefficient M of the study area, the age t of the source rock of the natural gas in the target gas reservoir is estimated by formula (3).

[0081]

[0082] S94, the estimated age t of the source rock of the natural gas in the target gas reservoir is compared with the actual stratigraphic age of the source rock existing in the study area, and the main source rock of the natural gas in the target gas reservoir is determined.

[0083] Based on the above natural gas source comparison method, the application also provides a natural gas source comparison system, comprising: a sampling container, a quadrupole mass spectrometer, a noble gas isotope mass spectrometer, a first calculation module, a second calculation module and a third calculation module.

[0084] Among them, the sampling container is used to obtain the natural gas sample of the target gas reservoir of the deep-ultra deep high evolution layer system research area; the quadrupole mass spectrometer is used to determine the measurement value of the Ar component abundance in the air standard sample and the natural gas sample; the noble gas isotope mass spectrometer is used to determine the measurement value of the Ar isotope in the air standard sample and the natural gas sample; the first calculation module is used to determine the correction value of the Ar component abundance in the natural gas sample by the Ar component abundance correction coefficient of the air standard sample and the measurement value of the Ar component abundance in the natural gas sample; the second calculation module is used to determine the radioactivity cause 40 Ar / 36 of the natural gas sample according to the actual value of the Ar isotope and the correction value of the Ar component abundance in the natural gas sample; the third calculation module is used to determine the age of the source rock of the natural gas in the target gas reservoir and the main source rock according to the radioactivity cause 40 Ar / 36 of the natural gas sample and the known geological data of the target gas reservoir. 40 40

[0085] For example, by the natural gas source comparison method and system of the application, the natural gas genesis and source of several target gas reservoirs in a certain basin are determined, the main source rock layer is determined, which provides important technical support for guiding deep-ultra deep natural gas exploration, and finally the data as shown in Table 1 and Figure 2 Table 2 are obtained.

[0086] Table 1

[0087] ​​

[0088] The natural gas source comparison method and system of the present application enriches and perfects the existing natural gas origin identification and source comparison technical index system based on the abundance of radioactive origin Ar component in natural gas, and can provide effective method and technical means for deep-ultra deep high evolution complex natural gas origin identification and source comparison, and has technical advantages in aspects of determining deep-ultra deep high evolution complex natural gas origin and source, determining main hydrocarbon source rock layer and its contribution, etc.

[0089] Although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those skilled in the art that the technical solutions recorded in the foregoing embodiments can still be modified, or some technical features can be replaced by equivalents; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method of natural gas gas source correlation, the method comprising: The method comprises the following steps: Obtaining a natural gas sample of a target gas reservoir in a deep-ultra deep high evolution layer series research area; determining a measured value of the abundance of the Ar component in the air standard and 40 Ar / 36 Ar isotope measured value; determining an Ar component abundance measurement in a natural gas sample and 40 Ar / 36 Ar isotope measurement; Determine the Ar component abundance correction value in the natural gas sample by the Ar component abundance correction coefficient of the air standard sample and the Ar component abundance measurement value in the natural gas sample; According to the actual value of Ar isotope and the correction value of Ar component abundance in the natural gas sample 40 Ar 36 The radioactivity origin of the natural gas sample is determined according to the actual value of Ar isotope and the correction value of Ar component abundance 40 The abundance value of Ar is as follows: C 40 Ar*=C 40 Ar 校 ×( 40 Ar / 36 Ar 实 -295.5) / ( 40 Ar / 36 Ar 实 +1.188) where C 40 Ar* represents the radioactivity origin of the natural gas sample 40 Ar abundance value, C 40 Ar 校 represents the Ar component abundance correction value in the natural gas sample, 40 Ar / 36 Ar 实 represents the Ar component abundance correction value in the natural gas sample, 40 Ar / 36 Ar isotope actual value; Determination of the radioactivity genesis of natural gas samples 40 The Ar abundance value and the known geological data of the target gas reservoir are used to determine the age of the hydrocarbon source rock and the main hydrocarbon source rock of the natural gas in the target gas reservoir, including the following steps: according to the obtained average value of the porosity of the known gas reservoir in the same research area as the target gas reservoir, the average value of the potassium content of the known hydrocarbon source rock system in the known gas reservoir, and the radioactivity genesis of multiple natural gas samples of the known gas reservoir 40 The average value of the Ar abundance is used to determine the balance compensation coefficient of the research area, specifically as follows: C 40 Ar* 平均 =[K]×X K ×Ø×M×λ K / λ(e λt -1) where C 40 Ar 平均 indicates the known radioactivity content of natural gas in known gas reservoirs 40 Ar abundance values average values, [K] indicates the average potassium content of known source rock series in known gas reservoirs, indicates the average porosity of known gas reservoirs, M indicates the equilibrium compensation factor of the study area, X K indicates the average potassium content of known source rock series in known gas reservoirs 40 K abundance, λ K indicates the average potassium content of known source rock series in known gas reservoirs 40 K decays into 40 Ar decay constant, λ indicates 40 K decays into 40 Ar and 40 Ca total decay constant, t indicates the age of known source rocks in known gas reservoirs in the study area; Based on the radiogenic origin of natural gas in the target gas reservoir 40 Ar abundance value, average potassium content of known source rocks in the target gas reservoir, and balance compensation coefficient of the study area are used to estimate the age of the source rocks of natural gas in the target gas reservoir. Compare the estimated age of the hydrocarbon source rock of the natural gas in the target gas reservoir with the actual stratigraphic age of the hydrocarbon source rock existing in the research area to determine the main hydrocarbon source rock of the natural gas in the target gas reservoir.

2. The natural gas gas source comparison method of claim 1, wherein, determining a measured value of the abundance of the Ar component in the air standard and 40 Ar / 36 Ar isotope measurement, comprising the steps of: Purify the air standard sample to obtain rare gases in the air standard sample; Measure the Ar component abundance in the rare gases in the air standard sample multiple times and take the average value to obtain the Ar component abundance measurement value in the air standard sample; The isotopic values of Ar in the air standard sample were measured multiple times and averaged to obtain the isotopic measurement value of Ar in the air standard sample. 40 Ar / 36 The isotopic values of Ar in the air standard sample were measured multiple times and averaged to obtain the isotopic measurement value of Ar in the air standard sample. 40 Ar / 36 The isotopic values of Ar in the air standard sample were measured multiple times and averaged to obtain the isotopic measurement value of Ar in the air standard sample.

3. The method of claim 1, wherein, The method further comprises the following steps: Determine the Ar component abundance correction coefficient of the air standard sample according to the average value of the Ar component abundance measurement values in the air standard sample and the Ar component abundance correction value of the standard gas.

4. The method of claim 1, wherein, The method further comprises the following steps: The average of the Ar isotope measurements and the standard gas 40 Ar 36 The average of the Ar isotope measurements and the standard gas 40 Ar 36 The average of the Ar isotope measurements and the standard gas 40 Ar 36 The average of the Ar isotope measurements and the standard gas 5. The method of claim 1, wherein, determining an Ar component abundance measurement in a natural gas sample and 40 Ar / 36 Ar isotope measurement, comprising the steps of: Purify the natural gas sample to obtain rare gases in the natural gas sample; Measure the Ar component abundance in the rare gases in the natural gas sample multiple times and take the average value to obtain the Ar component abundance measurement value in the natural gas sample; Multiple measurements of noble gases in a natural gas sample 40 Ar 36 Ar isotope values, and average to obtain an Ar isotope measurement value in the natural gas sample 40 Ar 36 Ar isotope measurement value.

6. The method of claim 1-5, wherein, The known geological data of the target gas reservoir include the average porosity of the target gas reservoir, the average potassium content of the possible hydrocarbon source rock series of the target gas reservoir, the average potassium content of the known hydrocarbon source rock series in the known gas reservoir in the research area, and the actual stratigraphic age of the known hydrocarbon source rock of the known gas reservoir.

7. A natural gas gas source correlation system characterized by, The method comprises the following steps: A sampling container for obtaining a natural gas sample of a target gas reservoir in a deep-ultra deep high evolution layer series research area; A quadrupole mass spectrometer for determining the Ar component abundance measurement value in the air standard sample and the natural gas sample; A noble gas isotope mass spectrometer for determining the isotopic composition of air standards and natural gas samples 40 Ar 36 Ar isotope measurements; A first calculation module for determining the Ar component abundance correction value in the natural gas sample by the Ar component abundance correction coefficient of the air standard sample and the Ar component abundance measurement value in the natural gas sample; The second calculation module is configured to determine the radioactivity origin of the natural gas sample according to the actual value of Ar isotope and the Ar component abundance correction value in the natural gas sample. 40 Ar 36 The actual value of Ar isotope and the Ar component abundance correction value are used to determine the radioactivity origin of the natural gas sample. 40 The abundance value of Ar is determined as follows: C 40 Ar*=C 40 Ar 校 ×( 40 Ar / 36 Ar 实 -295.5) / ( 40 Ar / 36 Ar 实 +1.188) where C 40 Ar* represents the radioactivity origin of the natural gas sample 40 Ar abundance value, C 40 Ar 校 represents the Ar component abundance correction value in the natural gas sample, 40 Ar / 36 Ar 实 represents the Ar component abundance correction value in the natural gas sample, 40 Ar / 36 Ar isotope actual value; The third calculation module is configured to determine the age of the source rock and the main source rock of the natural gas in the target gas reservoir according to the radioactivity genesis of the natural gas sample, the Ar abundance value and the known geological data of the target gas reservoir. 40 The Ar abundance value, the known geological data of the target gas reservoir and the radioactivity genesis of the natural gas sample are used to determine the age of the source rock and the main source rock of the natural gas in the target gas reservoir, including the following steps: according to the obtained average porosity value of the known gas reservoir in the same research area as the target gas reservoir, the average potassium content value of the known source rock system in the known gas reservoir and the radioactivity genesis of the plurality of natural gas samples in the known gas reservoir 40 The average value of the Ar abundance is used to determine the balance compensation coefficient of the research area, specifically as follows: C 40 Ar* 平均 =[K]×X K ×Ø×M×λ K / λ(e λt -1) where C 40 Ar 平均 indicates the known radioactivity content of natural gas in known gas reservoirs 40 Ar abundance values average values, [K] indicates the average potassium content of known source rock series in known gas reservoirs, indicates the average porosity of known gas reservoirs, M indicates the equilibrium compensation factor of the study area, X K indicates the average potassium content of known source rock series in known gas reservoirs 40 K abundance, λ K indicates the average potassium content of known source rock series in known gas reservoirs 40 K decays into 40 Ar decay constant, λ indicates 40 K decays into 40 Ar and 40 Ca total decay constant, t indicates the age of known source rocks in known gas reservoirs in the study area; Based on the radiogenic origin of natural gas in the target gas reservoir 40 Ar abundance value, average potassium content of known source rocks in the target gas reservoir, and balance compensation coefficient of the study area are used to estimate the age of the source rocks of natural gas in the target gas reservoir. Compare the estimated age of the hydrocarbon source rock of the natural gas in the target gas reservoir with the actual stratigraphic age of the hydrocarbon source rock existing in the research area to determine the main hydrocarbon source rock of the natural gas in the target gas reservoir.

8. The natural gas source comparison system of claim 7, wherein, The known geological data of the target gas reservoir include the average porosity of the target gas reservoir, the average potassium content of the possible hydrocarbon source rock series of the target gas reservoir, the average potassium content of the known hydrocarbon source rock series in the known gas reservoir in the research area, and the actual stratigraphic age of the known hydrocarbon source rock of the known gas reservoir.

Images