CN112382346A - Method and device for determining oil saturation, electronic device and storage medium - Google Patents

Method and device for determining oil saturation, electronic device and storage medium Download PDF

Info

Publication number
CN112382346A
CN112382346A CN202011418520.6A CN202011418520A CN112382346A CN 112382346 A CN112382346 A CN 112382346A CN 202011418520 A CN202011418520 A CN 202011418520A CN 112382346 A CN112382346 A CN 112382346A
Authority
CN
China
Prior art keywords
shale
oil
gas
determining
certain temperature
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011418520.6A
Other languages
Chinese (zh)
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 National Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging Co Ltd
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 National Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN202011418520.6A priority Critical patent/CN112382346A/en
Publication of CN112382346A publication Critical patent/CN112382346A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/30Prediction of properties of chemical compounds, compositions or mixtures

Landscapes

  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)

Abstract

The disclosure relates to a method and a device for determining oil saturation, an electronic device and a storage medium. The method for determining the oil saturation comprises the following steps: acquiring a preset shale oil-gas saturation model and measurement parameters thereof; and determining the oil saturation degree based on the measurement parameters and the shale oil and gas saturation degree model. Embodiments of the present disclosure may determine that the shale is saturated with oil.

Description

Method and device for determining oil saturation, electronic device and storage medium
Technical Field
The disclosure relates to the technical field of oil and gas exploration and development, and in particular to a method and a device for determining oil saturation, an electronic device and a storage medium.
Background
In the exploration and development of the oil field in Daqing at present, the shale oil-gas field gradually becomes an important field for oil-gas succession, the electrical relation of the reservoir rock is complex due to the problems of complex framework minerals, complex pore structure, complex pore wettability and the like of the reservoir, a non-Archie phenomenon occurs in the reservoir, the water saturation and the resistivity are not unique corresponding relation any more, factors such as the shale content, the pore structure, the wettability and the like can influence the water saturation, and therefore the Archie saturation model is not applicable any more. The evaluation of the oil-gas saturation is also the key of the evaluation of the shale oil-gas reservoir, so that the accurate quantitative calculation of the oil-gas saturation of the shale is an important link for assisting oil-gas discovery. The prior art can not realize the evaluation of the oil-gas saturation of the shale.
Disclosure of Invention
The disclosure provides a method and a device for determining oil saturation, electronic equipment and a storage medium technical scheme, which can determine the oil and gas saturation of shale and can indicate the oil saturation of the shale to a certain extent.
According to an aspect of the present disclosure, there is provided a method of determining oil saturation, including:
acquiring a preset shale oil-gas saturation model and measurement parameters thereof;
and determining the oil saturation degree based on the measurement parameters and the shale oil and gas saturation degree model.
Preferably, before the preset shale hydrocarbon saturation model is obtained, the preset shale hydrocarbon saturation model needs to be established, and the method for establishing the preset shale hydrocarbon saturation model includes:
obtaining the mass of oil gas in the shale, the volume of the oil gas in the shale and the effective porosity;
and establishing the preset shale oil-gas saturation model according to the quality of the oil gas in the shale, the volume of the oil gas in the shale and the effective porosity.
Preferably, before the obtaining of the mass of the hydrocarbons in the shale and the volume of the hydrocarbons in the shale, the mass of the hydrocarbons in the shale and the volume of the hydrocarbons in the shale need to be determined;
the method for determining the quality of hydrocarbons in shale comprises the following steps:
acquiring logging information, the quality of shale and scale conversion coefficients;
determining the content volume of hydrocarbon in the shale at different temperatures and/or different temperature intervals according to the logging information;
determining the mass of oil and gas in the shale based on the content volume of hydrocarbon in the shale at different temperatures and/or different temperature intervals, the mass of the shale and a scale conversion coefficient;
and/or the presence of a gas in the interior of the container,
the method for determining the volume of the oil and gas in the shale comprises the following steps:
acquiring the mass and density of oil and gas in the shale;
and determining the volume of the oil and gas in the shale according to the mass of the oil and gas in the shale and the density of the oil and gas.
Preferably, the method for acquiring the quality of the shale comprises the following steps:
acquiring the density of a rock skeleton, the bulk density of shale, the density of fluid in rock, the mass of the rock skeleton and the mass of fluid;
determining the volume ratio of shale pores according to the density of the rock skeleton, the bulk density of the shale and the density of fluid in the rock;
and determining the mass of the shale according to the shale pore volume ratio, the mass of the rock skeleton and the mass of the fluid.
Preferably, the method for determining oil saturation further comprises: determining the content volume ratio of hydrocarbon in the shale at a certain temperature and/or a certain temperature interval;
and determining a shale oil-gas saturation model corresponding to a certain temperature and/or a certain temperature interval by using the content volume ratio of the hydrocarbons in the shale at the certain temperature and/or the certain temperature interval to represent the content volumes of the hydrocarbons in the shale at other temperatures and/or other temperature intervals.
Preferably, the measuring parameters include: a given scale conversion coefficient, a measured content volume and a ratio of hydrocarbons in the shale at a certain temperature and/or a certain temperature interval, a measured bulk density of the shale, a measured oil and gas density and a measured effective porosity;
and/or the presence of a gas in the interior of the container,
the method for measuring the content volume of the hydrocarbon in the shale at a certain temperature and/or a certain temperature interval comprises the following steps:
acquiring the relation of the logging sensitivity curves of the hydrocarbon content volume, the deep lateral resistivity, the compensation density, the compensation sound wave and the compensation neutron in the shale at a certain temperature and/or a certain temperature interval;
and obtaining the content volume of the hydrocarbons in the shale at a certain temperature and/or a certain temperature interval based on the current deep lateral resistivity, the compensation density, the compensation sound wave and the compensation neutrons and the relation of the logging sensitivity curve.
7. The method for determining oil saturation according to claim 6, characterized by:
before obtaining a log sensitivity curve relation of hydrocarbon content volume, deep lateral resistivity, compensation density, compensation sound wave and compensation neutrons in shale at a certain temperature and/or a certain temperature interval, the method needs to establish the log sensitivity curve relation, and comprises the following steps:
obtaining a volume of hydrocarbon content in shale at a temperature and/or temperature interval for a plurality of samples;
and establishing the well logging sensitivity curve relation according to the content volume of the hydrocarbons in the shale of the multiple samples at a certain temperature and/or a certain temperature interval and the corresponding deep lateral resistivity, compensation density, compensation sound wave and compensation neutrons of the hydrocarbons.
According to an aspect of the present disclosure, there is provided an oil saturation determination device, including:
the acquiring unit is used for acquiring a preset shale oil-gas saturation model and measurement parameters thereof;
and the determining unit is used for determining the oil saturation based on the measurement parameters and the shale oil-gas saturation model.
According to an aspect of the present disclosure, there is provided an electronic device including:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to: the above-described determination method is performed.
According to an aspect of the present disclosure, there is provided a computer-readable storage medium having stored thereon computer program instructions which, when executed by a processor, implement the above-described determination method.
In the embodiment of the disclosure, the shale oil-gas saturation can be determined, and the shale oil-gas saturation can be indicated to a certain extent, so that the problem that an existing oil-gas saturation model is not applicable due to the existence and complex wettability (organic pore oil is wet and inorganic pore water is wet) of organic matters, pyrite and other special minerals in a shale reservoir is solved.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.
Other features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments, which proceeds with reference to the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and, together with the description, serve to explain the principles of the disclosure.
FIG. 1 shows a flow chart of a method of determining oil saturation according to an embodiment of the present disclosure;
FIG. 2 is a schematic diagram illustrating the effect of comparing S1 of well logging calculation with measured data according to an embodiment of the disclosure;
FIG. 3 is a schematic diagram illustrating a method for calculating the c value of the S1 ratio in (S0+ S1+ S2) according to an embodiment of the present disclosure;
FIG. 4 is a schematic diagram illustrating an actual well data processing effect according to an embodiment of the disclosure;
FIG. 5 is a block diagram illustrating an electronic device 800 in accordance with an exemplary embodiment;
fig. 6 is a block diagram illustrating an electronic device 1900 according to an example embodiment.
Detailed Description
Various exemplary embodiments, features and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. In the drawings, like reference numbers can indicate functionally identical or similar elements. While the various aspects of the embodiments are presented in drawings, the drawings are not necessarily drawn to scale unless specifically indicated.
The word "exemplary" is used exclusively herein to mean "serving as an example, embodiment, or illustration. Any embodiment described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other embodiments.
The term "and/or" herein is merely an association describing an associated object, meaning that three relationships may exist, e.g., a and/or B, may mean: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the term "at least one" herein means any one of a plurality or any combination of at least two of a plurality, for example, including at least one of A, B, C, and may mean including any one or more elements selected from the group consisting of A, B and C.
Furthermore, in the following detailed description, numerous specific details are set forth in order to provide a better understanding of the present disclosure. It will be understood by those skilled in the art that the present disclosure may be practiced without some of these specific details. In some instances, methods, means, elements and circuits that are well known to those skilled in the art have not been described in detail so as not to obscure the present disclosure.
Fig. 1 shows a flow chart of a method of determining oil saturation according to an embodiment of the present disclosure. As shown in fig. 1, the method for determining oil saturation includes: step S101: acquiring a preset shale oil-gas saturation model and measurement parameters thereof; step S102: and determining the oil saturation degree based on the measurement parameters and the shale oil and gas saturation degree model. The shale oil-gas saturation can be determined, the shale oil-gas saturation can be indicated to a certain extent, and the problem that an existing oil-gas saturation model is not applicable due to the existence and complex wettability (organic pore oil is wet and inorganic pore water is wet) of organic matters, pyrite and other special minerals in a shale reservoir is solved.
The present disclosure is illustrated by using the example of the gulong shale rock mass consisting of the rock skeleton mass, the oil gas mass in the pores, and the adsorbed oil gas mass, and it is believed that the oil-containing volume of the gulong shale is mainly composed of free hydrocarbons and organic matters in the pores and adsorbed hydrocarbons on the clay surface, and the free hydrocarbons are related to S0 and S1 in the pyrolysis data, and the adsorbed hydrocarbons are related to S2. Wherein S0 is the content of hydrocarbon in rock detected at 90 ℃, is a gaseous hydrocarbon component of C1-C7, and exists in liquid hydrocarbon in a pure oil layer in a dissolved gas mode; s1 is the content of hydrocarbon in rock detected at 300 ℃, and is light and medium liquid hydrocarbon components of C8-C29; s2 is the content of pyrolytic hydrocarbon in the rock detected at 300-600 ℃, and is a heavy crude oil component and colloid and asphaltene which are more than C29. The volumes of S0, S1 and S2 in the rock pyrolysis are obtained through logging information, and are combined with the pore volume to establish a non-electric oil saturation degree calculation model, so that an oil-containing property discrimination method based on the rock pyrolysis information is formed.
In the gulonic shale reservoir, the rock mass is composed of the rock skeleton mass, the organic matter mass and the pore fluid mass, and the mesochite or mudstone does not contain oil gas; an oil saturation model (shale oil and gas saturation model) can be established by utilizing logging information such as neutrons, density, sound waves, resistivity, effective porosity and the like.
Wherein the clay pores, matrix pores, and lamellar gaps are hydrophilic, water-wet, the organic pores are oleophilic, oil-wet, and the pores are 100% oleophilic, i.e., have 100% oleophilic saturation. The porosity of the gulonic shale reservoir matrix is mainly residual native porosity and unstable mineral dissolution pores.
Step S101: and acquiring a preset shale oil-gas saturation model and measurement parameters thereof.
In this disclosure, before the preset shale oil and gas saturation model is obtained, the preset shale oil and gas saturation model needs to be established, and the method for establishing the preset shale oil and gas saturation model includes: obtaining the mass of oil gas in the shale, the volume of the oil gas in the shale and the effective porosity; and establishing the preset shale oil-gas saturation model according to the quality of the oil gas in the shale, the volume of the oil gas in the shale and the effective porosity.
The following description is provided to enable those skilled in the art to understand and implement the present disclosure.
In the present disclosure, before the obtaining of the mass of hydrocarbons in the shale and the volume of hydrocarbons in the shale, the mass of hydrocarbons in the shale and the volume of hydrocarbons in the shale need to be determined.
In the present disclosure, the method of determining the quality of hydrocarbons in the shale comprises: acquiring logging information, the quality of shale and scale conversion coefficients; determining the content volume of hydrocarbon in the shale at different temperatures and/or different temperature intervals according to the logging information; and determining the mass of the oil and gas in the shale based on the content volume of the hydrocarbon in the shale at different temperatures and/or different temperature intervals, the mass of the shale and the scale conversion coefficient.
In the present disclosure, the method of determining the volume of hydrocarbons in the shale comprises: acquiring the mass and density of oil and gas in the shale; and determining the volume of the oil and gas in the shale according to the mass of the oil and gas in the shale and the density of the oil and gas.
In the present disclosure, the method for obtaining the quality of the shale includes: acquiring the density of a rock skeleton, the bulk density of shale, the density of fluid in rock, the mass of the rock skeleton and the mass of fluid; determining the volume ratio of shale pores according to the density of the rock skeleton, the bulk density of the shale and the density of fluid in the rock; and determining the mass of the shale according to the shale pore volume ratio, the mass of the rock skeleton and the mass of the fluid.
The mass M of the shale is the sum of the mass of the rock skeleton and the mass of the fluid:
M=ρma*(Vsum-Vf)+ρf*Vfformula 1
And determining the void volume ratio of the shale according to the density of the rock skeleton, the bulk density of the shale and the density of fluid in the rock. Where ρ isbIs the bulk density of the shale, and has the unit of g/cm3;ρfIs the density of the fluid in the rock, in g/cm3
The shale pore volume ratio is approximately:
Figure BDA0002819245930000081
substituting formula 2 into formula 1 to obtain the mass M of the shale:
M=ρbVsumformula 3
In the present disclosure, the volume of the hydrocarbon content in the shale at different temperatures and/or different temperature intervals is determined according to the logging data as follows: the content S0 of the hydrocarbon in the rock detected at 90 ℃, the content S1 of the hydrocarbon in the rock detected at 300 ℃ and the content S2 of the pyrolysis hydrocarbon in the rock detected at 300-600 ℃ are the volume of the hydrocarbon in the pyrolysis of the rock; wherein the content of hydrocarbons in the rock, S1, detected at 300 ℃, is taken as the volume of the content of hydrocarbons in the shale below a certain temperature and/or a certain temperature interval.
Specifically, the mass of oil and gas in shale is as follows:
Mog=(S0*M+S1*M+S2m) k formula 4
The volume of oil and gas in the shale is
Figure BDA0002819245930000082
Wherein M is the total mass of the rock and is given in g; k is a scale conversion coefficient and is dimensionless; mogIs the total mass of hydrocarbons in the rock in g; vogIs the hydrocarbon volume in the rock, and the unit is cm3;ρhcIs the density of oil gas in g/cm3
In particular, according to the rock skeleton density ρmaAnd the volume of the rock skeleton to obtain the mass of the rock skeleton rhoma*(Vsum-Vf). Where ρ ismaIs the density of the rock skeleton in g/cm3;VsumIs the total volume of rock in cm3;VfIs the volume of fluid in the rock in cm3
Obtaining the fluid mass rho according to the fluid density and the fluid volume in the rockf*Vf. Where ρ isfIs the density of the fluid in the rock, in g/cm3;VfIs the volume of fluid in the rock in cm3
The oil saturation is defined as: the ratio of the oil volume in the reservoir effective pore to the rock effective pore volume.
The hydrocarbon saturation in the shale is
Figure BDA0002819245930000091
Bringing formula 3 into formula 6 to obtain rock petroleum gas saturation
Figure BDA0002819245930000092
In the present disclosure, the method for determining oil saturation further includes: determining the content volume ratio of hydrocarbon in the shale at a certain temperature and/or a certain temperature interval; and determining a shale oil-gas saturation model corresponding to a certain temperature and/or a certain temperature interval by using the content volume ratio of the hydrocarbons in the shale at the certain temperature and/or the certain temperature interval to represent the content volumes of the hydrocarbons in the shale at other temperatures and/or other temperature intervals.
Specifically, one temperature and/or one temperature interval may be 300 ℃, and the other temperatures and/or other temperature intervals may be 90 ℃ and 300-600 ℃ respectively. And calculating the hydrocarbon content S0 in the rock detected at 90 ℃, the hydrocarbon content S1 in the rock detected at 300 ℃ and the pyrolysis hydrocarbon content S2 in the rock detected at 300-600 ℃ in the rock pyrolysis through logging data, and combining the volume with the pore volume to finally obtain the gulong shale oil saturation calculation model (shale oil and gas saturation model).
Establishing different depth sections (or vitrinite reflectivity) S through enclosed coring rock pyrolysis data analysis1Ratio c in pyrolysis hydrocarbon component:
Figure BDA0002819245930000101
substituting the formula 8 into the formula 7 to obtain the gulonic shale oil and gas saturation model:
Figure BDA0002819245930000102
in the formulas 1-9, M is the total mass of the rock and is expressed in g; vsumIs the total volume of rock in cm3;VfIs the volume of fluid in the rock in cm3;ρmaIs the density of the rock skeleton in g/cm3;ρfIs the density of the fluid in the rock, in g/cm3;MogIs the total mass of hydrocarbons in the rock in g; vogIs the hydrocarbon volume in the rock, and the unit is cm3;SogThe saturation of the oil and gas is shown in unit; s0、S1、S2In order to pyrolyze different hydrocarbon components, the unit is mg/g, which is calculated by logging information in practical application; rhobIs the bulk density of the shale, and has the unit of g/cm3;ρhcIs the density of oil gas in g/cm3;φeEffective porosity in v/v; k is a scale conversion coefficient and is dimensionless;
step S102: and determining the oil saturation degree based on the measurement parameters and the shale oil and gas saturation degree model.
Based on the above, the measuring parameters include: the method comprises the following steps of determining a given scale conversion coefficient, measuring the content volume and the proportion of hydrocarbons in the shale at a certain temperature and/or a certain temperature interval, measuring the bulk density of the shale, measuring the oil gas density and measuring the effective porosity. And substituting the measurement parameters into a preset shale oil-gas saturation model to obtain the oil-gas saturation.
Wherein the given scale conversion coefficient k can be 0.001, and the measured bulk density of the shale is obtained by well logging acquisition for removing a carbonate rock interlayerb(ii) a Measured effective porosity phieThe effective porosity of nuclear magnetic resonance logging can be applied, and when the nuclear magnetic resonance logging is not carried out, the effective porosity processed by optimization software can be applied, and the multivariate regression effective porosity calibrated by core analysis data can also be applied; measuringOil and gas density of volume ρhcThe oil density measurement is taken for the oil density in the test oil sample of the gulong pit.
In the present disclosure, the method for measuring the content volume of hydrocarbons in shale at a certain temperature and/or a certain temperature interval comprises: acquiring the relation of the logging sensitivity curves of the hydrocarbon content volume, the deep lateral resistivity, the compensation density, the compensation sound wave and the compensation neutron in the shale at a certain temperature and/or a certain temperature interval; and obtaining the content volume of the hydrocarbons in the shale at a certain temperature and/or a certain temperature interval based on the current deep lateral resistivity, the compensation density, the compensation sound wave and the compensation neutrons and the relation of the logging sensitivity curve.
In the present disclosure, before obtaining a log sensitivity curve relationship of a hydrocarbon content volume, a deep lateral resistivity, a compensation density, a compensation sound wave and a compensation neutron in a shale at a certain temperature and/or a certain temperature interval, the log sensitivity curve relationship needs to be established, and the method includes: obtaining a volume of hydrocarbon content in shale at a temperature and/or temperature interval for a plurality of samples; and establishing the well logging sensitivity curve relation according to the content volume of the hydrocarbons in the shale of the multiple samples at a certain temperature and/or a certain temperature interval and the corresponding deep lateral resistivity, compensation density, compensation sound wave and compensation neutrons of the hydrocarbons.
For example, pyrolysis S can be calculated by in situ analysis of core pyrolysis data for multiple samples of the gulon recess 2000, by multiple regression to establish a log sensitivity curve relationship (equation 10)1Content, from measured pyrolysis S1The consistency is better by comparison analysis, R is20.8775, as shown in fig. 2. FIG. 2 is a schematic diagram illustrating the effect of comparing the measured data with S1 of the well logging calculation according to the embodiment of the disclosure.
S1=f(RtbnDt) formula 10
In the formula, RtIs the deep lateral resistivity, in units of Ω/m; rhobTo compensate for density, the unit is g/cm3(ii) a dt is the compensating sound wave, and the unit is mus/ft; phi is anTo compensate for neutrons, the unit is P.U.
Meanwhile, the content volume of the hydrocarbon in the shale in a certain temperature and/or a certain temperature interval accounts for c value to be calculated: and establishing the proportion of each block in the pyrolysis hydrocarbons by combining the data of the regional closed coring analysis and the data of the vitrinite reflectivity (maturity).
Fig. 3 illustrates a schematic diagram of a calculation method of the S1 fraction c value in (S0+ S1+ S2) according to an embodiment of the present disclosure. Is approximated to obtain (S)0+S1+S2) As can be seen from FIG. 3, X-axis data represents vitrinite reflectance Ro, and Y-axis data represents core analysis S1The proportion c in pyrolysis hydrocarbons, the data intersection form and the trend line can be known, the correlation coefficient reaches 0.68, the c value has a good corresponding relation with the vitrinite reflectivity Ro (maturity), and a corresponding obtaining formula can be obtained through regional closed coring data statistics.
Let k be 0.001,. rhohc=0.788g/cm3(ancient page 1 well region), c takes the value as the attached figure 3, the oil saturation of the gulonic shale is calculated by the following formula,
Figure BDA0002819245930000121
compared with the oil saturation result of field one-dimensional nuclear magnetic experiment analysis, the method has the advantages that the comparison result is integrally well consistent, the calculation result in the lamellar shale section with high vitrinite reflectivity (maturity) is higher than the experiment result value and is possibly higher than the S0The components are volatilized and related, and the oil saturation of the gulong shale is calculated, drawn and displayed by utilizing the data and the macro editing function of Sinololog software.
In order to verify the applicability of the gulong shale oil and gas saturation model, the data of a gulong shale reservoir well in a research area are processed, fig. 4 is a schematic diagram of the processing effect of the data of an actual well, 1 st and 2 nd paths are geological stratification paths, 3 rd path is a lithologic logging curve path which comprises natural gamma, natural potential and well diameter curves, 4 th path is a depth path, 5 th path is deep medium shallow resistivity, 6 th path is a three-porosity logging curve path which comprises compensation density, compensation neutrons and acoustic wave time difference curves, and 7 th path is a logging calculation curvec value, S calculated for well log in lane 81And core analysis S1And comparing the oil saturation calculated for the logging curve and the one-dimensional nuclear magnetic saturation of the field core analysis in the 9 th path, wherein the integral comparison result has good coincidence, and the calculation result of the layered shale section with high vitrinite reflectivity (maturity) in the lower 2360-2376-meter layer section is higher than the experimental result value and possibly has better consistency with the S0The component is volatilized to cause the c value to be higher, and the analysis of the result proves the correctness of the method provided by the invention and the adaptability in the gulonic shale reservoir.
As prior art (1): sun jiangmeng a computational model of shale oil and gas saturation, university of petroleum in china (east china), 2015. Prior art (2): the key parameters and the solving method for Wangming shale oil evaluation are researched, the deposition is reported, 2 months in 2014, No. 32, No. 1, and pages 174-181. Sun's dream proposes that after parameters such as a, b, m, n and the like are determined in a rock-electricity experiment, the oil saturation is calculated according to an improved Archie formula, and the key parameters for Wangming shale oil-containing evaluation are TOC, free hydrocarbon and the like, and do not relate to saturation calculation.
The invention aims to solve the technical problem that an existing oil saturation model is not applicable due to the existence of organic matters, pyrite and other special minerals and complex wettability (organic pore oil is wet and inorganic pore water is wet) in a gulonic shale reservoir, and provides a method for pyrolyzing S in a large number of rock cores1On the basis, a non-electrical method hydrocarbon saturation model is established by combining a conventional logging sensitivity curve. Specifically, the disclosure calculates S from well log data1On the basis of S1The larger the value, the better the oil content of the gulong shale. The innovation of the invention is also the pyrolysis of S by rock1Can indicate the oil content of the gulong shale to a certain extent, namely S1Higher indicates better oil-bearing properties of the formation.
The method can solve the problem that an Archie model of the gulong shale reservoir is not suitable any more due to factors such as high shale content, complex pore structure and double wettability, accurately calculates the oil and gas saturation of the gulong shale reservoir, provides help for well logging and reservoir evaluation and oil field development, has incomparable advantages of other saturation models in the aspect of calculating the oil saturation of the gulong shale reservoir, and has obvious practical application effect, so the method has great popularization value.
It is understood that the above-mentioned method embodiments of the present disclosure can be combined with each other to form a combined embodiment without departing from the logic of the principle, which is limited by the space, and the detailed description of the present disclosure is omitted.
In addition, the present disclosure also provides a device, an electronic device, a computer-readable storage medium, and a program for determining oil saturation, which can be used to implement any method for determining oil saturation provided by the present disclosure, and the corresponding technical solutions and descriptions and corresponding descriptions in the methods section are not repeated.
The main body of the determination method of the oil saturation may be a signal processing apparatus, for example, the determination method of the oil saturation may be performed by a terminal device or a server or other processing device, wherein the terminal device may be a User Equipment (UE), a mobile device, a User terminal, a cellular phone, a cordless phone, a Personal Digital Assistant (PDA), a handheld device, a computing device, a vehicle-mounted device, a wearable device, or the like. In some possible implementations, the method for determining oil saturation may be implemented by a processor calling computer readable instructions stored in a memory.
It will be understood by those skilled in the art that in the method of the present invention, the order of writing the steps does not imply a strict order of execution and any limitations on the implementation, and the specific order of execution of the steps should be determined by their function and possible inherent logic.
The present disclosure provides a determination device for oil saturation, including: the acquiring unit is used for acquiring a preset shale oil-gas saturation model and measurement parameters thereof; and the determining unit is used for determining the oil saturation based on the measurement parameters and the shale oil-gas saturation model.
Simultaneously, this disclosure has still proposed an electronic equipment, includes: a processor; a memory for storing processor-executable instructions; wherein the processor is configured to invoke the memory-stored instructions to perform the above-described determination method.
The present disclosure also provides a computer-readable storage medium having computer program instructions stored thereon, wherein the computer program instructions are executed by a processor to perform the above-mentioned determination method.
In some embodiments, functions of or modules included in the apparatus provided in the embodiments of the present disclosure may be used to execute the method described in the above method embodiments, and specific implementation thereof may refer to the description of the above method embodiments, and for brevity, will not be described again here.
Embodiments of the present disclosure also provide a computer-readable storage medium having stored thereon computer program instructions, which when executed by a processor, implement the above-mentioned method. The computer readable storage medium may be a non-volatile computer readable storage medium.
An embodiment of the present disclosure further provides an electronic device, including: a processor; a memory for storing processor-executable instructions; wherein the processor is configured as the above method. The electronic device may be provided as a terminal, server, or other form of device.
Fig. 5 is a block diagram illustrating an electronic device 800 in accordance with an example embodiment. For example, the electronic device 800 may be a mobile phone, a computer, a digital broadcast terminal, a messaging device, a game console, a tablet device, a medical device, a fitness device, a personal digital assistant, or the like terminal.
Referring to fig. 5, electronic device 800 may include one or more of the following components: processing component 802, memory 804, power component 806, multimedia component 808, audio component 810, input/output (I/O) interface 812, sensor component 814, and communication component 816.
The processing component 802 generally controls overall operation of the electronic device 800, such as operations associated with display, telephone calls, data communications, camera operations, and recording operations. The processing components 802 may include one or more processors 820 to execute instructions to perform all or a portion of the steps of the methods described above. Further, the processing component 802 can include one or more modules that facilitate interaction between the processing component 802 and other components. For example, the processing component 802 can include a multimedia module to facilitate interaction between the multimedia component 808 and the processing component 802.
The memory 804 is configured to store various types of data to support operations at the electronic device 800. Examples of such data include instructions for any application or method operating on the electronic device 800, contact data, phonebook data, messages, pictures, videos, and so forth. The memory 804 may be implemented by any type or combination of volatile or non-volatile memory devices such as Static Random Access Memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic or optical disks.
The power supply component 806 provides power to the various components of the electronic device 800. The power components 806 may include a power management system, one or more power supplies, and other components associated with generating, managing, and distributing power for the electronic device 800.
The multimedia component 808 includes a screen that provides an output interface between the electronic device 800 and a user. In some embodiments, the screen may include a Liquid Crystal Display (LCD) and a Touch Panel (TP). If the screen includes a touch panel, the screen may be implemented as a touch screen to receive an input signal from a user. The touch panel includes one or more touch sensors to sense touch, slide, and gestures on the touch panel. The touch sensor may not only sense the boundary of a touch or slide action, but also detect the duration and pressure associated with the touch or slide operation. In some embodiments, the multimedia component 808 includes a front facing camera and/or a rear facing camera. The front camera and/or the rear camera may receive external multimedia data when the electronic device 800 is in an operation mode, such as a shooting mode or a video mode. Each front camera and rear camera may be a fixed optical lens system or have a focal length and optical zoom capability.
The audio component 810 is configured to output and/or input audio signals. For example, the audio component 810 includes a Microphone (MIC) configured to receive external audio signals when the electronic device 800 is in an operational mode, such as a call mode, a recording mode, and a voice recognition mode. The received audio signals may further be stored in the memory 804 or transmitted via the communication component 816. In some embodiments, audio component 810 also includes a speaker for outputting audio signals.
The I/O interface 812 provides an interface between the processing component 802 and peripheral interface modules, which may be keyboards, click wheels, buttons, etc. These buttons may include, but are not limited to: a home button, a volume button, a start button, and a lock button.
The sensor assembly 814 includes one or more sensors for providing various aspects of state assessment for the electronic device 800. For example, the sensor assembly 814 may detect an open/closed state of the electronic device 800, the relative positioning of components, such as a display and keypad of the electronic device 800, the sensor assembly 814 may also detect a change in the position of the electronic device 800 or a component of the electronic device 800, the presence or absence of user contact with the electronic device 800, orientation or acceleration/deceleration of the electronic device 800, and a change in the temperature of the electronic device 800. Sensor assembly 814 may include a proximity sensor configured to detect the presence of a nearby object without any physical contact. The sensor assembly 814 may also include a light sensor, such as a CMOS or CCD image sensor, for use in imaging applications. In some embodiments, the sensor assembly 814 may also include an acceleration sensor, a gyroscope sensor, a magnetic sensor, a pressure sensor, or a temperature sensor.
The communication component 816 is configured to facilitate wired or wireless communication between the electronic device 800 and other devices. The electronic device 800 may access a wireless network based on a communication standard, such as WiFi, 2G or 3G, or a combination thereof. In an exemplary embodiment, the communication component 816 receives a broadcast signal or broadcast related information from an external broadcast management system via a broadcast channel. In an exemplary embodiment, the communication component 816 further includes a Near Field Communication (NFC) module to facilitate short-range communications. For example, the NFC module may be implemented based on Radio Frequency Identification (RFID) technology, infrared data association (IrDA) technology, Ultra Wideband (UWB) technology, Bluetooth (BT) technology, and other technologies.
In an exemplary embodiment, the electronic device 800 may be implemented by one or more Application Specific Integrated Circuits (ASICs), Digital Signal Processors (DSPs), Digital Signal Processing Devices (DSPDs), Programmable Logic Devices (PLDs), Field Programmable Gate Arrays (FPGAs), controllers, micro-controllers, microprocessors or other electronic components for performing the above-described methods.
In an exemplary embodiment, a non-transitory computer-readable storage medium, such as the memory 804, is also provided that includes computer program instructions executable by the processor 820 of the electronic device 800 to perform the above-described methods.
Fig. 6 is a block diagram illustrating an electronic device 1900 according to an example embodiment. For example, the electronic device 1900 may be provided as a server. Referring to fig. 6, electronic device 1900 includes a processing component 1922 further including one or more processors and memory resources, represented by memory 1932, for storing instructions, e.g., applications, executable by processing component 1922. The application programs stored in memory 1932 may include one or more modules that each correspond to a set of instructions. Further, the processing component 1922 is configured to execute instructions to perform the above-described method.
The electronic device 1900 may also include a power component 1926 configured to perform power management of the electronic device 1900, a wired or wireless network interface 1950 configured to connect the electronic device 1900 to a network, and an input/output (I/O) interface 1958. The electronic device 1900 may operate based on an operating system stored in memory 1932, such as Windows Server, Mac OS XTM, UnixTM, LinuxTM, FreeBSDTM, or the like.
In an exemplary embodiment, a non-transitory computer readable storage medium, such as the memory 1932, is also provided that includes computer program instructions executable by the processing component 1922 of the electronic device 1900 to perform the above-described methods.
The present disclosure may be systems, methods, and/or computer program products. The computer program product may include a computer-readable storage medium having computer-readable program instructions embodied thereon for causing a processor to implement various aspects of the present disclosure.
The computer readable storage medium may be a tangible device that can hold and store the instructions for use by the instruction execution device. The computer readable storage medium may be, for example, but not limited to, an electronic memory device, a magnetic memory device, an optical memory device, an electromagnetic memory device, a semiconductor memory device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following: a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a Static Random Access Memory (SRAM), a portable compact disc read-only memory (CD-ROM), a Digital Versatile Disc (DVD), a memory stick, a floppy disk, a mechanical coding device, such as punch cards or in-groove projection structures having instructions stored thereon, and any suitable combination of the foregoing. Computer-readable storage media as used herein is not to be construed as transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission medium (e.g., optical pulses through a fiber optic cable), or electrical signals transmitted through electrical wires.
The computer-readable program instructions described herein may be downloaded from a computer-readable storage medium to a respective computing/processing device, or to an external computer or external storage device via a network, such as the internet, a local area network, a wide area network, and/or a wireless network. The network may include copper transmission cables, fiber optic transmission, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. The network adapter card or network interface in each computing/processing device receives computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in a computer-readable storage medium in the respective computing/processing device.
The computer program instructions for carrying out operations of the present disclosure may be assembler instructions, Instruction Set Architecture (ISA) instructions, machine-related instructions, microcode, firmware instructions, state setting data, or source or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C + + or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The computer-readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider). In some embodiments, the electronic circuitry that can execute the computer-readable program instructions implements aspects of the present disclosure by utilizing the state information of the computer-readable program instructions to personalize the electronic circuitry, such as a programmable logic circuit, a Field Programmable Gate Array (FPGA), or a Programmable Logic Array (PLA).
Various aspects of the present disclosure are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the disclosure. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer-readable program instructions.
These computer-readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer-readable program instructions may also be stored in a computer-readable storage medium that can direct a computer, programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer-readable medium storing the instructions comprises an article of manufacture including instructions which implement the function/act specified in the flowchart and/or block diagram block or blocks.
The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other devices to cause a series of operational steps to be performed on the computer, other programmable apparatus or other devices to produce a computer implemented process such that the instructions which execute on the computer, other programmable apparatus or other devices implement the functions/acts specified in the flowchart and/or block diagram block or blocks.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical function(s). In some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
Having described embodiments of the present disclosure, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The terminology used herein is chosen in order to best explain the principles of the embodiments, the practical application, or technical improvements to the technology in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.

Claims (10)

1. A method for determining oil saturation, comprising:
acquiring a preset shale oil-gas saturation model and measurement parameters thereof;
and determining the oil saturation degree based on the measurement parameters and the shale oil and gas saturation degree model.
2. The method for determining oil saturation according to claim 1, wherein before the obtaining of the preset shale oil and gas saturation model, the preset shale oil and gas saturation model needs to be established, and the method for establishing the preset shale oil and gas saturation model includes:
obtaining the mass of oil gas in the shale, the volume of the oil gas in the shale and the effective porosity;
and establishing the preset shale oil-gas saturation model according to the quality of the oil gas in the shale, the volume of the oil gas in the shale and the effective porosity.
3. The method of determining oil saturation according to claim 2, wherein the mass of hydrocarbons in the shale and the volume of hydrocarbons in the shale are determined before the obtaining of the mass of hydrocarbons in the shale and the volume of hydrocarbons in the shale;
the method for determining the quality of hydrocarbons in shale comprises the following steps:
acquiring logging information, the quality of shale and scale conversion coefficients;
determining the content volume of hydrocarbon in the shale at different temperatures and/or different temperature intervals according to the logging information;
determining the mass of oil and gas in the shale based on the content volume of hydrocarbon in the shale at different temperatures and/or different temperature intervals, the mass of the shale and a scale conversion coefficient;
and/or the presence of a gas in the interior of the container,
the method for determining the volume of the oil and gas in the shale comprises the following steps:
acquiring the mass and density of oil and gas in the shale;
and determining the volume of the oil and gas in the shale according to the mass of the oil and gas in the shale and the density of the oil and gas.
4. The method for determining oil saturation according to claim 3, wherein the method for obtaining the quality of the shale comprises:
acquiring the density of a rock skeleton, the bulk density of shale, the density of fluid in rock, the mass of the rock skeleton and the mass of fluid;
determining the volume ratio of shale pores according to the density of the rock skeleton, the bulk density of the shale and the density of fluid in the rock;
and determining the mass of the shale according to the shale pore volume ratio, the mass of the rock skeleton and the mass of the fluid.
5. The method for determining oil saturation according to any one of claims 1 to 4, further comprising: determining the content volume ratio of hydrocarbon in the shale at a certain temperature and/or a certain temperature interval;
and determining a shale oil-gas saturation model corresponding to a certain temperature and/or a certain temperature interval by using the content volume ratio of the hydrocarbons in the shale at the certain temperature and/or the certain temperature interval to represent the content volumes of the hydrocarbons in the shale at other temperatures and/or other temperature intervals.
6. The method for determining oil saturation according to any one of claims 1 to 4, wherein the measuring parameters include: a given scale conversion coefficient, a measured content volume and a ratio of hydrocarbons in the shale at a certain temperature and/or a certain temperature interval, a measured bulk density of the shale, a measured oil and gas density and a measured effective porosity;
and/or the presence of a gas in the interior of the container,
the method for measuring the content volume of the hydrocarbon in the shale at a certain temperature and/or a certain temperature interval comprises the following steps:
acquiring the relation of the logging sensitivity curves of the hydrocarbon content volume, the deep lateral resistivity, the compensation density, the compensation sound wave and the compensation neutron in the shale at a certain temperature and/or a certain temperature interval;
and obtaining the content volume of the hydrocarbons in the shale at a certain temperature and/or a certain temperature interval based on the current deep lateral resistivity, the compensation density, the compensation sound wave and the compensation neutrons and the relation of the logging sensitivity curve.
7. The method for determining oil saturation according to claim 6, characterized by:
before obtaining a log sensitivity curve relation of hydrocarbon content volume, deep lateral resistivity, compensation density, compensation sound wave and compensation neutrons in shale at a certain temperature and/or a certain temperature interval, the method needs to establish the log sensitivity curve relation, and comprises the following steps:
obtaining a volume of hydrocarbon content in shale at a temperature and/or temperature interval for a plurality of samples;
and establishing the well logging sensitivity curve relation according to the content volume of the hydrocarbons in the shale of the multiple samples at a certain temperature and/or a certain temperature interval and the corresponding deep lateral resistivity, compensation density, compensation sound wave and compensation neutrons of the hydrocarbons.
8. An apparatus for determining oil saturation, comprising:
the acquiring unit is used for acquiring a preset shale oil-gas saturation model and measurement parameters thereof;
and the determining unit is used for determining the oil saturation based on the measurement parameters and the shale oil-gas saturation model.
9. An electronic device, comprising:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to invoke the memory-stored instructions to perform the determination method of any one of claims 1 to 7.
10. A computer-readable storage medium having computer program instructions stored thereon, which, when executed by a processor, implement the determination method of any one of claims 1 to 7.
CN202011418520.6A 2020-12-06 2020-12-06 Method and device for determining oil saturation, electronic device and storage medium Pending CN112382346A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011418520.6A CN112382346A (en) 2020-12-06 2020-12-06 Method and device for determining oil saturation, electronic device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011418520.6A CN112382346A (en) 2020-12-06 2020-12-06 Method and device for determining oil saturation, electronic device and storage medium

Publications (1)

Publication Number Publication Date
CN112382346A true CN112382346A (en) 2021-02-19

Family

ID=74590655

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011418520.6A Pending CN112382346A (en) 2020-12-06 2020-12-06 Method and device for determining oil saturation, electronic device and storage medium

Country Status (1)

Country Link
CN (1) CN112382346A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113379209A (en) * 2021-05-31 2021-09-10 中石化经纬有限公司胜利测井公司 Shale oil saturation evaluation method based on geochemical pyrolysis correction
CN113781452A (en) * 2021-09-15 2021-12-10 大庆油田有限责任公司 Method and device for determining shale oil saturation, electronic equipment and storage medium

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113379209A (en) * 2021-05-31 2021-09-10 中石化经纬有限公司胜利测井公司 Shale oil saturation evaluation method based on geochemical pyrolysis correction
CN113781452A (en) * 2021-09-15 2021-12-10 大庆油田有限责任公司 Method and device for determining shale oil saturation, electronic equipment and storage medium
CN113781452B (en) * 2021-09-15 2022-09-27 大庆油田有限责任公司 Method and device for determining shale oil saturation, electronic equipment and storage medium

Similar Documents

Publication Publication Date Title
CN112382346A (en) Method and device for determining oil saturation, electronic device and storage medium
CN108708695B (en) Method and device for determining dominant seepage channel and storage medium
CN110837114B (en) Method and device for identifying rough-surface pyroclastic rock and electronic equipment
CN105372166A (en) Method and device for obtaining permeability of argillaceous sandstone
Yan et al. Prediction of seismic wave dispersion and attenuation from ultrasonic velocity measurements
Al-Ghamdi et al. Integrated high-resolution chemostratigraphy and facies-based stratigraphic architecture of the Lower Cretaceous (Aptian), Shu’aiba Formation, Saudi Arabia
Khodabakhshnejad et al. A fast method to forecast shale pressure depletion and well performance using geomechanical constraints-application to poro-elasticity modeling to predict mid and far field frac hits at an Eagle Ford and Wolfcamp Well
Akai et al. Pressure dependent permeability of tight rocks
CN115495865A (en) Clay mineral content determination method and device, electronic equipment and storage medium
CN116796503A (en) Shale organic carbon content calculation model establishment method and device
CN113781452B (en) Method and device for determining shale oil saturation, electronic equipment and storage medium
CN115935816A (en) Drilling parameter determination method, device, equipment and storage medium
Kukal et al. Critical problems hindering accurate log interpretation of tight gas sand reservoirs
CN112211624A (en) Reservoir fracture pressure determination method and device, electronic device and storage medium
Hursan et al. Slimhole NMR T1 logging while drilling enhances real-time petrophysics
CN111949924B (en) Determination method and apparatus for permeability, electronic device and storage medium
CN113123784B (en) Method and device for acquiring fracture permeability of low-hole fractured sandstone reservoir
CN112883606B (en) Permeability determination method, characterization method and device, electronic equipment and storage medium
CN115496311A (en) Method and device for determining favorable area of tight oil reservoir, electronic equipment and medium
CN115704309A (en) Method and device for determining productivity of complex conglomerate reservoir and electronic equipment
CN112922592B (en) Method and device for predicting flow channel and speed, electronic equipment and storage medium
CN113406697B (en) Method and device for acquiring transverse wave quality factor of reservoir
CN117852913A (en) Method and device for predicting daily gas production of single well of tight sandstone gas reservoir
CN111814879A (en) Method and device for classifying physical phases of reservoir rock, electronic equipment and storage medium
CN113791103B (en) Shale oil-containing property determination method and device, electronic equipment and storage medium

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