CN113175321B

CN113175321B - Method and device for determining fluid saturation parameters and computer equipment

Info

Publication number: CN113175321B
Application number: CN202110345901.4A
Authority: CN
Inventors: 谢伟彪; 殷秋丽; 司兆伟
Original assignee: Petrochina Co Ltd
Current assignee: Petrochina Co Ltd
Priority date: 2021-03-31
Filing date: 2021-03-31
Publication date: 2023-06-30
Anticipated expiration: 2041-03-31
Also published as: CN113175321A

Abstract

The application provides a method and a device for determining fluid saturation parameters and computer equipment, and belongs to the technical field of reservoir evaluation. The method comprises the following steps: determining a target reservoir of a region to be studied; for each depth of the target reservoir, determining a first map of the target reservoir at the depth based on one-dimensional nuclear magnetic resonance log data of a target well section of a standard well of the region to be studied; acquiring nuclear magnetic resonance experimental data of a plurality of rock samples of a target reservoir, wherein the nuclear magnetic resonance experimental data are obtained through nuclear magnetic resonance experiments of the plurality of rock samples; for each rock sample, determining a second profile of the rock sample based on nuclear magnetic resonance experimental data of the rock sample; determining a first parameter matrix of the target reservoir at depth based on the second patterns of the plurality of rock samples; based on the first map and the first parameter matrix, a fluid saturation parameter of the target aperture is determined, which reduces the cost of determining the fluid saturation parameter of the target aperture.

Description

Method and device for determining fluid saturation parameters and computer equipment

Technical Field

The present disclosure relates to the field of reservoir evaluation technologies, and in particular, to a method and an apparatus for determining a fluid saturation parameter, and a computer device.

Background

In the reservoir development process, the fluid saturation of pores in the reservoir is the core content of the evaluation reservoir, and is related to reservoir parameter calculation, reservoir reserve calculation, reservoir description and reservoir development scheme formulation; thus, determining the fluid saturation parameters of pores in a reservoir is important for efficient reservoir development.

In the related art, the fluid saturation of pores in a reservoir is determined mainly by using a two-dimensional nuclear magnetic resonance logging method, the method comprises the steps of firstly carrying out a two-dimensional nuclear magnetic resonance experiment on a rock sample, determining regional laws of bound water, oil gas and the like on a two-dimensional nuclear magnetic resonance spectrum, then carrying out regional division on the two-dimensional nuclear magnetic resonance spectrum according to the regional laws, and further determining the fluid saturation of the total pores of the reservoir based on the divided regions. Because the two-dimensional nuclear magnetic resonance experiment and the two-dimensional nuclear magnetic resonance logging are performed in two dimensions respectively, and the experimental instrument and the logging instrument for performing the experiment and the logging in the two dimensions respectively are high in price, the two-dimensional nuclear magnetic resonance logging has higher cost, and the cost for determining the fluid saturation of the pores in the reservoir is higher.

Disclosure of Invention

The embodiment of the application provides a method, a device and computer equipment for determining a fluid saturation parameter, which can reduce the cost for determining the fluid saturation parameter of a pore in a reservoir. The technical scheme is as follows:

In one aspect, a method for determining a fluid saturation parameter is provided, the method comprising:

determining a target reservoir of a region to be studied;

for each depth of the target reservoir, determining a first map of the target reservoir at the depth based on one-dimensional nuclear magnetic resonance logging information of a target well section of a standard well of the region to be studied, wherein the one-dimensional nuclear magnetic resonance logging information is obtained through one-dimensional nuclear magnetic resonance logging engineering of the target well section, the target well section is a well section of the standard well corresponding to the depth of the target reservoir, and the first map is a relation curve between pore parameters and time;

acquiring nuclear magnetic resonance experimental data of a plurality of rock samples of the target reservoir, wherein the nuclear magnetic resonance experimental data are obtained through nuclear magnetic resonance experiments of the plurality of rock samples;

for each rock sample, determining a second map of the rock sample based on nuclear magnetic resonance experimental data of the rock sample, wherein the second map is a relation curve between pore parameters and time;

determining a first parameter matrix of the target reservoir at the depth based on a second profile of the plurality of rock samples;

and determining a fluid saturation parameter of a target pore based on the first map and the first parameter matrix, wherein the target pore is a pore of the target reservoir at the depth.

In one possible implementation, the nuclear magnetic resonance experimental data of the rock sample includes saturated parameter water nuclear magnetic resonance experimental data, centrifugal nuclear magnetic resonance experimental data and saturated parameter oil nuclear magnetic resonance experimental data of the rock sample at different echo intervals;

the determining a second map of the rock sample based on nuclear magnetic resonance experimental data of the rock sample comprises:

determining a first saturated parameter water map and a second saturated parameter water map of the rock sample based on saturated parameter water nuclear magnetic resonance experimental data of the rock sample at a first echo interval time and a second echo interval time respectively;

determining a first centrifugal map and a second centrifugal map of the rock sample based on centrifugal nuclear magnetic resonance experimental data of the rock sample at the first echo interval time and the second echo interval time respectively;

determining a first saturated parameter oil spectrum and a second saturated parameter oil spectrum of the rock sample based on saturated parameter oil nuclear magnetic resonance experimental data of the rock sample at the first echo interval time and the second echo interval time respectively;

and forming the first saturated parameter water map, the second saturated parameter water map, the first centrifugal map, the second centrifugal map, the first saturated parameter oil map and the second saturated parameter oil map into the second map.

In one possible implementation, the determining a first parameter matrix of the target reservoir at the depth based on the second patterns of the plurality of rock samples includes:

for each rock sample, determining a first time point and a second time point of the rock sample, wherein the first time point is a time point corresponding to the maximum peak value of a first centrifugal map of the rock sample, and the second time point is a time point corresponding to the maximum peak value of a second centrifugal map of the rock sample;

determining a first bound water parameter and a second bound water parameter based on an average of a first time point and an average of a second time point of the plurality of rock samples, respectively;

determining a first movable water parameter based on a first saturation parameter water profile and a first centrifugal profile of the plurality of rock samples;

determining a second movable water parameter based on a second saturation parameter water profile and a second centrifugal profile of the plurality of rock samples;

determining a first oil parameter and a second oil parameter based on a first saturated parameter oil profile and a second saturated parameter oil profile of the plurality of rock samples, respectively;

generating the first parameter matrix comprising the first bound water parameter, the second bound water parameter, the first movable water parameter, the second movable water parameter, the first oil parameter, and the second oil parameter.

In one possible implementation, the determining a first movable water parameter based on the first saturation parameter water profile and the first centrifugal profile of the plurality of rock samples includes:

for each rock sample, determining a first difference spectrum of the rock sample based on a first saturation parameter water spectrum and a first centrifugal spectrum of the rock sample, wherein the first difference spectrum is a spectrum consisting of differences between first saturation water parameters and second centrifugal parameters at the same time in the first saturation parameter water spectrum and the first centrifugal spectrum;

determining a third time point corresponding to the maximum peak value of the first difference map;

the first mobile water parameter is determined based on an average of third time points of the plurality of rock samples.

In one possible implementation, the determining a second movable water parameter based on a second saturation parameter water profile and a second centrifugal profile of the plurality of rock samples includes:

for each rock sample, determining a second difference spectrum of the rock sample based on a second saturation parameter water spectrum and a second centrifugal spectrum of the rock sample, wherein the second difference spectrum is a spectrum consisting of differences between second saturation water parameters and second centrifugal parameters at the same time in the second saturation parameter water spectrum and the second centrifugal spectrum;

Determining a fourth time point corresponding to the maximum peak value of the second difference value map;

the second mobile water parameter is determined based on an average of a fourth time point of the plurality of rock samples.

In one possible implementation, the determining the fluid saturation parameter of the target aperture based on the first map and the first parameter matrix includes:

determining a second parameter matrix based on the first map and the first parameter matrix;

determining, based on the first map and the second parameter matrix, relationship data for the target reservoir at the depth, the relationship data being between the second parameter matrix and pore fluid volume;

based on the relationship data, a fluid saturation parameter of the target aperture is determined.

In one possible implementation, the first pattern includes a third target pattern and a fourth target pattern, the third target pattern and the fourth target pattern being patterns of the formation at a first echo interval time and a second echo interval time, respectively;

the determining a second parameter matrix based on the first map and the first parameter matrix includes:

determining a third tie water parameter, a third movable water parameter, and a third oil parameter based on the third target profile and the first parameter matrix;

Determining a fourth tie water parameter, a fourth movable water parameter, and a fourth oil parameter based on the fourth target profile and the first parameter matrix;

generating a second parameter matrix comprising the third bound water parameter, the third movable water parameter, the third oil parameter, the fourth bound water parameter, the fourth movable water parameter, and the fourth oil parameter.

In one possible implementation, the pore fluid comprises bound water, mobile water, and oil, and the relationship data comprises a second parameter matrix, a dependent variable matrix, and an independent variable matrix;

the determining, based on the first map and the second parameter matrix, relationship data for the target reservoir at the depth includes:

determining a first median time parameter and a second median time parameter based on the third target profile and the fourth target profile, respectively;

generating a dependent variable matrix comprising the first median time parameter and the second median time parameter;

the relationship data including the second parameter matrix and the independent variable matrix, and having the saturation of the bound water, the saturation of the movable water, and the saturation of the oil as independent variable matrices is generated.

In another aspect, there is provided a device for determining a fluid saturation parameter, the device comprising:

the first determining module is used for determining a target reservoir of the area to be researched;

the second determining module is used for determining a first map of the target reservoir at the depth based on one-dimensional nuclear magnetic resonance logging information of a target well section of a standard well of the region to be researched, wherein the one-dimensional nuclear magnetic resonance logging information is obtained through one-dimensional nuclear magnetic resonance logging engineering of the target well section, the target well section is a well section of the standard well corresponding to the depth of the target reservoir, and the first map is a relation curve between pore parameters and time;

the acquisition module is used for acquiring nuclear magnetic resonance experimental data of a plurality of rock samples of the target reservoir, wherein the nuclear magnetic resonance experimental data are obtained through nuclear magnetic resonance experiments of the plurality of rock samples;

a third determining module, configured to determine, for each rock sample, a second spectrum of the rock sample based on nuclear magnetic resonance experimental data of the rock sample, where the second spectrum is a relationship between pore parameters and time;

a fourth determination module for determining a first parameter matrix of the target reservoir at the depth based on a second profile of the plurality of rock samples;

And a fifth determining module, configured to determine a fluid saturation parameter of a target pore based on the first map and the first parameter matrix, where the target pore is a pore of the target reservoir at the depth.

the third determining module includes:

the first determining unit is used for determining a first saturated parameter water spectrum and a second saturated parameter water spectrum of the rock sample based on saturated parameter water nuclear magnetic resonance experimental data of the rock sample at a first echo interval time and a second echo interval time respectively;

a second determining unit, configured to determine a first centrifugal spectrum and a second centrifugal spectrum of the rock sample based on centrifugal nuclear magnetic resonance experimental data of the rock sample at the first echo interval time and the second echo interval time, respectively;

the third determining unit is used for determining a first saturated parameter oil spectrum and a second saturated parameter oil spectrum of the rock sample based on saturated parameter oil nuclear magnetic resonance experimental data of the rock sample at the first echo interval time and the second echo interval time respectively;

The composition unit is used for composing the first saturated parameter water spectrum, the second saturated parameter water spectrum, the first centrifugal spectrum, the second centrifugal spectrum, the first saturated parameter oil spectrum and the second saturated parameter oil spectrum into the second spectrum.

In one possible implementation manner, the fourth determining module includes:

a fourth determining unit, configured to determine, for each rock sample, a first time point and a second time point of the rock sample, where the first time point is a time point corresponding to a maximum peak value of a first centrifugal spectrum of the rock sample, and the second time point is a time point corresponding to a maximum peak value of a second centrifugal spectrum of the rock sample;

a fifth determining unit for determining a first bound water parameter and a second bound water parameter based on an average value of a first time point and an average value of a second time point of the plurality of rock samples, respectively;

a sixth determining unit for determining a first movable water parameter based on a first saturation parameter water profile and a first centrifugal profile of the plurality of rock samples;

a seventh determining unit for determining a second movable water parameter based on a second saturation parameter water profile and a second centrifugal profile of the plurality of rock samples;

An eighth determining unit for determining a first oil parameter and a second oil parameter based on a first saturation parameter oil map and a second saturation parameter oil map of the plurality of rock samples, respectively;

a generation unit for generating the first parameter matrix comprising the first bound water parameter, the second bound water parameter, the first movable water parameter, the second movable water parameter, the first oil parameter and the second oil parameter.

In one possible implementation manner, the sixth determining unit includes:

a first determining subunit, configured to determine, for each rock sample, a first difference spectrum of the rock sample based on a first saturation parameter water spectrum and a first centrifugal spectrum of the rock sample, where the first difference spectrum is a spectrum composed of differences between first saturation parameter water spectrum and second centrifugal parameter at the same time in the first saturation parameter water spectrum and the first centrifugal spectrum;

a second determining subunit, configured to determine a third time point corresponding to a maximum peak of the first difference map;

a third determination subunit for determining the first movable water parameter based on an average of third time points of the plurality of rock samples.

In one possible implementation manner, the seventh determining unit includes:

a fourth determination subunit, configured to determine, for each rock sample, a second difference spectrum of the rock sample based on a second saturation parameter water spectrum and a second centrifugal spectrum of the rock sample, where the second difference spectrum is a spectrum composed of differences between second saturation parameter water spectrum and second centrifugal parameter at the same time in the second saturation parameter water spectrum and the second centrifugal spectrum;

a fifth determining subunit, configured to determine a fourth time point corresponding to a maximum peak of the second difference map;

a sixth determination subunit for determining the second movable water parameter based on an average of fourth time points of the plurality of rock samples.

In one possible implementation manner, the fifth determining module includes:

a ninth determining unit configured to determine a second parameter matrix based on the first map and the first parameter matrix;

a tenth determination unit configured to determine, based on the first map and the second parameter matrix, relationship data of the target reservoir at the depth, the relationship data being relationship data between the second parameter matrix and a pore fluid volume;

An eleventh determining unit for determining a fluid saturation parameter of the target aperture based on the relationship data.

the ninth determination unit includes:

a seventh determining subunit configured to determine a third tie water parameter, a third movable water parameter, and a third oil parameter based on the third target map and the first parameter matrix;

an eighth determination subunit configured to determine a fourth bound water parameter, a fourth movable water parameter, and a fourth oil parameter based on the fourth target map and the first parameter matrix;

a first generation subunit for generating a second parameter matrix comprising the third bound water parameter, the third movable water parameter, the third oil parameter, the fourth bound water parameter, the fourth movable water parameter, and the fourth oil parameter.

The tenth determination unit includes:

a ninth determining subunit configured to determine a first median time parameter and a second median time parameter based on the third target profile and the fourth target profile, respectively;

a second generation subunit configured to generate a dependent variable matrix including the first median time parameter and the second median time parameter;

a third generation subunit configured to generate the relationship data including the second parameter matrix and the argument matrix, and having the saturation of the bound water, the saturation of the movable water, and the saturation of the oil as argument matrices.

In another aspect, a computer device is provided, the computer device including a processor and a memory, the memory storing at least one program code, the at least one program code loaded and executed by the processor to implement instructions of a method for determining a fluid saturation parameter according to any one of the implementations described above.

In another aspect, a computer readable storage medium having stored therein at least one program code loaded and executed by a processor to implement the steps in the method of determining a fluid saturation parameter according to any one of the above implementations is provided.

The beneficial effects of the technical scheme provided by the embodiment of the application at least comprise:

the embodiment of the application provides a method for determining fluid saturation parameters, wherein a first map and a second map are determined according to one-dimensional nuclear magnetic resonance logging data and nuclear magnetic resonance experimental data obtained through nuclear magnetic resonance experiments of one-dimensional nuclear magnetic resonance logging engineering and a rock sample, and the cost for determining the first map and the second map is reduced due to low measurement cost of the nuclear magnetic resonance experiments of the one-dimensional nuclear magnetic resonance logging engineering and the rock sample.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.

FIG. 1 is a flow chart of a method for determining a fluid saturation parameter provided in an embodiment of the present application;

FIG. 2 is a T2 spectrum provided by an embodiment of the present application;

FIG. 3 is a graph comparing the results of determining a fluid saturation parameter according to an embodiment of the present application;

FIG. 4 is a block diagram of a device for determining a fluid saturation parameter provided in an embodiment of the present application;

fig. 5 is a block diagram of a computer device provided in 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 below with reference to the accompanying drawings.

The terms "first," "second," "third," and "fourth" and the like in the description and in the claims of this application and in the drawings, are used for distinguishing between different objects and not for describing a particular sequential order. Furthermore, the terms "comprising," "including," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not limited to only those listed steps or elements but may include other steps or elements not listed or inherent to such process, method, article, or apparatus.

An embodiment of the present application provides a method for determining a fluid saturation parameter, referring to fig. 1, the method includes:

step 101: the computer device determines a target reservoir for the region under investigation.

The region to be researched is an oil-gas field region, and the target reservoir is a stratum rich in oil gas in the region to be researched.

Step 102: for each depth of the target reservoir, the computer device determines a first profile of the target reservoir at that depth based on one-dimensional nuclear magnetic resonance log data of a target well section of a standard well of the region to be studied.

The one-dimensional nuclear magnetic resonance logging information is obtained through one-dimensional nuclear magnetic resonance logging engineering of a target well section, the target well section is a well section of a standard well corresponding to the depth of a target reservoir, the first map is a relation curve between a pore parameter and time, the ordinate is the pore parameter, the pore parameter is a pore component, the abscissa is time, the time is relaxation time, and the relaxation time is time when transverse magnetization intensity disappears during nuclear magnetic resonance logging.

The one-dimensional nuclear magnetic resonance logging data are P-type nuclear magnetic resonance logging data, the data are obtained by logging engineering through a P-type nuclear magnetic resonance instrument, the first spectrum is a T2 spectrum obtained by logging through the P-type nuclear magnetic resonance instrument, and the T2 spectrum is an array spectrum which is obtained according to echo string signal inversion in the nuclear magnetic resonance logging data and reflects pore fluid distribution.

The P-type nuclear magnetic resonance logging data processing is standardized according to the procedure specified in the standard of the nuclear magnetic resonance logging data processing and interpretation specification SY/T6617-2016.

The first map comprises a third target map and a fourth target map, the third target map and the fourth target map are maps of the stratum under a first echo interval time and a second echo interval time respectively, and the first echo interval time is larger than the second echo interval time.

In the embodiment of the application, because the P-type nuclear magnetic resonance logging is only performed in one dimension, the logging cost is low, and the cost for determining the first map through one-dimensional nuclear magnetic resonance logging data is further reduced.

Step 103: the computer equipment acquires nuclear magnetic resonance experimental data of a plurality of rock samples of the target reservoir, wherein the nuclear magnetic resonance experimental data are obtained through nuclear magnetic resonance experiments of the plurality of rock samples.

The rock samples are rock samples of stratum at the same depth as the first map, the number of the rock samples is at least two, and the rock samples are selected according to a process specified by the SY/T6028-94 standard of the exploratory well test project sampling and achievement requirement. The nuclear magnetic resonance experiment of the rock sample is carried out according to the procedure specified in the standard of the rock sample nuclear magnetic resonance parameter laboratory measurement specification SY/T6490-2014.

The nuclear magnetic resonance experimental data of the rock sample comprise saturated parameter water nuclear magnetic resonance experimental data, centrifugal nuclear magnetic resonance experimental data and saturated parameter oil nuclear magnetic resonance experimental data of the rock sample at different echo interval times.

The saturated parameter water nuclear magnetic resonance experimental data are obtained by carrying out nuclear magnetic resonance experiments on the rock sample after the rock sample is saturated by water. The centrifugal nuclear magnetic resonance experimental data are data obtained by centrifuging a rock sample saturated with water and then performing a nuclear magnetic resonance experiment, and the centrifuged rock sample does not contain movable water. The saturated parameter oil nuclear magnetic resonance experimental data is data obtained by centrifuging a rock sample saturated with water, then saturating the rock sample with oil, and performing nuclear magnetic resonance experiments on the rock sample saturated with oil.

Before the saturation oil nuclear magnetic resonance experiment is performed, an oil viscosity experiment is performed on target oil in a target reservoir to obtain the viscosity of the target oil, and then the oil adopted in the saturation parameter oil nuclear magnetic resonance experiment can be the target oil or kerosene with the same viscosity as the target oil. The viscosity test of the oil was carried out according to the procedure specified in the standard of the crude oil viscosity determination rotational viscometer balance method SY/T0520-2008.

Step 104: the computer device determines, for each rock sample, a second profile of the rock sample based on nuclear magnetic resonance experimental data of the rock sample.

The second spectrum is a relation curve between the pore parameters and time, and the second spectrum is a T2 spectrum obtained through nuclear magnetic resonance experimental data of the rock sample.

This step can be achieved by the following steps (1) - (4):

(1) The computer equipment determines a first saturated parameter water map and a second saturated parameter water map of the rock sample based on saturated parameter water nuclear magnetic resonance experimental data of the rock sample at a first echo interval time and a second echo interval time respectively.

Wherein, the first saturation parameter water spectrum and the second saturation parameter water spectrum are both T2 spectrums. The saturated parameter water nuclear magnetic resonance experimental data under the first echo interval time comprises a plurality of pore parameters under the first echo interval time and time points corresponding to the pore parameters, and the computer equipment draws a first saturated parameter water map based on the pore parameters and the time points corresponding to the pore parameters. The saturated parameter water nuclear magnetic resonance experimental data at the second echo interval time comprises a plurality of pore parameters and time points corresponding to the pore parameters at the second echo interval time, and the computer equipment draws a second saturated parameter water map based on the pore parameters and the time points corresponding to the pore parameters.

(2) The computer device determines a first centrifugal map and a second centrifugal map of the rock sample based on centrifugal nuclear magnetic resonance experimental data of the rock sample at a first echo interval time and a second echo interval time, respectively.

Wherein, the first centrifugal spectrum and the second centrifugal spectrum are both T2 spectrums. The centrifugal nuclear magnetic resonance experimental data at the first echo interval time includes a plurality of pore parameters at the first echo interval time and time points corresponding to the plurality of pore parameters, and the computer device draws a first centrifugal map based on the plurality of pore parameters and the time points corresponding to the plurality of pore parameters. The centrifugal nuclear magnetic resonance experimental data at the second echo interval time includes a plurality of pore parameters at the second echo interval time and time points corresponding to the plurality of pore parameters, and the computer device draws a second centrifugal map based on the plurality of pore parameters and the time points corresponding to the plurality of pore parameters.

(3) The computer equipment determines a first saturated parameter oil spectrum and a second saturated parameter oil spectrum of the rock sample based on saturated parameter oil nuclear magnetic resonance experimental data of the rock sample at a first echo interval time and a second echo interval time respectively.

Wherein, the first saturation parameter oil spectrum and the second saturation parameter oil spectrum are both T2 spectrums. The saturated parameter oil nuclear magnetic resonance experimental data at the first echo interval time comprises a plurality of pore parameters and time points corresponding to the pore parameters at the first echo interval time, and the computer equipment draws a first saturated parameter oil map based on the pore parameters and the time points corresponding to the pore parameters. The saturated parameter oil nuclear magnetic resonance experimental data at the second echo interval time comprises a plurality of pore parameters and time points corresponding to the pore parameters at the second echo interval time, and the computer equipment draws a second saturated parameter oil map based on the pore parameters and the time points corresponding to the pore parameters.

It should be noted that, steps (1) to (3) are not strictly time-sequential, and may be performed in any order, or may be performed simultaneously.

(4) The computer equipment forms a first saturated parameter water map, a second saturated parameter water map, a first centrifugal map, a second centrifugal map, a first saturated parameter oil map and a second saturated parameter oil map into a second map.

Step 105: the computer device determines a first parameter matrix for the target reservoir at depth based on a second map of the plurality of rock samples.

Wherein the first parameter matrix is a matrix of 2 rows and 3 columns. This step can be achieved by the following steps (1) - (6):

(1) For each rock sample, the computer device determines a first time point and a second time point of the rock sample, wherein the first time point is a time point corresponding to the maximum peak value of a first centrifugal spectrum of the rock sample, and the second time point is a time point corresponding to the maximum peak value of a second centrifugal spectrum of the rock sample.

Wherein, the first centrifugal spectrum and the second centrifugal spectrum are both in a unimodal form. The computer device determines a target time point corresponding to a maximum peak in the map from the first centrifugal map, determines the target time point as a first time point, determines a target time point corresponding to a maximum peak in the map from the second centrifugal map, and determines the target time point as a second time point.

The maximum peak value in the first centrifugal map represents the maximum pore component in the rock sample, the maximum pore component corresponds to the peak volume of the bound water at the first echo interval time, and the first time point is the peak volume time point of the bound water. The maximum peak in the second centrifugal profile represents the maximum pore component in the rock sample, which corresponds to the peak volume of bound water at the second echo interval time, i.e. the second point in time is the peak volume of bound water point in time.

(2) The computer device determines a first bound water parameter and a second bound water parameter based on an average of a first time point and an average of a second time point of the plurality of rock samples, respectively.

The computer device averages a first time point of the plurality of rock samples, takes the average value as a first bound water parameter of the first parameter matrix, averages a third time point of the plurality of rock samples, and takes the average value as a second bound water parameter of the first parameter matrix.

Wherein the first bound water parameter represents an average peak volume time point of bound water at a first echo interval time and the second bound water parameter represents an average peak volume time point of bound water at a second echo interval time.

The first bound water parameter is an element A11 of the 1 st row and the 1 st column in the first parameter matrix, and the second bound water parameter is an element A21 of the 2 nd row and the 1 st column in the first parameter matrix.

(3) The computer device determines a first movable water parameter based on a first saturation parameter water profile and a first centrifugal profile of the plurality of rock samples.

This step is achieved by the following steps A1-A3:

a1: the computer device determines, for each rock sample, a first difference profile of the rock sample based on a first saturation parameter water profile and a first centrifugal profile of the rock sample.

The first difference value map is a map formed by differences between the first saturated water parameter and the first centrifugal parameter at the same time in the first saturated water parameter map and the first centrifugal map.

Wherein the first saturated water parameter comprises a pore component for each time point in the first saturated water parameter water profile and the first centrifugal parameter comprises a pore component for each time point in the first centrifugal profile. The computer equipment subtracts the pore component of the same time point of the first centrifugal spectrum from the pore component of each time point in the first saturated parameter water spectrum to obtain a pore component difference value of each time point under the first echo interval time, and draws a first difference value spectrum based on the pore component difference value and each time point.

The pore component of each time point of the first saturation parameter water map is subtracted by the pore component of the same time point of the first centrifugal map, namely the pore component of saturated water at each time point is subtracted by the pore component of bound water, so that the pore component of movable water is obtained, namely the first difference value map represents the T2 spectrum of the movable water in the rock sample at the first echo interval time.

A2: the computer device determines a third point in time corresponding to a maximum peak of the first difference map.

The computer equipment determines a target time point corresponding to the maximum peak value of the first difference value spectrum, and takes the target time point as a third time point.

The maximum peak value in the first difference value map represents the maximum pore component in the rock sample, the maximum pore component corresponds to the peak volume of the movable water at the first echo interval time, and the third time point is the peak volume time point of the movable water.

A3: the computer device determines a first movable water parameter based on an average of a third time point of the plurality of rock samples.

Wherein the computer device averages a third time point of the plurality of rock samples and takes the average as the first movable water parameter.

Wherein the first movable water parameter represents an average peak volume time point of the movable water at the first echo interval time. The first movable water parameter is element a12 of row 1 and column 2 in the first parameter matrix.

(4) The computer device determines a second movable water parameter based on a second saturation parameter water profile and a second centrifugal profile of the plurality of rock samples.

This step is achieved by the following steps A1-A3:

a1: the computer equipment determines a second difference value map of the rock sample based on a second saturated parameter water map and a second centrifugal map of the rock sample for each rock sample, wherein the second difference value map is a map formed by differences between second saturated water parameters and second centrifugal parameters at the same time in the second saturated parameter water map and the second centrifugal map.

Wherein the second saturated water parameter comprises a pore component for each time point in the second saturated water parameter water profile and the second centrifugal parameter comprises a pore component for each time point in the second centrifugal profile. The computer equipment subtracts the pore component of the second centrifugal spectrum at the same time point from the pore component of each time point in the second saturated parameter water spectrum to obtain a pore component difference value of each time point under the second echo interval time, and draws a second difference value spectrum based on the pore component difference value and each time point.

The pore component of each time point in the second saturation parameter water map is subtracted by the pore component of the same time point of the second centrifugal map, namely, the pore component of the saturated water at each time point is subtracted by the pore component of the bound water, so that the pore component of the movable water is obtained, namely, the second difference value map represents the T2 spectrum of the movable water in the rock sample at the second echo interval time.

A2: the computer device determines a fourth point in time corresponding to a maximum peak of the second difference map.

The computer equipment determines a target time point corresponding to the maximum peak value of the second difference value spectrum, and takes the target time point as a fourth time point.

The maximum peak value in the second difference value map represents the maximum pore component in the rock sample, the maximum pore component corresponds to the peak volume of the movable water at the second echo interval time, and the fourth time point is the peak volume time point of the movable water.

A3: the computer device determines a second movable water parameter based on an average of a fourth time point of the plurality of rock samples.

Wherein the computer device averages a fourth point in time of the plurality of rock samples and takes the average as the second movable water parameter.

Wherein the second movable water parameter represents an average peak volume time point of the movable water at the second echo interval time. The second movable water parameter is element a22 of row 2 and column 2 in the first parameter matrix.

(5) The computer device determines a first oil parameter and a second oil parameter based on the first saturated parameter oil spectrum and the second saturated parameter oil spectrum, respectively, of the plurality of rock samples.

The first oil parameter is an element A13 of the 1 st row and the 3 rd column in the first parameter matrix, and the computer equipment determines the first oil parameter based on first saturated parameter oil patterns of a plurality of rock samples.

If the first saturated parameter oil patterns of the plurality of rock samples are in a unimodal form, the computer equipment determines an average value of a fifth time point corresponding to the maximum peak value of the plurality of first saturated parameter oil patterns, and takes the average value of the fifth time point as the first oil parameter.

If the first saturated parameter oil patterns of the plurality of rock samples are in a bimodal form, the computer equipment determines an average value of a fifth time point corresponding to the peak value of the second peak of the plurality of first saturated parameter oil patterns, and takes the average value of the fifth time point as the first oil parameter.

If the first saturated parameter oil patterns of the plurality of rock samples are in a three-peak form, the computer equipment determines an average value of a fifth time point corresponding to the peak value of the third peak of the plurality of first saturated parameter oil patterns, and takes the average value of the fifth time point as the first oil parameter.

The fifth time point is a peak volume time point of oil in the oil map of the first saturation parameter, and the first oil parameter represents an average peak volume time point of oil at the first echo interval time.

Wherein the second oil parameter is element A23 of the 2 nd row and 3 rd column in the first parameter matrix, and the computer equipment determines the second oil parameter based on second saturated parameter oil patterns of the plurality of rock samples.

If the second saturated parameter oil patterns of the plurality of rock samples are in a unimodal form, the computer equipment determines an average value of a sixth time point corresponding to the maximum peak value of the plurality of second saturated parameter oil patterns, and takes the average value of the sixth time point as the second oil parameter.

If the second saturated parameter oil patterns of the plurality of rock samples are in a bimodal form, the computer equipment determines a sixth time point corresponding to the peak value of the second peaks of the plurality of second saturated parameter oil patterns, and takes the average value of the sixth time point as the second oil parameter.

If the second saturated parameter oil spectrum of the plurality of rock samples is in a three-peak form, the computer equipment determines a sixth time point corresponding to the peak value of the third peak of the second saturated parameter oil spectrum, and the average value of the sixth time point is taken as the second oil parameter.

The sixth time point is a peak volume time point of oil in the oil map of the second saturation parameter, and the second oil parameter represents an average peak volume time point of oil under the first echo interval time.

It should be noted that, steps (2) to (5) are not strictly time-sequential, and may be performed in any order, or may be performed simultaneously.

(6) The computer device generates a first parameter matrix comprising a first bound water parameter, a second bound water parameter, a first movable water parameter, a second movable water parameter, a first oil parameter, and a second oil parameter.

The first parameter matrix is a matrix of 2 rows and 3 columns, and the first bound water parameter, the second bound water parameter, the first movable water parameter, the second movable water parameter, the first oil parameter and the second oil parameter are respectively used as elements of the 1 st row, the 1 st column, the 2 nd row, the 2 nd column, the 1 st row, the 3 rd column and the 2 nd row and the 3 rd column in the first parameter matrix.

Step 106: the computer device determines a fluid saturation parameter of the target aperture based on the first map and the first parameter matrix.

Wherein the target pore is a pore of the target reservoir at the depth, the fluid in the target pore comprises bound water, mobile water and oil, and the fluid saturation parameters of the target pore comprise saturation of the bound water, saturation of the mobile water and saturation of the oil.

This step can be achieved by the following steps (1) - (4):

(1) The computer device determines a second parameter matrix based on the first map and the first parameter matrix.

Wherein, this step is realized by the following steps A1-A4:

a1: the computer device determines a third tie water parameter, a third movable water parameter, and a third oil parameter based on the third target profile and the first parameter matrix.

The third constraint water parameter is an element C11 of the 1 st row and the 1 st column in the second parameter matrix, the third movable water parameter is an element C12 of the 1 st row and the 2 nd column in the second parameter matrix, and the third oil parameter is an element C13 of the 1 st row and the 3 rd column in the second parameter matrix.

If the third target map is in a unimodal state, determining a first peak time point corresponding to the maximum peak of the third target map, marking as F1, and under the condition that the absolute value of F1-A11 is less than or equal to the absolute value of F1-A12, the first peak time point is close to the peak volume time point of bound water, indicating that the first peak time point is the peak volume time point of bound water, assigning the first peak time point F1 to C11, assigning A12 to C12, and assigning A13 to C13. In the case of |F1-A11| > |F1-A12| and |F1-A13| > |F1-A12|, the first peak time point is close to the peak volume time point of the movable water, indicating that the first peak time point is the peak volume time point of the movable water, assigning A11 to C11, F1 to C12, and A13 to C13. In the case of |F1-A11| > |F1-A12| > |F1-A13|, the first peak time point is close to the peak volume time point of the oil, indicating that the first peak time point is the peak volume time point of the oil, A11 is assigned to C11, A12 is assigned to C12, and F1 is assigned to C13.

If the third target map is in a bimodal state, determining a first peak time point corresponding to a peak value of a first peak of the third target map, marking as F1, determining a second peak time point corresponding to a peak value of a second peak of the third target map, marking as F2, and when |F1-A11| is less than or equal to |F1-A12|, the first peak time point is close to the peak time point of the bound water, which means that the first peak time point is the peak volume time point of the bound water, the second peak time point is the peak volume time point of the movable water, assigning F1 to C11, F2 to C12, and A13 to C13. In the case of |F1-A11| > |F1-A12|, the first peak time point is close to the peak time point of the movable water, which means that the first peak time point is the peak volume time point of the movable water, the second peak time point is the peak volume time point of the oil, A11 is assigned to C11, F1 is assigned to C12, and F2 is assigned to C13.

If the third target map is in a three-peak state, determining a first peak time point corresponding to a peak value of a first peak of the third target map, marking as F1, determining a second peak time point corresponding to a peak value of a second peak of the third target map, marking as F2, determining a third peak time point corresponding to a peak value of a third peak of the third target map, marking as F3, wherein the first peak time point, the second peak time point and the third peak time point are respectively a peak volume time point of bound water, a peak volume time point of movable water and a peak volume time point of oil, assigning F1 to C11, F2 to C12 and F3 to C13.

Referring to fig. 2, the solid line pattern in fig. 2 is a third target pattern, which is in a three-peak state, and the first peak time point, the second peak time point, and the third peak time point are assigned to C11, C12, and C13, respectively.

A2: the computer device determines a fourth bound water parameter, a fourth movable water parameter, and a fourth oil parameter based on the fourth target profile and the first parameter matrix.

The fourth bound water parameter is an element C21 of the 2 nd row and the 1 st column in the second parameter matrix, the fourth movable water parameter is an element C22 of the 2 nd row and the 2 nd column in the second parameter matrix, and the fourth oil parameter is an element C23 of the 2 nd row and the 3 rd column in the second parameter matrix.

If the fourth target map is in a unimodal state, determining a first peak time point corresponding to the maximum peak of the fourth target map, marking as E1, and under the condition that the absolute value of E1-A21 is less than or equal to the absolute value of E1-A22, enabling the first peak time point to be close to the peak volume time point of bound water, indicating that the first peak time point is the peak volume time point of bound water, assigning E1 to C21, assigning A22 to C22, and assigning A23 to C23. In the case of |E1-A21| > E1-A22| and |E1-A23| > E1-A22|, the first peak time point is close to the peak volume time point of the movable water, indicating that the first peak time point is the peak volume time point of the movable water, A11 is assigned to C11, A21 is assigned to C21, E1 is assigned to C22, and A23 is assigned to C23. In the case of |E1-A21| > E1-A22| and |E1-A22| > E1-A23|, the first peak time point is close to the peak volume time point of the oil, indicating that the first peak time point is the peak volume time point of the oil, assigning A21 to C21, A22 to C22, and E1 to C23.

If the fourth target map is in a bimodal state, determining a first peak time point corresponding to a peak value of a first peak of the fourth target map, marking as E1, determining a second peak time point corresponding to a peak value of a second peak of the fourth target map, marking as E2, and when the absolute value of E1-A21 is less than or equal to the absolute value of E1-A22, the first peak time point is close to the peak time point of bound water, which means that the first peak time point is the peak volume time point of bound water, the second peak time point is the peak volume time point of movable water, assigning E1 to C21, assigning E2 to C22, and assigning A23 to C23. In the case of |E1-A21| > |E1-A22|, the first peak time point is close to the peak time point of the movable water, which means that the first peak time point is the peak volume time point of the movable water, the second peak time point is the peak volume time point of the oil, A21 is assigned to C21, E1 is assigned to C22, and E2 is assigned to C23.

With continued reference to fig. 2, the solid line pattern in fig. 2 is a fourth target pattern that is bimodal, the first peak time point is the peak volume time point of bound water, the second peak time point is the peak volume time point of mobile water, E1 is assigned to C21, and E2 is assigned to C22.

If the fourth target map is in a three-peak state, determining a first peak time point corresponding to a peak value of a first peak of the fourth target map, marking as E1, determining a second peak time point corresponding to a peak value of a second peak of the fourth target map, marking as E2, determining a third peak time point corresponding to a peak value of a third peak of the fourth target map, marking as E3, wherein the first peak time point, the second peak time point and the third peak time point are respectively a peak volume time point of bound water, a peak volume time point of movable water and a peak volume time point of oil, assigning E1 to C21, E2 to C22, and E3 to C23.

A3: the computer device generates a second parameter matrix comprising a third bound water parameter, a third movable water parameter, a third oil parameter, a fourth bound water parameter, a fourth movable water parameter, and a fourth oil parameter.

The second parameter matrix is a 3-row 3-column matrix, the third bound water parameter, the fourth bound water parameter, the third movable water parameter, the fourth movable water parameter, the third oil parameter and the fourth oil parameter are respectively used as elements of a 1 st row, a 1 st column, a 2 nd row, a 2 nd column, a 2 nd row, a3 rd column, a 1 st row and a3 rd column in the second parameter matrix, and the elements of the 3 rd row, the 1 st column, the 3 rd row, the 2 nd column and the 3 rd row and the 3 rd column in the second parameter matrix are respectively assigned to be 10.

In the embodiment of the present application, since the first spectrum is a T2 spectrum of an actual underground stratum obtained through a P-type nmr logging process, the P-type nmr logging process may not measure partial peak data due to errors and other reasons in the measurement process, so that partial peaks are missing on the first spectrum, and the second spectrum is a T2 spectrum of a rock sample obtained by taking out the rock sample and performing a nmr experiment on the ground, the T2 spectrum includes T2 spectrums of a plurality of rock samples, peak data of the obtained T2 spectrum is complete, and further the missing peaks in the first spectrum are supplemented through the second spectrum, so that accuracy of a second parameter matrix determined based on the actual logging first spectrum can be ensured.

(2) The computer device determines, based on the first map and the second parameter matrix, relationship data for the target reservoir at depth, the relationship data being between the second parameter matrix and pore fluid saturation parameters.

Wherein the pore fluid comprises bound water, mobile water and oil, and the saturation parameters of the pore fluid comprise saturation of the bound water, saturation of the mobile water and saturation of the oil.

This step can be achieved by the following steps A1-A2:

a1: the computer device determines a first median time parameter and a second median time parameter based on the third target profile and the fourth target profile, respectively, and generates a dependent variable matrix comprising the first median time parameter and the second median time parameter.

The first median time parameter is a relaxation time point corresponding to 50% of the total porosity of the pore components in the third target spectrum, and the second median time parameter is a relaxation time point corresponding to 50% of the total porosity of the pore components in the fourth target spectrum. The first median time parameter is the element H11 of row 1 and column 1 of the dependent variable matrix, and the second median time parameter is the element H21 of row 2 and column 1 of the dependent variable matrix. The total porosity is obtained by integrating the pore components over the T2 spectrum.

The dependent variable matrix is a matrix of 3 rows and 1 columns, the first median time parameter and the second median time parameter are respectively used as elements of the 1 st row, the 1 st column and the 2 nd row and the 1 st column in the dependent variable matrix, and the elements of the 3 rd row and the 1 st column are assigned to be 10.

The computer equipment determines an integral value of each time point on the third target map through a target integral relation based on the third target map, determines a target integral value and a target time point corresponding to the target integral value from the integral values of each time point and each time point, and takes the target time point as a first element H11 of the dependent variable matrix; wherein the target integrated value is 0.5, and the target time point is a time point at which the integrated value is 0.5.

Wherein the target integral relation is an integral relation between the relaxation time and the pore component, and the target relation may be:

wherein F (T) is the integral value at time point T, T2_min is the minimum value of the relaxation time on the third target map, T2_max is the maximum value of the relaxation time on the third target map,

and phi is the total porosity of the third target map, which is the pore component at the time point t on the third target map.

The computer equipment determines an integral value of each time point on the fourth target map through a target integral relation based on the fourth target map, determines a target integral value and a target time point corresponding to the target integral value from the integral values of each time point and each time point, and takes the target time point as a second element H21 of the dependent variable matrix; wherein the target integrated value is 0.5, and the target time point is a time point at which the integrated value is 0.5.

wherein F (T) is an integral value at a time point of T, T2_min is a minimum value of the relaxation time on the fourth preset map, T2_max is a maximum value of the relaxation time on the fourth preset map,

And (2) the pore component of the t time point on the fourth preset map, wherein phi is the total porosity of the fourth preset map. />

A2: the computer device generates relational data comprising the second parameter matrix and the dependent variable matrix with the saturation of the tie water, the saturation of the movable water, and the saturation of the oil as the independent variable matrices.

Wherein the relationship data may be the following set of linear equations:

wherein,,

for the second parameter matrix->

Is an independent variable matrix>

Is a dependent variable matrix.

The saturation of the bound water, the saturation of the movable water and the saturation of the oil form an independent variable matrix, wherein the independent variable matrix is a matrix of 3 rows and 1 columns, and an element S1 of the 1 st row and the 1 st column in the independent variable matrix represents the saturation of the bound water; element S2 of row 2, column 1 in the argument matrix represents the saturation of movable water; element S3 in row 3, column 1 of the argument matrix represents the saturation of oil.

(3) The computer device determines a fluid saturation parameter of the target aperture based on the relationship data.

The second parameter matrix and the dependent variable matrix in the relation data are of known quantity, the independent variable matrix is of unknown quantity, the second parameter matrix and the dependent variable matrix are substituted into the relation data, the independent variable matrix is obtained through solving, and then the saturation of bound water, the saturation of movable water and the saturation of oil are obtained.

In the embodiment of the application, the fluid saturation parameters of each depth of the target reservoir are determined one by the methods from step 101 to step 106, so as to obtain the fluid saturation parameters of different depths of the whole target reservoir.

Referring to fig. 3, fig. 3 is a comparison graph of the results of fluid saturation parameters determined by the methods provided in the examples of the present application. The first trace in FIG. 3 is a lithology curve, where GR is a natural gamma curve in API, SP is a natural potential curve in mv; the second channel is a depth channel, and the depth unit is m; the third is a resistivity trace, wherein RLLD is a deep lateral resistivity curve in Ω.m, RLLS is a shallow lateral resistivity curve in Ω.m; the fourth is a porosity tract, wherein AC is an acoustic curve in μs/ft, CNL is a compensated neutron curve in% and DEN is a compensated density curve in g/cm3; the fifth lane is a standard T2 spectrum (T2 spectrum of the second echo interval time), and TASPEC is a T2 spectrum curve; the sixth track is a T2 spectrum trace of the first echo interval time, wherein TBSPEC is a T2 spectrum trace of the first echo interval time; the seventh lane is the irreducible water saturation lane, in units of; the eighth lane is the pore fluid saturation lane, where SW is free water saturation, SWI is irreducible water saturation, SO is oil saturation, in v/v. As can be seen from the graph, the result of the fluid saturation parameter determined by the embodiment of the application is matched with the logging response characteristics of each curve in the graph, which shows that the result of the fluid saturation parameter determined by the embodiment of the application has high accuracy, thereby being capable of effectively guiding the identification and evaluation of the reservoir oil content and providing guiding basis for the effective development of the reservoir

In one aspect, an embodiment of the present application provides a device for determining a fluid saturation parameter, referring to fig. 4, the device includes:

a first determining module 401 for determining a target reservoir of a region to be studied;

a second determining module 402, configured to determine, for each depth of the target reservoir, a first spectrum of the target reservoir at the depth based on one-dimensional nmr logging data of a target well section of a standard well in the region to be studied, where the one-dimensional nmr logging data is obtained by one-dimensional nmr logging engineering of the target well section, the target well section is a well section of the standard well corresponding to the depth of the target reservoir, and the first spectrum is a relationship curve between pore parameters and time;

The acquiring module 403 is configured to acquire nuclear magnetic resonance experimental data of a plurality of rock samples of a target reservoir, where the nuclear magnetic resonance experimental data is obtained through nuclear magnetic resonance experiments of the plurality of rock samples;

a third determining module 404, configured to determine, for each rock sample, a second spectrum of the rock sample based on nuclear magnetic resonance experimental data of the rock sample, where the second spectrum is a relationship between pore parameters and time;

a fourth determination module 405 for determining a first parameter matrix for the target reservoir at depth based on the second patterns of the plurality of rock samples;

a fifth determination module 406 is configured to determine a fluid saturation parameter of a target pore based on the first map and the first parameter matrix, the target pore being a pore of the target reservoir at the depth.

In one possible implementation, the nuclear magnetic resonance experimental data of the rock sample include saturated parameter water nuclear magnetic resonance experimental data, centrifugal nuclear magnetic resonance experimental data and saturated parameter oil nuclear magnetic resonance experimental data of the rock sample at different echo intervals;

the third determination module 404 includes:

the first determining unit is used for determining a first saturated parameter water map and a second saturated parameter water map of the rock sample based on saturated parameter water nuclear magnetic resonance experimental data of the rock sample at a first echo interval time and a second echo interval time respectively;

The second determining unit is used for determining a first centrifugal spectrum and a second centrifugal spectrum of the rock sample based on centrifugal nuclear magnetic resonance experimental data of the rock sample at the first echo interval time and the second echo interval time respectively;

the composition unit is used for forming a first saturated parameter water spectrum, a second saturated parameter water spectrum, a first centrifugal spectrum, a second centrifugal spectrum, a first saturated parameter oil spectrum and a second saturated parameter oil spectrum into a second spectrum.

In one possible implementation, the fourth determining module 405 includes:

an eighth determining unit for determining a first oil parameter and a second oil parameter based on the first saturation parameter oil map and the second saturation parameter oil map of the plurality of rock samples, respectively;

a generation unit for generating a first parameter matrix comprising a first bound water parameter, a second bound water parameter, a first movable water parameter, a second movable water parameter, a first oil parameter and a second oil parameter.

In one possible implementation manner, the sixth determining unit includes:

a first determining subunit, configured to determine, for each rock sample, a first difference spectrum of the rock sample based on a first saturation parameter water spectrum and a first centrifugal spectrum of the rock sample, where the first difference spectrum is a spectrum composed of differences between first saturation parameter water spectrum and first saturation parameter and second centrifugal parameter at the same time in the first saturation parameter water spectrum and the first centrifugal spectrum;

a second determining subunit, configured to determine a third time point corresponding to the maximum peak of the first difference map;

In one possible implementation manner, the seventh determining unit includes:

a fourth determination subunit, configured to determine, for each rock sample, a second difference spectrum of the rock sample based on a second saturation parameter water spectrum and a second centrifugal spectrum of the rock sample, where the second difference spectrum is a spectrum formed by differences between second saturation parameter water spectrum and second centrifugal parameter at the same time in the second saturation parameter water spectrum and the second centrifugal spectrum;

a sixth determination subunit for determining a second movable water parameter based on an average of fourth time points of the plurality of rock samples.

In one possible implementation, the fifth determining module 406 includes:

a tenth determining unit, configured to determine, based on the first map and the second parameter matrix, relationship data of the target reservoir at depth, where the relationship data is relationship data between the second parameter matrix and the pore fluid volume;

In one possible implementation, the first pattern includes a third target pattern and a fourth target pattern, the third target pattern and the fourth target pattern being patterns of the formation at a first echo interval and a second echo interval, respectively;

a ninth determination unit including:

a first generation subunit for generating a second parameter matrix comprising a third bound water parameter, a third movable water parameter, a third oil parameter, a fourth bound water parameter, a fourth movable water parameter, and a fourth oil parameter.

a tenth determination unit including:

a second generation subunit for generating a dependent variable matrix comprising a first median time parameter and a second median time parameter;

a third generation subunit for generating relationship data including a second parameter matrix and a dependent variable matrix, and using the saturation of the bound water, the saturation of the movable water and the saturation of the oil as the independent variable matrix

The embodiment of the application provides a computer device, which comprises a processor and a memory, wherein at least one program code is stored in the memory, and the at least one program code is loaded and executed by the processor to realize the instructions of the method for determining the fluid saturation parameter in any implementation manner.

Fig. 5 shows a block diagram of a computer device 500 provided in an exemplary embodiment of the present application. The computer device 500 may be a portable mobile computer device such as: a smart phone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III, motion picture expert compression standard audio plane 3), an MP4 (Moving Picture Experts Group Audio Layer IV, motion picture expert compression standard audio plane 4) player, a notebook computer, or a desktop computer. Computer device 500 may also be referred to by other names of user devices, portable computer devices, laptop computer devices, desktop computer devices, and the like.

In general, the computer device 500 includes: a processor 501 and a memory 502.

Processor 501 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and the like. The processor 501 may be implemented in at least one hardware form of DSP (Digital Signal Processing ), FPGA (Field-Programmable Gate Array, field programmable gate array), PLA (Programmable Logic Array ). The processor 501 may also include a main processor and a coprocessor, the main processor being a processor for processing data in an awake state, also referred to as a CPU (Central Processing Unit ); a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 501 may be integrated with a GPU (Graphics Processing Unit, image processor) for taking care of rendering and rendering of content that the display screen is required to display. In some embodiments, the processor 501 may also include an AI (Artificial Intelligence ) processor for processing computing operations related to machine learning.

Memory 502 may include one or more computer-readable storage media, which may be non-transitory. Memory 502 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 502 is used to store at least one instruction for execution by processor 501 to implement the method of determining a fluid saturation parameter provided by a method embodiment in the present application.

In some embodiments, the computer device 500 may further optionally include: a peripheral interface 503 and at least one peripheral. The processor 501, memory 502, and peripheral interface 503 may be connected by buses or signal lines. The individual peripheral devices may be connected to the peripheral device interface 503 by buses, signal lines or circuit boards. Specifically, the peripheral device includes: at least one of radio frequency circuitry 504, a display 505, a camera assembly 506, audio circuitry 507, a positioning assembly 508, and a power supply 509.

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

The Radio Frequency circuit 504 is configured to receive and transmit RF (Radio Frequency) signals, also known as electromagnetic signals. The radio frequency circuitry 504 communicates with a communication network and other communication devices via electromagnetic signals. The radio frequency circuit 504 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 504 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 504 may communicate with other computer devices 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 circuitry 504 may also include NFC (Near Field Communication ) related circuitry, which is not limited in this application.

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

The camera assembly 506 is used to capture images or video. Optionally, the camera assembly 506 includes a front camera and a rear camera. Typically, the front camera is disposed on a front panel of the computer device and the rear camera is disposed on a rear surface of the computer device. In some embodiments, the at least two rear cameras are any one of a main camera, a depth camera, a wide-angle camera and a tele camera, so as to realize that the main camera and the depth camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize a panoramic shooting and Virtual Reality (VR) shooting function or other fusion shooting functions. In some embodiments, camera assembly 506 may also include a flash. The flash lamp can be a single-color temperature flash lamp or a double-color temperature flash lamp. The dual-color temperature flash lamp refers to a combination of a warm light flash lamp and a cold light flash lamp, and can be used for light compensation under different color temperatures.

The audio circuitry 507 may include a microphone and a speaker. The microphone is used for collecting sound waves of users and environments, converting the sound waves into electric signals, and inputting the electric signals to the processor 501 for processing, or inputting the electric signals to the radio frequency circuit 504 for voice communication. The microphone may be provided in a plurality of different locations of the computer device 500 for stereo acquisition or noise reduction purposes. The microphone may also be an array microphone or an omni-directional pickup microphone. The speaker is used to convert electrical signals from the processor 501 or the radio frequency circuit 504 into sound waves. The speaker may be a conventional thin film speaker or a piezoelectric ceramic speaker. When the speaker is a piezoelectric ceramic speaker, not only the electric signal can be converted into a sound wave audible to humans, but also the electric signal can be converted into a sound wave inaudible to humans for ranging and other purposes. In some embodiments, audio circuitry 507 may also include a headphone jack.

The location component 508 is used to locate the current geographic location of the computer device 500 to enable navigation or LBS (Location Based Service, location-based services). The positioning component 508 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 509 is used to power the various components in the computer device 500. The power supply 509 may be an alternating current, a direct current, a disposable battery, or a rechargeable battery. When the power supply 509 comprises 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.

In some embodiments, the computer device 500 further includes one or more sensors 510. The one or more sensors 510 include, but are not limited to: an acceleration sensor 511, a gyro sensor 512, a pressure sensor 513, a fingerprint sensor 514, an optical sensor 515, and a proximity sensor 516.

The acceleration sensor 511 can detect the magnitudes of accelerations on three coordinate axes of the coordinate system established with the computer device 500. For example, the acceleration sensor 511 may be used to detect components of gravitational acceleration on three coordinate axes. The processor 501 may control the display 505 to display a user interface in a landscape view or a portrait view according to a gravitational acceleration signal acquired by the acceleration sensor 511. The acceleration sensor 511 may also be used for acquisition of motion data of a game or a user.

The gyro sensor 512 may detect a body direction and a rotation angle of the computer device 500, and the gyro sensor 512 may collect a 3D motion of the user on the computer device 500 in cooperation with the acceleration sensor 511. The processor 501 may implement the following functions based on the data collected by the gyro sensor 512: motion sensing (e.g., changing UI according to a tilting operation by a user), image stabilization at shooting, game control, and inertial navigation.

The pressure sensor 513 may be disposed on a side frame of the computer device 500 and/or on an underlying layer of the display 505. When the pressure sensor 513 is disposed on the side frame of the computer device 500, a grip signal of the computer device 500 by a user may be detected, and the processor 501 performs left-right hand recognition or quick operation according to the grip signal collected by the pressure sensor 513. When the pressure sensor 513 is disposed at the lower layer of the display screen 505, the processor 501 controls the operability control on the UI interface according to the pressure operation of the user on the display screen 505. The operability controls include at least one of a button control, a scroll bar control, an icon control, and a menu control.

The fingerprint sensor 514 is used for collecting the fingerprint of the user, and the processor 501 identifies the identity of the user according to the fingerprint collected by the fingerprint sensor 514, or the fingerprint sensor 514 identifies the identity of the user according to the collected fingerprint. Upon recognizing that the user's identity is a trusted identity, the user is authorized by the processor 501 to perform relevant sensitive operations including unlocking the screen, viewing encrypted information, downloading software, paying for and changing settings, etc. The fingerprint sensor 514 may be disposed on the front, back, or side of the computer device 500. When a physical key or vendor Logo is provided on the computer device 500, the fingerprint sensor 514 may be integrated with the physical key or vendor Logo.

The optical sensor 515 is used to collect the ambient light intensity. In one embodiment, the processor 501 may control the display brightness of the display screen 505 based on the intensity of ambient light collected by the optical sensor 515. Specifically, when the intensity of the ambient light is high, the display brightness of the display screen 505 is turned up; when the ambient light intensity is low, the display brightness of the display screen 505 is turned down. In another embodiment, the processor 501 may also dynamically adjust the shooting parameters of the camera assembly 506 based on the ambient light intensity collected by the optical sensor 515.

A proximity sensor 516, also referred to as a distance sensor, is typically provided on the front panel of the computer device 500. The proximity sensor 516 is used to collect the distance between the user and the front of the computer device 500. In one embodiment, when the proximity sensor 516 detects a gradual decrease in the distance between the user and the front of the computer device 500, the processor 501 controls the display 505 to switch from the bright screen state to the off screen state; when the proximity sensor 516 detects that the distance between the user and the front of the computer device 500 gradually increases, the display 505 is controlled by the processor 501 to switch from the off-screen state to the on-screen state.

Those skilled in the art will appreciate that the architecture shown in fig. 5 is not limiting as to the computer device 500, and may include more or fewer components than shown, or may combine certain components, or employ a different arrangement of components.

Embodiments of the present application provide a computer readable storage medium having at least one program code stored therein, the at least one program code being loaded and executed by a processor to implement the steps in the method for determining a fluid saturation parameter of any of the above implementations.

The foregoing description of the preferred embodiments is merely exemplary in nature and is in no way intended to limit the invention, since it is intended that all modifications, equivalents, improvements, etc. that fall within the spirit and scope of the invention.

Claims

1. A method of determining a fluid saturation parameter, the method comprising:

determining a target reservoir of a region to be studied;

2. The method of determining fluid saturation parameters according to claim 1, wherein the nuclear magnetic resonance experimental data of the rock sample includes saturated parameter water nuclear magnetic resonance experimental data, centrifugal nuclear magnetic resonance experimental data and saturated parameter oil nuclear magnetic resonance experimental data of the rock sample at different echo intervals;

3. The method of determining fluid saturation parameters of claim 2, wherein the determining a first parameter matrix for the target reservoir at the depth based on the second patterns of the plurality of rock samples comprises:

4. A method of determining a fluid saturation parameter according to claim 3, wherein the determining a first movable water parameter based on a first saturation parameter water profile of the plurality of rock samples and the first centrifugal profile comprises:

5. A method of determining a fluid saturation parameter according to claim 3, wherein the determining a second movable water parameter based on a second saturation parameter water profile and a second centrifugal profile of the plurality of rock samples comprises:

6. The method of determining a fluid saturation parameter according to claim 1, wherein the determining a fluid saturation parameter of a target aperture based on the first map and the first parameter matrix comprises:

7. The method of determining a fluid saturation parameter of claim 6, wherein the first map includes a third target map and a fourth target map, the third target map and the fourth target map being maps of the formation at a first echo interval time and a second echo interval time, respectively;

8. The method of determining a fluid saturation parameter of claim 7, wherein the pore fluid includes bound water, mobile water, and oil, and the relationship data includes a second parameter matrix, a dependent variable matrix, and an independent variable matrix;

The relationship data including the second parameter matrix and the dependent variable matrix and having the saturation of the tie water, the saturation of the movable water, and the saturation of the oil as independent variable matrices is generated.

9. A device for determining a fluid saturation parameter, the device comprising:

10. A computer device comprising a processor and a memory, the memory storing at least one program code, the at least one program code loaded and executed by the processor to implement instructions of a method of determining a fluid saturation parameter according to any one of claims 1 to 8.