CN106324008B - Method for correcting nuclear magnetic resonance logging spectrum - Google Patents

Method for correcting nuclear magnetic resonance logging spectrum Download PDF

Info

Publication number
CN106324008B
CN106324008B CN201510383206.1A CN201510383206A CN106324008B CN 106324008 B CN106324008 B CN 106324008B CN 201510383206 A CN201510383206 A CN 201510383206A CN 106324008 B CN106324008 B CN 106324008B
Authority
CN
China
Prior art keywords
oil
spectrum
water
formula
sample
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510383206.1A
Other languages
Chinese (zh)
Other versions
CN106324008A (en
Inventor
王振林
孙中春
齐红岩
王志维
于静
王伟
黄亚雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201510383206.1A priority Critical patent/CN106324008B/en
Publication of CN106324008A publication Critical patent/CN106324008A/en
Application granted granted Critical
Publication of CN106324008B publication Critical patent/CN106324008B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The invention provides a method for correcting a nuclear magnetic resonance logging spectrum. The correction method comprises the following steps: step S1, dividing the nuclear magnetic resonance spectrum of the neutral reservoir into a water-containing spectrum and an oil-containing spectrum by adopting a nuclear magnetic resonance spectrum separation model; step S2, calculating the water surface relaxation rate and the oil surface relaxation rate of the pores of the neutral reservoir by adopting a surface relaxation rate calculation model; step S3, calculating the pore throat radius distribution of the oil-containing part of the neutral reservoir by using the oil surface relaxation rate; step S4, correcting the pore throat radius distribution of the oil-containing part according to the water surface relaxation rate to obtain a nuclear magnetic resonance spectrum of saturated water of the oil-containing part of the corrected neutral reservoir; and step S5, superposing the water spectrum and the corrected nuclear magnetic resonance spectrum of the saturated water of the oil-containing part to obtain the corrected nuclear magnetic resonance spectrum of the saturated water of the neutral reservoir. The correction method can obtain accurate nuclear magnetic resonance spectrum of saturated water of the neutral reservoir, and further can accurately analyze the neutral reservoir.

Description

The bearing calibration of nuclear magnetic resonance log wave spectrum
Technical field
The present invention relates to oil fields, in particular to a kind of bearing calibration of nuclear magnetic resonance log wave spectrum.
Background technology
Petroleum industry passes through the development of over half a century, and the target of exploration and development becomes increasingly complex, the infiltration of reservoir Rate is lower and lower, hyposmosis at present with fine and close oily reservoir be or will be exploration and development major heading, such reservoir comments One of core content of valence is exactly pore structure evaluation, and existing pore structure assessment technique mainly applies nuclear magnetic resonance log, When the NMR relaxation mechanism of fluid in pore media is Volume Relaxation, the Response Mechanism of NMR spectrum is as follows:
Wherein, T1,2For NMR relaxation time, ρ1,2For the surface relaxivity of rock, S is the surface area of hole, and V is The volume of hole.The model is the evaluation model of most physical significance, but it is reservoir containing only single to apply the premise of the model Pore-fluid.Subsequent countries in the world well log interpretation expert proposes the evaluation model of some series, before the application of these models It is for hydrophilic reservoir to carry, and nuclear magnetic resonance log is influenced to want small by oily.
Most of conventional reservoir is hydrophilic rock at present, and the reservoir in nuclear magnetic resonance log investigative range is because by mud The influence of intrusion is starched, residual oil saturation is low, and nuclear magnetic resonance log at this time is influenced very faint, response spy by oiliness Sign is single fluid, i.e. the nuclear magnetic resonance of water, and at this moment above-mentioned model can be used for evaluating the pore throat radius distribution of reservoir.
But for the reservoir of neutral wetting, such as fine and close oily reservoir, oil saturation is higher, and pore throat radius It is small, mud be difficult to reservoir carry out displacement, nuclear magnetic resonance log is big by the effect of signals of crude oil, at this moment using existing method come The pore throat radius distribution of the fine and close oily reservoir of analysis calculating has the great error of generation.
Analysis method in the prior art does not fully take into account influence of the neutral reservoir oiliness to nuclear magnetic resonance log, If being difficult to correctly characterize the pore throat radius of neutral reservoir without carrying out oiliness correction to the nuclear magnetic resonance log wave spectrum of reservoir Distribution, therefore, it is necessary to a kind of methods of nuclear magnetic resonance log wave spectrum oil gas correction.
Invention content
The main purpose of the present invention is to provide a kind of bearing calibrations of nuclear magnetic resonance log wave spectrum, to solve the prior art In the NMR spectrum of water the problem of cannot accurately analyzing the distribution of neutral reservoir pore throat radius.
To achieve the goals above, according to an aspect of the invention, there is provided a kind of school of nuclear magnetic resonance log wave spectrum Correction method, above-mentioned bearing calibration include:Step S1, using NMR spectrum disjunctive model by the nuclear magnetic resonance wave of neutral reservoir Spectrum is divided into aqueous spectrum and is composed with oil-containing;Step S2 calculates the water of the hole of above-mentioned neutral reservoir using surface relaxivity computation model Surface relaxivity and oil meter face relaxation rate;Step S3 calculates the oil-containing portion of above-mentioned neutral reservoir using above-mentioned oil meter face relaxation rate The pore throat radius distribution divided;Step S4 is distributed the pore throat radius of above-mentioned oil-containing part according to above-mentioned water surface relaxation rate and carries out Correction, the NMR spectrum for being full of water of the oil-containing part of the above-mentioned neutral reservoir after being corrected;And step S5, to upper It states aqueous spectrum and the NMR spectrum for being full of water of the oil-containing part of above-mentioned correction is overlapped, after being corrected among the above The NMR spectrum for being full of water of property reservoir.
Further, above-mentioned NMR spectrum disjunctive model is:T1w,2w∪T1o,2o=T1,2(1), in formula (1):T1,2For The NMR spectrum of above-mentioned neutrality reservoir;T1w,2wFor the aqueous spectrum of above-mentioned neutral reservoir;T1o,2oFor containing for above-mentioned neutral reservoir Oil spectrum.
Further, above-mentioned steps S2 includes:
Step A chooses coring sample, measures the nuclear-magnetism for being full of water after the NMR spectrum of above-mentioned coring sample, washing oil Resonance wave spectrum and total pore throat radius of above-mentioned neutral reservoir are distributed;
Step B is counted according to the NMR spectrum and capillary pressure curve for being full of water after above-mentioned washing oil using formula (2) It counts in stating water surface relaxation rate, above-mentioned formula (2) is:In above-mentioned formula (2):ρwFor above-mentioned water surface Relaxation rate, θ are the contact angle of mercury, and γ is the surface tension of mercury, PtotalFor capillary pressure, T1tw,2twFor being full of after above-mentioned washing oil The nuclear magnetic resonance lateral relaxation time T of water2twOr longitudinal relaxation time T1tw, unit ms;
The NMR spectrum for being full of grease of above-mentioned coring sample is divided into the aqueous spectrum of coring sample and coring by step C Sample oil-containing is composed;
Step D calculates above-mentioned coring according to above-mentioned water surface relaxation rate and the above-mentioned aqueous spectrum of coring sample using formula (3) The pore throat radius distribution containing water section of sample, above-mentioned formula (3) are:rw=a × ρw×T1w,2wIn above-mentioned formula (3):rwIt is taken to be above-mentioned The pore throat radius distribution containing water section of heart sample, ρwFor above-mentioned water surface relaxation rate, a is constant, when above-mentioned pore shape point Not Wei cylindrical, spherical, shape channel-shaped when, a is respectively 2,3,1, T1w,2wFor the lateral relaxation time of the above-mentioned aqueous spectrum of coring sample T2wOr longitudinal relaxation time T1w
Step E is distributed according to above-mentioned total pore throat radius distribution and the pore throat radius containing water section of above-mentioned coring sample, adopts The pore throat radius distribution of the oil-containing part of above-mentioned coring sample is calculated with formula (4), above-mentioned formula (4) isAbove-mentioned formula (4) in, roFor the pore throat radius distribution of the oil-containing part of above-mentioned coring sample, rwFor the hole containing water section of above-mentioned coring sample Larynx radius distribution, rtotalIt is distributed for above-mentioned total pore throat radius;
Step F is composed according to the distribution of the pore throat radius of the oil-containing part of above-mentioned coring sample and the oil-containing of above-mentioned coring sample, Above-mentioned oil meter face relaxation rate is calculated using formula (5), above-mentioned formula (5) isIn above-mentioned formula (5), ρoIt is above-mentioned Oil meter face relaxation rate;θ is the contact angle of above-mentioned mercury, and γ is the surface tension of above-mentioned mercury, and a is constant, when above-mentioned pore shape point Not Wei cylindrical, spherical, shape channel-shaped when, a is respectively 2,3,1, PoFor the corresponding hair in aperture of above-mentioned neutral reservoir oil-containing part Pipe pressure is calculated, T according to the distribution of the pore throat radius of the oil-containing part of above-mentioned coring sample1o,2oFor the cross of above-mentioned oil-containing spectrum To relaxation time T2oOr longitudinal relaxation time T1o
Further, above-mentioned steps A includes:Sealing core drilling sample is chosen, the nuclear-magnetism for measuring above-mentioned sealing core drilling sample is total Vibration wave is composed;Washing oil is carried out to above-mentioned sealing core drilling sample, measures the nuclear magnetic resonance for being full of water of the above-mentioned coring sample after washing oil Wave spectrum;And pressure mercury experiment is carried out to above-mentioned sealing core drilling sample, obtain above-mentioned total pore throat radius distribution.
Further, above-mentioned steps A includes:Non- sealing core drilling sample is chosen, above-mentioned non-sealing core drilling sample is washed Oil measures the NMR spectrum for being full of water of the above-mentioned non-sealing core drilling sample after washing oil;To above-mentioned non-after above-mentioned washing oil Sealing core drilling sample carries out pressure mercury experiment, obtains above-mentioned total pore throat radius distribution;And above-mentioned non-sealing core drilling sample is carried out Core Location, using practical logging NMR spectrum as the NMR spectrum of above-mentioned non-sealing core drilling sample.
Further, above-mentioned core Location is the core Location of Centimeter Level.
It further, the nuclear magnetic resonance time corresponding to based on crude preservation lower limiting value, will be above-mentioned in above-mentioned steps C The NMR spectrum of coring sample is divided into the aqueous spectrum of coring sample and is composed with coring sample oil-containing.
Further, above-mentioned steps S3 includes according to the above-mentioned oil-containing in above-mentioned oil meter face relaxation rate and above-mentioned steps S1 Spectrum, the pore throat radius that above-mentioned neutral reservoir oil-containing part is calculated using formula (6) are distributed, and above-mentioned formula (6) is:Ro=a × ρo× T1o,2o, in above-mentioned formula (6), RoFor the pore throat radius distribution of above-mentioned neutral reservoir oil-containing part;ρoFor above-mentioned oil meter face relaxation rate; A is constant, and 2 are taken when above-mentioned pore shape is cylindrical, takes 3 when being spherical, 1 is taken when being channel-shaped;T1o,2oIt is composed for above-mentioned oil-containing Lateral relaxation time T2oOr longitudinal relaxation time T1o
Further, in above-mentioned steps S4, according to above-mentioned water surface relaxation rate, using formula (7) to above-mentioned neutral reservoir The pore throat radius distribution of oil-containing part is corrected, and obtains the nuclear magnetic resonance wave for being full of water of above-mentioned neutral reservoir oil-containing part Spectrum, above-mentioned formula (7) is,In above-mentioned formula (7), T1jw,2jwIt is full of water for above-mentioned neutral reservoir oil-containing part NMR spectrum;RoFor the pore throat radius distribution of the oil-containing part of above-mentioned coring sample;ρwFor above-mentioned water surface relaxation rate;a For constant, when above-mentioned pore shape is respectively cylinder, spherical shape, shape channel-shaped, a is respectively 2,3,1.
Further, in above-mentioned steps S5, using formula (8) by above-mentioned steps S1 it is above-mentioned it is aqueous spectrum with above-mentioned neutrality The NMR spectrum for being full of water of reservoir oil-containing part is overlapped, and obtains the nuclear magnetic resonance wave that above-mentioned neutral reservoir is full of water Spectrum, above-mentioned formula (8) is T1tw,2tw=T1w,2w∪T1jw,2jw, in above-mentioned formula (8), T1tw,2twIt is full of the core of water for above-mentioned neutral reservoir Magnetic Resonance Spectrum;T1w,2wFor above-mentioned aqueous spectrum;T1jw,2jwFor the nuclear magnetic resonance wave for being full of water of above-mentioned neutral reservoir oil-containing part Spectrum.
Using the bearing calibration of the present invention, influence of the oiliness to nuclear magnetic resonance log of neutral reservoir has been fully considered, Hole using the oil-containing part of neutral reservoir is full of water NMR spectrum, is corrected to the aqueous spectrum of neutral reservoir, school Wave spectrum after just is the NMR spectrum for being full of water of neutral reservoir, which can correctly characterize the pore throat half of neutral reservoir Diameter is distributed so that those skilled in the art can accurately analyze neutral reservoir, and then carry out rational exploration and development, have Conducive to the development of petroleum industry.
Description of the drawings
The accompanying drawings which form a part of this application are used to provide further understanding of the present invention, and of the invention shows Meaning property embodiment and its explanation are not constituted improper limitations of the present invention for explaining the present invention.In the accompanying drawings:
Fig. 1 shows a kind of flow signal of the bearing calibration for the nuclear magnetic resonance log wave spectrum that exemplary embodiment provides Figure;
Fig. 2 shows a kind of nuclear magnetic resonance T for being full of water-oil phase of preferred embodiment reservoir2It is full of after wave spectrum and washing oil The NMR spectrum comparison diagram of water;
Fig. 3 shows that a well sealing core drilling well section that a kind of preferred embodiment of the application provides changes using different AT values The Analysis of Mean Square Error figure that generation calculates;
Fig. 4 a show the NMR spectrum for being full of water-oil phase that a kind of preferred embodiment provides be separated into water spectrum with The schematic diagram of oil spectrum;
Fig. 4 b show the schematic diagram for the aqueous spectrum that the NMR spectrum from Fig. 4 a is isolated;
Fig. 4 c show the schematic diagram for the oil-containing spectrum that the NMR spectrum from Fig. 4 a is isolated;
Fig. 5 a show a kind of capillary pressure curve that preferred embodiment provides;
The NMR spectrum capillary pressure curve of water is provided after the washing oil provided Figure 5b shows that a kind of preferred embodiment;
Fig. 5 c show the NMR spectrum capillary pressure curve for being full of water after the washing oil that a kind of preferred embodiment provides;
Fig. 6 a show a kind of aqueous spectrum for the coring sample that preferred embodiment provides;
Fig. 6 b show a kind of pore throat radius distribution containing water section for the coring sample that preferred embodiment provides;
Fig. 7 a show a kind of total pore throat radius distribution schematic diagram that preferred embodiment provides;
Fig. 7 b show that a kind of preferred embodiment provides aqueous pore throat radius distribution schematic diagram;
Fig. 7 c show a kind of pore throat radius for the oil-containing part obtained according to Fig. 7 a and Fig. 7 b that preferred embodiment provides Distribution schematic diagram;
Fig. 8 shows the test result figure that a kind of preferred embodiment of the application provides;And
Fig. 9 shows the 5th some depth point of region in Fig. 8, corrects to obtain completely aqueous spectrum meter using by oiliness The pore throat radius distribution curve comparison diagram that the throat distribution feature of calculation is obtained with experiment.
Specific implementation mode
It should be noted that in the absence of conflict, the features in the embodiments and the embodiments of the present application can phase Mutually combination.The present invention will be described in detail below with reference to the accompanying drawings and embodiments.
This application provides a kind of bearing calibrations of nuclear magnetic resonance log wave spectrum, as shown in Figure 1, the bearing calibration includes: The NMR spectrum of neutral reservoir is divided into aqueous spectrum using NMR spectrum disjunctive model and is composed with oil-containing by step S1;Step Rapid S2 calculates the water surface relaxation rate and oil meter face relaxation of the hole of above-mentioned neutral reservoir using surface relaxivity computation model Rate;Step S3, the pore throat radius that the oil-containing part of above-mentioned neutral reservoir is calculated using above-mentioned oil meter face relaxation rate are distributed;Step S4 is corrected the pore throat radius distribution of above-mentioned oil-containing part according to above-mentioned water surface relaxation rate, above-mentioned after being corrected The NMR spectrum for being full of water of the oil-containing part of neutral reservoir;And step S5, to above-mentioned aqueous spectrum and above-mentioned correction The NMR spectrum for being full of water of oil-containing part is overlapped, the nuclear-magnetism for being full of water of the above-mentioned neutral reservoir after being corrected Resonate wave spectrum.All steps are the core samples for coring section in above-mentioned steps S2, and above-mentioned steps S3, S4 and S5 are to be directed to The practical NMR spectrum measured in non-coring section underground.
Above-mentioned bearing calibration has fully considered influence of the oiliness to nuclear magnetic resonance log of neutral reservoir, using neutrality The hole of the oil-containing part of reservoir is full of water NMR spectrum, is corrected to the aqueous spectrum of neutral reservoir, the wave after correction Spectrum is the NMR spectrum for being full of water of neutral reservoir, which can correctly characterize the pore throat radius distribution of neutral reservoir, Those skilled in the art accurately analyze neutral reservoir, and then carry out rational exploration and development, are conducive to stone Oily industrial expansion.
It is composed with oil-containing in order to which the NMR spectrum of neutral reservoir is more accurately divided into aqueous spectrum, the application is preferred Above-mentioned NMR spectrum disjunctive model is:
In formula (1), T1,2For the NMR spectrum of above-mentioned neutral reservoir;T1w,2wFor the aqueous spectrum of above-mentioned neutral reservoir; T1o,2oFor the oil-containing spectrum of above-mentioned neutral reservoir.
NMR spectrum can be longitudinal relaxation time (T1) spectrum, can also be lateral relaxation time (T2) spectrum.For oil The disjunctive model of water spectrum is exactly by finding a NMR relaxation time point (T1With T2), which corresponds to former The lower limiting value of oily preservation.The wave spectrum on the nuclear magnetic resonance time point left side (being less than the time) is aqueous spectrum, and rock is mainly Hydrophilic rock;The wave spectrum at the nuclear magnetic resonance time point (being more than the time) is composed for oil-containing, and rock is oil-wet rock, on the whole Rock shows as neutral wetting.The NMR relaxation time point can be calculated using grease characteristic nmr difference, The corresponding NMR relaxation time value of reservoir crude oil minimum filling radius determines when can also be applied to Tibetan, can also be by two kinds of sides Method, which combines, to be determined.
In another preferred embodiment of the application, above-mentioned steps S2 includes:
Step A chooses coring sample, measures the nuclear-magnetism for being full of water after the NMR spectrum of above-mentioned coring sample, washing oil Resonance wave spectrum and total pore throat radius of above-mentioned neutral reservoir are distributed;
Step B is counted according to the NMR spectrum and capillary pressure curve for being full of water after above-mentioned washing oil using formula (2) It counts in stating water surface relaxation rate, above-mentioned formula (2) is:
In formula (2), ρwFor above-mentioned water surface relaxation rate, unit 10-3μm/ms, θ is the contact angle of mercury, and γ is the table of mercury Face tension, unit psicm, PtotalFor capillary pressure, unit MPa, T1tw,2twFor the core for being full of water after above-mentioned washing oil Magnetic resonance lateral relaxation time (T2tw) or longitudinal relaxation time (T1tw), unit ms;
The NMR spectrum for being full of grease of above-mentioned coring sample is divided into the aqueous spectrum of coring sample and coring by step C Sample oil-containing is composed;
Step D calculates above-mentioned coring according to above-mentioned water surface relaxation rate and the above-mentioned aqueous spectrum of coring sample using formula (3) The pore throat radius distribution containing water section of sample, above-mentioned formula (3) are:
rw=a × ρw×T1w,2w……………….(3)
In above-mentioned formula (3), rwThe pore throat radius distribution containing water section that is above-mentioned coring sample, unit are μm ρwIt is above-mentioned Water surface relaxation rate, unit 10-3μm/ms, and a is constant, when above-mentioned pore shape is respectively cylinder, spherical shape, shape channel-shaped, A is respectively 2,3,1, T1w,2wFor the lateral relaxation time (T of the above-mentioned aqueous spectrum of coring sample2w) or longitudinal relaxation time (T1w), it is single Position is ms;
Step E is distributed according to above-mentioned total pore throat radius distribution and the pore throat radius containing water section of above-mentioned coring sample, adopts The pore throat radius distribution of the oil-containing part of above-mentioned coring sample is calculated with formula (4), above-mentioned formula (4) is:
In above-mentioned formula (4):roFor the pore throat radius distribution of the oil-containing part of above-mentioned coring sample, rwFor above-mentioned coring sample The pore throat radius containing water section point, rtotalIt is distributed for above-mentioned total pore throat radius;
Step F is composed according to the distribution of the pore throat radius of the oil-containing part of above-mentioned coring sample and the oil-containing of above-mentioned coring sample, Above-mentioned oil meter face relaxation rate is calculated using formula (5), calculation formula is:
In above-mentioned formula (5):ρoFor above-mentioned oil meter face relaxation rate, unit 10-3μm/ms;θ is the contact angle of above-mentioned mercury, single Position for °, γ is the surface tension of above-mentioned mercury, and unit psicm, a are constant, when above-mentioned pore shape be respectively it is cylindrical, When spherical shape, shape channel-shaped, a is respectively 2,3,1, PoFor the corresponding capillary pressure in aperture of above-mentioned neutral reservoir oil-containing part, according to The pore throat radius distribution of the oil-containing part of above-mentioned coring sample is calculated, unit MPa, T1o,2oFor the transverse direction of above-mentioned oil-containing spectrum Relaxation time (T2o) or longitudinal relaxation time (T1o), unit ms.
It adopts with the aforedescribed process, is capable of the water surface relaxation rate and oil meter of the more accurately hole of calculated neutral reservoir Face relaxation rate further ensures the accuracy of the bearing calibration.
In order to further obtain being full of the nuclear magnetic resonance of water after the NMR spectrum of accurate above-mentioned coring sample, washing oil Total pore throat radius of wave spectrum and above-mentioned neutral reservoir is distributed, and those skilled in the art can choose sealing core drilling according to actual conditions Sample is tested with non-sealing core drilling sample.
When sample is sealing core drilling sample, above-mentioned steps A includes:Sealing core drilling sample is chosen, above-mentioned closed take is measured The NMR spectrum of heart sample;Washing oil is carried out to above-mentioned sealing core drilling sample, measures the full of the above-mentioned coring sample after washing oil Aqueous NMR spectrum;And pressure mercury experiment is carried out to above-mentioned sealing core drilling sample, obtain above-mentioned total pore throat radius distribution.
When sample is non-sealing core drilling sample, above-mentioned steps A includes:Non- sealing core drilling sample is chosen, to above-mentioned non-close It closes coring sample and carries out washing oil, measure the NMR spectrum for being full of water of the above-mentioned non-sealing core drilling sample after washing oil;To upper It states the above-mentioned non-sealing core drilling sample after washing oil and carries out pressure mercury experiment, obtain above-mentioned total pore throat radius distribution;And to above-mentioned non- Sealing core drilling sample carries out core Location, is total to practical logging NMR spectrum as the nuclear-magnetism of above-mentioned non-sealing core drilling sample Vibration wave is composed.
In another preferred embodiment of the application, above-mentioned core Location is the core Location of Centimeter Level, in this way can The work of true layer into where restoring rock core, and then further can accurately obtain the nuclear magnetic resonance for being full of grease of coring sample Wave spectrum.
It, can be more quickly by the core of coring sample in order to improve the efficiency for the NMR spectrum for detaching coring sample The efficiency of Magnetic Resonance Spectrum is divided into the aqueous spectrum of coring sample and is composed with coring sample oil-containing, preferably in above-mentioned steps C, based on crude The NMR spectrum of above-mentioned coring sample is divided into the aqueous spectrum of coring sample by the nuclear magnetic resonance time corresponding to preservation lower limiting value It is composed with coring sample oil-containing.
In another preferred embodiment of the application, above-mentioned steps S3 includes:
It is composed according to the above-mentioned oil-containing in above-mentioned oil meter face relaxation rate and above-mentioned steps S1, above-mentioned neutrality is calculated using formula (6) The pore throat radius of reservoir oil-containing part is distributed, and above-mentioned formula (6) is:
Ro=a × ρo×T1o,2o……………….(6)
In above-mentioned formula (6):RoFor the pore throat radius distribution of above-mentioned neutral reservoir oil-containing part, unit is μm;ρoFor above-mentioned oil Surface relaxivity, unit 10-3μm/ms;A is constant, and 2 are taken when above-mentioned pore shape is cylindrical, takes 3 when being spherical, is 1 is taken when channel-shaped;T1o,2oFor the lateral relaxation time (T of above-mentioned oil-containing spectrum2o) or longitudinal relaxation time (T1o), unit ms.
In order to further accurately obtain the NMR spectrum for being full of water of above-mentioned neutral reservoir oil-containing part, in above-mentioned step In rapid S4, according to above-mentioned water surface relaxation rate, the pore throat radius of above-mentioned neutral reservoir oil-containing part is distributed using formula (7) and is carried out Correction, obtains the NMR spectrum for being full of water of above-mentioned neutral reservoir oil-containing part, and above-mentioned formula (7) is:
In above-mentioned formula (7), T1jw,2jwFor the NMR spectrum for being full of water of above-mentioned neutral reservoir oil-containing part;RoIt is upper The pore throat radius distribution of the oil-containing part of coring sample is stated, unit is μm;ρwFor above-mentioned water surface relaxation rate, unit 10-3μm/ ms;A is constant, and when above-mentioned pore shape is respectively cylinder, spherical shape, shape channel-shaped, a is respectively 2,3,1.
It, will be upper in above-mentioned steps S1 using formula (8) in above-mentioned steps S5 in another preferred embodiment of the application It states aqueous spectrum and the NMR spectrum for being full of water of above-mentioned neutral reservoir oil-containing part is overlapped, obtain above-mentioned neutral reservoir It is full of the NMR spectrum of water, above-mentioned formula (8) is:
T1tw,2tw=T1w,2w∪T1jw,2jw……………….(8)
In above-mentioned formula (8), T1tw,2twIt is full of the NMR spectrum of water for above-mentioned neutral reservoir;T1w,2wIt is above-mentioned aqueous Spectrum;T1jw,2jwFor the NMR spectrum for being full of water of above-mentioned neutral reservoir oil-containing part.
In order to enable those skilled in the art clearly to understand the technical solution of the application, below with reference to implementation Example is described in detail with attached drawing.
Some fine and close oilfield of the selection Junggar Basin Permian System, the reservoir of the densification oilfield is salified lacustrine environment Fine-grained sediment, principal lithologic are clastic rock and carbonate rock transitionality rock class.The porosity of reservoir covers pressure between 6~16% Permeability is generally less than 0.1md, and for the diameter of reservoir pore space based on nanoscale, the oily saturation degree of reservoir is higher, average 70% Left and right, wetability the experimental results showed that:Reservoir wettability is mainly typical fine and close oily reservoir based on neutral wetting.Exploration Initial stage does not consider influence of the crude oil to nuclear magnetic resonance log, is directly evaluated using routine evaluations method, result show hole Larynx radius is bigger than normal, and distribution characteristics is similar to conventional reservoir, is unable to meet production the technology requirement of fine and close oily pore throat radius characterization.
Influence of the crude oil to nuclear magnetic resonance log has been fully considered according to technical solution provided by the invention, using oil-containing portion The hole divided is full of water NMR spectrum and is corrected to the aqueous spectrum of reservoir, and specific flow is divided into five steps, respectively For:
The first step, the petrophysics experiment of application system establish the water-oil separating mould of the NMR spectrum of area's reservoir The NMR spectrum of reservoir is divided into aqueous spectrum and is composed with oil-containing by type.
Choose polylith sealing core drilling sample, carry out reservoir be full of water-oil phase nuclear magnetic resonance (i.e. reservoir nuclear-magnetism be total to Vibration wave is composed) experiment, after measuring its oil saturation, then sample is cleaned, measures the NMR spectrum that sample is full of water.
The lateral relaxation time spectrum that this research area is all made of, so the NMR spectrum in this example is T2Spectrum.
The experimental results showed that as shown in Fig. 2, the nuclear magnetic resonance T for being full of water-oil phase of reservoir2It is full of after wave spectrum and washing oil The NMR spectrum of water is compared, and reservoir is full of the nuclear magnetic resonance T of water-oil phase2Wave spectrum is longer, and then demonstrates crude oil to T2Spectrum Large effect can be generated.
Utilize oil saturation analysis data and nuclear magnetic resonance T2Spectral data binding analysis analyzes different time to count meters Relative error between the oil saturation that the oil saturation of calculation and experiment obtain, time to count when relative error minimum is just For the disengaging time of NMR spectrum water-oil separating model, i.e. lateral relaxation time threshold value.
Specifically determining the disengaging time of NMR spectrum water-oil separating model using iterative method, root mean square is calculated T2 values when error minimum are threshold value A T2, and calculation formula is as follows:
Wherein, WT2 (j) is that the side of j-th of iteration T2 threshold value calculates error, and n is of oil saturation experimental data Number, SOi are the saturation degree measurement data of i-th of sampling point, SSOjiFor i-th of calculating saturation degree of j-th of iteration T2 threshold value, meter It is as follows to calculate formula:
Wherein, AT (j) is that oil saturation starts nuclear magnetic resonance lateral relaxation time, unit ms;ATS is active porosity The nuclear magnetic resonance transverse relaxation time to count of degree;ATD is that the nuclear magnetic resonance transverse relaxation of effecive porosity terminates the time.
Error analysis the result shows that, as shown in figure 3, when time to count is 6ms, oil-containing oil saturation degree error is minimum, because This, using 6ms as the disengaging time of oil spectrum and water spectrum in reservoir NMR spectrum as shown in fig. 4 a, i.e., under crude oil preservation Time value is limited, the separation of aqueous spectrum and oil-containing spectrum can be completed with the value, obtains Fig. 4 b aqueous spectrum and Fig. 4 c institutes with shown in The oil-containing spectrum shown.
Second, the petrophysics experiment of application system establishes the computation model of surface relaxivity respectively, calculates neutral reservoir Hole water surface relaxation rate and surface relaxivity of the oil meter face relaxation rate Liao Gai area's reservoirs containing water voids and oil-containing hole table Face relaxation rate computation model.
Step A obtains the T for being full of water-oil phase of sealing core drilling sample2It composes, the nuclear magnetic resonance T of water is full of after washing oil2Spectrum It is distributed with total pore throat radius.
Water surface relaxation rate is the surface relaxivity containing water voids, and oil meter face relaxation rate is oil-containing pore surface relaxation rate. In order to establish water surface relaxation rate and computation model containing surface relaxivity, the whole district has selected the representational sealing core drilling sample of polylith Product, successively carried out sealing core drilling sample is full of the nuclear magnetic resonance wave for being full of water after water-oil phase Nuclear Magnetic Resonance Measurement, washing oil Spectrometry and high-pressure mercury experiment, obtain the T for being full of water-oil phase of sealing core drilling sample2It composes, the core of water is full of after washing oil Magnetic resonance T2Spectrum is distributed with total pore throat radius.
For the depth segment of not sealing core drilling, washing oil first is carried out to non-closed sample, measures being full of after obtaining washing oil The nuclear magnetic resonance T of water2Spectrum;Then, after carrying out pressure mercury translocation experiment, total pore throat radius distribution is obtained;Finally to non-sealing core drilling Sample carries out the core Location of Centimeter Level, obtains real well and corresponds to the nuclear magnetic resonance log of depth and compose, as being full of grease The T of two-phase2Spectrum.
Step B calculates above-mentioned water meter according to the NMR spectrum and capillary pressure curve for being full of water after above-mentioned washing oil Face relaxation rate.
Calculating process after the washing oil of the capillary pressure curve and Fig. 5 b of Fig. 5 a as shown in figure 5, be full of the nuclear magnetic resonance of water When wave spectrum overlaps, as shown in Figure 5 c, water surface relaxation rate is obtained using following formula:
In formula, ρwFor water surface relaxation rate, unit is μm/ms;θ is the contact angle of mercury, is 140 ° in this example;γ is mercury Surface tension, be 7.15 × 10 in this example-3psi·cm;A is constant, is 3 in this example;PtotalIt is single for capillary pressure Position is Mpa;T2twFor the nuclear magnetic resonance lateral relaxation time for being full of water after washing oil, unit ms.
Surface relaxivity need to divide the water surface relaxation rate that lithology calculates different lithology, and it is thin that this example calculation obtains dolomite dust The water surface relaxation rate ρ of sandstonew=2.5 × 10-3μm/ms。
Coring sample is full of grease two by step C, the water-oil separating model using the NMR spectrum of area's reservoir The NMR spectrum of phase is divided into the aqueous spectrum of coring sample and is composed with coring sample oil-containing.The cross that 6ms is detached as water-oil phase To relaxation time threshold value, by T2Part of the time more than 6ms is that the oil-containing of coring sample is composed, by T2Part of the time less than 6ms is The aqueous spectrum of coring sample.
Water surface relaxation rate is applied to step C and obtained in aqueous spectrum as shown in Figure 6 a, is calculated above-mentioned by step D The pore throat radius distribution containing water section of coring sample, as shown in Figure 6 b, calculation formula rw=a × ρw×T2w.In this example, For white clouds matter Extra-fine sand rock, ρw=2.5 × 10-3μm/ms, a=3.
Step E is distributed, meter according to above-mentioned total pore throat radius distribution and the pore throat radius containing water section of above-mentioned coring sample Count stating the pore throat radius distribution of the oil-containing part of coring sample in, calculation formula isWherein, roFor above-mentioned coring The pore throat radius of the oil-containing part of sample is distributed, rwFor the pore throat radius containing water section point of above-mentioned coring sample, rtotalIt is upper Total pore throat radius distribution is stated, total pore throat radius distribution shown in Fig. 7 a, which is subtracted aqueous pore throat radius shown in Fig. 7 b, to be distributed, and is obtained Pore throat radius to oil-containing part shown in Fig. 7 c is distributed;
Step F is composed according to the distribution of the pore throat radius of the oil-containing part of above-mentioned coring sample and the oil-containing of above-mentioned coring sample, Above-mentioned oil meter face relaxation rate is calculated using formula (5), above-mentioned formula (5) is
In above-mentioned formula, ρoFor above-mentioned oil meter face relaxation rate, θ is the contact angle of above-mentioned mercury, is 140 ° in this example, γ is above-mentioned The surface tension of mercury, this example are 7.15 × 10-3Psicm, a 3, PoAperture for above-mentioned neutral reservoir oil-containing part is corresponding Capillary pressure, T2oFor the lateral relaxation time of above-mentioned oil-containing spectrum.In this example, white clouds matter Extra-fine sand rock is calculated, ρo It is 0.25 × 10-3μm/ms。
Above-mentioned is the core experiment for cored interval from step A to step D, and water surface is obtained with core experiment data The parameter, is then applied to the core of the underground survey for the non-cored interval that subsequent step is directed to by relaxation rate and oil meter face relaxation rate Magnetic Resonance Spectrum.
Third uses above-mentioned oil meter face relaxation rate ρoTo above-mentioned T2Spectrum oil-containing peak is handled to obtain above-mentioned reservoir oil-containing Partial pore throat radius distribution, specific formula for calculation are as follows
Ro=a × ρo×T2o
R in above formulaoFor the pore diameter of reservoir oil-containing part, unit is μm;A is 3;T2oLateral for above-mentioned oil-containing spectrum is relaxed Henan time, unit ms, ρo=0.25 × 10-3 μm/ms.
4th, the pore throat radius distribution of above-mentioned reservoir oil-containing part is corrected using water surface relaxation rate, is obtained The NMR spectrum for being full of water of the oil-containing part of above-mentioned neutral reservoir after correction, specific formula for calculation are as follows:
In above formula, T1jw,2jwThe aqueous spectrum obtained after oil-containing corrects is composed for reservoir oil-containing;RoFor the hole of reservoir oil-containing part Larynx radius, unit are μm;ρw=2.5 × 10-3μm/ms;A is 3.
5th, using formula T2tw=T2w∪T2jw, by T2jwWith T2wIt is overlapped, obtains the nuclear-magnetism that neutral reservoir is full of water Resonate wave spectrum.
Fig. 8 is finally obtained test result figure, is divided into seven regions in figure, is gamma ray curve in first area (GR), second area is the coordinate of depth, and third region is original nuclear magnetic resonance T2It composes (i.e. original well logging spectrum), in the fourth region It is the nuclear magnetic resonance T for being full of water2It composes (i.e. aqueous spectrum completely), is that calculated pore throat radius distribution curve is (thin in the 5th region Line) the analysis pore throat radius distribution curve (thick line) that obtains with experimental analysis, the 6th region is to calculate capillary pressure curve (filament) It is duty pressure and duty pressure in analysis in calculating in the 7th region with the obtained capillary pressure curve (thick line) of pressure mercury experimental analysis (dot).
Fig. 9 shows that the corresponding calculated pore throat radius distribution curve of a depth value in the 5th region in Fig. 8 is (empty Line) the analysis pore throat radius distribution curve (solid line) that obtains with experimental analysis.
As shown in Figure 9, the neutral reservoir after the correction obtained using the application is full of the NMR spectrum calculating of water Pore throat radius, which is distributed the pore throat radius distribution obtained with experiment, has preferable consistency.
It can be seen from the above description that the above embodiments of the present invention realize following technique effect:
The bearing calibration of the application has fully considered influence of the oiliness to nuclear magnetic resonance log of neutral reservoir, in The hole of the oil-containing part of property reservoir is full of water NMR spectrum, is corrected to the aqueous spectrum of neutral reservoir, after correction Wave spectrum is the NMR spectrum for being full of water of neutral reservoir, which can correctly characterize the pore throat radius point of neutral reservoir Cloth so that those skilled in the art can accurately analyze neutral reservoir, and then carry out rational exploration and development, be conducive to The development of petroleum industry.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, any made by repair Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.

Claims (7)

1. a kind of bearing calibration of nuclear magnetic resonance log wave spectrum, which is characterized in that the bearing calibration includes:
The NMR spectrum of neutral reservoir is divided into aqueous spectrum and oil-containing by step S1 using NMR spectrum disjunctive model Spectrum;
Step S2 calculates the water surface relaxation rate and oil meter face of the hole of the neutral reservoir using surface relaxivity computation model Relaxation rate;
Step S3, the pore throat radius that the oil-containing part of the neutral reservoir is calculated using the oil meter face relaxation rate are distributed;
Step S4 is corrected the pore throat radius distribution of the oil-containing part according to the water surface relaxation rate, is corrected The NMR spectrum for being full of water of the oil-containing part of the neutral reservoir afterwards;And
Step S5 is overlapped the aqueous spectrum and the NMR spectrum for being full of water of the oil-containing part of the correction, obtains The NMR spectrum for being full of water of the neutral reservoir after to correction,
The step S2 includes:
Step A chooses coring sample, measures the nuclear magnetic resonance for being full of water after the NMR spectrum of the coring sample, washing oil Total pore throat radius of wave spectrum and the neutral reservoir is distributed;
Step B calculates institute according to the NMR spectrum and capillary pressure curve for being full of water after the washing oil using formula (2) Water surface relaxation rate is stated, the formula (2) is:
In the formula (2):
ρwFor the water surface relaxation rate,
θ is the contact angle of mercury,
γ is the surface tension of mercury,
PtotalFor capillary pressure,
T1tw,2twFor the nuclear magnetic resonance lateral relaxation time T for being full of water after the washing oil2twOr longitudinal relaxation time T1tw, unit For ms;
The NMR spectrum for being full of grease of the coring sample is divided into the aqueous spectrum of coring sample and coring sample by step C Oil-containing is composed;
Step D calculates the coring sample according to the water surface relaxation rate and the aqueous spectrum of coring sample using formula (3) The pore throat radius distribution containing water section, the formula (3) is:
rw=a × ρw×T1w,2w……………….(3)
In the formula (3):
rwIt is distributed for the pore throat radius containing water section of the coring sample,
ρwFor the water surface relaxation rate,
A is constant, and when the pore shape is respectively cylinder, spherical shape, shape channel-shaped, a is respectively 2,3,1,
T1w,2wFor the lateral relaxation time T of the aqueous spectrum of coring sample2wOr longitudinal relaxation time T1w
Step E is distributed according to total pore throat radius distribution and the pore throat radius containing water section of the coring sample, using formula (4) the pore throat radius distribution of the oil-containing part of the coring sample is calculated, the formula (4) is:
In the formula (4):
roPore throat radius for the oil-containing part of the coring sample is distributed,
rwIt is distributed for the pore throat radius containing water section of the coring sample,
rtotalIt is distributed for total pore throat radius;
Step F is composed according to the distribution of the pore throat radius of the oil-containing part of the coring sample and the oil-containing of the coring sample, is used Formula (5) calculates the oil meter face relaxation rate, and the formula (5) is
In the formula (5):
ρoFor the oil meter face relaxation rate;
θ is the contact angle of the mercury,
γ is the surface tension of the mercury,
A is constant, and when the pore shape is respectively cylinder, spherical shape, shape channel-shaped, a is respectively 2,3,1,
PoFor the corresponding capillary pressure in aperture of the neutral reservoir oil-containing part, according to the oil-containing part of the coring sample Pore throat radius distribution is calculated,
T1o,2oFor the lateral relaxation time T of oil-containing spectrum2oOr longitudinal relaxation time T1o,
The step S3 includes:
It is composed according to the oil-containing in the oil meter face relaxation rate and the step S1, the neutral reservoir is calculated using formula (6) The pore throat radius of oil-containing part is distributed, and the formula (6) is:
Ro=a × ρo×T1o,2o……………….(6)
In the formula (6):
RoFor the pore throat radius distribution of the neutral reservoir oil-containing part;
ρoFor the oil meter face relaxation rate;
A is constant, and 2 are taken when the pore shape is cylindrical, takes 3 when being spherical, 1 is taken when being channel-shaped;
T1o,2oFor the lateral relaxation time T of oil-containing spectrum2oOr longitudinal relaxation time T1o,
In the step S4, according to the water surface relaxation rate, using formula (7) to the pore throat of the neutral reservoir oil-containing part Radius distribution is corrected, and obtains the NMR spectrum for being full of water of the neutral reservoir oil-containing part, and the formula (7) is:
In the formula (7):
T1jw,2jwFor the NMR spectrum for being full of water of the neutral reservoir oil-containing part;
RoFor the pore throat radius distribution of the oil-containing part of the coring sample;
ρwFor the water surface relaxation rate;
A is constant, and when the pore shape is respectively cylinder, spherical shape, shape channel-shaped, a is respectively 2,3,1.
2. bearing calibration according to claim 1, which is characterized in that the NMR spectrum disjunctive model is:
T1w,2w∪T1o,2o=T1,2……………….(1)
In formula (1):
T1,2For the NMR spectrum of the neutral reservoir;
T1w,2wFor the aqueous spectrum of the neutral reservoir;
T1o,2oFor the oil-containing spectrum of the neutral reservoir.
3. bearing calibration according to claim 1, which is characterized in that the step A includes:
Sealing core drilling sample is chosen, the NMR spectrum of the sealing core drilling sample is measured;
Washing oil is carried out to the sealing core drilling sample, measures the nuclear magnetic resonance wave for being full of water of the coring sample after washing oil Spectrum;And
Pressure mercury experiment is carried out to the sealing core drilling sample, obtains total pore throat radius distribution.
4. bearing calibration according to claim 1, which is characterized in that the step A includes:
Non- sealing core drilling sample is chosen, washing oil is carried out to the non-sealing core drilling sample, measure after washing oil described non-closed takes The NMR spectrum for being full of water of heart sample;
Pressure mercury experiment is carried out to the non-sealing core drilling sample after the washing oil, obtains total pore throat radius distribution;And
Core Location is carried out to the non-sealing core drilling sample, using practical logging NMR spectrum as the non-sealing core drilling The NMR spectrum of sample.
5. bearing calibration according to claim 4, which is characterized in that the core Location is the core Location of Centimeter Level.
6. bearing calibration according to claim 1, which is characterized in that in the step C, based on crude preservation lower limiting value The NMR spectrum of the coring sample is divided into the aqueous spectrum of coring sample and coring sample by the corresponding nuclear magnetic resonance time Oil-containing is composed.
7. bearing calibration according to claim 1, which is characterized in that in the step S5, using formula (8) by the step The aqueous spectrum and the NMR spectrum for being full of water of the neutral reservoir oil-containing part in rapid S1 are overlapped, and obtain institute The NMR spectrum that neutral reservoir is full of water is stated, the formula (8) is:
T1tw,2tw=T1w,2w∪T1jw,2jw……………….(8)
In the formula (8):
T1tw,2twIt is full of the NMR spectrum of water for the neutral reservoir;
T1w,2wFor the aqueous spectrum;
T1jw,2jwFor the NMR spectrum for being full of water of the neutral reservoir oil-containing part.
CN201510383206.1A 2015-07-02 2015-07-02 Method for correcting nuclear magnetic resonance logging spectrum Active CN106324008B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510383206.1A CN106324008B (en) 2015-07-02 2015-07-02 Method for correcting nuclear magnetic resonance logging spectrum

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510383206.1A CN106324008B (en) 2015-07-02 2015-07-02 Method for correcting nuclear magnetic resonance logging spectrum

Publications (2)

Publication Number Publication Date
CN106324008A CN106324008A (en) 2017-01-11
CN106324008B true CN106324008B (en) 2018-10-16

Family

ID=57726343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510383206.1A Active CN106324008B (en) 2015-07-02 2015-07-02 Method for correcting nuclear magnetic resonance logging spectrum

Country Status (1)

Country Link
CN (1) CN106324008B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108593697A (en) * 2018-05-09 2018-09-28 中国地质大学(北京) A kind of method and apparatus measuring Coal Under High Pressure wetability using low-field nuclear magnetic resonance
CN110082382A (en) * 2019-05-08 2019-08-02 中国科学院武汉物理与数学研究所 A method of the separation of water-oil phase NMR signal is carried out using dynamical nuclear polarization
CN110410058B (en) * 2019-06-20 2023-01-31 中国石油化工集团有限公司 Method for correcting core experiment result scale two-dimensional nuclear magnetic logging
CN113027421B (en) * 2021-04-30 2024-04-02 中海石油(中国)有限公司上海分公司 Nuclear magnetic resonance logging T2 spectrum morphology correction method in oil-based drilling fluid environment

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5696448A (en) * 1995-06-26 1997-12-09 Numar Corporation NMR system and method for formation evaluation using diffusion and relaxation log measurements
EP1073916A1 (en) * 1998-04-22 2001-02-07 Southwest Research Institute Porosity and permeability measurement of underground formations containing crude oil, using epr response data
ES2327912T3 (en) * 2002-03-07 2009-11-05 Carnegie Mellon University CONTRACT AGENTS FOR THE FORMATION OF MAGNETIC RESONANCE IMAGES AND METHODS RELATED TO THEMSELVES.
WO2004017095A2 (en) * 2002-08-14 2004-02-26 Jp Laboratories, Inc. Thick radiation sensitive devices
GB0618514D0 (en) * 2006-09-20 2006-11-01 Univ Nottingham Trent Method of detecting interactions on a microarray using nuclear magnetic resonance
CN103257151B (en) * 2013-04-24 2016-03-23 西安石油大学 In a kind of quantitative evaluation oil and gas secondary migration process, pore throat employs the method for rule
CN103513285B (en) * 2013-09-27 2016-01-06 中国石油天然气股份有限公司 Method and device for determining transverse surface relaxation rate
CN103926267B (en) * 2014-04-28 2016-05-04 西安石油大学 A kind of method of pore throat intensity of variation in quantitative assessment stress sensitive process

Also Published As

Publication number Publication date
CN106324008A (en) 2017-01-11

Similar Documents

Publication Publication Date Title
Meng et al. A critical review on fundamental mechanisms of spontaneous imbibition and the impact of boundary condition, fluid viscosity and wettability
CN104847340B (en) A kind of logging in water flooded layer quantitative evaluation method
CN106324008B (en) Method for correcting nuclear magnetic resonance logging spectrum
CN105804734B (en) A method of identifying thick oil reservoir using nuclear magnetic resonance log
Nazemi et al. The effect of carbonate reservoir heterogeneity on Archie's exponents (a and m), an example from Kangan and Dalan gas formations in the central Persian Gulf
CN110231272A (en) Tight sand aperture and nuclear magnetic resonance T2It is worth the determination method and system of transformational relation
CN106050225A (en) Method for determining 100% pure water spectrum through nuclear magnetic resonance logging spectrum T2
CN106154343B (en) Method for calculating oil saturation of tight oil reservoir
CN104358565A (en) Determining method of gas saturation in gas reservoir
CN106897531B (en) Quantitative evaluation method for permeability of low-permeability limestone reservoir
CN104213899A (en) Logging identifying method of stratum rock framework
CN109459795B (en) Method for evaluating sealing property of porous sandstone fracture zone
CN103630939A (en) Air layer identification and evaluation method
Askari et al. A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir
CN112363242A (en) Reservoir fluid identification method and device based on logging fusion
CN111812736A (en) Method for evaluating gas content of compact sandstone anhydrous gas reservoir
Kumar et al. Mapping 3D pore scale fluid distributions: how rock resistivity is influenced by wettability and saturation history
CN112145165A (en) Dynamic and static permeability conversion method for micro-crack-pore type reservoir
CN105628580B (en) Reservoir baste identifies and physical property bearing calibration
CN106351652A (en) Shape correcting method for nuclear magnetic resonance logging T2 spectrum containing hydrocarbon reservoir layers
CN106468169B (en) The method for calculating oil reservoir oil saturation
CN205297569U (en) Device for determining saturation index of tight sandstone
CN112115592A (en) Large-inclination/horizontal well resistivity correction method based on rock physics experiment
CN104748906B (en) carbonate rock capillary pressure determination method and device
CN114060015B (en) Method and device for evaluating gas content of compact sandstone

Legal Events

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