CN117150207B - Method, equipment and storage medium for evaluating associated helium resource quantity - Google Patents

Method, equipment and storage medium for evaluating associated helium resource quantity Download PDF

Info

Publication number
CN117150207B
CN117150207B CN202310967793.3A CN202310967793A CN117150207B CN 117150207 B CN117150207 B CN 117150207B CN 202310967793 A CN202310967793 A CN 202310967793A CN 117150207 B CN117150207 B CN 117150207B
Authority
CN
China
Prior art keywords
helium
percentage
natural gas
block
evaluated
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.)
Active
Application number
CN202310967793.3A
Other languages
Chinese (zh)
Other versions
CN117150207A (en
Inventor
李振
张金川
莫宣学
丁江辉
仝忠正
李兴起
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN202310967793.3A priority Critical patent/CN117150207B/en
Publication of CN117150207A publication Critical patent/CN117150207A/en
Application granted granted Critical
Publication of CN117150207B publication Critical patent/CN117150207B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Geology (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Business, Economics & Management (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Agronomy & Crop Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Algebra (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The invention provides an evaluation method, equipment and storage medium of associated helium resource amount, wherein the method comprises the following steps: acquiring sampling data of a plurality of natural gas wells in a block to be evaluated; determining the gas-containing area of the block to be evaluated based on the reservoir thickness and the reservoir gas content of each natural gas well and the bottom plate contour map; determining a probability density function of natural gas resource abundance in the block to be evaluated based on reservoir rock density, reservoir gas content and reservoir thickness of the same reservoir section of each natural gas well; determining the average value and standard deviation of the total helium percentage content of the block to be evaluated based on the total helium percentage content of each sampling sample, and determining a probability density function of the total helium percentage content in the block to be evaluated; and determining the associated helium resource amount of the block to be evaluated based on the probability density function of the gas-containing area of the reservoir, the abundance of the natural gas resource and the probability density function of the total helium percentage in the block to be evaluated. The method can improve the accuracy of the evaluation of the associated helium resource quantity.

Description

Method, equipment and storage medium for evaluating associated helium resource quantity
Technical Field
The invention relates to the technical field of helium resource exploration, in particular to an evaluation method, equipment and storage medium of associated helium resource quantity.
Background
Helium is an extremely important natural resource, and is widely applied to the fields of medicine, aerospace, military, nuclear industry and the like, and particularly has special application in many high-precision scientific researches.
Helium is a non-renewable resource and is produced in very thin yields. Helium in the stratum is mainly endowed in a natural gas reservoir in a natural gas associated state, and the calculation of associated helium resources is also carried out by means of a natural gas resource evaluation method, but the formation, the source, the migration, the aggregation and other aspects have obvious differences, so that the existing natural gas resource evaluation method is not suitable for the associated helium.
How to accurately calculate the natural gas associated helium resource amount becomes a technical problem to be solved at present.
Disclosure of Invention
The embodiment of the invention provides an evaluation method, equipment and a storage medium of associated helium resource quantity, which are used for solving the problem of inaccurate calculation of the current associated helium resource quantity.
In a first aspect, an embodiment of the present invention provides a method for evaluating an associated helium resource amount, including:
acquiring sampling data of a plurality of natural gas wells in a block to be evaluated; the sampling data comprise reservoir rock density, reservoir gas content and reservoir thickness of each natural gas well, and total helium percentage content of each sampling sample;
Determining the gas-containing area of the block to be evaluated based on the reservoir thickness and the reservoir gas content of each natural gas well and the bottom plate contour map;
determining a probability density function of natural gas resource abundance in the block to be evaluated based on reservoir rock density, reservoir gas content and reservoir thickness of the same reservoir section of each natural gas well;
determining an average value and a standard deviation of the total helium percentage of the block to be evaluated based on the total helium percentage of each sampling sample, and analyzing the normal distribution of the total helium percentage of the block to be evaluated based on the average value and the standard deviation of the total helium percentage of the block to be evaluated to determine a probability density function of the total helium percentage in the block to be evaluated;
and determining the associated helium resource quantity of the block to be evaluated under a set probability based on the probability density function of the gas-containing area of the reservoir, the abundance of the natural gas resource and the probability density function of the total helium percentage in the block to be evaluated.
In one possible implementation, the sampled data further includes helium isotope ratios for each sampled sample;
the evaluation method further comprises the following steps:
determining the percentage of shell source helium and/or the percentage of veil source helium for each sample based on the total helium percentage and helium isotope ratio for each sample;
Determining an average value and a standard deviation of the percentage content of the shell source helium of the block to be evaluated based on the percentage content of the shell source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the shell source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the shell source helium of the block to be evaluated to determine a probability density function of the percentage content of the shell source helium in the block to be evaluated; and/or
Determining an average value and a standard deviation of the percentage content of the curtain source helium of the block to be evaluated based on the percentage content of the curtain source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the curtain source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the curtain source helium of the block to be evaluated to determine a probability density function of the percentage content of the curtain source helium in the block to be evaluated.
In one possible implementation, determining the percentage of shell source helium and/or the percentage of veil source helium for each sample based on the total helium percentage and helium isotope ratio for each sample comprises:
determining a shell source helium percentage content of each sampling sample based on a total helium percentage content, a helium isotope ratio and a shell source helium percentage content determination formula of each sampling sample;
The percentage of veil source helium for each sample is determined based on the total percentage helium for each sample, the helium isotope ratio, and the veil source helium percentage determination formula.
In one possible implementation, the shell source helium percentage determination formula χ (4 He ) The method comprises the following steps:
determination formula χ of percentage of valance source helium (3 He ) The method comprises the following steps:
the probability density function f (χ (He)) of the total helium percentage in the block to be evaluated is:
probability density function f (χ (4) He ) Is) is:
probability density function f (χ (3) He ) Is) is:
wherein,and->The average value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage is respectively; sigma (sigma) χ(He) 、σ χ(4He) 、σ χ(3He) The standard deviation of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; χ (He) i The value range of (C) is [ χ (He) min ，χ(He) max ]，χ(4 He ) i The value range of (X) is% 4 He) min ，χ( 4 He) max ]，χ(3 He ) i The value range of (X) is% 3 He) min ，χ( 3 He) max ]，χ(He) min 、χ(He) max 、χ( 4 He) min 、χ( 4 He) max 、χ( 3 He) min 、χ( 3 He) max The minimum value and the maximum value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; chi (He) is the percentage of total helium He; />Is helium isotope 3 He/ 4 He ratio. ,
in one possible implementation, determining the associated helium resource amount for the block under evaluation at a set probability based on a probability density function of reservoir gas area, natural gas resource abundance, and a probability density function of total helium percentage within the block under evaluation, comprising:
Integrating the probability density function of the natural gas resource abundance and the probability density function of the total helium percentage content respectively, and determining the natural gas resource abundance value and the total helium percentage content value under the set probability;
based on the natural gas resource abundance value and the total helium percentage content value under the set probability, determining the associated helium resource quantity of the block to be evaluated under the set probability according to the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability;
preset resource quantity calculating formula Q (X) P The method comprises the following steps:
Q(X) P ＝χ(X) P ×R P ×S；
wherein X is gas, χ (X) P R is the percentage of X gas under the probability P P The probability is the abundance of natural gas resources under P, S is the gas-containing area of the reservoir, P is the probability value, and the value range is [0%,100%]。
In one possible implementation, the evaluation method further includes:
integrating the probability density function of the shell source helium percentage and/or the probability density function of the veil source helium percentage to determine the shell source helium percentage value and/or the veil source helium percentage value under the set probability;
determining the shell source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and the shell source helium percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated, and a preset resource quantity calculation formula under the set probability; and/or
And determining the valance source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and valance source helium percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability.
In one possible implementation, determining the gas bearing area of the zone to be evaluated based on the reservoir thickness and reservoir gas content for each gas well, and the floor contour map, comprises:
determining a bottom plate contour map of a block to be evaluated according to the elevation of a drilling hole of each natural gas well and the depth data of each natural gas well;
drawing a reservoir thickness contour map of a block to be evaluated based on the reservoir thickness of each natural gas well;
drawing a gas content contour map of a block to be evaluated based on the gas content of a reservoir of each natural gas well;
and superposing the reservoir thickness contour map and the air content contour map on the bottom plate contour map, and determining the air content area of the block to be evaluated.
In one possible implementation, determining a probability density function of natural gas resource abundance within a zone to be evaluated based on reservoir rock density, reservoir gas content, and reservoir thickness for the same reservoir interval for each natural gas well, comprises:
Determining natural gas resource abundance of each natural gas well based on reservoir rock density, reservoir gas content and reservoir thickness of the same reservoir section of each natural gas well;
determining natural gas resource abundance average values and natural gas resource abundance standard deviations of the to-be-evaluated blocks based on the natural gas resource abundance of all natural gas wells;
analyzing the natural gas resource abundance sampling normal distribution of the block to be evaluated based on the natural gas resource abundance average value and the natural gas resource abundance standard deviation of the block to be evaluated so as to determine a probability density function of the natural gas resource abundance in the block to be evaluated;
wherein, natural gas resource abundance R of ith natural gas well i The method comprises the following steps:
the probability density function f (R) of the abundance of natural gas resources is:
wherein m is the total number of sections of the ith natural gas well, ρ ij For the rock density of the j th section of the reservoir in the i th natural gas well, H ij C is the thickness of the j section of the reservoir in the ith natural gas well ij The gas content of the jth section of the reservoir in the ith natural gas well; sigma (sigma) R Is the standard deviation of the abundance of natural gas resources,the value range of R is [ R ] which is the average value of the abundance of natural gas resources min ，R max ]，R min Is the minimum value of natural gas resource abundance, R max Is the maximum value of natural gas resource abundance.
In a second aspect, an embodiment of the present invention provides an apparatus for evaluating an associated helium gas resource amount, including:
the acquisition module is used for acquiring sampling data of a plurality of natural gas wells in the block to be evaluated; the sampling data comprise reservoir rock density, reservoir gas content and reservoir thickness of each natural gas well, and total helium percentage content of each sampling sample;
the first determining module is used for determining the gas containing area of the block to be evaluated based on the reservoir thickness and the reservoir gas content of each natural gas well and the bottom plate contour map;
the second determining module is used for determining a probability density function of the abundance of the natural gas resources in the block to be evaluated based on the reservoir rock density, the reservoir gas content and the reservoir thickness of the same reservoir section of each natural gas well;
the third determining module is used for determining the average value and the standard deviation of the total helium percentage of the block to be evaluated based on the total helium percentage of each sampling sample, and analyzing the normal distribution of the total helium percentage of the block to be evaluated based on the average value and the standard deviation of the total helium percentage of the block to be evaluated so as to determine the probability density function of the total helium percentage in the block to be evaluated;
And the fourth determining module is used for determining the accompanying helium resource quantity of the block to be evaluated under a set probability based on the probability density function of the gas-containing area of the reservoir layer, the abundance of the natural gas resource and the probability density function of the total helium percentage in the block to be evaluated.
In one possible implementation, the sampled data further includes helium isotope ratios for each sampled sample;
the evaluation device further includes: a third determination module for determining the percentage of shell source helium and/or the percentage of veil source helium for each sample based on the total helium percentage and helium isotope ratio for each sample;
determining an average value and a standard deviation of the percentage content of the shell source helium of the block to be evaluated based on the percentage content of the shell source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the shell source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the shell source helium of the block to be evaluated to determine a probability density function of the percentage content of the shell source helium in the block to be evaluated; and/or
Determining an average value and a standard deviation of the percentage content of the curtain source helium of the block to be evaluated based on the percentage content of the curtain source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the curtain source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the curtain source helium of the block to be evaluated to determine a probability density function of the percentage content of the curtain source helium in the block to be evaluated.
In one possible implementation, the third determining module is configured to determine a shell source helium percentage for each sample based on a total helium percentage, a helium isotope ratio, and a shell source helium percentage determination formula for each sample;
the percentage of veil source helium for each sample is determined based on the total percentage helium for each sample, the helium isotope ratio, and the veil source helium percentage determination formula.
In one possible implementation, the shell source helium percentage determination formula χ (4 He ) The method comprises the following steps:
determination formula χ of percentage of valance source helium (3 He ) The method comprises the following steps:
the probability density function f (χ (He)) of the total helium percentage in the block to be evaluated is:
probability density function f (χ (4) He ) Is) is:
probability density function f (χ (3) He ) Is) is:
wherein,and->The average value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage is respectively; sigma (sigma) χ(He) 、σ χ(4He) 、σ χ(3He) The standard deviation of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; χ (He) i The value range of (C) is [ χ (He) min ，χ(He) max ]，χ(4 He ) i The value range of (X) is% 4 He) min ，χ( 4 He) max ]，χ(3 He ) i The value range of (X) is% 3 He) min ，χ( 3 He) max ]，χ(He) min 、χ(He) max 、χ( 4 He) min 、χ( 4 He) max 、χ( 3 He) min 、χ( 3 He) max The minimum value and the maximum value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; chi (He) is the percentage of total helium He; />Is helium isotope 3 He/ 4 He ratio.
In one possible implementation manner, the fourth determining module is configured to integrate the probability density function of the abundance of the natural gas resource and the probability density function of the total helium percentage, respectively, to determine the abundance value of the natural gas resource and the total helium percentage value under the set probability;
based on the natural gas resource abundance value and the total helium percentage content value under the set probability, determining the associated helium resource quantity of the block to be evaluated under the set probability according to the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability;
calculation formula of preset resource quantityQ(X) P The method comprises the following steps:
Q(X) P ＝χ(X) P ×R P ×S；
wherein X is gas, χ (X) P R is the percentage of X gas under the probability P P The probability is the abundance of natural gas resources under P, S is the gas-containing area of the reservoir, P is the probability value, and the value range is [0%,100%]。
In one possible implementation, the evaluation device further includes:
the fourth determining module is used for integrating the probability density function of the percentage content of the shell source helium and/or the probability density function of the percentage content of the mantle source helium, and determining the percentage content value of the shell source helium and/or the percentage content value of the mantle source helium under the set probability;
Determining the shell source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and the shell source helium percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated, and a preset resource quantity calculation formula under the set probability; and/or
And determining the valance source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and valance source helium percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability.
In one possible implementation manner, a first determining module is used for determining a bottom plate contour map of a block to be evaluated according to the elevation of a drilling hole of each natural gas well and the depth data of each natural gas well;
drawing a reservoir thickness contour map of a block to be evaluated based on the reservoir thickness of each natural gas well;
drawing a gas content contour map of a block to be evaluated based on the gas content of a reservoir of each natural gas well;
and superposing the reservoir thickness contour map and the air content contour map on the bottom plate contour map, and determining the air content area of the block to be evaluated.
In one possible implementation, a second determination module is configured to determine natural gas resource abundance of each natural gas well based on reservoir rock density, reservoir gas content, and reservoir thickness of the same reservoir section of each natural gas well;
Determining natural gas resource abundance average values and natural gas resource abundance standard deviations of the to-be-evaluated blocks based on the natural gas resource abundance of all natural gas wells;
analyzing the natural gas resource abundance sampling normal distribution of the block to be evaluated based on the natural gas resource abundance average value and the natural gas resource abundance standard deviation of the block to be evaluated so as to determine a probability density function of the natural gas resource abundance in the block to be evaluated;
wherein, natural gas resource abundance R of ith natural gas well i The method comprises the following steps:
the probability density function f (R) of the abundance of natural gas resources is:
wherein m is the total number of sections of the ith natural gas well, ρ ij For the rock density of the j th section of the reservoir in the i th natural gas well, H ij C is the thickness of the j section of the reservoir in the ith natural gas well ij The gas content of the jth section of the reservoir in the ith natural gas well; sigma (sigma) R Is the standard deviation of the abundance of natural gas resources,the value range of R is [ R ] which is the average value of the abundance of natural gas resources min ，R max ]，R min Is the minimum value of natural gas resource abundance, R max Is the maximum value of natural gas resource abundance.
In a third aspect, an embodiment of the present invention provides an electronic device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method according to the first aspect or any one of the possible implementations of the first aspect, when the computer program is executed by the processor.
In a fourth aspect, embodiments of the present invention provide a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method as described above in the first aspect or any one of the possible implementations of the first aspect.
The embodiment of the invention provides an evaluation method, equipment and a storage medium for associated helium resource quantity. Next, the gas bearing area of the block to be evaluated is determined based on the reservoir thickness and reservoir gas content of each gas well, and the floor contour map. Next, a probability density function of the abundance of natural gas resources within the zone to be evaluated is determined based on the reservoir rock density, reservoir gas content, and reservoir thickness of the same reservoir section for each natural gas well. Then, based on the total helium percentage content of each sampling sample, determining an average value and a standard deviation of the total helium percentage content of the block to be evaluated, and based on the average value and the standard deviation of the total helium percentage content of the block to be evaluated, analyzing the normal distribution of the total helium percentage content of the block to be evaluated to determine a probability density function of the total helium percentage content in the block to be evaluated. And finally, determining the associated helium resource quantity of the block to be evaluated under a set probability based on the probability density function of the gas-containing area of the reservoir layer, the abundance of the natural gas resource and the probability density function of the total helium percentage in the block to be evaluated.
According to the invention, on the basis of acquiring the sampling data of the plurality of natural gas wells in the block to be evaluated, the natural gas resource abundance, the gas-containing area of the reservoir and the percentage content of helium are subjected to statistical analysis, so that the helium resource quantity calculation result under the set probability can be calculated, and the accuracy of the helium resource quantity evaluation result is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments or the description of the prior art will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of an implementation of a method for evaluating the amount of associated helium resources provided by an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of an apparatus for evaluating the amount of associated helium resources according to an embodiment of the present invention;
fig. 3 is a schematic diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth such as the particular system architecture, techniques, etc., in order to provide a thorough understanding of the embodiments of the present invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the following description will be made by way of specific embodiments with reference to the accompanying drawings.
As described in the background, helium is increasingly widely used and in greater demand with the continued development of socioeconomic and scientific technologies. Helium is a non-renewable resource and the yield is very scarce, and we need to find new helium resources in order to meet the development needs.
Helium in the stratum is mainly endowed in a natural gas reservoir in a natural gas associated state, and the natural gas associated helium is found out to be a premise for improving the efficient utilization of helium resources. Therefore, how to perform a scientific and reasonable natural gas associated helium resource amount calculation method has become a research direction for helium exploration knowledge and development planning.
The calculation of helium resources is mostly carried out by means of a natural gas resource evaluation method, and currently, common methods for evaluating the amount of natural gas resources include: the volumetric method is based on the determination of the average value or probability of parameters such as gas-containing area, porosity, gas-containing saturation, volume coefficient and the like, so as to calculate the natural gas resource quantity. The causative method comprises the following steps: through the research of hydrocarbon generation and drainage and related geological processes, the resource quantity of the natural gas is determined mainly by determining the retention coefficient, but the difficulty of obtaining the retention coefficient is relatively high. The analogy method takes a standard region with similar geological conditions as an analogy object, and obtains analogy coefficients by selecting area, volume, gas content, deposition rate, EUR and other parameters which can be subjected to analogy, and the natural gas resource quantity of a region to be evaluated. And the evaluation unit dividing method divides the evaluation units according to geological conditions, exploration degrees and the like, and performs summation calculation on the resource quantity. The yield curve is subtracted, and the yield decreasing rule of the production well is researched, so that the natural gas resource quantity in the stratum can be calculated in a recursive manner. The material balance method is used for calculating the natural gas yield and the resource amount based on the material balance principle. And accumulating natural gas output even in the history by using a history statistical method to obtain related resource quantity. The terfei method is a comprehensive weighting method, and the resource quantity is obtained by weighting and summing the evaluation results of different evaluation methods.
Overall, the existing natural gas resource evaluation method mainly has the following problems and disadvantages when applied to helium gas resource amount calculation: the first and existing methods are not aimed at the calculation of the resource amount of the natural gas associated helium, so that the applicability and scientific effectiveness of the natural gas associated helium in the calculation of the resource amount of the associated helium cannot be ensured. Secondly, the existing method lacks constraint on helium geological awareness and data grasping degree, and cannot reflect the influence of the helium data grasping degree on a resource amount calculation result. Third, the calculation of the natural gas associated helium resource amount by the prior method is only considered from the voluntary aspect of helium, so that the relationship between the natural gas and the associated helium resource amount is ambiguous. Fourth, the existing method is difficult to consider the non-uniformity of helium distribution in the calculation of the natural gas associated helium resource amount, so that the result of the helium resource amount is not reasonable and scientific enough. Fifth, the existing method is mainly to calculate the total helium resource amount associated with natural gas, and is not capable of distinguishing the resource amounts of helium gas due to the shell source and the veil source. Sixthly, the existing method has few consideration on the aspects of data rationality, reliability, sources and the like, and further reduces the credibility of the resource calculation result to a certain extent.
It can be seen that the existing resource amount calculation method is not suitable for resource amount evaluation of natural gas associated helium.
In order to solve the problems in the prior art, the embodiment of the invention provides an evaluation method, equipment and a storage medium of associated helium resource quantity. The method provided by the embodiment of the invention is first described below.
Referring to fig. 1, a flowchart of an implementation of an evaluation method for an associated helium resource amount according to an embodiment of the present invention is shown, and details are as follows:
and step S110, acquiring sampling data of a plurality of natural gas wells in the block to be evaluated.
The sampled data includes reservoir rock density, reservoir gas content, and reservoir thickness for each gas well, as well as the total helium percentage for each sampled sample.
The reasonable calculation of the associated helium resource quantity depends on reliable data materials, so that before the associated helium resource quantity is evaluated, materials in aspects of geological exploration, geophysics, drilling, logging, experimental testing, development production and the like in a block to be evaluated are required to be collected, and the obtained materials are further subjected to rationality, representativeness and representativeness analysis, and related data materials with low reliability, unreasonable and invalid are removed. The data is analyzed and processed based on the subsequent need to use statistical knowledge, so that at least 30 effective natural gas well data are required to be counted.
After collecting the data of the natural gas well, the data of each well needs to be tidied. The data set of the 30 natural gas wells is w= (W) 1 ，W 2 ，W 3 ，……，W n ) Where n is the number of statistically valid wells. Ith natural gas well dataset W i ＝(ρ i ，C i ，H i ) Reservoir thickness H i ＝H i1 +H i2 +H i3 +……+H ip P is the reservoir rock density, reservoir gas content and reservoir in the ith natural gas wellNumber of layer segments of layer thickness. Ith natural gas well reservoir rock density dataset ρ i ＝(ρ i1 ，ρ i2 ，ρ i3 ，……，ρ ip ) (units: t/m 3 ) The method comprises the steps of carrying out a first treatment on the surface of the Gas content data set C of ith natural gas well i ＝(C i1 ，C i2 ，C i3 ，……，C ip ) (units: m is m 3 T); gas well reservoir thickness dataset H i ＝(H i1 ，H i2 ，H i3 ，……，H ip ) (units: m).
Because helium gas has a lower research degree than natural gas, data are less, and statistical carding is difficult to carry out by taking a well as a unit, the data of the total helium percentage content detected by each sampling sample in the area to be evaluated is used as the data of helium gas. Data set χ (He) = (χ (He) for total helium percentage 1 ，χ(He) 2 ，χ(He) 3 ，……，χ(He) q ) (units: % of the total number of samples taken) q.
And step S120, determining the gas-containing area of the block to be evaluated based on the reservoir thickness and the reservoir gas content of each natural gas well and the bottom plate contour map.
In some embodiments, to determine the gas bearing area of the region to be evaluated, the gas bearing area of the region to be evaluated may be determined in combination with the geological related data, reservoir thickness, and reservoir gas content of the region to be evaluated. The specific steps are as follows:
And step 1201, determining a bottom plate contour map of the block to be evaluated according to the elevation of the drilling hole of each natural gas well and the depth data of each natural gas well.
The bottom plate contour map can be comprehensively compiled according to the data of field observation, drilling control, seismic interpretation and the like. The floor contour map is a projection map on a plan view based on the intersection line of the floor surface of the block to be evaluated and the level surfaces of the respective elevations. According to the elevation of the bottom plate surface of the block to be evaluated, a horizontal projection method is adopted, and a projection drawing piece, also called a construction drawing of the block to be evaluated, is represented by using contour lines.
And step 1202, drawing a reservoir thickness contour map of the block to be evaluated based on the reservoir thickness of each natural gas well.
Based on the obtained reservoir thickness of each natural gas well, and combining the fault distribution of the block to be evaluated, the pinch-out of the reservoir, the erosion of magma, the formation denudation, the burial depth, the wind oxidation zone, the lower limit of the reservoir thickness and the like, and drawing a reservoir thickness contour map of the block to be evaluated under the condition of considering the deposition environment.
And step S1203, drawing a gas content contour map of the block to be evaluated based on the gas content of the reservoir of each natural gas well.
Based on the reservoir gas content of each natural gas well obtained in the above, on the basis of analysis of the reservoir conditions and mechanisms, a gas content contour map of a block to be evaluated is drawn according to main control factors affecting natural gas enrichment.
And step S1204, superposing the reservoir thickness contour map and the air content contour map on the bottom plate contour map, and determining the air content area of the block to be evaluated.
Simultaneously drawing an air content contour map and a reservoir thickness contour map on a bottom plate contour map, and further obtaining an air content area S (unit: m) of a block to be evaluated by a method of overlapping the two contour maps 2 ) Is a value of (2).
Through the steps, the gas-containing area of the block to be evaluated can be accurately determined according to the geological related data, the reservoir thickness and the reservoir gas content.
And step S130, determining a probability density function of the abundance of the natural gas resources in the block to be evaluated based on the reservoir rock density, the reservoir gas content and the reservoir thickness of the same reservoir section of each natural gas well.
Considering that the research degree of natural gas is higher than that of helium at present, the invention performs statistical analysis on the related parameters of the abundance of the natural gas resources by taking a well as a unit, and performs statistical analysis on the abundance of the natural gas resources of the whole block to be evaluated on the basis. The method comprises the following specific steps:
And step S1301, determining the natural gas resource abundance of each natural gas well based on the reservoir rock density, the reservoir gas content and the reservoir thickness of the same reservoir section of each natural gas well.
Specifically, natural gas resource abundance R of ith natural gas well i The method comprises the following steps:
wherein R is i The resource abundance for the ith natural gas well, unit: m. ρ ij Rock density in the j th section of the reservoir in the i th natural gas well, in units: t/m 3 。H ij The thickness of the j section of the reservoir in the ith natural gas well is given in units of: m. C (C) ij The unit is the gas content of the j th section of the reservoir in the i th natural gas well: m is m 3 /t。
Step S1302, determining the natural gas resource abundance average value and the natural gas resource abundance standard deviation of the to-be-evaluated block based on the natural gas resource abundance of all the natural gas wells.
Acquiring a natural gas resource abundance data database R= (R) of all wells 1 ，R 2 ，R 3 ，……，R n ) Further calculating an average of the natural gas resource abundance of these wells(unit: m) and natural gas resource abundance standard deviation sigma R (unit: m).
Natural gas resource abundance average(unit: m) is:
standard deviation sigma of natural gas resource abundance R (unit: m) is:
step S1303, based on the natural gas resource abundance average value and the natural gas resource abundance standard deviation of the block to be evaluated, analyzing the natural gas resource abundance sampling normal distribution of the block to be evaluated to determine a probability density function of the natural gas resource abundance in the block to be evaluated.
The natural gas resource abundance is analyzed by adopting normal distribution, and the probability density function is expressed as follows:
step S140, determining an average value and a standard deviation of the total helium percentage of the block to be evaluated based on the total helium percentage of each sampling sample, and analyzing the normal distribution of the total helium percentage of the block to be evaluated based on the average value and the standard deviation of the total helium percentage of the block to be evaluated to determine a probability density function of the total helium percentage in the block to be evaluated.
Data base χ (He) = (χ (He)) for total helium percentage in all natural gas for all sampled samples of the block to be evaluated 1 ，χ(He) 2 ，χ(He) 3 ，……，χ(He) q ) On the basis of (a), further calculation of the average value of the percentage of total helium in the natural gas of these wellsAnd standard deviation sigma of total helium percentage χ(He) 。
Average value of percentage of total helium (He) associated with natural gasThe method comprises the following steps:
standard deviation sigma of percentage content of total helium associated with natural gas χ(He) The method comprises the following steps:
the natural gas associated total helium percentage content is analyzed by adopting normal distribution, and the probability density function is expressed as follows:
step S150, determining the associated helium resource quantity of the block to be evaluated under a set probability based on the probability density function of the gas-containing area of the reservoir, the abundance of the natural gas resource and the probability density function of the total helium percentage in the block to be evaluated.
Through the steps, the associated helium resource quantity under the set probability can be calculated on the basis of determining the probability density function of the gas-containing area of the reservoir, the abundance of the natural gas resource and the probability density function of the total helium percentage in the block to be evaluated. The specific steps are as follows:
and step S1501, integrating the probability density function of the natural gas resource abundance and the probability density function of the total helium percentage content respectively, and determining the natural gas resource abundance value and the total helium percentage content value under the set probability.
Here, the probability is assumed to be P, and the value range is [0%,100% ].
Enabling the integrated value of the natural gas resource abundance probability density function to be equal to P, and obtaining the corresponding natural gas resource abundance R under the probability P condition P . In the following formula, the value range of R is [ R ] min ，R max ]。
/>
The corresponding natural gas resource abundance R under the probability P can be determined through calculation P 。
When the integral value of the probability density function of the total helium percentage is equal to P, the corresponding total helium percentage χ (He) under the probability P condition can be obtained P . In the following formula, the value range of χ (He) is χ (He) min ，χ(He) max ]。
The corresponding total helium percentage χ (He) under the probability P can be determined by calculation P 。
Step S1502, determining the associated helium resource amount of the block to be evaluated under the set probability based on the natural gas resource abundance value and the total helium percentage content value under the set probability, the gas-containing area of the reservoir in the block to be evaluated, and a preset resource amount calculation formula under the set probability.
Preset resource quantity calculating formula Q (X) P The method comprises the following steps:
Q(He) P ＝χ(He) P ×R P ×S；
wherein, χ (X) P R is the percentage of He gas under probability P P The probability is the abundance of natural gas resources under P, and S is the gas-containing area of the reservoir.
And (3) carrying the natural gas resource abundance value and the total helium percentage content value under the set probability P into a preset resource quantity calculation formula to determine the associated helium resource quantity of the to-be-evaluated block under the set probability.
In some embodiments, the natural gas accompanies helium with helium primarily from the veil source 3 Helium of He and Shell Source 4 He causes, so that data calculation can be performed on helium with two causes, and understanding of potential of natural gas associated helium can be deepened. The resource quantity of the valance source helium and the resource quantity of the shell source helium under the condition that the parameter probability is P can be calculated respectively, and the specific calculation process is as follows:
in step S210, in addition to the data in step S110, the helium isotope ratio of each sample needs to be obtained at the same time when the sample data is obtained.
Helium isotope ratio dataset(dimensionless).
Step S220, determining the percentage of shell source helium and/or the percentage of mantle source helium of each sampling sample based on the total helium percentage and the helium isotope ratio of each sampling sample.
Whether the resource amount of the curtain source helium and the resource amount of the shell source helium need to be calculated at the same time can be determined according to actual application scenes, and the calculation mode of the resource amount of the curtain source helium and the resource amount of the shell source helium can be introduced together.
The percent shell source helium for each sample may be first determined based on a total helium percentage, helium isotope ratio, and shell source helium percentage determination formula for each sample.
Shell source helium percentage determination formula χ (4 he ) The method comprises the following steps:
chi (He) is the percentage of total helium He,is helium isotope 3 He/ 4 He ratio.
The veil source helium percentage for each sample is then determined based on the total helium percentage, helium isotope ratio, and veil source helium percentage determination formula for each sample.
Determination formula χ of percentage of valance source helium (3 He ) The method comprises the following steps:
step S230, determining the average value and the standard deviation of the percentage content of the shell source helium of the block to be evaluated based on the percentage content of the shell source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the shell source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the shell source helium of the block to be evaluated to determine the probability density function of the percentage content of the shell source helium in the block to be evaluated.
Obtaining and researching a data base χ of percentage content of shell source helium in all natural gas 4 He)＝(χ( 4 He) 1 ，χ( 4 He) 2 ，χ( 4 He) 3 ，……，χ( 4 He) q ) On the basis of (1) further calculating the average value of the percentage of shell source helium in the natural gas of the wellsAnd shell source helium percentage standard deviation sigma χ(4He) 。
Helium with shell source 4 Average value of He percentageThe method comprises the following steps:
standard deviation sigma of percentage content of shell source helium χ(4He) The method comprises the following steps:
the shell source helium percentage content is analyzed by adopting normal distribution, and the probability density function is expressed as follows:
step S240, determining the average value and the standard deviation of the percentage content of the curtain source helium of the block to be evaluated based on the percentage content of the curtain source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the curtain source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the curtain source helium of the block to be evaluated to determine the probability density function of the percentage content of the curtain source helium in the block to be evaluated.
Obtaining and researching X of a valance source helium percentage content data base in all natural gases 3 He)＝(χ( 3 He) 1 ，χ( 3 He) 2 ，χ( 3 He) 3 ，……，χ( 3 He) q ) On the basis of (1) further calculating the average value of the percentage of the curtain source helium in the natural gas of the wellsStandard deviation sigma of sum valance source helium percentage χ(3He) 。
Curtain source helium 3 Average value of He percentageThe method comprises the following steps:
standard deviation sigma of valance source helium percentage χ(3He) The method comprises the following steps:
the natural gas associated valance source helium percentage content is analyzed by adopting normal distribution, and the probability density function is expressed as follows:
and S250, integrating the probability density function of the shell source helium percentage and the probability density function of the mantle source helium percentage respectively, and determining the shell source helium percentage value and the mantle source helium percentage value under the set probability.
The probability density function integral value of the percentage content of the shell source helium is equal to P, and the probability P condition can be obtainedCorresponding shell source helium percentage χ% 4 He) P . In the following formula, χ is% 4 He) is within the range of [ χ ] 4 He) min ，χ( 4 He) max ]。
Solving the above method to determine the percent χ of the shell source helium under the probability P 4 He) P 。
Enabling the probability density function integral value of the percentage content of the valance source helium to be equal to P, and obtaining the corresponding valance source helium percentage content χ under the probability P condition 3 He) P . In the following formula, χ is% 3 He) is within the range of [ χ ] 3 He) min ，χ( 3 He) max ]。
Solving the above method to determine the X percentage content of valance source helium under the probability P 3 He) P 。
Step S260, determining the shell source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and the shell source helium percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability.
The probability of the parameters is the helium of the shell source under the condition of P 4 He) resource amount calculation formula is:
Q P (4 He )＝χ(4 He ) P *R P *S；
p is probability value, its value range is [0%,100%]；Q( 4 He) P Is shell source helium under the condition of probability P 4 He) amount of resources, unit: m is m 3 ；χ( 4 He) P Is the percentage of shell source helium with probability P; r is R P For the day with probability PAbundance of natural gas resources, unit: m; s is the gas-containing area of the reservoir, unit: m is m 2 。
And inputting the natural gas resource abundance value and the shell source helium percentage content value under the set probability P into a preset resource quantity calculation formula to determine the shell source helium resource quantity of the block to be evaluated under the set probability.
Step S270, determining the valance source helium gas resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and the valance source helium gas percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated, and a preset resource quantity calculation formula under the set probability.
Curtain source helium under the condition of P parameter probability 3 He) resource amount calculation formula is:
Q P (3 He )＝χ(3 He ) P *R P *S；
in the formula, P is a probability value, and the value range is [0%,100%]；Q( 3 He) P Is valance source helium under the condition of probability P 3 He) amount of resources, unit: m is m 3 ；χ( 3 He) P Is the percentage of valance source helium with probability P; r is R P The probability is the abundance of natural gas resources under the condition of P, and the unit is: m; s is the gas-containing area of the reservoir, unit: m is m 2 。
And inputting the natural gas resource abundance value and the valance source helium percentage content value under the set probability P into a preset resource quantity calculation formula to determine the valance source helium resource quantity of the block to be evaluated under the set probability.
Through the steps, the total helium (He) and valance source helium gas under different probability conditions can be obtained 3 He) and helium from shell source 4 He) resource amount, and further has more comprehensive and scientific knowledge on natural gas associated helium resources in a research area.
The invention provides an evaluation method, equipment and a storage medium. Next, the gas bearing area of the block to be evaluated is determined based on the reservoir thickness and reservoir gas content of each gas well, and the floor contour map. Next, a probability density function of the abundance of natural gas resources within the zone to be evaluated is determined based on the reservoir rock density, reservoir gas content, and reservoir thickness of the same reservoir section for each natural gas well. Then, based on the total helium percentage content of each sampling sample, determining an average value and a standard deviation of the total helium percentage content of the block to be evaluated, and based on the average value and the standard deviation of the total helium percentage content of the block to be evaluated, analyzing the normal distribution of the total helium percentage content of the block to be evaluated to determine a probability density function of the total helium percentage content in the block to be evaluated. And finally, determining the associated helium resource quantity of the block to be evaluated under a set probability based on the probability density function of the gas-containing area of the reservoir layer, the abundance of the natural gas resource and the probability density function of the total helium percentage in the block to be evaluated.
According to the invention, on the basis of acquiring the sampling data of the plurality of natural gas wells in the block to be evaluated, the natural gas resource abundance, the gas-containing area of the reservoir and the percentage content of helium are subjected to statistical analysis, so that the helium resource quantity calculation result under the set probability can be calculated, and the accuracy of the helium resource quantity evaluation result is improved.
According to the invention, the evaluation parameters are determined on the basis of acquiring the statistical results of different parameters and combining with the geological rules, and the resource amount calculation work can be carried out on the basis of knowing the geological and experimental result rules, so that the calculation result has strong geological and statistical basis and basis. By considering the difference of the grasping degree of the natural gas data and the helium gas data, the natural gas is statistically analyzed by taking the well as a unit, and the helium gas is counted by focusing on the test value of the sampling data in the whole region to be evaluated, so that the evaluation result has stronger pertinence.
In addition, the invention can further obtain the resource quantity parameters of natural gas and associated helium gas simultaneously through calculation and statistics of the abundance of natural gas resources and statistical analysis of total helium, shell source helium gas and valance source helium gas. The helium resources of two causes of the shell source and the veil source are distinguished, so that the knowledge of the helium resource amount is more comprehensive.
It should be understood that the sequence number of each step in the foregoing embodiment does not mean that the execution sequence of each process should be determined by the function and the internal logic, and should not limit the implementation process of the embodiment of the present invention.
Based on the method for evaluating the associated helium gas resource amount provided by the embodiment, correspondingly, the invention further provides a specific implementation mode of the device for evaluating the associated helium gas resource amount, which is applied to the method for evaluating the associated helium gas resource amount. Please refer to the following examples.
As shown in fig. 2, there is provided an apparatus 200 for evaluating an associated helium gas resource amount, the apparatus comprising:
an acquisition module 210, configured to acquire sampling data of a plurality of natural gas wells in a block to be evaluated; the sampling data comprise reservoir rock density, reservoir gas content and reservoir thickness of each natural gas well, and total helium percentage content of each sampling sample;
a first determining module 220, configured to determine a gas-containing area of the block to be evaluated based on a reservoir thickness and a reservoir gas content of each natural gas well, and a bottom plate contour map;
a second determination module 230 for determining a probability density function of natural gas resource abundance within the block to be evaluated based on reservoir rock density, reservoir gas content, and reservoir thickness for the same reservoir interval for each natural gas well;
A third determining module 240, configured to determine an average value and a standard deviation of the total helium percentage of the block to be evaluated based on the total helium percentage of each sampling sample, and analyze a normal distribution of the total helium percentage of the block to be evaluated based on the average value and the standard deviation of the total helium percentage of the block to be evaluated, so as to determine a probability density function of the total helium percentage in the block to be evaluated;
a fourth determining module 250, configured to determine an associated helium resource amount of the block to be evaluated under a set probability based on a probability density function of a reservoir gas area, a natural gas resource abundance, and a probability density function of a total helium percentage in the block to be evaluated.
In one possible implementation, the sampled data further includes helium isotope ratios for each sampled sample;
the evaluation device further includes: a third determination module 240 for determining a percentage of shell source helium and/or a percentage of veil source helium for each sample based on the total helium percentage and helium isotope ratio for each sample;
determining an average value and a standard deviation of the percentage content of the shell source helium of the block to be evaluated based on the percentage content of the shell source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the shell source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the shell source helium of the block to be evaluated to determine a probability density function of the percentage content of the shell source helium in the block to be evaluated; and/or
Determining an average value and a standard deviation of the percentage content of the curtain source helium of the block to be evaluated based on the percentage content of the curtain source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the curtain source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the curtain source helium of the block to be evaluated to determine a probability density function of the percentage content of the curtain source helium in the block to be evaluated.
In one possible implementation, the third determining module 240 is configured to determine the shell source helium percentage for each sample based on a total helium percentage, helium isotope ratio, and shell source helium percentage determination formula for each sample;
the percentage of veil source helium for each sample is determined based on the total percentage helium for each sample, the helium isotope ratio, and the veil source helium percentage determination formula.
In one possible implementation, the shell source helium percentage determination formula χ (4 He ) The method comprises the following steps:
determination formula χ of percentage of valance source helium (3 He ) The method comprises the following steps:
the probability density function f (χ (He)) of the total helium percentage in the block to be evaluated is:
probability density function f (χ (4) He ) Is) is:
probability density function f (χ (3) He ) Is) is:
wherein,and->The average value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage is respectively; sigma (sigma) χ(He) 、σ χ(4He) 、σ χ(3He) The standard deviation of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; χ (He) i The value range of (C) is [ χ (He) min ，χ(He) max ]，χ(4 He ) i The value range of (X) is% 4 He) min ，χ( 4 He) max ]，χ(3 He ) i The value range of (X) is% 3 He) min ，χ( 3 He) max ]，χ(He) min 、χ(He) max 、χ( 4 He) min 、χ( 4 He) max 、χ( 3 He) min 、χ( 3 He) max The minimum value and the maximum value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; chi (He) is the percentage of total helium He; />Is helium isotope 3 He/ 4 He ratio.
In one possible implementation, the fourth determining module 250 is configured to integrate the probability density function of the abundance of the natural gas resource and the probability density function of the total helium percentage, respectively, to determine the abundance value of the natural gas resource and the total helium percentage value at the set probability;
based on the natural gas resource abundance value and the total helium percentage content value under the set probability, determining the associated helium resource quantity of the block to be evaluated under the set probability according to the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability;
Preset resource quantity calculating formula Q (X) P The method comprises the following steps:
Q(X) P ＝χ(X) P ×R P ×S；
wherein X is gas, χ (X) P R is the percentage of X gas under the probability P P The probability is the abundance of natural gas resources under P, S is the gas-containing area of the reservoir, P is the probability value, and the value range is [0%,100%]。
In one possible implementation, the evaluation device further includes:
a fourth determining module 250, configured to integrate the probability density function of the percentage content of shell source helium and/or the probability density function of the percentage content of mantle source helium, and determine a shell source helium percentage content value and/or a mantle source helium percentage content value under a set probability;
determining the shell source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and the shell source helium percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated, and a preset resource quantity calculation formula under the set probability; and/or
And determining the valance source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and valance source helium percentage content value under the set probability, the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability.
In one possible implementation, the first determining module 220 is configured to determine a bottom plate contour map of the block to be evaluated according to the drilling hole elevation of each natural gas well and the depth data of each natural gas well;
Drawing a reservoir thickness contour map of a block to be evaluated based on the reservoir thickness of each natural gas well;
drawing a gas content contour map of a block to be evaluated based on the gas content of a reservoir of each natural gas well;
and superposing the reservoir thickness contour map and the air content contour map on the bottom plate contour map, and determining the air content area of the block to be evaluated.
In one possible implementation, the second determining module 230 is configured to determine the natural gas resource abundance of each natural gas well based on the reservoir rock density, the reservoir gas content, and the reservoir thickness of the same reservoir interval of each natural gas well;
determining natural gas resource abundance average values and natural gas resource abundance standard deviations of the to-be-evaluated blocks based on the natural gas resource abundance of all natural gas wells;
analyzing the natural gas resource abundance sampling normal distribution of the block to be evaluated based on the natural gas resource abundance average value and the natural gas resource abundance standard deviation of the block to be evaluated so as to determine a probability density function of the natural gas resource abundance in the block to be evaluated;
wherein, natural gas resource abundance R of ith natural gas well i The method comprises the following steps:
the probability density function f (R) of the abundance of natural gas resources is:
wherein m is the total number of sections of the ith natural gas well, ρ ij For the rock density of the j th section of the reservoir in the i th natural gas well, H ij C is the thickness of the j section of the reservoir in the ith natural gas well ij The gas content of the jth section of the reservoir in the ith natural gas well; sigma (sigma) R Is the standard deviation of the abundance of natural gas resources,the value range of R is [ R ] which is the average value of the abundance of natural gas resources min ，R max ]，R min Is the minimum value of natural gas resource abundance, R max Is the maximum value of natural gas resource abundance.
Fig. 3 is a schematic diagram of an electronic device according to an embodiment of the present invention. As shown in fig. 3, the electronic apparatus 3 of this embodiment includes: a processor 30, a memory 31 and a computer program 32 stored in said memory 31 and executable on said processor 30. The processor 30, when executing the computer program 32, performs the steps of the above-described embodiments of the associated helium gas resource amount evaluation method, such as steps 110 through 150 shown in fig. 1. Alternatively, the processor 30 may perform the functions of the modules of the apparatus embodiments described above, such as the functions of the modules 210-250 of fig. 2, when executing the computer program 32.
Illustratively, the computer program 32 may be partitioned into one or more modules that are stored in the memory 31 and executed by the processor 30 to complete the present invention. The one or more modules may be a series of computer program instruction segments capable of performing the specified functions for describing the execution of the computer program 32 in the electronic device 3. For example, the computer program 32 may be partitioned into modules 210 through 250 shown in FIG. 2.
The electronic device 3 may include, but is not limited to, a processor 30, a memory 31. It will be appreciated by those skilled in the art that fig. 3 is merely an example of the electronic device 3 and does not constitute a limitation of the electronic device 3, and may include more or fewer components than shown, or may combine certain components, or different components, e.g., the electronic device may further include an input-output device, a network access device, a bus, etc.
The processor 30 may be a central processing unit (Central Processing Unit, CPU), other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field-programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 31 may be an internal storage unit of the electronic device 3, such as a hard disk or a memory of the electronic device 3. The memory 31 may be an external storage device of the electronic device 3, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card) or the like, which are provided on the electronic device 3. Further, the memory 31 may also include both an internal storage unit and an external storage device of the electronic device 3. The memory 31 is used for storing the computer program and other programs and data required by the electronic device. The memory 31 may also be used for temporarily storing data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, specific names of the functional units and modules are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and in part, not described or illustrated in any particular embodiment, reference is made to the related descriptions of other embodiments.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus/electronic device and method may be implemented in other manners. For example, the apparatus/electronic device embodiments described above are merely illustrative, e.g., the division of the modules or units is merely a logical function division, and there may be additional divisions in actual implementation, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection via interfaces, devices or units, which may be in electrical, mechanical or other forms.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated modules/units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present invention may be implemented in whole or in part by a computer program for instructing related hardware to perform the steps of the above-described method embodiment of the associated helium gas resource amount evaluation method, wherein the computer program may be stored in a computer readable storage medium, and the computer program may be executed by a processor. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth.
The above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention, and are intended to be included in the scope of the present invention.

Claims (9)

1. An associated helium gas resource quantity evaluation method is characterized by comprising the following steps:
acquiring sampling data of a plurality of natural gas wells in a block to be evaluated; wherein the sampled data includes reservoir rock density, reservoir gas content, and reservoir thickness for each gas well, and total helium percentage for each sampled sample;
determining the gas-containing area of the block to be evaluated based on the reservoir thickness and the reservoir gas content of each natural gas well and a bottom plate contour map;
determining a probability density function of natural gas resource abundance in the block to be evaluated based on reservoir rock density, reservoir gas content and reservoir thickness of the same reservoir section of each natural gas well;
Determining an average value and a standard deviation of the total helium percentage of the block to be evaluated based on the total helium percentage of each sampling sample, and analyzing the normal distribution of the total helium percentage of the block to be evaluated based on the average value and the standard deviation of the total helium percentage of the block to be evaluated to determine a probability density function of the total helium percentage in the block to be evaluated;
determining the associated helium resource quantity of the block to be evaluated under a set probability based on the reservoir gas-containing area in the block to be evaluated, the probability density function of the natural gas resource abundance and the probability density function of the total helium percentage;
wherein the determining the associated helium resource amount of the block to be evaluated under a set probability based on the probability density function of the reservoir gas-containing area, the natural gas resource abundance, and the probability density function of the total helium percentage content in the block to be evaluated comprises:
integrating the probability density function of the natural gas resource abundance and the probability density function of the total helium percentage content respectively, and determining the natural gas resource abundance value and the total helium percentage content value under the set probability; based on the natural gas resource abundance value and the total helium percentage content value under the set probability, determining the associated helium resource quantity of the block to be evaluated under the set probability according to the reservoir gas-containing area in the block to be evaluated and a preset resource quantity calculation formula under the set probability;
The preset resource amount calculation formula Q (X) P The method comprises the following steps:
Q(X) P ＝χ(X) P ×R P ×S；
wherein X is gas, χ (X) P R is the percentage of X gas under the probability P P The probability is the abundance of natural gas resources under P, S is the gas-containing area of the reservoir, P is the probability value, and the value range is [0%,100%]。
2. The method of evaluating of claim 1, wherein the sampled data further comprises a helium isotope ratio for each sampled sample;
the evaluation method further includes:
determining the percentage of shell source helium and/or the percentage of mantle source helium of each sampling sample based on the total helium percentage and helium isotope ratio of each sampling sample;
determining an average value and a standard deviation of the percentage content of the shell source helium of the block to be evaluated based on the percentage content of the shell source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the shell source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the shell source helium of the block to be evaluated to determine a probability density function of the percentage content of the shell source helium in the block to be evaluated; and/or
And determining the average value and the standard deviation of the percentage content of the curtain source helium of the block to be evaluated based on the percentage content of the curtain source helium of each sampling sample, and analyzing the normal distribution of the percentage content of the curtain source helium of the block to be evaluated based on the average value and the standard deviation of the percentage content of the curtain source helium of the block to be evaluated so as to determine the probability density function of the percentage content of the curtain source helium in the block to be evaluated.
3. The method of evaluating of claim 2, wherein determining the percentage of shell source helium and/or the percentage of veil source helium for each sample based on the total helium percentage and helium isotope ratio for each sample comprises:
determining the percentage content of the shell source helium of each sampling sample based on the total percentage content of helium, the helium isotope ratio and the percentage content of the shell source helium of each sampling sample;
determining the percentage of the mantle source helium for each sample based on the formula for determining the percentage of total helium, the ratio of helium isotopes and the percentage of mantle source helium for each sample.
4. The method of claim 3, wherein said shell source helium percentage determination formula χ (4 He ) The method comprises the following steps:
determination formula χ of percentage of veil source helium (3 He ) The method comprises the following steps:
the probability density function f (χ (He)) of the total helium percentage in the block to be evaluated is:
probability density of shell source helium percentage in the block to be evaluatedFunction f (χ (4) He ) Is) is:
probability density function f (χ (3) Hr ) Is) is:
Wherein,and->The average value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage is respectively; sigma (sigma) χ(He) 、σ χ(4He) 、σ χ(3He) The standard deviation of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; χ (He) i The value range of (C) is [ χ (He) min ，χ(He) max ]，χ(4 He ) i The value range of (X) is% 4 He) min ，χ( 4 He) max ]，χ(3 He ) i The value range of (X) is% 3 He) min ，χ( 3 He) max ]，χ(He) min 、χ(He) max 、χ( 4 He) min 、χ( 4 He) max 、χ( 3 He) min 、χ( 3 He) max The minimum value and the maximum value of the total helium percentage, the shell source helium percentage and the mantle source helium percentage are respectively; chi (He) is the percentage of total helium He; />Is helium isotope 3 He/ 4 He ratio.
5. The evaluation method according to claim 2, characterized in that the evaluation method further comprises:
integrating the probability density function of the shell source helium percentage and/or the probability density function of the mantle source helium percentage to determine a shell source helium percentage value and/or a mantle source helium percentage value under the set probability;
determining the shell source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and the shell source helium percentage value under the set probability, the reservoir gas-containing area in the block to be evaluated, and the preset resource quantity calculation formula under the set probability; and/or
And determining the curtain source helium resource quantity of the block to be evaluated under the set probability based on the natural gas resource abundance value and the curtain source helium percentage value under the set probability, and the reservoir gas-containing area in the block to be evaluated and the preset resource quantity calculation formula under the set probability.
6. The method of evaluating any one of claims 1-5, wherein the determining the gas bearing area of the zone to be evaluated based on the reservoir thickness and reservoir gas content of each natural gas well and a floor contour map comprises:
determining a bottom plate contour map of the block to be evaluated according to the elevation of the drilling hole of each natural gas well and the depth data of each natural gas well;
drawing a reservoir thickness contour map of the block to be evaluated based on the reservoir thickness of each natural gas well;
drawing a gas content contour map of the block to be evaluated based on the reservoir gas content of each natural gas well;
and superposing the reservoir thickness contour map and the gas content contour map on the bottom plate contour map, and determining the gas content area of the block to be evaluated.
7. The method of evaluating any of claims 1-5, wherein the determining a probability density function of natural gas resource abundance within the zone to be evaluated based on reservoir rock density, reservoir gas content, and reservoir thickness for the same reservoir interval for each of the natural gas wells comprises:
Determining natural gas resource abundance of each natural gas well based on reservoir rock density, reservoir gas content and reservoir thickness of the same reservoir section of each natural gas well;
determining natural gas resource abundance averages and natural gas resource abundance standard deviations of the to-be-evaluated blocks based on natural gas resource abundance of all natural gas wells;
analyzing the natural gas resource abundance sampling normal distribution of the block to be evaluated based on the natural gas resource abundance average value and the natural gas resource abundance standard deviation of the block to be evaluated so as to determine a probability density function of the natural gas resource abundance in the block to be evaluated;
wherein, natural gas resource abundance R of ith natural gas well i The method comprises the following steps:
the probability density function f (R) of the abundance of natural gas resources is:
wherein m is the total number of sections of the ith natural gas well, ρ ij For the rock density of the j th section of the reservoir in the i th natural gas well, H ij C is the thickness of the j section of the reservoir in the ith natural gas well ij The gas content of the jth section of the reservoir in the ith natural gas well; sigma (sigma) R Is the standard deviation of the abundance of natural gas resources,the value range of R is [ R ] which is the average value of the abundance of natural gas resources min ，R max ]，R min Is the minimum value of natural gas resource abundance, R max Is the maximum value of natural gas resource abundance.
8. An electronic device comprising a memory for storing a computer program and a processor for invoking and running the computer program stored in the memory to perform the method of any of claims 1 to 7.
9. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the method according to any one of claims 1 to 7.
CN202310967793.3A 2023-08-02 2023-08-02 Method, equipment and storage medium for evaluating associated helium resource quantity Active CN117150207B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310967793.3A CN117150207B (en) 2023-08-02 2023-08-02 Method, equipment and storage medium for evaluating associated helium resource quantity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310967793.3A CN117150207B (en) 2023-08-02 2023-08-02 Method, equipment and storage medium for evaluating associated helium resource quantity

Publications (2)

Publication Number Publication Date
CN117150207A CN117150207A (en) 2023-12-01
CN117150207B true CN117150207B (en) 2024-01-26

Family

ID=88907036

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310967793.3A Active CN117150207B (en) 2023-08-02 2023-08-02 Method, equipment and storage medium for evaluating associated helium resource quantity

Country Status (1)

Country Link
CN (1) CN117150207B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105528522A (en) * 2015-12-15 2016-04-27 中国石油天然气股份有限公司 Method and device for calculating quantity of resources of continuous oil and gas reservoir based on space grids
CN110322094A (en) * 2018-03-30 2019-10-11 中国石油化工股份有限公司 A method of evaluation shale petroleum resources
CN110991933A (en) * 2019-12-19 2020-04-10 西南石油大学 Mountain shale gas resource evaluation method and system
CN114910976A (en) * 2022-04-18 2022-08-16 中国科学院西北生态环境资源研究院 Geological evaluation method for helium resource potential in low-exploration-degree area
CN115221675A (en) * 2022-05-18 2022-10-21 中国石油天然气集团有限公司 Helium gas resource scale sequence determination method, device and equipment
CN116047624A (en) * 2022-12-29 2023-05-02 中国地质调查局西安地质调查中心(西北地质科技创新中心) Helium resource amount evaluation method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105528522A (en) * 2015-12-15 2016-04-27 中国石油天然气股份有限公司 Method and device for calculating quantity of resources of continuous oil and gas reservoir based on space grids
CN110322094A (en) * 2018-03-30 2019-10-11 中国石油化工股份有限公司 A method of evaluation shale petroleum resources
CN110991933A (en) * 2019-12-19 2020-04-10 西南石油大学 Mountain shale gas resource evaluation method and system
CN114910976A (en) * 2022-04-18 2022-08-16 中国科学院西北生态环境资源研究院 Geological evaluation method for helium resource potential in low-exploration-degree area
CN115221675A (en) * 2022-05-18 2022-10-21 中国石油天然气集团有限公司 Helium gas resource scale sequence determination method, device and equipment
CN116047624A (en) * 2022-12-29 2023-05-02 中国地质调查局西安地质调查中心(西北地质科技创新中心) Helium resource amount evaluation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
致密气储量计算方法研究;陈晓智;;中国矿业(第06期);全文 *

Also Published As

Publication number Publication date
CN117150207A (en) 2023-12-01

Similar Documents

Publication Publication Date Title
CN103345566B (en) Based on the geochemical anomaly discrimination and evaluation method of Geological Connotation
CN106355011A (en) Geochemical data element sequence structure analysis method and device
CN109113729B (en) Lithology identification method and device based on well logging curve
Schiffelbein et al. Direct assessment of stable isotope variability in planktonic foraminifera populations
CN107764974A (en) A kind of granite type U-ore is into ore deposit incident heat dating method
CN106649339A (en) Method and device for mining outlier
CN111413741B (en) Sandstone-type uranium ore resource amount calculation method and device
CN114240212A (en) Method and equipment for determining influence weight of geological parameters on resource quantity
CN117150207B (en) Method, equipment and storage medium for evaluating associated helium resource quantity
WO2024012222A1 (en) Method and apparatus for selecting key geological parameter of block to be surveyed, and device
CN116071191A (en) Quantitative characterization method for contribution share of helium source rock in stabilized Keraton basin helium-rich field
CN110895704A (en) Microorganism dune complex reservoir type identification method and device and storage medium
CN115146976A (en) Method and device for selecting main geological parameters influencing evaluation of block to be surveyed
CN111710157B (en) Method for extracting hot spot area of taxi
CN114186879A (en) Method and equipment for evaluating influence of geological parameters on resource quantity calculation errors
CN114970135A (en) Sandstone proportion determination method and system, electronic device and storage medium
CN116842327B (en) Method, device and equipment for processing abnormal data in resource quantity evaluation
RU2639649C2 (en) Computer system for processing and analysis of geophysical data
CN112630839A (en) Well logging curve standardization method and system
CN116108368B (en) Deposition microphase identification method and device based on deep learning mixed model
CN116842327A (en) Method, device and equipment for processing abnormal data in resource quantity evaluation
CN109143334B (en) Level interpolation method, device and system
CN117575106B (en) Method, system, electronic equipment and medium for predicting gas production profile of coal-bed gas well
Xu et al. A quadtree-based batch method for local singularity mapping and application for geochemical anomaly identification
CN118039029A (en) Method and system for identifying granite type based on machine learning and zircon component

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant