CN113495046B

CN113495046B - Method and device for determining reservoir type and storage medium

Info

Publication number: CN113495046B
Application number: CN202010258280.1A
Authority: CN
Inventors: 王俊杰; 胡勇; 张连进; 兰雪梅; 何溥为; 杨东凡; 徐昌海
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2020-04-03
Filing date: 2020-04-03
Publication date: 2024-06-25
Anticipated expiration: 2040-04-03
Also published as: CN113495046A

Abstract

The application discloses a method and a device for determining reservoir types and a storage medium, and belongs to the technical field of oil and gas reservoir exploration. The method comprises the following steps: acquiring the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal and nuclear magnetic conversion coefficient of a core of a target sample; based on the length, radius, nuclear magnetic T2 spectrum and total nuclear magnetic signal of the target core, respectively determining the nuclear magnetic porosity and nuclear magnetic permeability of the target sample core; determining a plurality of pore structure parameters of the target sample core based on the nuclear magnetic T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core; and determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the pore structure parameters of the target sample core. According to the method, the reservoir type corresponding to the target sample core is determined, so that technical guidance is provided for subsequent reservoir evaluation and development.

Description

Method and device for determining reservoir type and storage medium

Technical Field

The application relates to the technical field of oil and gas reservoir exploration, in particular to a method and a device for determining reservoir types and a storage medium.

Background

Because the pores, holes and slots of carbonate reservoirs are important sites for oil and gas enrichment, carbonate reservoirs are very closely related to oil and gas. Most of carbonate reservoirs are formed in multi-rotation superimposed basins, and because the carbonate reservoirs have the characteristic of easy change, the carbonate reservoirs are more easily subjected to strong diagenetic transformation, so that the types of carbonate reservoirs in the carbonate reservoirs are more and the heterogeneity is strong. The type of the carbonate reservoir is a key for evaluating the carbonate reservoir, and is a basis for guiding subsequent development. Thus, the reservoir type of the carbonate reservoir may be predetermined prior to the development of the carbonate reservoir.

Disclosure of Invention

The application provides a method and a device for determining a reservoir type and a storage medium, which can solve the problem of reservoir type determination of a carbonate reservoir. The technical scheme is as follows:

in a first aspect, there is provided a method of determining a reservoir type, the method comprising:

Acquiring the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal and nuclear magnetic conversion coefficient of a core of a target sample;

Based on the length, radius, nuclear magnetic T2 spectrum and total nuclear magnetic signal of the target core, respectively determining the nuclear magnetic porosity and nuclear magnetic permeability of the target sample core;

determining a plurality of pore structure parameters of the target sample core based on a nuclear magnetic T2 spectrum and a nuclear magnetic conversion coefficient of the target sample core;

and determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the plurality of pore structure parameters of the target sample core.

Optionally, the determining the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core based on the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target core includes:

determining the nuclear magnetic porosity of the target sample core based on the length, radius and nuclear magnetic signal total amount of the target sample core;

determining the geometric relaxation time of the target sample core based on a plurality of relaxation times and a plurality of nuclear magnetic signal amplitude values which are in one-to-one correspondence on a nuclear magnetic T2 spectrum of the target sample core;

And determining the nuclear magnetic permeability of the target sample core based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core.

Optionally, the determining the nuclear magnetic porosity of the target sample core according to the first conversion formula based on the length, the radius and the total nuclear magnetic signal of the target sample core includes:

Determining the nuclear magnetic porosity of the target sample core through a first conversion formula based on the length, the radius and the total nuclear magnetic signal amount of the target sample core;

the first conversion formula:

Wherein, in the first conversion formula, Φ _nmr refers to the nuclear magnetic porosity of the sample core; m refers to the total nuclear magnetic signal of the sample core; l is the length of the sample core; r is the radius of the sample core; a refers to a first parameter; b refers to a second parameter;

The determining the nuclear magnetic permeability of the target sample core based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core comprises the following steps:

Determining the nuclear magnetic permeability of the target sample core through a second conversion formula based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core;

The second conversion formula: k _nmr＝E·Φ_nmr ⁴·T_2g ²

Wherein, in the first conversion formula and the second conversion formula, K _nmr refers to the nuclear magnetic permeability of the sample core; e refers to the average osmotic conversion coefficient of the sample core; phi _nmr refers to the nuclear magnetic porosity of the sample core; t _2g refers to the geometric relaxation time of the sample core.

Optionally, before determining the nuclear magnetic porosity of the target sample core according to the first conversion formula based on the length, the radius and the total nuclear magnetic signal of the target sample core, the method further includes:

acquiring the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal and theoretical porosity of each sample core in a plurality of sample cores;

Determining nuclear magnetic signal densities of the plurality of sample cores based on the lengths, the radii and the total nuclear magnetic signal amounts of the plurality of sample cores, respectively;

And performing linear fitting on the nuclear magnetic signal densities and the theoretical porosities of the plurality of sample cores to obtain the first conversion formula.

Optionally, before determining the nuclear magnetic permeability of the target sample core according to the second conversion formula based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core, the method further includes:

obtaining theoretical permeability of the plurality of sample cores;

Based on the lengths, the radiuses and the total nuclear magnetic signal amounts of the plurality of sample cores, respectively determining the nuclear magnetic porosity of each sample core through the first conversion formula;

Based on the nuclear magnetic T2 spectra of the plurality of sample cores, respectively determining the geometric relaxation time of each sample core;

Determining an average permeability conversion coefficient based on the core magnetic porosity, the geometric relaxation time, and the theoretical permeability of each of the plurality of sample cores;

the second conversion formula is generated based on the average osmotic conversion coefficient, the nuclear magnetic porosity, and the geometric relaxation time.

Optionally, the determining a plurality of pore structure parameters of the target sample core based on the nuclear magnetic T2 spectrum and the nuclear magnetic transformation coefficient of the target sample core includes:

Converting the nuclear magnetic T2 spectrum of the target sample core based on a nuclear magnetic conversion coefficient to obtain a nuclear magnetic pore size distribution curve of the target sample core;

determining a plurality of nuclear magnetic radii corresponding to a plurality of accumulation frequencies one by one based on the distribution frequency of the nuclear magnetic radii on the nuclear magnetic pore diameter distribution curve of the target sample core;

respectively performing geometric transformation on the nuclear magnetic radii to obtain a plurality of geometric series in one-to-one correspondence;

A plurality of pore structure parameters of the target sample core are determined based on the plurality of geometric progression.

Optionally, the determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the plurality of pore structure parameters of the target sample core includes:

Determining a characteristic constant of the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and a plurality of pore structure parameters of the target sample core;

determining a characteristic range in which a characteristic constant of the target sample core is located from a plurality of characteristic ranges;

And determining the reservoir type corresponding to the target sample core from the corresponding relation between the pre-stored characteristic range and the reservoir type.

Optionally, the plurality of pore structure parameters includes a sorting coefficient, kurtosis, skewness, radius average, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient, and relative sorting coefficient of the target sample core;

The determining the characteristic constant of the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the plurality of pore structure parameters of the target sample core comprises:

Determining a characteristic constant of the target sample core through a type discrimination formula based on the nuclear magnetic porosity, the nuclear magnetic permeability and a plurality of pore structure parameters of the target sample core;

The type discrimination formula:

Wherein, in the type discrimination formula, T refers to a characteristic constant of the sample core; s _p is the sorting coefficient of the sample core; k _p refers to the kurtosis of the sample core; s _kp refers to the skewness of the sample core; d _M refers to the radius average of the sample core; r ₅₀ is the median nuclear magnetic radius of the sample core; mean core magnetic radius of the sample core; r ₁₀₀ is the maximum nuclear magnetic radius of the sample core; alpha refers to the homogeneity coefficient of the sample core; d refers to the relative sorting coefficient of the sample core; k _nmr refers to the nuclear magnetic permeability of the sample core.

In a second aspect, there is provided an apparatus for determining a reservoir type, the apparatus comprising:

the acquisition module is used for acquiring the length, the radius, the nuclear magnetic T2 spectrum, the total nuclear magnetic signal amount and the nuclear magnetic conversion coefficient of the core of the target sample;

the first determining module is used for respectively determining the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core based on the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target core;

The second determining module is used for determining a plurality of pore structure parameters of the target sample core based on the nuclear magnetic T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core;

and the third determining module is used for determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the pore structure parameters of the target sample core.

Optionally, the first determining module includes:

A sixth determining unit, configured to determine a nuclear magnetic porosity of the target sample core based on a length, a radius, and a total nuclear magnetic signal of the target sample core;

A seventh determining unit, configured to determine a geometric relaxation time of the target sample core based on a plurality of relaxation times on a nuclear magnetic T2 spectrum of the target sample core and a plurality of nuclear magnetic signal amplitude values that are in one-to-one correspondence;

and an eighth determining unit, configured to determine a nuclear magnetic permeability of the target sample core based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core.

Optionally, the sixth determining unit is mainly configured to:

the first conversion formula:

The eighth determination unit is mainly configured to:

The second conversion formula: k _nmr＝E·Φ_nmr ⁴·T_2g ²

Optionally, the sixth determining unit is further configured to:

Optionally, the eighth determining unit is further configured to:

obtaining theoretical permeability of the plurality of sample cores;

Optionally, the second determining module includes:

The first conversion unit is used for converting the nuclear magnetic T2 spectrum of the target sample core based on the nuclear magnetic conversion coefficient to obtain a nuclear magnetic pore diameter distribution curve of the target sample core;

the first determining unit is used for determining a plurality of nuclear magnetic radii corresponding to a plurality of accumulated frequencies one by one based on the distribution frequency of the nuclear magnetic radii on the nuclear magnetic pore diameter distribution curve of the target sample core;

the second conversion unit is used for respectively carrying out geometric conversion on the plurality of nuclear magnetic radii to obtain a plurality of geometric series in one-to-one correspondence;

And the second determining unit is used for determining a plurality of pore structure parameters of the target sample core based on the geometric series.

Optionally, the third determining module includes:

A third determining unit, configured to determine a characteristic constant of the target sample core based on a nuclear magnetic porosity, a nuclear magnetic permeability, and a plurality of pore structure parameters of the target sample core;

a fourth determining unit, configured to determine a feature range in which a feature constant of the target sample core is located from a plurality of feature ranges;

and a fifth determining unit, configured to determine a reservoir type corresponding to the target sample core from a pre-stored correspondence between a feature range and a reservoir type.

the third determining unit is mainly configured to:

The type discrimination formula:

In a third aspect, there is provided a computer readable storage medium having stored therein a computer program which, when executed by a processor, implements the method of any of the first aspects provided above.

In a fourth aspect, there is provided a computer program product containing instructions which, when run on a computer, cause the computer to perform any of the methods provided in the first aspect.

The technical scheme provided by the application has the beneficial effects that at least the following steps are included:

And determining the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core through the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target sample core, and determining a plurality of pore structure parameters of the target sample core through the nuclear magnetic T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core. And then determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the pore structure parameters of the target sample core, so that technical guidance is provided for the evaluation of subsequent reservoirs, and the development of the subsequent reservoirs is facilitated.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for the description of the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present application, and 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 a method for determining a reservoir type according to an embodiment of the present application;

FIG. 2 is a flow chart of another method for determining reservoir type according to an embodiment of the present application;

FIG. 3 is a schematic structural diagram of a reservoir type determining device according to an embodiment of the present application;

fig. 4 is a block diagram of a terminal according to an embodiment of the present application.

Detailed Description

For the purpose of making the objects, technical solutions and advantages of the present application more apparent, the embodiments of the present application will be described in further detail with reference to the accompanying drawings.

Fig. 1 is a flow chart of a method for determining a reservoir type according to an embodiment of the present application. Referring to fig. 1, the method includes the following steps.

Step 101: and acquiring the length, radius, nuclear magnetic T2 spectrum, nuclear magnetic signal total quantity and nuclear magnetic conversion coefficient of the core of the target sample.

Step 102: and respectively determining the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core based on the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target core.

Step 103: and determining a plurality of pore structure parameters of the target sample core based on the nuclear magnetic T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core.

Step 104: and determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the pore structure parameters of the target sample core.

In the embodiment of the application, the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core are determined according to the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target sample core, and a plurality of pore structure parameters of the target sample core are determined according to the greetings T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core. And then determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the pore structure parameters of the target sample core, so that technical guidance is provided for the evaluation of subsequent reservoirs, and the development of the subsequent reservoirs is facilitated.

Optionally, determining a plurality of pore structure parameters of the target sample core based on the nuclear magnetic T2 spectrum and the nuclear magnetic transformation coefficient of the target sample core includes:

converting the nuclear magnetic T2 spectrum of the target sample core based on the nuclear magnetic conversion coefficient to obtain a nuclear magnetic pore size distribution curve of the target sample core;

determining a plurality of nuclear magnetic radii corresponding to the plurality of accumulation frequencies one by one based on the distribution frequency of the nuclear magnetic radii on the nuclear magnetic pore diameter distribution curve of the target sample core;

Respectively performing geometric transformation on the plurality of nuclear magnetic radii to obtain a plurality of geometric series in one-to-one correspondence;

Optionally, determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability, and the plurality of pore structure parameters of the target sample core includes:

Determining a characteristic constant of the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the plurality of pore structure parameters of the target sample core;

Determining a characteristic constant of the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability, and the plurality of pore structure parameters of the target sample core, comprising:

Type discrimination formula:

Optionally, determining the core magnetic porosity and the core magnetic permeability of the target sample core based on the length, the radius, the core magnetic T2 spectrum and the total amount of the core magnetic signals, respectively, includes:

determining the nuclear magnetic porosity of the target sample core based on the length, the radius and the total nuclear magnetic signal amount of the target sample core;

Determining the geometric relaxation time of the target sample core based on a plurality of relaxation times on a nuclear magnetic T2 spectrum of the target sample core and a plurality of nuclear magnetic signal amplitude values in one-to-one correspondence;

Optionally, determining the nuclear magnetic porosity of the target sample core by the first conversion formula based on the length, the radius, and the total nuclear magnetic signal of the target sample core includes:

A first conversion formula:

determining the nuclear magnetic permeability of the target sample core based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core, comprising:

A second conversion formula: k _nmr＝E·Φ_nmr ⁴·T_2g ²

Optionally, before determining the nuclear magnetic porosity of the target sample core according to the first conversion formula based on the length, the radius and the total nuclear magnetic signal of the target sample core, the method further comprises:

Respectively determining the nuclear magnetic signal densities of the plurality of sample cores based on the lengths, the radiuses and the total nuclear magnetic signal amounts of the plurality of sample cores;

And performing linear fitting on the nuclear magnetic signal densities and the theoretical porosities of the plurality of sample cores to obtain a first conversion formula.

obtaining theoretical permeability of a plurality of sample cores;

Based on the length, the radius and the total nuclear magnetic signal amount of the plurality of sample cores, respectively determining the nuclear magnetic porosity of each sample core through a first conversion formula;

Based on the nuclear magnetic T2 spectra of a plurality of sample cores, respectively determining the geometric relaxation time of each sample core;

determining an average permeability conversion coefficient based on the nuclear magnetic porosity, the geometric relaxation time, and the theoretical permeability of each of the plurality of sample cores;

a second conversion formula is generated based on the average osmotic conversion coefficient, the nuclear magnetic porosity, and the geometric relaxation time.

All the above optional technical solutions may be combined according to any choice to form an optional embodiment of the present application, and the embodiments of the present application will not be described in detail.

Fig. 2 is a flow chart of a method for determining a reservoir type according to an embodiment of the present application. The method is applied to a reservoir type determining device which is arranged in a terminal, wherein the terminal can be a computer, a mobile phone, a palm computer, a tablet personal computer and the like. Referring to fig. 2, the method includes the following steps.

In order to ensure accuracy of reservoir type determination, a corresponding sample core is generally obtained from a reservoir, and then the type of the reservoir corresponding to the sample core is determined through analysis of the sample core. In the embodiment of the present application, the reservoir type of the reservoir corresponding to the target sample core may be determined through the following steps 201-207.

Step 201: and acquiring the length, radius, nuclear magnetic T2 spectrum, nuclear magnetic signal total quantity and nuclear magnetic conversion coefficient of the core of the target sample.

For the target sample core, the terminal can display a first parameter acquisition interface to acquire the length, the radius, the nuclear magnetic T2 spectrum, the total nuclear magnetic signal and the nuclear magnetic conversion coefficient of the target sample core input by a user at the first parameter acquisition interface. That is, the user may input the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal, and nuclear magnetic conversion coefficient of the target sample core at the first parameter acquisition interface, so that the terminal may acquire these parameters from the first parameter acquisition interface. Of course, the terminal may also communicate with the data storage device to obtain the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal and nuclear magnetic conversion coefficient of the target sample core from the storage device.

Wherein, the abscissa of the nuclear magnetic T2 spectrum can be relaxation time, and the ordinate can be nuclear magnetic signal amplitude value.

It should be noted that, before the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal and nuclear magnetic conversion coefficient of the target sample core are obtained, parameters of the target sample core may be tested by various test methods, respectively.

Illustratively, the length and radius of the target sample core may be tested using a vernier caliper; the first sample core can be placed in the center of a nuclear magnetic resonance coil of a nuclear magnetic resonance device, and the nuclear magnetic T2 spectrum and the total nuclear magnetic signal amount of the first sample core are tested through the nuclear magnetic resonance device.

In some embodiments, a theoretical pore size distribution curve and a nuclear magnetic T2 spectrum of a plurality of sample cores may be obtained, and then based on the theoretical pore size distribution curve and the nuclear magnetic T2 spectrum of each sample core, the nuclear magnetic conversion coefficient of each sample core is determined by a maximum error method, and then an average value of the nuclear magnetic conversion coefficients of the plurality of sample cores is used as the nuclear magnetic conversion coefficient of the target sample core.

The determination of the nuclear magnetic conversion coefficient of each sample core by the maximum error method may refer to related technologies, and the embodiments of the present application are not described herein.

Step 202: and respectively determining the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core based on the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target core.

In some embodiments, the core magnetic porosity and core magnetic permeability of the target sample core may be determined based on the length, radius, core magnetic T2 spectrum, and core magnetic signal total amount, respectively, according to the following steps (1) - (3).

(1) And determining the nuclear magnetic porosity of the target sample core based on the length, the radius and the total nuclear magnetic signal amount of the target sample core.

In some embodiments, the nuclear magnetic porosity of the target sample core may be determined by a first transformation formula based on the length, radius, and nuclear magnetic signal total amount of the target sample core.

Specifically, the length, the radius and the total nuclear magnetic signal amount of the target sample core can be brought into a first conversion formula, and the nuclear magnetic void fraction of the target sample core can be obtained.

Wherein, the first conversion formula:

In the first conversion formula, Φ _nmr refers to the nuclear magnetic porosity of the sample core; m refers to the total nuclear magnetic signal of the sample core; l is the length of the sample core; r is the radius of the sample core; pi refers to the circumference ratio, taking the constant 3.14; a refers to a first parameter; b refers to a second parameter. The first parameter a may be the slope of the fitted line and the second parameter b may be the intercept of the fitted line.

(2) And determining the geometric relaxation time of the target sample core based on the relaxation times on the nuclear magnetic T2 spectrum of the target sample core and the nuclear magnetic signal amplitude values in one-to-one correspondence.

It should be noted that, the method for determining the geometric relaxation time of the target sample core and the method for determining the geometric relaxation time of the first sample core are the same or similar, which is not described in detail in the embodiment of the present application.

(3) And determining the nuclear magnetic permeability of the target sample core based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core.

In some embodiments, the nuclear magnetic permeability of the target sample core may be determined by a second conversion formula based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core.

Specifically, the nuclear magnetic porosity and the geometric relaxation time of the target sample core can be brought into a second conversion formula, and the nuclear magnetic permeability of the target sample core can be obtained.

Wherein, the second conversion formula: k _nmr＝E·Φ_nmr ⁴·T_2g ²

In the second conversion formula, K _nmr refers to the nuclear magnetic permeability of the sample core; e refers to the average osmotic conversion coefficient of the sample core; phi _nmr refers to the nuclear magnetic porosity of the sample core; t _2g refers to the geometric relaxation time of the sample core.

It should be noted that, before the core magnetic porosity of the target sample core is determined by the first conversion formula in the step (1), in some embodiments, the first conversion formula and the second conversion formula may be determined by the following steps (1) - (7).

(1) And acquiring the length, the radius, the nuclear magnetic T2 spectrum, the total nuclear magnetic signal, the theoretical porosity and the theoretical permeability of each of the plurality of sample cores.

The plurality of sample cores may be all sample cores obtained from a reservoir, or may be selected from all sample cores obtained. The screening may be based on the theoretical permeability of each sample core, although other means may be used.

For a first sample core of the plurality of sample cores, the implementation manner of acquiring the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal, theoretical porosity and theoretical permeability of the first sample core may be the same as or similar to the implementation manner of acquiring the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal and nuclear magnetic conversion coefficient of the target sample core described above, which is not limited in the embodiment of the present application.

The first sample core refers to any one of a plurality of sample cores.

Illustratively, the length, radius, total nuclear magnetic signal, theoretical porosity, and theoretical permeability of the 6 sample cores obtained may be as shown in table 1 below.

TABLE 1

It should be noted that, before the length, radius, nuclear magnetic T2 spectrum, total nuclear magnetic signal, theoretical porosity and theoretical permeability of the first sample core are obtained, the parameters of the first sample core may be tested by various test methods, respectively.

Illustratively, the length and radius of the first sample core may be tested using a vernier caliper; the theoretical porosity of the first sample core can be tested by a GRI test method, the theoretical permeability of the first sample core can be tested by a pulse attenuation method, the first sample core can be placed at the center of a nuclear magnetic resonance coil of nuclear magnetic resonance equipment, and the total nuclear magnetic T2 spectrum and nuclear magnetic signal of the first sample core can be tested by the nuclear magnetic resonance equipment.

(2) And determining the nuclear magnetic signal densities of the plurality of sample cores respectively based on the lengths, the radii and the total nuclear magnetic signal amounts of the plurality of sample cores.

For a first sample core of the plurality of sample cores, a volume of the first sample core may be determined based on a length and a radius of the first sample core, and then a ratio between a total amount of nuclear magnetic signals of the first sample core and the volume may be determined as a nuclear magnetic signal density of the first sample core. That is, the nuclear magnetic signal density of the first sample core may be determined according to the following formula (1) based on the length, radius, and total nuclear magnetic signal amount of the first sample core.

Wherein, in the above formula (1),The nuclear magnetic signal density of the first sample core; m _i refers to the total amount of nuclear magnetic signals of the first sample core; l _i is the length of the first sample core; r _i refers to the radius of the first sample core; pi refers to the circumference ratio, taking the constant 3.14.

Since the determination manner of the nuclear magnetic signal density of each sample core in the plurality of sample cores may be the same, the determination of the nuclear magnetic signal density of other sample cores in the plurality of sample cores except the first sample core is not described herein in detail.

Continuing with the above example, the nuclear magnetic signal densities of the 6 sample cores determined by the above equation (1) were 175.5061, 207.1986, 259.0871, 259.172, 325.4913, 326.8172, respectively, in order.

(3) And performing linear fitting on the nuclear magnetic signal densities and the theoretical porosities of the plurality of sample cores to obtain a first conversion formula.

In some embodiments, the nuclear magnetic signal density and the theoretical porosity of the plurality of sample cores may be linearly fitted by a least square method with the nuclear magnetic signal density as a dependent variable and the theoretical porosity as an independent variable, and an expression of the fitted straight line may be a first conversion formula:

continuing with the above example, after straight line fitting is performed on the nuclear magnetic signal densities and the theoretical porosities of the 6 sample cores, the expression of the obtained fitting straight line may be: That is, the first parameter is 31.04 and the second parameter is 100.57.

Of course, the line fitting may be performed by other manners besides the least square method, which is not limited in the embodiment of the present application.

(4) And determining the nuclear magnetic porosity of each sample core respectively through a first conversion formula based on the lengths, the radiuses and the total nuclear magnetic signal amount of the plurality of sample cores.

After the first conversion formula is determined in the step (3), for the first sample cores in the plurality of sample cores, the total nuclear magnetic signal of the first sample core may be brought into the first conversion formula, so as to obtain the nuclear magnetic porosity of the first sample core. That is, the product of the nuclear magnetic signal density of the first sample core and the first parameter may be determined first, and then the sum of the product and the second parameter may be determined as the nuclear magnetic porosity of the first sample core.

Since the determination manner of the core magnetic porosity of each of the plurality of sample cores may be the same, the determination of the core magnetic porosity of the other sample cores of the plurality of sample cores except the first sample core is not described herein.

Continuing with the above example, the determined core magnetic porosities of the 6 sample cores were 3.44, 5.11, 2.41, 7.29, 7.25 in that order.

(5) And respectively determining the geometric relaxation time of each sample core based on the nuclear magnetic T2 spectrums of the plurality of sample cores.

For a first sample core of the plurality of sample cores, in some embodiments, a plurality of relaxation times may be first determined from a nuclear magnetic T2 spectrum of the first sample core, and a plurality of nuclear magnetic signal amplitude values corresponding one-to-one to the plurality of relaxation times, followed by determining a geometric relaxation time of the first sample core based on the plurality of relaxation times and the plurality of nuclear magnetic signal amplitude values corresponding one-to-one.

For example, the geometric relaxation time of the first sample core may be determined according to the following equation (2) based on a plurality of relaxation times and a plurality of nuclear magnetic signal amplitude values that are in one-to-one correspondence.

Wherein, in the above formula (2), T _2gi refers to the geometric relaxation time of the first sample core; t _1i refers to the first relaxation time of the first sample core; t _2i refers to the second relaxation time of the first sample core; t _100i refers to the 100 th relaxation time of the first sample core; n _1i is the nuclear magnetic signal amplitude value corresponding to the 1 st relaxation time of the first sample core; n _2i is the nuclear magnetic signal amplitude value corresponding to the 2 nd relaxation time of the first sample core; n _100i refers to the magnitude of the nuclear magnetic signal corresponding to the 100 th relaxation time of the first sample core.

It should be noted that, since the abscissa of the nuclear magnetic T2 spectrum of the first sample core may include innumerable relaxation times, in order to simplify the operation, a plurality of relaxation times may be selected from the innumerable relaxation times according to a certain step, that is, the plurality of relaxation times is selected on the abscissa of the nuclear magnetic T2 spectrum of the first sample core according to a certain step.

Since the determination manner of the geometric relaxation time of each sample core in the plurality of sample cores may be the same, the determination process of the geometric relaxation time of other sample cores in the plurality of sample cores except the first sample core is not repeated herein.

Continuing with the above example, the geometric relaxation times of the 6 sample cores were 31.30, 35.14, 37.86, 23.68, 41.39, 37.27 in that order.

(6) An average permeability conversion coefficient is determined based on the core magnetic porosity, the geometric relaxation time, and the theoretical porosity of each of the plurality of sample cores.

For a first sample core of the plurality of sample cores, in some embodiments, the permeability conversion coefficient of the first sample core may be determined according to the following equation (3) based on the nuclear magnetic porosity, the theoretical permeability, and the geometric relaxation time of the first sample core.

Wherein, in the above formula (3), E _i refers to the osmotic conversion coefficient of the first sample core; phi _nmri refers to the nuclear magnetic porosity of the first sample core; t _2gi refers to the geometric relaxation time of the first sample core; k _gi refers to the theoretical permeability of the first sample core.

Since the determination manner of the osmotic conversion coefficient of each sample core in the plurality of sample cores may be the same, the determination process of the osmotic conversion coefficient of the other sample cores in the plurality of sample cores except the first sample core is not described herein.

And then, averaging the osmotic conversion coefficients of the plurality of sample cores to obtain the average osmotic conversion coefficient of the plurality of sample cores.

Continuing with the above example, the osmotic conversion coefficients of the 6 sample cores were 36.49, 68.83, 44.44, 73.51, 41.16, 61.63,6 sample cores in that order, and the average osmotic conversion coefficient was 54.23.

(7) A second transformation formula is generated based on the average osmotic transformation coefficient, the nuclear magnetic porosity, and the geometric relaxation time.

Specifically, the second conversion formula K _nmr＝E·Φ_nmr ⁴·T_2g ² is obtained by taking the nuclear magnetic porosity as a parameter, taking the nuclear magnetic porosity and the geometric relaxation time as a dependent variable, taking the average osmotic conversion coefficient as a slope, and taking 0 as an intercept.

Continuing with the above example, since the average osmotic conversion coefficient of the 6 sample cores is 54.23, the second conversion formula may be: k _nmr＝57.23·Φ_nmr ⁴·T_2g ².

And determining a plurality of pore structure parameters of the target sample core based on the nuclear magnetic T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core. In some embodiments, a plurality of pore structure parameters of the target sample core may be determined as follows steps 203-206.

Step 203: and converting the nuclear magnetic T2 spectrum of the target sample core based on the nuclear magnetic conversion coefficient to obtain a nuclear magnetic pore size distribution curve of the target sample core.

Firstly, converting a plurality of relaxation times included in a nuclear magnetic T2 spectrum of a target sample core according to a nuclear magnetic conversion coefficient to obtain a plurality of nuclear magnetic radii of the target sample core, wherein the nuclear magnetic radii correspond to the relaxation times one by one.

It should be noted that, since the abscissa of the nuclear magnetic T2 spectrum of the target pattern core may include innumerable relaxation times, in order to simplify the operation, a plurality of relaxation times may be selected from the innumerable relaxation times, and then, for the plurality of relaxation times, a ratio between each relaxation time and the nuclear magnetic transformation coefficient is determined as a corresponding one of the nuclear magnetic radii, thereby obtaining a plurality of nuclear magnetic radii.

And secondly, converting a plurality of nuclear magnetic signal amplitude values included in the nuclear magnetic T2 spectrum of the target sample core to respectively obtain a plurality of distribution frequencies corresponding to the target sample core.

In some embodiments, the nuclear magnetic T2 spectrum of the target sample core may be integrated to determine the area enclosed by the nuclear magnetic T2 spectrum, that is, determine the total amplitude value of the nuclear magnetic signal of the nuclear magnetic T2 spectrum. Then, since the ordinate of the nuclear magnetic T2 spectrum of the target pattern core may include innumerable nuclear magnetic signal amplitude values, in order to simplify the operation, a plurality of nuclear magnetic signal amplitude values corresponding one to one may be selected from the innumerable nuclear magnetic signal amplitude values based on a plurality of relaxation times. And then, determining the ratio of the amplitude value of each nuclear magnetic signal to the total amplitude value of the nuclear magnetic signal as a corresponding distribution frequency, and further obtaining a plurality of distribution frequencies.

And finally, generating a nuclear magnetic pore size distribution curve of the target sample core based on a plurality of nuclear magnetic radii and a plurality of distribution frequencies of the target sample core.

Step 204: and determining a plurality of nuclear magnetic radii corresponding to the plurality of accumulation frequencies one by one based on the distribution frequency of the nuclear magnetic radii on the nuclear magnetic pore diameter distribution curve of the target sample core.

In some embodiments, for each nuclear magnetic radius included on the nuclear magnetic pore size distribution curve of the target sample core, a section of curve on the nuclear magnetic pore size distribution curve smaller than the corresponding nuclear magnetic radius may be integrated and summed to obtain an area of an area surrounded by the section of curve and the coordinate axis, and then the area of the surrounding area may be used as the cumulative frequency of the corresponding nuclear magnetic radius. Then, based on the plurality of accumulation frequencies, a plurality of nuclear magnetic radii corresponding to the plurality of accumulation frequencies one by one are respectively determined from the correspondence between the nuclear magnetic radii and the accumulation frequencies.

The specific implementation process of integrating and summing the curves can refer to the related technology, and the embodiment of the application is not limited to this. For example, the plurality of cumulative frequencies may be: 5%, 16%, 25%, 50%, 75%, 84%, 95% and 100%.

Step 205: and respectively performing geometric transformation on the plurality of nuclear magnetic radii to obtain a plurality of geometric series in one-to-one correspondence.

After determining the plurality of nuclear magnetic radii through step 205, geometric transformation is performed on each nuclear magnetic radius according to the following formula (4) to obtain a geometric progression corresponding to the corresponding nuclear magnetic radius.

Ψ_j＝log₂R_j (4)

Wherein, in the above formula (4), ψ _j is the geometric progression when the cumulative frequency is j; r _j is the nuclear magnetic radius at the cumulative frequency j; j refers to any one of a plurality of cumulative frequencies.

Step 206: a plurality of pore structure parameters of the target sample core are determined based on the plurality of geometric progression.

Wherein the plurality of pore structure parameters may include at least one of a sorting coefficient, kurtosis, skewness, radius average, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient, and relative sorting coefficient.

In some embodiments, the sorting coefficient may be determined according to the following equation (5), the kurtosis according to the following equation (6), the skewness according to the following equation (7), the radius average according to the following equation (8), the median nuclear magnetic radius according to the following equation (9), the average nuclear magnetic radius according to the following equation (10), the maximum nuclear magnetic radius according to the following equation (11), the homogeneity coefficient according to the following equation (12), and the relative sorting coefficient according to the following equation (13) based on the resulting plurality of geometric progression.

/>

Wherein, in the formulas (5) - (13), S _p refers to the sorting coefficient of the target sample core; k _p refers to the kurtosis of the target sample core; s _kp refers to the skewness of the target sample core; d _M refers to the radius average of the target sample core; r ₅₀ is the median nuclear magnetic radius of the target sample core; Mean core magnetic radius of the target sample core; r _i is the ith nuclear magnetic radius of the target sample core; omega _i refers to the distribution frequency corresponding to the ith nuclear magnetic radius of the target sample core; n refers to the number of the nuclear magnetic radii of the core of the target sample; r ₁₀₀ is the maximum nuclear magnetic radius of the target sample core; alpha refers to the homogeneity coefficient of the core of the target sample; d refers to the relative sorting coefficient of the target sample core; ψ ₅ refers to the geometric progression of the target sample core at an accumulated frequency of 5%; psi ₁₆ is the target sample

Geometric progression when the cumulative frequency of the core is 16%; psi ₂₅ is the geometric progression of the target sample core at an accumulated frequency of 25%; psi ₅₀ is the geometric progression of the target sample core at an accumulated frequency of 50%; psi ₇₅ is the geometric progression of the target sample core at 75% cumulative frequency; psi ₈₄ is the geometric progression of the target sample core at 84% cumulative frequency; psi ₉₅ is the geometric progression of the target sample core at an accumulated frequency of 95%; psi ₁₀₀ refers to the geometric progression at which the cumulative frequency of the target sample core is 100%.

Step 207: and determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the pore structure parameters of the target sample core.

In some embodiments, the reservoir type corresponding to the target sample core may be determined based on the core magnetic porosity, core magnetic permeability, and a plurality of structural parameters of the target sample core according to steps (1) - (3) below.

(1) And determining a characteristic constant of the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the plurality of pore structure parameters of the target sample core.

In some embodiments, the plurality of pore structure parameters includes a sorting coefficient, kurtosis, skewness, radius average, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient, and relative sorting coefficient of the target sample core.

The characteristic constant of the target sample core can be determined by a type discrimination formula based on the nuclear magnetic porosity, the nuclear magnetic permeability and a plurality of pore structure parameters of the target sample core.

Specifically, the core magnetic permeability, the sorting coefficient, the kurtosis, the skewness, the radius average value, the median core magnetic radius, the average core magnetic radius, the maximum core magnetic radius, the homogeneity coefficient and the relative sorting coefficient of the target sample core can be substituted into a type discrimination formula to determine the characteristic constant of the target sample core.

Wherein, the category judgment formula:

In the type discrimination formula, T refers to a characteristic constant of a sample core; s _p is the sorting coefficient of the sample core; k _p refers to the kurtosis of the sample core; s _kp refers to the skewness of the sample core; d _M refers to the radius average of the sample core; r ₅₀ is the median nuclear magnetic radius of the sample core; mean core magnetic radius of the sample core; r ₁₀₀ is the maximum nuclear magnetic radius of the sample core; alpha refers to the homogeneity coefficient of the sample core; d refers to the relative sorting coefficient of the sample core; k _nmr refers to the nuclear magnetic permeability of the sample core.

(2) And determining the characteristic range of the characteristic constant of the target sample core from the characteristic ranges.

The characteristic ranges are preset according to the types of the reservoir, and the characteristic ranges correspond to the types of the reservoir one by one. Therefore, after the characteristic constant of the target sample core is determined, the characteristic range in which the characteristic constant is located can be determined from a plurality of characteristic ranges.

Illustratively, when the theoretical porosity is less than 2%, the corresponding reservoir is a type of reservoir; when the theoretical porosity is between 2 and 6 percent, the corresponding reservoir is a second-class reservoir; when the theoretical porosity is between 6 and 12 percent, the corresponding reservoirs are three types of reservoirs; when the theoretical porosity is more than 12%, the corresponding reservoirs are four types of reservoirs; and for a type of reservoir, the preset characteristic range can be 0.5-1; for the second-class reservoir, the preset characteristic range can be 1-1.5; for three types of reservoirs, the preset characteristic range can be 1.5-2; for four types of reservoirs, the predetermined characteristic range may be 2 to 2.5.

(3) And determining the reservoir type corresponding to the target sample core from the corresponding relation between the pre-stored characteristic range and the reservoir type.

For example, when the characteristic constant is 1.816, the characteristic constant is in the range of 1.5-2, and the reservoir corresponding to the target sample core can be determined to be three types of reservoirs.

In the embodiment of the present application, before determining the characteristic constant of the target sample core through the type discrimination formula in the step (1), the type discrimination formula may be determined as follows.

In some embodiments, the core magnetic porosity and core magnetic permeability of each sample core may be determined based on the length, radius, core magnetic T2 spectrum, and core magnetic signal total amount of the plurality of sample cores obtained above, and the plurality of pore structure parameters of each sample core may be determined based on the core magnetic T2 spectrum and core magnetic conversion coefficient of each sample core. The plurality of pore structure parameters may include a sorting coefficient, kurtosis, skewness, radius average, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient, and relative sorting coefficient. For the nuclear magnetic porosity, the nuclear magnetic permeability, the sorting coefficient, the kurtosis, the skewness, the radius average value, the median nuclear magnetic radius, the average nuclear magnetic radius, the maximum nuclear magnetic radius, the homogeneity coefficient and the relative sorting coefficient of a plurality of sample cores, the reference constant of each sample core can be preset based on the nuclear magnetic porosity of each sample core, then the reference constant is taken as a dependent variable, the nuclear magnetic permeability, the sorting coefficient, the kurtosis, the skewness, the radius average value, the median nuclear magnetic radius, the average nuclear magnetic radius, the maximum nuclear magnetic radius, the homogeneity coefficient and the relative sorting coefficient are taken as independent variables, and the determined category judgment formula is as follows by a multiple regression analysis method:

The specific implementation process of the determined category determination formula by the multiple regression analysis method may refer to related technologies, and the embodiments of the present application are not described herein.

For the specific implementation process of determining the core magnetic porosity and core magnetic permeability of the target sample core and the specific implementation process of determining the pore structure parameters of the target sample core based on the length, radius, core magnetic T2 spectrum and core magnetic signal total amount of the first sample core and the core magnetic T2 spectrum and core magnetic conversion coefficient of the first sample core, the embodiment of the present application will not be repeated herein.

Continuing with the above example, the number of geometric progression for each of the 6 sample cores may be as shown in table 2 below, and the pore structure parameters for each of the 6 sample cores may be as shown in table 3 below.

TABLE 2

TABLE 3 Table 3

Assuming a nuclear magnetic porosity of less than 2%, the corresponding reference constant is 0.75, the nuclear magnetic porosity is greater than 2% and less than 6%, the corresponding reference constant is 1.25; the nuclear magnetic porosity is more than 6% and less than 12%, and the corresponding reference constant is 1.75; the core magnetic porosity is greater than 12%, and the corresponding reference constant is 2.25, so that the reference constants of the 6 sample cores can be determined to be 1.25, 1.75 and 1.75 in sequence.

And determining a multiple regression equation between the characteristic constant and the pore structure parameter through multiple regression analysis based on the reference constant, the nuclear magnetic permeability, the sorting coefficient, the kurtosis, the skewness, the radius average value, the median nuclear magnetic radius, the average nuclear magnetic radius, the maximum nuclear magnetic radius, the homogeneity coefficient and the relative sorting coefficient of the 6 sample cores.

Further, after determining the reservoir type corresponding to the target sample core, a suitable development strategy can be selected to develop the reservoir corresponding to the target sample core based on the determined reservoir type, so as to improve the development efficiency of the reservoir corresponding to the target sample core. The development strategy mainly comprises wellhead position, drilling depth, fracturing mode and the like. For example, when the level of the determined reservoir type is higher, a wellhead may be set in the region of the reservoir corresponding to the target pattern core with the highest level, and a drilling depth of the oil-gas well may be determined based on the depth of the reservoir with the highest level (the drilling depth is greater than the depth of the reservoir with the highest level).

In the embodiment of the application, the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core are determined according to the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target sample core, and a plurality of pore structure parameters of the target sample core are determined according to the greetings T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core. And then determining the reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability and the pore structure parameters of the target sample core, so that the determination process of the reservoir type is simplified, technical guidance is provided for the evaluation of the subsequent reservoir, and the development of the subsequent reservoir is facilitated. In addition, (1) the method adopts a small amount of core to carry out porosity and permeability calibration, establishes a nuclear magnetic resonance T2 spectrum conversion nuclear magnetic porosity and nuclear magnetic permeability conversion method, develops a single experiment type, and realizes reservoir classification evaluation. (2) The method is mainly based on nuclear magnetic resonance test, can be repeated for a plurality of times, and further completes other analysis tests of the sample, thereby increasing the utilization degree of the sample. (3) According to the method, nuclear magnetic pore size distribution is adopted to obtain nuclear magnetic porosity, nuclear magnetic permeability and pore structure parameters, so that various utilization modes of nuclear magnetic resonance data are realized, and the mode of acquiring the parameters is simpler and more convenient. (4) The method has certain reference significance for reservoir classification evaluation of sandstone oil and gas reservoirs.

Fig. 3 is a schematic structural diagram of a device for determining a reservoir type according to an embodiment of the present application. The device is integrated with a terminal, which may be a computer, a notebook computer, etc. Referring to fig. 3, the apparatus includes:

The acquisition module 301 is configured to acquire a length, a radius, a nuclear magnetic T2 spectrum, a total nuclear magnetic signal amount, and a nuclear magnetic conversion coefficient of a core of the target sample.

The first determining module 302 is configured to determine the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core based on the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target core, respectively.

The second determining module 303 is configured to determine a plurality of pore structure parameters of the target sample core based on the nuclear magnetic T2 spectrum and the nuclear magnetic conversion coefficient of the target sample core.

The third determining module 304 is configured to determine a reservoir type corresponding to the target sample core based on the nuclear magnetic porosity, the nuclear magnetic permeability, and the plurality of pore structure parameters of the target sample core.

Optionally, the first determining module 302 includes:

And an eighth determination unit for determining the nuclear magnetic permeability of the target sample core based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core.

Optionally, the sixth determining unit is mainly configured to:

A first conversion formula:

The eighth determination unit is mainly used for:

A second conversion formula: k _nmr＝E·Φ_nmr ⁴·T_2g ²

Optionally, the sixth determining unit is further configured to:

Optionally, the eighth determining unit is further configured to:

obtaining theoretical permeability of a plurality of sample cores;

Optionally, the second determining module includes:

Optionally, the third determining module includes:

the third determining unit is mainly configured to:

The type discrimination formula:

It should be noted that: in the determining device for a reservoir type provided in the foregoing embodiment, when determining a reservoir type corresponding to a target sample core, only the division of the foregoing functional modules is used for illustrating, in practical application, the foregoing functional allocation may be performed by different functional modules according to needs, that is, the internal structure of the apparatus is divided into different functional modules, so as to complete all or part of the functions described above. In addition, the device for determining the reservoir type provided in the foregoing embodiment belongs to the same concept as the method embodiment for determining the reservoir type, and detailed implementation processes of the device are shown in the method embodiment, which is not repeated herein.

Fig. 4 illustrates a block diagram of a terminal 400 according to an exemplary embodiment of the present application. Referring to fig. 4, the terminal 400 may be: smart phones, tablet computers, notebook computers or desktop computers. The terminal 400 may also be referred to by other names as user equipment, portable terminal, laptop terminal, desktop terminal, etc. Referring to fig. 4, the terminal 400 may include a processor 401 and a memory 402.

Processor 401 may include one or more processing cores such as a 4-core processor, an 8-core processor, etc. The processor 401 may be implemented in at least one hardware form of DSP (DIGITAL SIGNAL Processing), FPGA (Field-Programmable gate array), PLA (Programmable Logic Array ). Processor 401 may also include a main processor, which is a processor for processing data in an awake state, also called a CPU (Central Processing Unit ), and a coprocessor; a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 401 may integrate a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content required to be displayed by the display screen. In some embodiments, the processor 401 may also include an AI (ARTIFICIAL INTELLIGENCE ) processor for processing computing operations related to machine learning.

Memory 402 may include one or more computer-readable storage media, which may be non-transitory. Memory 402 may also include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 402 is used to store at least one instruction for execution by processor 401 to implement a method of determining a reservoir type provided by an embodiment of the method of the present application.

In some embodiments, the terminal 400 may further optionally include: a peripheral interface 403 and at least one peripheral. The processor 401, memory 402, and peripheral interface 403 may be connected by a bus or signal line. The individual peripheral devices may be connected to the peripheral device interface 403 via buses, signal lines or a circuit board. Specifically, the peripheral device includes: at least one of radio frequency circuitry 404, a display screen 405, a positioning component 406, and a power supply 407.

Peripheral interface 403 may be used to connect at least one Input/Output (I/O) related peripheral to processor 401 and memory 402. In some embodiments, processor 401, memory 402, and peripheral interface 403 are integrated on the same chip or circuit board; in some other embodiments, either or both of the processor 401, memory 402, and peripheral interface 403 may be implemented on separate chips or circuit boards, which is not limited in this embodiment.

The Radio Frequency circuit 404 is configured to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuitry 404 communicates with a communication network and other communication devices via electromagnetic signals. The radio frequency circuit 404 converts an electrical signal into an electromagnetic signal for transmission, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 404 includes: antenna systems, RF transceivers, one or more amplifiers, tuners, oscillators, digital signal processors, codec chipsets, subscriber identity module cards, and so forth. The radio frequency circuitry 404 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocol includes, but is not limited to: the world wide web, metropolitan area networks, intranets, generation mobile communication networks (2G, 3G, 4G, and 5G), wireless local area networks, and/or WiFi (WIRELESS FIDELITY ) networks. In some embodiments, the radio frequency circuit 404 may further include NFC (NEAR FIELD Communication) related circuits, which is not limited by the present application.

The display screen 405 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display screen 405 is a display screen, the display screen 405 also has the ability to collect touch signals at or above the surface of the display screen 405. The touch signal may be input as a control signal to the processor 401 for processing. At this time, the display screen 405 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, the display 405 may be one, providing a front panel of the terminal 400; in other embodiments, the display 405 may be at least two, and disposed on different surfaces of the terminal 400 or in a folded design; in still other embodiments, the display 405 may be a flexible display disposed on a curved surface or a folded surface of the terminal 400. Even more, the display screen 405 may be arranged in an irregular pattern that is not rectangular, i.e. a shaped screen. The display screen 405 may be made of materials such as an LCD (Liquid CRYSTAL DISPLAY) and an OLED (Organic Light-Emitting Diode).

The locating component 406 is used to locate the current geographic location of the terminal 400 to enable navigation or LBS (Location Based Service, location-based services). The positioning component 406 may be a positioning component based on the United states GPS (Global Positioning System ), the Beidou system of China, or the Galileo system of Russia.

The power supply 407 is used to power the various components in the terminal 400. The power supply 407 may be an alternating current, a direct current, a disposable battery, or a rechargeable battery. When the power supply 407 includes a rechargeable battery, the rechargeable battery may be a wired rechargeable battery or a wireless rechargeable battery. The wired rechargeable battery is a battery charged through a wired line, and the wireless rechargeable battery is a battery charged through a wireless coil. The rechargeable battery may also be used to support fast charge technology.

Those skilled in the art will appreciate that the structure shown in fig. 4 is not limiting of the terminal 400 and may include more or fewer components than shown, or may combine certain components, or may employ a different arrangement of components.

In the above embodiments, there is also provided a non-transitory computer readable storage medium comprising instructions for storing at least one instruction for execution by a processor to implement the method provided by the above embodiments shown in fig. 1 or fig. 2.

Embodiments of the present application also provide a computer program product comprising instructions which, when run on a computer, cause the computer to perform the method provided by the embodiments shown in fig. 1 or fig. 2 described above.

It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program for instructing relevant hardware, where the program may be stored in a computer readable storage medium, and the storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.

The foregoing description of the preferred embodiments of the application is not intended to limit the application to the precise form disclosed, and any such modifications, equivalents, and alternatives falling within the spirit and scope of the application are intended to be included within the scope of the application.

Claims

1. A method of determining a reservoir type, the method comprising:

Based on the length, radius, nuclear magnetic T2 spectrum and total nuclear magnetic signal of the target sample core, respectively determining the nuclear magnetic porosity and nuclear magnetic permeability of the target sample core;

converting the nuclear magnetic T2 spectrum of the target sample core based on the nuclear magnetic conversion coefficient of the target sample core to obtain a nuclear magnetic pore size distribution curve of the target sample core;

Determining a plurality of pore structure parameters of the target sample core based on the plurality of geometric progression, the plurality of pore structure parameters including a sorting coefficient, kurtosis, skewness, radius average, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient, and relative sorting coefficient of the target sample core;

The type discrimination formula:

Wherein, in the type discrimination formula, T refers to a characteristic constant of the sample core; s _p is the sorting coefficient of the sample core; k _p refers to the kurtosis of the sample core; s _kp refers to the skewness of the sample core; d _M refers to the radius average of the sample core; r ₅₀ is the median nuclear magnetic radius of the sample core; Mean core magnetic radius of the sample core; r ₁₀₀ is the maximum nuclear magnetic radius of the sample core; alpha refers to the homogeneity coefficient of the sample core; d refers to the relative sorting coefficient of the sample core; k _nmr refers to the nuclear magnetic permeability of the sample core; the determining mode of the type discrimination formula comprises the following steps: setting a reference constant of each sample core based on the nuclear magnetic porosity of each sample core in a plurality of sample cores, and determining the type discrimination formula through multiple regression analysis by taking the reference constant as a dependent variable and taking the nuclear magnetic permeability, sorting coefficient, kurtosis, skewness, radius average value, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient and relative sorting coefficient of each sample core as independent variables;

2. The method of claim 1, wherein the determining the core magnetic porosity and core magnetic permeability of the target sample core based on the length, radius, core magnetic T2 spectrum, and core magnetic signal total amount, respectively, comprises:

3. The method of claim 2, wherein the determining the nuclear magnetic porosity of the target sample core by a first transformation formula based on the length, radius, and nuclear magnetic signal total amount of the target sample core comprises:

the first conversion formula:

the second conversion formula: k _nmr＝E·Φ_nmr ⁴·T_2g ²;

4. The method of claim 3, wherein prior to determining the nuclear magnetic porosity of the target sample core by the first transformation formula based on the length, radius, and nuclear magnetic signal total amount of the target sample core, further comprising:

5. The method of claim 3, wherein prior to determining the nuclear magnetic permeability of the target sample core by the second conversion formula based on the nuclear magnetic porosity and the geometric relaxation time of the target sample core, further comprising:

obtaining theoretical permeability of the plurality of sample cores;

6. A reservoir type determination apparatus, the apparatus comprising:

The first determining module is used for respectively determining the nuclear magnetic porosity and the nuclear magnetic permeability of the target sample core based on the length, the radius, the nuclear magnetic T2 spectrum and the total nuclear magnetic signal of the target sample core;

The second determining module is used for converting the nuclear magnetic T2 spectrum of the target sample core based on the nuclear magnetic conversion coefficient of the target sample core to obtain a nuclear magnetic pore diameter distribution curve of the target sample core; determining a plurality of nuclear magnetic radii corresponding to a plurality of accumulation frequencies one by one based on the distribution frequency of the nuclear magnetic radii on the nuclear magnetic pore diameter distribution curve of the target sample core; respectively performing geometric transformation on the nuclear magnetic radii to obtain a plurality of geometric series in one-to-one correspondence; determining a plurality of pore structure parameters of the target sample core based on the plurality of geometric progression, the plurality of pore structure parameters including a sorting coefficient, kurtosis, skewness, radius average, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient, and relative sorting coefficient of the target sample core;

The third determining module is used for determining a characteristic constant of the target sample core through a type discrimination formula based on the nuclear magnetic porosity, the nuclear magnetic permeability and a plurality of pore structure parameters of the target sample core; the type discrimination formula: Wherein, in the type discrimination formula, T refers to a characteristic constant of the sample core; s _p is the sorting coefficient of the sample core; k _p refers to the kurtosis of the sample core; s _kp refers to the skewness of the sample core; d _M refers to the radius average of the sample core; r ₅₀ is the median nuclear magnetic radius of the sample core; /(I) Mean core magnetic radius of the sample core; r ₁₀₀ is the maximum nuclear magnetic radius of the sample core; alpha refers to the homogeneity coefficient of the sample core; d refers to the relative sorting coefficient of the sample core; k _nmr refers to the nuclear magnetic permeability of the sample core; the determining mode of the type discrimination formula comprises the following steps: setting a reference constant of each sample core based on the nuclear magnetic porosity of each sample core in a plurality of sample cores, and determining the type discrimination formula through multiple regression analysis by taking the reference constant as a dependent variable and taking the nuclear magnetic permeability, sorting coefficient, kurtosis, skewness, radius average value, median nuclear magnetic radius, average nuclear magnetic radius, maximum nuclear magnetic radius, homogeneity coefficient and relative sorting coefficient of each sample core as independent variables; determining a characteristic range in which a characteristic constant of the target sample core is located from a plurality of characteristic ranges; and determining the reservoir type corresponding to the target sample core from the corresponding relation between the pre-stored characteristic range and the reservoir type.

7. A computer readable storage medium, characterized in that the storage medium has stored therein a computer program which, when executed by a processor, implements the method of any of claims 1-5.