CN106324008A - Correction method of nuclear magnetic resonance logging wave spectrum - Google Patents

Abstract

The invention provides a correction method of nuclear magnetic resonance logging wave spectrum. The correction method comprises the following steps: a step S1, separating the nuclear magnetic resonance wave spectrum in a neutral reservoir into a water-containing spectrum and an oil-containing spectrum by using a nuclear magnetic resonance wave spectrum separating model; a step S2, calculating a water surface relaxation rate and an oil surface relaxation rate of an aperture in the neutral reservoir by using a surface relaxation rate calculating model; a step S3, calculating the pore throat radius distribution of an oil-containing part in the neutral reservoir by using the oil surface relaxation rate; a step S4, correcting the pore throat radius distribution of the oil-containing part according to the water surface relaxation rate, in order to obtain the corrected water saturated nuclear magnetic resonance wave spectrum of the oil-containing part in the neutral reservoir; a step S5, carrying out superposition for the oil-containing spectrum and the corrected water saturated nuclear magnetic resonance wave spectrum of the oil-containing part, in order to obtain the corrected water saturated nuclear magnetic resonance wave spectrum of the neutral reservoir; the accurate water saturated nuclear magnetic resonance wave spectrum of the neutral reservoir is obtained by using the correction method, and accurate analysis is carried out for the neutral reservoir.

Description

The bearing calibration of nuclear magnetic resonance log wave spectrum

Technical field

The present invention relates to oil field, in particular to the bearing calibration of a kind of nuclear magnetic resonance log wave spectrum.

Background technology

Petroleum industry becomes increasingly complex through the development of over half a century, the target of its exploration and development, and the permeability of its reservoir is more come The lowest, current hyposmosis and densification oil reservoir have been will to be maybe the major heading of exploration and development, the core of such evaluating reservoir One of content is exactly pore structure evaluation, and existing pore structure assessment technique mainly applies nuclear magnetic resonance log, works as pore media When the NMR relaxation mechanism of middle fluid is Volume Relaxation, the Response Mechanism of NMR (Nuclear Magnetic Resonance) spectrum is as follows:

$T_{1,2}={\mathrm{\ρ}}_{1,2}\×{\left(\frac{S}{V}\right)}_{por}$

Wherein, T_1,2For NMR relaxation time, ρ_1,2For the surface relaxivity of rock, S is the surface area of hole, and V is hole The volume of gap.This model is the evaluation model of most physical significance, but applying the premise of this model is that reservoir is containing only single hole Clearance flow body.Well log interpretation expert in countries in the world proposes the evaluation model of some series subsequently, and the premise of these models application is all For hydrophilic reservoir, and nuclear magnetic resonance log is affected little by oily.

Current most conventional reservoir is that the reservoir in hydrophilic rock, and nuclear magnetic resonance log investigative range is because of by slurry compounding Impact, its residual oil saturation is low, and now nuclear magnetic resonance log is affected the faintest by oiliness, and its response characteristic is single The nuclear magnetic resonance, NMR of one fluid, i.e. water, the most above-mentioned model all can be used to evaluate the pore throat radius distribution of reservoir.

But, for the reservoir of neutral wetting, such as fine and close oil reservoir, its oil saturation is higher, and pore throat radius is little, Mud is difficult to reservoir is carried out displacement, and nuclear magnetic resonance log is big by the effect of signals of crude oil, at this moment uses existing method to analyze The pore throat radius distribution calculating fine and close oil reservoir has the great error of generation.

Analysis method of the prior art does not fully take into account the neutral reservoir oiliness impact on nuclear magnetic resonance log, if not The nuclear magnetic resonance log wave spectrum of reservoir is carried out oiliness correction and is just difficult to correctly characterize the pore throat radius distribution of neutral reservoir, therefore, Need a kind of method that nuclear magnetic resonance log wave spectrum oil gas corrects.

Summary of the invention

Present invention is primarily targeted at the bearing calibration that a kind of nuclear magnetic resonance log wave spectrum is provided, to solve water of the prior art The NMR (Nuclear Magnetic Resonance) spectrum problem that can not accurately analyze the distribution of neutral reservoir pore throat radius.

To achieve these goals, according to an aspect of the invention, it is provided the bearing calibration of a kind of nuclear magnetic resonance log wave spectrum, Above-mentioned bearing calibration includes: step S1, uses NMR (Nuclear Magnetic Resonance) spectrum disjunctive model the NMR (Nuclear Magnetic Resonance) spectrum of neutral reservoir to be divided into and contains Water spectrum is composed with oil-containing；Step S2, uses the water surface relaxation rate of the hole of the above-mentioned neutral reservoir of surface relaxivity computation model calculating With oil surface relaxivity；Step S3, utilizes the pore throat radius of the oil-containing part of the above-mentioned above-mentioned neutral reservoir of oil surface relaxivity calculating Distribution；Step S4, is corrected the pore throat radius distribution of above-mentioned oil-containing part according to above-mentioned water surface relaxation rate, is corrected After the NMR (Nuclear Magnetic Resonance) spectrum being full of water of oil-containing part of above-mentioned neutral reservoir；And step S5, to above-mentioned aqueous spectrum with above-mentioned The NMR (Nuclear Magnetic Resonance) spectrum being full of water of oil-containing part of correction is overlapped, above-mentioned neutrality reservoir after correct be full of water NMR (Nuclear Magnetic Resonance) spectrum.

Further, above-mentioned NMR (Nuclear Magnetic Resonance) spectrum disjunctive model is: T_1w,2w∪T_1o,2o=T_1,2(1), in formula (1): T_1,2For The NMR (Nuclear Magnetic Resonance) spectrum of above-mentioned neutral reservoir；T_1w,2wAqueous spectrum for above-mentioned neutral reservoir；T_1o,2oOil-containing for above-mentioned neutral reservoir Spectrum.

Further, above-mentioned steps S2 includes:

Step A, chooses sample of coring, and is full of the nuclear magnetic resonance, NMR of water after measuring the NMR (Nuclear Magnetic Resonance) spectrum of above-mentioned sample of coring, washing oil Wave spectrum is distributed with total pore throat radius of above-mentioned neutral reservoir；

Step B, according to the NMR (Nuclear Magnetic Resonance) spectrum being full of water after above-mentioned washing oil and capillary pressure curve, uses in formula (2) calculating Stating water surface relaxation rate, above-mentioned formula (2) is:In above-mentioned formula (2): ρ_wFor above-mentioned Water surface relaxation rate, θ is the contact angle of hydrargyrum, and γ is the surface tension of hydrargyrum, P_totalFor capillary pressure, T_1tw,2twAfter above-mentioned washing oil The nuclear magnetic resonance, NMR T2 T being full of water_2twOr longitudinal relaxation time T_1tw, unit is ms；

Step C, is divided into the aqueous spectrum of sample of coring and sample oil-containing of coring by the NMR (Nuclear Magnetic Resonance) spectrum being full of profit of above-mentioned sample of coring Spectrum；

Step D, according to above-mentioned water surface relaxation rate and the above-mentioned aqueous spectrum of sample of coring, uses formula (3) to calculate above-mentioned sample of coring The pore throat radius containing water section distribution, above-mentioned formula (3) is: r_w=a × ρ_w×T_1w,2wIn above-mentioned formula (3): r_wFor upper State the distribution of the pore throat radius containing water section of sample of coring, ρ_wFor above-mentioned water surface relaxation rate, a is constant, when above-mentioned pore shape During the most cylindrical, spherical, shape channel-shaped, a is respectively 2,3,1, T_1w,2wTransverse relaxation for the above-mentioned aqueous spectrum of sample of coring Time T_2wOr longitudinal relaxation time T_1w；

Step E, according to the distribution of the pore throat radius containing water section of the distribution of above-mentioned total pore throat radius with above-mentioned sample of coring, uses formula (4) Calculating the pore throat radius distribution of the oil-containing part of above-mentioned sample of coring, above-mentioned formula (4) isIn above-mentioned formula (4), r_oFor the pore throat radius distribution of the oil-containing part of above-mentioned sample of coring, r_wThe pore throat radius containing water section for above-mentioned sample of coring divides Cloth, r_totalIt is distributed for above-mentioned total pore throat radius；

Step F, composes with the oil-containing of above-mentioned sample of coring according to the pore throat radius distribution of oil-containing part of above-mentioned sample of coring, uses formula (5) calculating above-mentioned oil surface relaxivity, above-mentioned formula (5) isIn above-mentioned formula (5), ρ_oFor upper State oil surface relaxivity；θ is the contact angle of above-mentioned hydrargyrum, and γ is the surface tension of above-mentioned hydrargyrum, and a is constant, when above-mentioned pore shape During the most cylindrical, spherical, shape channel-shaped, a is respectively 2,3,1, P_oAperture for above-mentioned neutral reservoir oil-containing part is corresponding Capillary pressure, be calculated according to the pore throat radius distribution of the oil-containing part of above-mentioned sample of coring, T_1o,2oCompose for above-mentioned oil-containing T2 T_2oOr longitudinal relaxation time T_1o。

Further, above-mentioned steps A includes: choose sealing core drilling sample, measures the nuclear magnetic resonance, NMR ripple of above-mentioned sealing core drilling sample Spectrum；Above-mentioned sealing core drilling sample is carried out washing oil, measures the NMR (Nuclear Magnetic Resonance) spectrum being full of water of the above-mentioned sample of coring after washing oil； And carry out pressing hydrargyrum experiment to above-mentioned sealing core drilling sample, it is thus achieved that above-mentioned total pore throat radius is distributed.

Further, above-mentioned steps A includes: choose non-sealing core drilling sample, and above-mentioned non-sealing core drilling sample is carried out washing oil, Measure the NMR (Nuclear Magnetic Resonance) spectrum being full of water of the above-mentioned non-sealing core drilling sample after washing oil；Above-mentioned non-airtight after above-mentioned washing oil is taken Heart sample carries out pressing hydrargyrum experiment, it is thus achieved that above-mentioned total pore throat radius is distributed；And above-mentioned non-sealing core drilling sample is carried out core Location, Using practical logging NMR (Nuclear Magnetic Resonance) spectrum as the NMR (Nuclear Magnetic Resonance) spectrum of above-mentioned non-sealing core drilling sample.

Further, above-mentioned core Location is the core Location of Centimeter Level.

Further, in above-mentioned steps C, based on crude is composed and is left the nuclear magnetic resonance, NMR time corresponding to limit value, cores above-mentioned The NMR (Nuclear Magnetic Resonance) spectrum of sample is divided into the aqueous spectrum of sample of coring and sample oil-containing spectrum of coring.

Further, above-mentioned steps S3 includes, according to above-mentioned oil surface relaxivity and the above-mentioned oil-containing spectrum in above-mentioned steps S1, adopting Calculate the pore throat radius distribution of above-mentioned neutral reservoir oil-containing part by formula (6), above-mentioned formula (6) is: R_o=a × ρ_o×T_1o,2o, In above-mentioned formula (6), R_oPore throat radius distribution for above-mentioned neutral reservoir oil-containing part；ρ_oFor above-mentioned oil surface relaxivity；A is Constant, when above-mentioned pore shape is for taking 2 time cylindrical, for taking 3 time spherical, takes 1 during for channel-shaped；T_1o,2oCompose for above-mentioned oil-containing T2 T_2oOr longitudinal relaxation time T_1o。

Further, in above-mentioned steps S4, according to above-mentioned water surface relaxation rate, use formula (7) that above-mentioned neutral reservoir is contained The pore throat radius distribution of oil part is corrected, and obtains the NMR (Nuclear Magnetic Resonance) spectrum being full of water of above-mentioned neutral reservoir oil-containing part, on Stating formula (7) is,In above-mentioned formula (7), T_1jw,2jwThe water that is full of for above-mentioned neutral reservoir oil-containing part NMR (Nuclear Magnetic Resonance) spectrum；R_oPore throat radius distribution for the oil-containing part of above-mentioned sample of coring；ρ_wFor above-mentioned water surface relaxation rate；a For constant, when above-mentioned pore shape is the most cylindrical, spherical, shape channel-shaped time, a is respectively 2,3,1.

Further, in above-mentioned steps S5, use formula (8) by the above-mentioned aqueous spectrum in above-mentioned steps S1 and above-mentioned neutral storage The NMR (Nuclear Magnetic Resonance) spectrum being full of water of layer oil-containing part is overlapped, and obtains above-mentioned neutral reservoir and is full of the NMR (Nuclear Magnetic Resonance) spectrum of water, Above-mentioned formula (8) is, T_1tw,2tw=T_1w,2w∪T_1jw,2jw, in above-mentioned formula (8), T_1tw,2twIt is full of water for above-mentioned neutral reservoir NMR (Nuclear Magnetic Resonance) spectrum；T_1w,2wFor above-mentioned aqueous spectrum；T_1jw,2jwThe nuclear magnetic resonance, NMR ripple being full of water for above-mentioned neutral reservoir oil-containing part Spectrum.

The bearing calibration of the application present invention, has taken into full account the oiliness impact on nuclear magnetic resonance log of neutral reservoir, in employing The hole of the oil-containing part of property reservoir is full of water NMR (Nuclear Magnetic Resonance) spectrum, and the aqueous spectrum of centering reservoir is corrected, the ripple after correction Spectrum is the NMR (Nuclear Magnetic Resonance) spectrum being full of water of neutral reservoir, and this wave spectrum can correctly characterize the pore throat radius distribution of neutral reservoir, make Obtain those skilled in the art can accurately analyze by centering reservoir, and then carry out rational exploration and development, beneficially oil work The development of industry.

Accompanying drawing explanation

The Figure of description of the part constituting the application is used for providing a further understanding of the present invention, and the present invention's is schematic real Execute example and illustrate for explaining the present invention, being not intended that inappropriate limitation of the present invention.In the accompanying drawings:

Fig. 1 shows the schematic flow sheet of the bearing calibration of the nuclear magnetic resonance log wave spectrum that a kind of exemplary embodiment provides；

Fig. 2 shows the nuclear magnetic resonance, NMR T being full of water-oil phase of a kind of preferred embodiment reservoir₂The core of water it is full of after wave spectrum and washing oil Magnetic Resonance Spectrum comparison diagram；

Fig. 3 shows that a bite well sealed coring well section that a kind of preferred embodiment of the application provides applies different AT value iterative computation Analysis of Mean Square Error figure；

Fig. 4 a shows that the NMR (Nuclear Magnetic Resonance) spectrum being full of water-oil phase that a kind of preferred embodiment provides is separated into water spectrum and showing that oil is composed It is intended to；

Fig. 4 b shows the schematic diagram from the isolated aqueous spectrum of the NMR (Nuclear Magnetic Resonance) spectrum of Fig. 4 a；

Fig. 4 c shows the schematic diagram composed from the isolated oil-containing of NMR (Nuclear Magnetic Resonance) spectrum of Fig. 4 a；

Fig. 5 a shows the capillary pressure curve that a kind of preferred embodiment provides；

Fig. 5 b is full of the NMR (Nuclear Magnetic Resonance) spectrum capillary pressure curve of water after showing the washing oil that a kind of preferred embodiment provides；

Fig. 5 c is full of the NMR (Nuclear Magnetic Resonance) spectrum capillary pressure curve of water after showing the washing oil that a kind of preferred embodiment provides；

Fig. 6 a shows the aqueous spectrum of the sample of coring that a kind of preferred embodiment provides；

Fig. 6 b shows the distribution of the pore throat radius containing water section of the sample of coring that a kind of preferred embodiment provides；

Fig. 7 a shows total pore throat radius distribution schematic diagram that a kind of preferred embodiment provides；

Fig. 7 b shows that a kind of preferred embodiment provides aqueous pore throat radius distribution schematic diagram；

The pore throat radius distribution of the oil-containing part drawn according to Fig. 7 a and Fig. 7 b that Fig. 7 c shows that a kind of preferred embodiment provides is shown It is intended to；

Fig. 8 shows the test result figure that a kind of preferred embodiment of the application provides；And

Fig. 9 illustrates the 5th some depth point, region in Fig. 8, utilizes and obtains what complete aqueous spectrum calculated through oiliness correction The pore throat radius distribution curve comparison diagram that throat distribution feature obtains with experiment.

Detailed description of the invention

It should be noted that in the case of not conflicting, the embodiment in the application and the feature in embodiment can be mutually combined. Describe the present invention below with reference to the accompanying drawings and in conjunction with the embodiments in detail.

This application provides the bearing calibration of a kind of nuclear magnetic resonance log wave spectrum, as it is shown in figure 1, this bearing calibration includes: step S1, uses NMR (Nuclear Magnetic Resonance) spectrum disjunctive model to be divided into aqueous spectrum to compose with oil-containing the NMR (Nuclear Magnetic Resonance) spectrum of neutral reservoir；Step S2, Use water surface relaxation rate and the oil surface relaxivity of the hole of the above-mentioned neutral reservoir of surface relaxivity computation model calculating；Step S3, Utilize the pore throat radius distribution of the oil-containing part of the above-mentioned above-mentioned neutral reservoir of oil surface relaxivity calculating；Step S4, according to above-mentioned water The pore throat radius of above-mentioned oil-containing part is distributed and is corrected by surface relaxivity, the oil-containing portion of the above-mentioned neutral reservoir after being corrected The NMR (Nuclear Magnetic Resonance) spectrum being full of water divided；And step S5, the water that is full of to the oil-containing part of above-mentioned aqueous spectrum and above-mentioned correction NMR (Nuclear Magnetic Resonance) spectrum is overlapped, the NMR (Nuclear Magnetic Resonance) spectrum being full of water of the above-mentioned neutral reservoir after being corrected.Above-mentioned steps S2 Middle is the most all the core sample for the section of coring, and above-mentioned steps S3, S4 Yu S5 are to survey for the non-section of coring down-hole is actual The NMR (Nuclear Magnetic Resonance) spectrum of amount.

Above-mentioned bearing calibration has taken into full account the oiliness impact on nuclear magnetic resonance log of neutral reservoir, uses neutral reservoir The hole of oil-containing part is full of water NMR (Nuclear Magnetic Resonance) spectrum, and the aqueous spectrum of centering reservoir is corrected, and the wave spectrum after correction is neutral The NMR (Nuclear Magnetic Resonance) spectrum being full of water of reservoir, this wave spectrum can correctly characterize the pore throat radius distribution of neutral reservoir so that this area Technical staff can accurately analyze by centering reservoir, and then carries out the development of rational exploration and development, beneficially petroleum industry.

In order to be divided into aqueous spectrum to compose with oil-containing the NMR (Nuclear Magnetic Resonance) spectrum of neutral reservoir more accurately, the preferred above-mentioned core of the application Magnetic Resonance Spectrum disjunctive model is:

$T_{1w,2w}\∪T_{1o,2o}=T_{1,2}\mathrm{...}\left(1\right)$

In formula (1), T_1,2NMR (Nuclear Magnetic Resonance) spectrum for above-mentioned neutral reservoir；T_1w,2wAqueous spectrum for above-mentioned neutral reservoir；T_1o,2o Oil-containing spectrum for above-mentioned neutral reservoir.

NMR (Nuclear Magnetic Resonance) spectrum can be longitudinal relaxation time (T₁) spectrum, it is also possible to it is T2 (T₂) spectrum.For profit The disjunctive model of spectrum is through finding a NMR relaxation time point (T₁With T₂), this time point correspondence crude oil is composed The lower limit deposited.The wave spectrum on this nuclear magnetic resonance, NMR time point left side (less than this time) is aqueous spectrum, and rock is the most hydrophilic Rock；The wave spectrum of this nuclear magnetic resonance, NMR time point (more than this time) is oil-containing spectrum, and rock is oil-wet rock, on the whole rock Stone shows as neutral wetting.This NMR relaxation time point can apply profit characteristic nmr difference to calculate, it is possible to The NMR relaxation time value that when being applied to hide, reservoir crude oil minimum fills radius corresponding determines, it is also possible to combined by two kinds of methods Determine.

In another preferred embodiment of the application, above-mentioned steps S2 includes:

Step A, chooses sample of coring, and is full of the nuclear magnetic resonance, NMR of water after measuring the NMR (Nuclear Magnetic Resonance) spectrum of above-mentioned sample of coring, washing oil Wave spectrum is distributed with total pore throat radius of above-mentioned neutral reservoir；

Step B, according to the NMR (Nuclear Magnetic Resonance) spectrum being full of water after above-mentioned washing oil and capillary pressure curve, uses in formula (2) calculating Stating water surface relaxation rate, above-mentioned formula (2) is:

${\mathrm{\ρ}}_w=\frac{2\×\mathrm{\γ}\×cos\mathrm{\θ}}{a\×T_{1tw,2tw}\×P_{total}}\mathrm{...}\left(2\right)$

In formula (2), ρ_wFor above-mentioned water surface relaxation rate, unit is 10^-3μm/ms, θ are the contact angle of hydrargyrum, and γ is the surface tension of hydrargyrum, Unit is psi cm, P_totalFor capillary pressure, unit is MPa, T_1tw,2twHorizontal for the nuclear magnetic resonance, NMR being full of water after above-mentioned washing oil To relaxation time (T_2tw) or longitudinal relaxation time (T_1tw), unit is ms；

Step C, is divided into the aqueous spectrum of sample of coring and sample oil-containing of coring by the NMR (Nuclear Magnetic Resonance) spectrum being full of profit of above-mentioned sample of coring Spectrum；

Step D, according to above-mentioned water surface relaxation rate and the above-mentioned aqueous spectrum of sample of coring, uses formula (3) to calculate above-mentioned sample of coring The pore throat radius containing water section distribution, above-mentioned formula (3) is:

r_w=a × ρ_w×T_1w,2w……………….(3)

In above-mentioned formula (3), r_wThe pore throat radius containing water section for above-mentioned sample of coring is distributed, and unit is μm, ρ_wFor above-mentioned water meter Face relaxation rate, unit is 10^-3μm/ms, a are constant, when above-mentioned pore shape is the most cylindrical, spherical, shape channel-shaped time, a It is respectively 2,3,1, T_1w,2wT2 (T for the above-mentioned aqueous spectrum of sample of coring_2w) or longitudinal relaxation time (T_1w), Unit is ms；

Step E, according to the distribution of the pore throat radius containing water section of the distribution of above-mentioned total pore throat radius with above-mentioned sample of coring, uses formula (4) Calculating the pore throat radius distribution of the oil-containing part of above-mentioned sample of coring, above-mentioned formula (4) is:

$r_o=\stackrel{\‾}{r_w\∩r_{total}}\mathrm{...}\left(4\right)$

In above-mentioned formula (4): r_oFor the pore throat radius distribution of the oil-containing part of above-mentioned sample of coring, r_wAqueous portion for above-mentioned sample of coring The pore throat radius divided divides, r_totalIt is distributed for above-mentioned total pore throat radius；

Step F, composes with the oil-containing of above-mentioned sample of coring according to the pore throat radius distribution of oil-containing part of above-mentioned sample of coring, uses formula (5) calculating above-mentioned oil surface relaxivity, computing formula is:

${\mathrm{\ρ}}_o=\frac{2\×\mathrm{\γ}\×cos\mathrm{\θ}}{a\×T_{1o,2o}\×P_o}\mathrm{...}\left(5\right)$

In above-mentioned formula (5): ρ_oFor above-mentioned oil surface relaxivity, unit is 10^-3μm/ms；θ is the contact angle of above-mentioned hydrargyrum, unit For °, γ is the surface tension of above-mentioned hydrargyrum, and unit is psi cm, and a is constant, when above-mentioned pore shape is respectively cylindrical, ball When shape, shape channel-shaped, a is respectively 2,3,1, P_oFor the capillary pressure that the aperture of above-mentioned neutral reservoir oil-containing part is corresponding, according to The pore throat radius distribution of the oil-containing part of above-mentioned sample of coring is calculated, and unit is MPa, T_1o,2oFor above-mentioned oil-containing compose horizontal Relaxation time (T_2o) or longitudinal relaxation time (T_1o), unit is ms.

Adopt with the aforedescribed process, it is possible to the water surface relaxation rate of the hole of the neutral reservoir calculated more accurately and oil surface relaxation Rate, further ensures the accuracy of this bearing calibration.

In order to be full of after obtaining the NMR (Nuclear Magnetic Resonance) spectrum of the most above-mentioned sample of coring, washing oil further the NMR (Nuclear Magnetic Resonance) spectrum of water with Total pore throat radius distribution of above-mentioned neutral reservoir, those skilled in the art can choose sealing core drilling sample with non-according to practical situation Sealing core drilling sample is tested.

When sample is sealing core drilling sample, above-mentioned steps A includes: choose sealing core drilling sample, measures above-mentioned sealing core drilling sample The NMR (Nuclear Magnetic Resonance) spectrum of product；Above-mentioned sealing core drilling sample is carried out washing oil, measures the water that is full of of above-mentioned sample of coring after washing oil NMR (Nuclear Magnetic Resonance) spectrum；And carry out pressing hydrargyrum experiment to above-mentioned sealing core drilling sample, it is thus achieved that above-mentioned total pore throat radius is distributed.

When sample is non-sealing core drilling sample, above-mentioned steps A includes: choose non-sealing core drilling sample, non-airtight takes above-mentioned Heart sample carries out washing oil, measures the NMR (Nuclear Magnetic Resonance) spectrum being full of water of the above-mentioned non-sealing core drilling sample after washing oil；To above-mentioned washing oil After above-mentioned non-sealing core drilling sample carry out press hydrargyrum experiment, it is thus achieved that above-mentioned total pore throat radius be distributed；And to above-mentioned non-sealing core drilling Sample carries out core Location, using practical logging NMR (Nuclear Magnetic Resonance) spectrum as the NMR (Nuclear Magnetic Resonance) spectrum of above-mentioned non-sealing core drilling sample.

In another preferred embodiment of the application, above-mentioned core Location is the core Location of Centimeter Level, so can enter to recover The work of rock core place true layer, and then can the most accurately obtain the NMR (Nuclear Magnetic Resonance) spectrum being full of profit of sample of coring.

In order to improve the efficiency of the NMR (Nuclear Magnetic Resonance) spectrum separating sample of coring, it is possible to more quickly by the nuclear magnetic resonance, NMR of sample of coring The efficiency of wave spectrum is divided into the aqueous spectrum of sample of coring and sample oil-containing spectrum of coring, and preferably in above-mentioned steps C, based on crude is composed and left The nuclear magnetic resonance, NMR time corresponding to limit value, the NMR (Nuclear Magnetic Resonance) spectrum of above-mentioned sample of coring is divided into the aqueous spectrum of sample of coring and sample of coring Product oil-containing is composed.

In the another kind of preferred embodiment of the application, above-mentioned steps S3 includes:

According to above-mentioned oil surface relaxivity and the above-mentioned oil-containing spectrum in above-mentioned steps S1, formula (6) is used to calculate above-mentioned neutral reservoir The pore throat radius distribution of oil-containing part, above-mentioned formula (6) is:

R_o=a × ρ_o×T_1o,2o……………….(6)

In above-mentioned formula (6): R_oFor the pore throat radius distribution of above-mentioned neutral reservoir oil-containing part, unit is μm；ρ_oFor above-mentioned oil Surface relaxivity, unit is 10^-3μm/ms；A is constant, when above-mentioned pore shape is for taking 2 time cylindrical, for taking 3 time spherical, 1 is taken during for channel-shaped；T_1o,2oT2 (the T composed for above-mentioned oil-containing_2o) or longitudinal relaxation time (T_1o), unit is ms.

In order to the most accurately obtain the NMR (Nuclear Magnetic Resonance) spectrum being full of water of above-mentioned neutral reservoir oil-containing part, in above-mentioned steps S4, According to above-mentioned water surface relaxation rate, use formula (7) that the pore throat radius distribution of above-mentioned neutral reservoir oil-containing part is corrected, To the NMR (Nuclear Magnetic Resonance) spectrum being full of water of above-mentioned neutral reservoir oil-containing part, above-mentioned formula (7) is:

$T_{1jw,2jw}=\frac{R_o}{a\×{\mathrm{\ρ}}_w}\mathrm{...}\left(7\right)$

In above-mentioned formula (7), T_1jw,2jwThe NMR (Nuclear Magnetic Resonance) spectrum being full of water for above-mentioned neutral reservoir oil-containing part；R_oFor above-mentioned sample of coring The pore throat radius distribution of the oil-containing part of product, unit is μm；ρ_wFor above-mentioned water surface relaxation rate, unit is 10^-3μm/ms；a For constant, when above-mentioned pore shape is the most cylindrical, spherical, shape channel-shaped time, a is respectively 2,3,1.

In the another kind of preferred embodiment of the application, in above-mentioned steps S5, use above-mentioned by above-mentioned steps S1 of formula (8) Aqueous spectrum is overlapped with the NMR (Nuclear Magnetic Resonance) spectrum being full of water of above-mentioned neutral reservoir oil-containing part, obtains above-mentioned neutral reservoir and is full of The NMR (Nuclear Magnetic Resonance) spectrum of water, above-mentioned formula (8) is:

T_1tw,2tw=T_1w,2w∪T_1jw,2jw……………….(8)

In above-mentioned formula (8), T_1tw,2twThe NMR (Nuclear Magnetic Resonance) spectrum of water it is full of for above-mentioned neutral reservoir；T_1w,2wFor above-mentioned aqueous spectrum；T_1jw,2jw The NMR (Nuclear Magnetic Resonance) spectrum being full of water for above-mentioned neutral reservoir oil-containing part.

In order to enable those skilled in the art clearly to understand the technical scheme of the application, below with reference to embodiment with attached Figure is described in detail.

Choosing certain fine and close oilfield of the Junggar Basin Permian System, the reservoir of this densification oilfield is that the particulate of salified lacustrine environment sinks Long-pending, principal lithologic is clastic rock and carbonate rock transitionality rock class.The porosity of reservoir is between 6～16%, and overburden permeability is general All over less than 0.1md, the diameter of reservoir pore space is based on nanoscale, and the oily saturation of reservoir is higher, averagely about 70%, and moistening Property test result indicate that: reservoir wettability is mainly based on neutral wetting, for typical fine and close oil reservoir.The exploration initial stage does not considers The crude oil impact on nuclear magnetic resonance log, directly application routine evaluations method is evaluated, and its result display pore throat radius is bigger than normal, Distribution characteristics is similar to conventional reservoir, it is impossible to meets and produces the technology requirement that upper fine and close oil pore throat radius characterizes.

The technical scheme provided according to the present invention has taken into full account the crude oil impact on nuclear magnetic resonance log, uses the hole of oil-containing part Gap is full of water NMR (Nuclear Magnetic Resonance) spectrum and is corrected the aqueous spectrum of reservoir, and concrete flow process is divided into five steps, is respectively as follows:

The first step, the petrophysics experiment of application system, set up the oil-water separation model of the NMR (Nuclear Magnetic Resonance) spectrum of this district's reservoir, will The NMR (Nuclear Magnetic Resonance) spectrum of reservoir is divided into aqueous spectrum to compose with oil-containing.

Choose polylith sealing core drilling sample, carry out the nuclear magnetic resonance, NMR (i.e. the NMR (Nuclear Magnetic Resonance) spectrum of reservoir) being full of water-oil phase of reservoir Experiment, after measuring its oil saturation, then is carried out sample, and measuring samples is full of the NMR (Nuclear Magnetic Resonance) spectrum of water.

The T2 spectrum that this study area all uses, so the NMR (Nuclear Magnetic Resonance) spectrum in this example is T₂Spectrum.

Test result indicate that, as in figure 2 it is shown, the nuclear magnetic resonance, NMR T being full of water-oil phase of reservoir₂It is full of water after wave spectrum and washing oil NMR (Nuclear Magnetic Resonance) spectrum is compared, and reservoir is full of the nuclear magnetic resonance, NMR T of water-oil phase₂Wave spectrum is longer, and then demonstrates crude oil to T₂Compose permissible Produce large effect.

Utilize oil saturation analytical data and nuclear magnetic resonance, NMR T₂Spectral data binding analysis, analyzes containing of different time to count calculating Relative error between the oil saturation that oil saturation and experiment obtain, time to count during relative error minimum is just for nuclear-magnetism altogether Shake the disengaging time of wave spectrum oil-water separation model, i.e. T2 threshold value.

Specifically, it is to use iterative method to determine the disengaging time of NMR (Nuclear Magnetic Resonance) spectrum oil-water separation model, calculates root-mean-square error Hour T2 value be threshold value A T2, its computing formula is as follows:

$WT2\left(j\right)=\underset{j=1}{\overset{m}{\Σ}}\frac{1}{n}\underset{i=1}{\overset{n}{\Σ}}{({\mathrm{SO}}_i-{\mathrm{SSO}}_{ji})}^2$

Wherein, WT2 (j) is the mean square calculating error of jth iteration T2 threshold value, and n is the number of oil saturation experimental data, SOi is the saturation measurement data of i-th sampling point, SSO_jiI-th for jth iteration T2 threshold value calculates saturation, its meter Calculation formula is as follows:

${\mathrm{SSO}}_{ji}=1-\left(\underset{i=ATS}{\overset{AT\left(j\right)}{\Σ}}{\mathrm{\φ}}_i\right)/\underset{i=ATS}{\overset{ATD}{\Σ}}{\mathrm{\φ}}_i$

Wherein, AT (j) is that oil saturation starts at nuclear magnetic resonance, NMR T2, and unit is ms；ATS is effecive porosity Nuclear magnetic resonance, NMR transverse relaxation time to count；ATD is that the nuclear magnetic resonance, NMR transverse relaxation of effecive porosity terminates the time.

Error analysis result shows, as it is shown on figure 3, when time to count is 6ms, oil-containing oil saturation error is minimum, therefore, The disengaging time composed as oil spectrum in reservoir NMR (Nuclear Magnetic Resonance) spectrum as shown in fig. 4 a with water by 6ms, when i.e. the lower limit deposited composed by crude oil Between be worth, just can complete aqueous spectrum and oil-containing by this value and compose and separate, obtain Fig. 4 b and shown aqueous spectrum and containing shown in Fig. 4 c Oil spectrum.

Second, the petrophysics experiment of application system, set up the computation model of surface relaxivity respectively, calculate the hole of neutral reservoir The water surface relaxation rate of gap contains surface relaxivity and the oil-containing pore surface relaxation rate of water voids with oil surface relaxivity Liao Gai district reservoir Computation model.

Step A, obtains the T being full of water-oil phase of sealing core drilling sample₂Spectrum, is full of the nuclear magnetic resonance, NMR T of water after washing oil₂Spectrum is with total Pore throat radius is distributed.

Water surface relaxation rate is the surface relaxivity containing water voids, and oil surface relaxivity is oil-containing pore surface relaxation rate.In order to build Vertical water surface relaxation rate is with containing surface relaxivity computation model, and the whole district have selected polylith representational sealing core drilling sample, successively Carried out sealing core drilling sample be full of water-oil phase Nuclear Magnetic Resonance Measurement, washing oil after be full of water NMR (Nuclear Magnetic Resonance) spectrum measure and High-pressure mercury is tested, and obtains the T being full of water-oil phase of sealing core drilling sample₂Spectrum, is full of the nuclear magnetic resonance, NMR T of water after washing oil₂Spectrum with Total pore throat radius distribution.

For there is no the depth segment of sealing core drilling, first non-airtight sample is carried out washing oil, measure the core being full of water after obtaining washing oil Magnetic resonance T₂Spectrum；Then, after carrying out pressing hydrargyrum translocation experiment, it is thus achieved that total pore throat radius distribution；Finally non-sealing core drilling sample is entered The core Location of row Centimeter Level, it is thus achieved that the nuclear magnetic resonance log spectrum of the real well correspondence degree of depth, as the T being full of water-oil phase₂ Spectrum.

Step B, according to the NMR (Nuclear Magnetic Resonance) spectrum being full of water after above-mentioned washing oil and capillary pressure curve, calculates above-mentioned water surface relaxation Rate.

Calculating process is as it is shown in figure 5, when the NMR (Nuclear Magnetic Resonance) spectrum weight being full of water after the capillary pressure curve of Fig. 5 a and the washing oil of Fig. 5 b During conjunction, as shown in Figure 5 c, utilize following formula obtain water surface relaxation rate:

${\mathrm{\ρ}}_w=\frac{2\×\mathrm{\γ}\×cos\mathrm{\θ}}{a\×T_{2tw}\×P_{total}}$

In formula, ρ_wFor water surface relaxation rate, unit is μm/ms；θ is the contact angle of hydrargyrum, is 140 ° in this example；γ is hydrargyrum Surface tension, this example is 7.15 × 10^-3psi·cm；A is constant, is 3 in this example；P_totalFor capillary pressure, list Position is Mpa；T_2twFor the nuclear magnetic resonance, NMR T2 being full of water after washing oil, unit is ms.

Surface relaxivity need to divide lithology to calculate the water surface relaxation rate of different lithology, and this example calculation obtains dolomite dust packsand Water surface relaxation rate ρ_w=2.5 × 10^-3μm/ms。

Step C, utilizes the oil-water separation model of NMR (Nuclear Magnetic Resonance) spectrum of this district's reservoir by the nuclear-magnetism being full of water-oil phase of sample of coring Resonance wave spectrum is divided into the aqueous spectrum of sample of coring and sample oil-containing spectrum of coring.The T2 that 6ms is separated as water-oil phase Threshold value, by T₂Time is the oil-containing spectrum of sample of coring more than the part of 6ms, by T₂The time part less than 6ms is sample of coring The aqueous spectrum of product.

Step D, is applied to water surface relaxation rate in the aqueous spectrum that step C obtains as shown in Figure 6 a, calculates above-mentioned coring Pore throat radius containing the water section distribution of sample, as shown in Figure 6 b, computing formula is r_w=a × ρ_w×T_2w.In this example, For white clouds matter Extra-fine sand rock, ρ_w=2.5 × 10^-3μm/ms, a=3.

Step E, according to the distribution of the pore throat radius containing water section of the distribution of above-mentioned total pore throat radius with above-mentioned sample of coring, calculates above-mentioned Core sample oil-containing part pore throat radius distribution, computing formula isWherein, r_oFor above-mentioned sample of coring The pore throat radius distribution of oil-containing part, r_wThe pore throat radius containing water section for above-mentioned sample of coring divides, r_totalFor above-mentioned total pore throat half Footpath is distributed, and the total pore throat radius shown in Fig. 7 a is distributed the aqueous pore throat radius distribution deducted shown in Fig. 7 b, obtains shown in Fig. 7 c Oil-containing part pore throat radius distribution；

Step F, composes with the oil-containing of above-mentioned sample of coring according to the pore throat radius distribution of oil-containing part of above-mentioned sample of coring, uses formula (5) calculating above-mentioned oil surface relaxivity, above-mentioned formula (5) is

${\mathrm{\ρ}}_o=\frac{2\×\mathrm{\γ}\×cos\mathrm{\θ}}{a\×T_{2o}\×P_o}$

In above-mentioned formula, ρ_oFor above-mentioned oil surface relaxivity, θ is the contact angle of above-mentioned hydrargyrum, is 140 ° in this example, and γ is above-mentioned hydrargyrum Surface tension, this example is 7.15 × 10^-3Psi cm, a are 3, P_oFor the hair that the aperture of above-mentioned neutral reservoir oil-containing part is corresponding Pipe pressure, T_2oThe T2 composed for above-mentioned oil-containing.In this example, white clouds matter Extra-fine sand rock is calculated, ρ_oFor 0.25×10^-3μm/ms。

Above-mentioned is all the core experiment for cored interval from step A to step D, obtains water surface relaxation by core experiment data Rate and oil surface relaxivity, then this parameter is applied to subsequent step for the nuclear magnetic resonance, NMR of underground survey of non-cored interval Wave spectrum.

3rd, use above-mentioned oil surface relaxivity ρ_oTo above-mentioned T₂The carrying out at spectrum oil-containing peak processes and obtains above-mentioned reservoir oil-containing part Pore throat radius is distributed, and specific formula for calculation is as follows

R_o=a × ρ_o×T_2o

R in above formula_oFor the pore diameter of reservoir oil-containing part, unit is μm；A is 3；T_2oThe transverse relaxation composed for above-mentioned oil-containing Time, unit is ms, ρ_o=0.25 × 10-3 μm/ms.

4th, use water surface relaxation rate that the pore throat radius of above-mentioned reservoir oil-containing part is distributed and be corrected, after being corrected The NMR (Nuclear Magnetic Resonance) spectrum being full of water of oil-containing part of above-mentioned neutral reservoir, specific formula for calculation is as follows:

$T_{1jw,2jw}=\frac{R_o}{a\×{\mathrm{\ρ}}_w}$

In above formula, T_1jw,2jwThe aqueous spectrum obtained after oil-containing corrects for reservoir oil-containing spectrum；R_oPore throat half for reservoir oil-containing part Footpath, unit is μm；ρ_w=2.5 × 10^-3μm/ms；A is 3.

5th, use formula T_2tw=T_2w∪T_2jw, by T_2jwWith T_2wIt is overlapped, it is thus achieved that neutral reservoir is full of the nuclear-magnetism of water Resonance wave spectrum.

Fig. 8 is the test result figure finally given, and is divided into seven regions in figure, is gamma ray curve (GR) in first area, Second area is the coordinate of the degree of depth, and the 3rd region is parent magnetic resonance T₂Spectrum (the most original well logging spectrum), is to be full of in the 4th region The nuclear magnetic resonance, NMR T of water₂Spectrum (the most complete aqueous spectrum), is that the pore throat radius distribution curve (fine rule) calculated is with real in the 5th region Testing and analyze analysis pore throat radius distribution curve (thick line) obtained, the 6th region is for calculating capillary pressure curve (fine rule) and pressure hydrargyrum The capillary pressure curve (thick line) that experimental analysis obtains, is duty pressure and duty pressure (round dot) in analysis in calculating in the 7th region.

Fig. 9 show the pore throat radius distribution curve (dotted line) calculated that a depth value in the 5th region in Fig. 8 is corresponding with Analysis pore throat radius distribution curve (solid line) that experimental analysis obtains.

As shown in Figure 9, the pore throat half that the NMR (Nuclear Magnetic Resonance) spectrum utilizing the neutral reservoir after the correction that the application obtains to be full of water calculates The pore throat radius distribution that footpath distribution and test obtain has preferable concordance.

As can be seen from the above description, the above embodiments of the present invention achieve following technique effect:

The bearing calibration of the application has taken into full account the oiliness impact on nuclear magnetic resonance log of neutral reservoir, uses neutral reservoir The hole of oil-containing part be full of water NMR (Nuclear Magnetic Resonance) spectrum, the aqueous spectrum of centering reservoir is corrected, during the wave spectrum after correction is Property reservoir the NMR (Nuclear Magnetic Resonance) spectrum being full of water, this wave spectrum can correctly characterize the pore throat radius distribution of neutral reservoir so that ability Field technique personnel can accurately analyze by centering reservoir, and then carries out sending out of rational exploration and development, beneficially petroleum industry Exhibition.

The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for those skilled in the art For, the present invention can have various modifications and variations.All within the spirit and principles in the present invention, any amendment of being made, etc. With replacement, improvement etc., should be included within the scope of the present invention.

Claims (10)

1. the bearing calibration of a nuclear magnetic resonance log wave spectrum, it is characterised in that described bearing calibration includes:

Step S1, uses NMR (Nuclear Magnetic Resonance) spectrum disjunctive model that the NMR (Nuclear Magnetic Resonance) spectrum of neutral reservoir is divided into aqueous spectrum and oil-containing Spectrum；

Step S2, uses water surface relaxation rate and the oil meter of the hole of the described neutral reservoir of surface relaxivity computation model calculating Face relaxation rate；

Step S3, utilizes the pore throat radius distribution of the oil-containing part of the described described neutral reservoir of oil surface relaxivity calculating；

Step S4, is corrected the pore throat radius distribution of described oil-containing part according to described water surface relaxation rate, obtains school The NMR (Nuclear Magnetic Resonance) spectrum being full of water of the oil-containing part of the described neutral reservoir after just；And

Step S5, is overlapped the NMR (Nuclear Magnetic Resonance) spectrum being full of water of described aqueous spectrum with the oil-containing part of described correction, The NMR (Nuclear Magnetic Resonance) spectrum being full of water of the described neutral reservoir after being corrected.

Bearing calibration the most according to claim 1, it is characterised in that described NMR (Nuclear Magnetic Resonance) spectrum disjunctive model is:

T_1w,2w∪T_1o,2o=T_1,2...................(1)

In formula (1):

T_1,2NMR (Nuclear Magnetic Resonance) spectrum for described neutral reservoir；

T_1w,2wAqueous spectrum for described neutral reservoir；

T_1o,2oOil-containing spectrum for described neutral reservoir.

Bearing calibration the most according to claim 1, it is characterised in that described step S2 includes:

Step A, chooses sample of coring, and is full of the nuclear-magnetism of water described in measurement after the core NMR (Nuclear Magnetic Resonance) spectrum of sample, washing oil Resonance wave spectrum is distributed with total pore throat radius of described neutral reservoir；

Step B, according to the NMR (Nuclear Magnetic Resonance) spectrum being full of water after described washing oil and capillary pressure curve, uses formula (2) meter Calculating described water surface relaxation rate, described formula (2) is:

${\mathrm{\ρ}}_w=\frac{2\×\mathrm{\γ}\×cos\mathrm{\θ}}{a\×T_{1tw,2tw}\×P_{total}}\mathrm{...}\left(2\right)$

In described formula (2):

ρ_wFor described water surface relaxation rate,

θ is the contact angle of hydrargyrum,

γ is the surface tension of hydrargyrum,

P_totalFor capillary pressure,

T_1tw,2twFor the nuclear magnetic resonance, NMR T2 T being full of water after described washing oil_2twOr longitudinal relaxation time T_1tw, single Position is ms；

Step C, is divided into the aqueous spectrum of sample of coring and sample of coring by the NMR (Nuclear Magnetic Resonance) spectrum being full of profit of described sample of coring Product oil-containing is composed；

Step D, according to described water surface relaxation rate and the described aqueous spectrum of sample of coring, uses and cores described in formula (3) calculating Pore throat radius containing the water section distribution of sample, described formula (3) is:

r_w=a × ρ_w×T_1w,2w...................(3)

In described formula (3):

r_wThe pore throat radius containing water section for described sample of coring is distributed,

ρ_wFor described water surface relaxation rate,

A is constant, when described pore shape is the most cylindrical, spherical, shape channel-shaped time, a is respectively 2,3,1,

T_1w,2wT2 T for the described aqueous spectrum of sample of coring_2wOr longitudinal relaxation time T_1w；

Step E, according to the distribution of the pore throat radius containing water section of the distribution of described total pore throat radius with described sample of coring, uses Formula (4) calculate described in core sample oil-containing part pore throat radius distribution, described formula (4) is:

$r_o=\stackrel{\‾}{r_w\∩r_{total}}\mathrm{...}\left(4\right)$

In described formula (4):

r_oPore throat radius for the oil-containing part of described sample of coring is distributed,

r_wThe pore throat radius containing water section for described sample of coring is distributed,

r_totalIt is distributed for described total pore throat radius；

Step F, composes with the oil-containing of described sample of coring according to the pore throat radius distribution of oil-containing part of described sample of coring, adopts Calculating described oil surface relaxivity by formula (5), described formula (5) is

${\mathrm{\ρ}}_o=\frac{2\×\mathrm{\γ}\×cos\mathrm{\θ}}{a\×T_{1o,2o}\×P_o}\mathrm{...}\left(5\right)$

In described formula (5):

ρ_oFor described oil surface relaxivity；

θ is the contact angle of described hydrargyrum,

γ is the surface tension of described hydrargyrum,

A is constant, when described pore shape is the most cylindrical, spherical, shape channel-shaped time, a is respectively 2,3,1,

P_oFor the capillary pressure that the aperture of described neutral reservoir oil-containing part is corresponding, according to the oil-containing part of described sample of coring Pore throat radius distribution be calculated,

T_1o,2oThe T2 T composed for described oil-containing_2oOr longitudinal relaxation time T_1o。

Bearing calibration the most according to claim 3, it is characterised in that described step A includes:

Choose sealing core drilling sample, measure the NMR (Nuclear Magnetic Resonance) spectrum of described sealing core drilling sample；

Described sealing core drilling sample is carried out washing oil, the nuclear magnetic resonance, NMR ripple being full of water of sample of coring described in after measuring washing oil Spectrum；And

Carry out described sealing core drilling sample pressing hydrargyrum experiment, it is thus achieved that described total pore throat radius is distributed.

Bearing calibration the most according to claim 3, it is characterised in that described step A includes:

Choosing non-sealing core drilling sample, described non-sealing core drilling sample carries out washing oil, that measures after washing oil is described non-airtight Core the NMR (Nuclear Magnetic Resonance) spectrum being full of water of sample；

Carry out the described non-sealing core drilling sample after described washing oil pressing hydrargyrum experiment, it is thus achieved that described total pore throat radius is distributed；With And

Described non-sealing core drilling sample is carried out core Location, practical logging NMR (Nuclear Magnetic Resonance) spectrum non-airtight is taken as described The NMR (Nuclear Magnetic Resonance) spectrum of heart sample.

Bearing calibration the most according to claim 5, it is characterised in that described core Location is the core Location of Centimeter Level.

Bearing calibration the most according to claim 3, it is characterised in that in described step C, based on crude is composed and is left limit value The corresponding nuclear magnetic resonance, NMR time, the NMR (Nuclear Magnetic Resonance) spectrum of described sample of coring is divided into the aqueous spectrum of sample of coring and sample of coring Oil-containing is composed.

Bearing calibration the most according to claim 1, it is characterised in that described step S3 includes:

According to described oil surface relaxivity and the described oil-containing spectrum in described step S1, formula (6) is used to calculate described neutrality The pore throat radius distribution of reservoir oil-containing part, described formula (6) is:

R_o=a × ρ_o×T_1o,2o...................(6)

In described formula (6):

R_oPore throat radius distribution for described neutral reservoir oil-containing part；

ρ_oFor described oil surface relaxivity；

A is constant, when described pore shape is for taking 2 time cylindrical, for taking 3 time spherical, takes 1 during for channel-shaped；

T_1o,2oThe T2 T composed for described oil-containing_2oOr longitudinal relaxation time T_1o。

Bearing calibration the most according to claim 3, it is characterised in that in described step S4, according to described water surface relaxation Rate, uses formula (7) to be corrected the pore throat radius distribution of described neutral reservoir oil-containing part, obtains described neutral reservoir The NMR (Nuclear Magnetic Resonance) spectrum being full of water of oil-containing part, described formula (7) is:

$T_{1jw,2jw}=\frac{R_o}{a\×{\mathrm{\ρ}}_w}\mathrm{...}\left(7\right)$

In described formula (7):

T_1jw,2jwThe NMR (Nuclear Magnetic Resonance) spectrum being full of water for described neutral reservoir oil-containing part；

R_oPore throat radius distribution for the oil-containing part of described sample of coring；

ρ_wFor described water surface relaxation rate；

A is constant, when described pore shape is the most cylindrical, spherical, shape channel-shaped time, a is respectively 2,3,1.

Bearing calibration the most according to claim 1, it is characterised in that in described step S5, uses formula (8) by described step Described aqueous spectrum in rapid S1 is overlapped with the NMR (Nuclear Magnetic Resonance) spectrum being full of water of described neutral reservoir oil-containing part, obtains Described neutral reservoir is full of the NMR (Nuclear Magnetic Resonance) spectrum of water, and described formula (8) is:

T_1tw,2tw=T_1w,2w∪T_1jw,2jw...................(8)

In described formula (8):

T_1tw,2twThe NMR (Nuclear Magnetic Resonance) spectrum of water it is full of for described neutral reservoir；

T_1w,2wFor described aqueous spectrum；

T_1jw,2jwThe NMR (Nuclear Magnetic Resonance) spectrum being full of water for described neutral reservoir oil-containing part.