CN106154343A - Method for calculating oil saturation of tight oil reservoir - Google Patents
Method for calculating oil saturation of tight oil reservoir Download PDFInfo
- Publication number
- CN106154343A CN106154343A CN201510202590.0A CN201510202590A CN106154343A CN 106154343 A CN106154343 A CN 106154343A CN 201510202590 A CN201510202590 A CN 201510202590A CN 106154343 A CN106154343 A CN 106154343A
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- nuclear magnetic
- saturation
- oil
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 88
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 117
- 238000001228 spectrum Methods 0.000 claims abstract description 62
- 239000011148 porous material Substances 0.000 claims abstract description 15
- 230000007935 neutral effect Effects 0.000 claims abstract description 8
- 238000009736 wetting Methods 0.000 claims abstract description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 30
- 239000011435 rock Substances 0.000 claims description 21
- 238000004458 analytical method Methods 0.000 claims description 20
- 238000005553 drilling Methods 0.000 claims description 14
- 238000007789 sealing Methods 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 10
- 230000005311 nuclear magnetism Effects 0.000 claims description 7
- 238000003024 molecular redistribution determination method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 229920006395 saturated elastomer Polymers 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 2
- 239000003921 oil Substances 0.000 description 117
- 239000000523 sample Substances 0.000 description 21
- 238000005406 washing Methods 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 10
- 238000000280 densification Methods 0.000 description 7
- 230000035699 permeability Effects 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 239000006101 laboratory sample Substances 0.000 description 5
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 4
- 239000007789 gas Substances 0.000 description 4
- 239000001307 helium Substances 0.000 description 4
- 229910052734 helium Inorganic materials 0.000 description 4
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 239000010779 crude oil Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000008398 formation water Substances 0.000 description 3
- 230000005251 gamma ray Effects 0.000 description 3
- 238000004611 spectroscopical analysis Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005138 cryopreservation Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000001704 evaporation Methods 0.000 description 2
- 230000008020 evaporation Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 210000000867 larynx Anatomy 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000005945 translocation Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 229910052925 anhydrite Inorganic materials 0.000 description 1
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000706 filtrate Substances 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 238000009715 pressure infiltration Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for calculating oil saturation of a compact oil reservoir. The method comprises the following steps: collecting data by using a nuclear magnetic resonance logging instrument, and obtaining transverse relaxation time T through data processing2A spectrum; determining a transverse relaxation time threshold for the tight oil reservoir that is hydrophilic or neutrally wetted, corresponding to a lower limit of the pore throat radius; and continuously calculating the oil saturation of the compact oil reservoir according to the transverse relaxation time T2 spectrum and the transverse relaxation time threshold value. The method is based on nuclear magnetic resonance logging data, and forms a brand new transverse relaxation time T applying nuclear magnetic resonance by determining the nuclear magnetic resonance transverse relaxation minimum time threshold corresponding to the lower limit of the oil-containing pore throat radius of a compact oil reservoir under the condition of hydrophilic or neutral wetting2A method for directly calculating the oil saturation of a compact oil reservoir by a spectrum.
Description
Technical field
The present invention relates to petroleum exploration domain, in particular to a kind of method of oil saturation calculating fine and close oil reservoir.
Background technology
Oil saturation is the important parameter that reservoir oiliness characterizes.For conventional reservoir, obtain reservoir oil-containing by the method for well logging
The method of saturation generally has electrical log, nuclear magnetic resonance log and C/O to compare spectrometry logging.For fine and close oil reservoir, above-mentioned survey
The application of well method is all restricted.
Wherein, electrical log calculates saturation to be needed and porosity logging use in conjunction, application following formula calculating oil saturation:In formula: Φ is the porosity value that porosity logging obtains, decimal；A, b, m, n are rock electricity
Experiment parameter, dimensionless；RT is the formation resistivity that electrical log obtains, Ω .m；Rw is formation water resistivity, Ω .m.
But, for fine and close oil, water is not the most produced on stratum, it is impossible to obtain formation water resistivity data, thus cannot be the most direct
Above-mentioned formula is applied to calculate oil saturation.C/O is to have higher porosity than the condition that spectrometry logging is applied, and could obtain noise
The highest well-log information, fine and close oil reservoir porosity is the most relatively low, and the method is the most inapplicable.
Application NMR logging data calculates oil saturation generally two kinds of methods: a kind of is by the relaxation of nuclear magnetic resonance log
Wave spectrum is converted to pseudo-capillary pressure curve, and in the case of known oil-water interfaces and profit density contrast, the oil-containing calculating every bit is satisfied
And degree.The achievement in research of domestic this aspect is more, obvious in clastic reservoir rock and volcanic rock reservoir application effect.Fault Lai training etc.,
2003 (southwest petroleum journals) have studied application nuclear magnetic resonance T 2 spectrum and are converted to pseudo-capillary pressure curve, thus calculate saturation
Method；From the beginning of 2000, rectify the beginning of spring, Mao Zhiqiang, Sun Zhongchun, He Dan, Xiao Liang etc. combine Junggar Basin low-resistance, low hole,
Low permeability sandstone reservoirs and the feature of volcanic rock reservoir, carried out the method for application NMR logging data conversion capillary pressure curve
Research, it is proposed that the method for segmentation power exponent conversion, and obtain national inventing patent (2012), in clastic rock, volcanic rock storage
Layer saturation computation has seen preferable geological effect.Another kind of method assumes that reservoir height is sufficiently large, surveys with nuclear magnetic resonance, NMR
The irreducible water saturation estimation oil saturation that well obtains.
But, compact oil reservoir has source storage one, and non-buoyancy becomes the feature hidden, do not has clear and definite oil-water interfaces, first method
Cannot directly apply.Owing to fine and close oil reservoir is based on nano-pore, cannot flow under formation conditions, fluid almost all is
Constraint fluid, second method also cannot use.It is to say, under existing technical conditions, it is to cause that oil saturation calculates
The technical barrier of close oil reservoir characterization and the key issue needing solution badly.
Summary of the invention
Present invention is primarily targeted at a kind of method that oil saturation calculating fine and close oil reservoir is provided.
To achieve these goals, according to an aspect of the invention, it is provided a kind of calculating is hydrophilic or neutral wetting fine and close oil storage
The method of the oil saturation of layer, comprises the following steps: application nuclear magnetic resonance logging instrument gathers data, and obtains through data process
Obtain T2 T2Spectrum；Determine that hydrophilic or neutral wetting fine and close oil reservoir relaxes corresponding to the horizontal of oil-containing pore throat radius lower limit
Henan time threshold；According to T2 T2Spectrum and the oil saturation of T2 threshold value Continuous plus densification oil reservoir.
Further, it is thus achieved that T2 T2The method of spectrum includes: utilizes the minimum echo interval of nuclear magnetic resonance analyser to gather and causes
The log data of close oil reservoir, and utilize log data to obtain the T2 T of proton2Spectrum.
Further, selecting two kinds of methods to determine T2 threshold value, first method is the hydrogen-free desaturation of core sample
With nuclear magnetic resonance, NMR translocation method, second method is that sealing core drilling analyzes oil saturation data and nuclear magnetic resonance log T2Threshold value iterative method.
Further, the hydrogen-free desaturation of core sample includes with the step of nuclear magnetic resonance, NMR translocation method: water conservation hydrogen-free rock core is processed；
Nuclear Magnetic Resonance Measurement during water conservation hydrogen-free desaturation；The identification of water spectrum and T2The determination of threshold value.
Further, sealing core drilling analyzes oil saturation data and nuclear magnetic resonance log T2The step of threshold value iterative method includes: utilize every
Rice combines playback not less than what 3 data lacunarity analysis data carried out saturation analysis data；Utilize micro resistor
Data carries out playback fine setting to the rock core of fine and close oil reservoir；Use iterative method to determine T2 threshold value, calculate mean square error
T time minimum2Value is saturation computation threshold value A T2, and its computing formula is as follows:Wherein, WT2 (j) is jth iteration T2The mean square calculating error of threshold value, n
For the number of oil saturation experimental data, SOiFor the saturation measurement data of i-th sampling point, SSOjiFor jth iteration T2
The i-th of threshold value calculates saturation, and its computing formula is as follows:
Wherein AT (j) is that oil saturation starts at nuclear magnetic resonance, NMR T2,
ms；ATS is the nuclear magnetic resonance, NMR transverse relaxation time to count of effecive porosity；ATD is that the nuclear magnetic resonance, NMR of effecive porosity is horizontal
The time is terminated to relaxation.
Further, first method and second method use respectively；Or first method and second method apply simultaneously.
Further, according to T2 T2The oil-containing of spectrum and T2 threshold calculations densification oil reservoir is saturated
The formula of degree is:
Application technical scheme, the method, based on nuclear magnetic resonance log data, is determined by fine and close oil reservoir parent
In the case of water or neutral wetting, the nuclear magnetic resonance, NMR transverse relaxation minimum time threshold value that its oil-containing pore throat radius lower limit is corresponding, thus shape
Become the T2 T of a kind of brand-new application nuclear magnetic resonance, NMR2Spectrum directly calculates the method for fine and close oil reservoir oil saturation.
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:
The schematic flow sheet of the method for the oil saturation calculating fine and close oil reservoir that Fig. 1 is provided by the application embodiment；
Fig. 2 is the hydrogen-free desaturation experimental result with the acquisition of nuclear magnetic resonance, NMR translocation method of one piece of core sample；
Fig. 3 is the playback figure of a bite well sealing core drilling saturation analysis data well section；
Fig. 4 is the Analysis of Mean Square Error figure that a bite well sealed coring well section applies different AT2 value iterative computation；
Fig. 5 is that the AT2 and nuclear magnetic resonance log T determined applied by a bite well2The oil saturation figure that spectrum calculates；And
Fig. 6 is that 59 groups of 3 mouthfuls of wells calculate oil saturation and experimental analysis saturation error analysis figure.
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 application below with reference to the accompanying drawings and in conjunction with the embodiments in detail.
It should be noted that term used herein above merely to describe detailed description of the invention, and be not intended to restricted root according to this Shen
Illustrative embodiments please.As used herein, unless the context clearly indicates otherwise, otherwise singulative is also intended to
Including plural form, additionally, it should be understood that, when using term " to comprise " in this manual and/or time " including ",
It indicates existing characteristics, step, operation, device, assembly and/or combinations thereof.
The present invention, on the basis of rock core desaturation nuclear magnetic resonance experiment data and continuous print sealing core data are analyzed, establishes really
The nuclear magnetic resonance, NMR transverse relaxation minimum time threshold value that fixed fine and close oil reservoir minimum oil-containing pore throat radius is corresponding, thus define a kind of complete
New application nuclear magnetic resonance, NMR T2 wave spectrum directly calculates the method for fine and close oil reservoir oil saturation.As it is shown in figure 1,
The method comprises the following steps: application nuclear magnetic resonance logging instrument gathers data, and processes acquisition T2 T through data2
Spectrum；Determine that hydrophilic or neutral wetting fine and close oil reservoir corresponds to the T2 threshold value of oil-containing pore throat radius lower limit；According to
T2 T2Spectrum and the oil saturation of T2 threshold value Continuous plus densification oil reservoir.
The present invention is on the basis of the log data of fine and close oil reservoir, and the nuclear-magnetism that fine and close oil reservoir minimum oil-containing pore throat radius is corresponding is total to
Shake transverse relaxation minimum time threshold value, thus define the T2 T of a kind of brand-new application nuclear magnetic resonance, NMR2Spectrum directly meter
The method calculating fine and close oil reservoir oil saturation.
The exemplary reality of the method for the oil saturation that calculate fine and close oil reservoir according to present invention offer is described in more detail below
Execute mode.But, these illustrative embodiments can be implemented by multiple different form, and is not construed as only
It is limited to embodiments set forth herein.It should be appreciated that these embodiments are provided so that disclosure herein is thorough
The end and complete, and the design of these illustrative embodiments is fully conveyed to those of ordinary skill in the art.
First, application nuclear magnetic resonance logging instrument gathers data, and processes acquisition T2 T through data2Spectrum.This step
In, the feature grown according to fine and close oil reservoir nanoaperture, the general minimum echo using instrument to provide interval gathers, and
Control is tested the speed, and improves the signal to noise ratio of well-log information, obtains the T2 T of proton with well-log information2Spectrum (i.e. T2Wave spectrum).
Then, it is determined that hydrophilic or neutral wetting fine and close oil reservoir corresponds to the T2 threshold value of oil-containing pore throat radius lower limit.
The invention provides two kinds of methods: a kind of hydrogen-free desaturation for core sample and nuclear magnetic resonance, NMR translocation method.Another kind is close
Close core analysis oil saturation data and nuclear magnetic resonance log T2Threshold value iterative method.The former is applicable to the feelings without system sealing core data
Condition, cores can complete by routine, and the latter is applicable to the situation having system sealing core drilling.Two kinds of methods also can be applied simultaneously, phase
Confirmation mutually, to improve saturation T2The determination precision of threshold value.
Wherein, the hydrogen-free desaturation of core sample includes with nuclear magnetic resonance, NMR translocation method experimental procedure: 1., water conservation hydrogen-free rock core processing；
2., the Nuclear Magnetic Resonance Measurement during the desaturation of water conservation hydrogen-free；What 3., water was composed identifies the determination with T2 threshold value.The core of this experiment
The heart is that water conservation obtains substantially without the spectrum of the NMR water containing oil samples under the conditions of hydrocarbon.The feature of fine and close oil reservoir is with nanoaperture
Being main, oil saturation is higher, composes containing the water that content in oil samples is relatively low and is stored in the interstitial space that pore diameter is less, relatively
It is easier to preserve.Further, since reservoir samples is fine and close, during coring, mud filtrate invasion is shallower, easily obtains representational sample
Product.
Drilling through of laboratory sample: the full-hole core cryopreservation at a temperature of 30 DEG C obtained of coring, chooses representational rock core
Drill through the sample of 1 inch with liquid nitrogen, removing two ends has the part of certain pollution, and middle 4cm sample is as laboratory sample, it is ensured that
The anhydrous intrusion of specimen sample, evaporation.
Measure the NMR (Nuclear Magnetic Resonance) spectrum that the remaining oil before non-washing oil, water are biphase, spectral line in Fig. 2.What deserves to be explained is, due to
Sample is that routine is cored acquisition, and in this case, the light components of crude oil has scattered and disappeared, but complete aqueous spectrum is obtained by this
Must affect little.
Water conservation hydrogen-free desaturation and NMR (Nuclear Magnetic Resonance) spectrum are measured: fine and close oil reservoir samples permeability is extremely low, and overburden permeability is the least
In 0.1md, cannot realize with the method desaturation of displacement.Use the method desaturation of carbon dioxide washing oil, first wash away in big pore throat
Oil gas, wash away the oil gas in smaller hole larynx again, the corresponding NMR (Nuclear Magnetic Resonance) spectrum measuring core sample in this course, until
Substantially the hydro carbons in core sample is washed away, it is thus achieved that be substantially free of the residue water spectrum of hydro carbons, in Fig. 2 × spectral line of symbol.
The identification of water spectrum and T2The determination of threshold value: analyze the feature of residue water spectrum, determine aqueous volume, obtain oil saturation meter
The T calculated2Threshold value.T is obtained from remaining complete aqueous spectrum2Threshold value A T2For 6ms.
In the another kind of preferred implementation of the present invention, sealing core drilling is used to analyze oil saturation data and nuclear magnetic resonance log T2Threshold
1., the accurate playback of saturation data the step of value iterative method includes:；2. when, iterative method determines saturation computation transverse relaxation
Between T2Threshold value.
The accurate playback of saturation data: first, carry out saturation analysis with every meter not less than 3 data lacunarity analysis data
Data combine playback；Then, the fine setting of the playback of rock core is carried out by micro resistor data, it is ensured that playback error is not
Big 0.1m (Fig. 3).In Fig. 3, CALi, GR and SP are lithology curve, are used for dividing reservoir and non-reservoir, wherein CALi
For CAL, GR is Natural Gamma-ray Logging Curves, and SP is nutural potential logging curve；RI, RXO and RT are resistivity
Curve, is used for judging the oiliness of reservoir, and wherein RXO is flushed zone formation resistivity；DEN, CNL and AC are porosity
Curve, is used for calculating the size of reservoir porosity, and wherein DEN is density curve；T2 is NMR (Nuclear Magnetic Resonance) spectrum.
Sealing core drilling analyzes oil saturation data and nuclear magnetic resonance log T2The step of threshold value iterative method includes: utilize every meter not less than 3
What individual data lacunarity analysis data carried out saturation analysis data combines playback；Utilize micro resistor data to densification
The rock core of oil reservoir carries out playback fine setting；Use iterative method to determine T2 threshold value, calculate T during mean square error minimum2
Value is saturation computation threshold value A T2, and its computing formula is as follows: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, and SOi is
The saturation measurement data of i sampling point, SSOji is that the i-th of jth iteration T2 threshold value calculates saturation, and its computing formula is such as
Under:
Wherein AT (j) is that oil saturation starts at nuclear magnetic resonance, NMR T2,
ms；ATS is the nuclear magnetic resonance, NMR transverse relaxation time to count of effecive porosity；ATD is that the nuclear magnetic resonance, NMR of effecive porosity is horizontal
The time is terminated to relaxation.
Wherein, first method and second method use respectively；Or first method and second method apply simultaneously.
Calculate the T that mean square error is minimum2Value is the T determined2Saturation computation threshold value A T2.Fig. 4 is a bite well sealed coring well section
Apply different AT2The mean square error that value calculates, AT corresponding during mean square error minimum2For 6ms, complete with core experiment result
Complete consistent.
Finally, according to T2 wave spectrum and the oil saturation of T2 threshold value Continuous plus densification oil reservoir.
The AT that application determines2The continuous T obtained with nuclear magnetic resonance log2Wave spectrum is by the saturation of the following equation each measuring point of calculating:
Wherein, So is the oil saturation of described fine and close oil reservoir, φiIt is hole corresponding to i-th millisecond of NMR relaxation time
Relative volume, AT2 is described T2 threshold value, and ATS is the nuclear magnetic resonance, NMR transverse relaxation time to count of effecive porosity,
ms；ATD is that the nuclear magnetic resonance, NMR transverse relaxation of effecive porosity terminates time, ms.
Fig. 5 is the continuous T that a bite well applies that the AT2 determined and nuclear magnetic resonance log obtain2Wave spectrum calculates the example of saturation.Meter
Calculating result to contrast with measured result, the form of two kinds of results is consistent, and the relative error calculating saturation is relatively small.3 mouthfuls of wells are airtight
The saturation contrast that saturation of coring and the present invention calculate, within calculating error is mainly distributed on 10%, mean absolute error is
2.3%, computational accuracy is high, is fully able to meet the required precision that reserves calculate.In Fig. 5, CALi, GR and SP are that lithology is bent
Line, is used for dividing reservoir and non-reservoir, and wherein CALi is CAL, and GR is Natural Gamma-ray Logging Curves, and SP is nature electricity
Position log；RI, RXO and RT are resistivity curve, are used for judging the oiliness of reservoir, and wherein RXO is flushed zone ground
Layer resistivity；DEN, CNL and AC are porosity curve, are used for calculating the size of reservoir porosity, and wherein DEN is density
Curve；T2 is NMR (Nuclear Magnetic Resonance) spectrum.
As a example by a fine and close oilfield, the techniqueflow of the present invention will be described in detail below, but it is not intended that to this
Bright can the restriction of practical range.
Illustrate that block is one fine and close oilfield of the Junggar Basin Permian System.Fine and close oil reservoir is the fine-grained sediment of salified lacustrine environment,
Principal lithologic is clastic rock and carbonate rock transitionality rock class.The porosity of reservoir is mainly distributed between 6～16%, covers pressure infiltration
Rate is generally less than 0.1md, for typical fine and close oil reservoir.The oil saturation of reservoir is the highest, analytical data oil saturation
Mostly higher than 70%.Whole district's block many wells formation testing, without natural production capacity, has higher yield after volume fracturing.The all fine and close oil of the whole district
Block Pressure Curve in Oil Testing Well Duan Jun does not produce water, calculates saturation by resistivity method and cannot obtain formation water resistivity data.The all densifications of the whole district
Oil well has all carried out nuclear magnetic resonance log.The nuclear magnetic resonance log saturation computation method of exploration initial stage application conventional reservoir, calculates
Saturation error bigger.Cannot meet and produce the technology requirement that upper fine and close oil oiliness characterizes.The method provided according to the present invention
Technology preferably solves the technical barrier of saturation computation.
Detailed description of the invention and implementing procedure be:
1) feature, according to fine and close oil reservoir nanoaperture grown, uses minimum echo interval to gather, and controls to test the speed, it is thus achieved that
Qualified NMR logging data, Inversion Calculation, obtain the T2 T of proton2Wave spectrum.
2), choose representational rock core and carry out hydrogen-free desaturation and nuclear magnetic resonance, NMR translocation experiment
Drilling through of laboratory sample: the full-hole core cryopreservation at a temperature of 30 DEG C obtained of coring, chooses representational rock core
Drill through the sample of 1 inch with liquid nitrogen, removing two ends has the part of certain pollution, and middle 4cm sample is as laboratory sample, it is ensured that
The anhydrous intrusion of specimen sample, evaporation.First, the NMR (Nuclear Magnetic Resonance) spectrum that the remaining oil before non-washing oil, water are biphase is measured.It is worth saying
Bright, owing to sample is that routine is cored acquisition, in this case, the light components of crude oil has scattered and disappeared, but this is to complete
The acquisition impact of complete aqueous spectrum is little.
Water conservation hydrogen-free desaturation and NMR (Nuclear Magnetic Resonance) spectrum are measured: fine and close oil reservoir samples permeability is extremely low, and overburden permeability is the least
In 0.1md, cannot realize with the method desaturation of displacement.Use the method desaturation of carbon dioxide washing oil, first wash away in big pore throat
Oil gas, wash away the oil gas in smaller hole larynx again, the corresponding NMR (Nuclear Magnetic Resonance) spectrum measuring core sample in this course, until
Substantially the hydro carbons in core sample is washed away, it is thus achieved that be substantially free of the residue water spectrum of hydro carbons.
The identification of water spectrum and T2The determination of threshold value: analyze the feature of residue water spectrum, determine aqueous volume, obtain oil saturation meter
The T calculated2Threshold value.T is obtained from remaining complete aqueous spectrum2Threshold value A T2.
This densification oilfield have chosen 4 pieces of samples to carry out hydrogen-free desaturation and nuclear magnetic resonance, NMR translocation and tests, respectively sand formation cuttings cloud rock,
Cloud bits siltstone, Anhydrite landwaste siltstone and cloud matter siltstone.Now as a example by the experiment of cloud matter siltstone, introduce experiment and determine T2Threshold
The method of value AT2.The porosity of laboratory sample is 13.4%, and overburden permeability is 0.09md, and experimental result is shown in accompanying drawing 2.Figure
In: T2Spectrum is the T before washing oil2Spectrum；◇ is washing oil T after 2 days2Spectrum；× for the water spectrum obtained after washing oil 6 days.Use note
The method of the supporting measurement of helium porosity can detect light components volatile quantity and wash away the pore volume of crude oil.Before washing oil, note helium is surveyed
Porosity be 3.4%, nuclear magnetic resonance, NMR the porosity obtained is 10.04%, it is thus achieved that NMR (Nuclear Magnetic Resonance) spectrum the oil-bearing characteristics obvious.
Hereafter, every washing oil one day, carry out nuclear-magnetism, the note supporting measurement of helium, washing oil measurement result two days later and magnetic resonance spectroscopy show,
Sample still has the oil-bearing characteristics, until washing oil is after 6 days, it is thus achieved that the most aqueous nuclear-magnetism spectrum.Now, sample is basic
Unstressed configuration, nuclear-magnetism spectrum water spectrum signature is obvious, and the measured value of note helium porosity is 10.94%, it is thus achieved that the oil saturation of sample is 82%.
The AT2 threshold value obtained from aqueous spectrum is 6ms.
3), this district systematic sealing core drilling saturation survey data, application sealing core drilling analyzes oil saturation data and nuclear-magnetism
Well logging T2Threshold value iterative method verifies core experiment result further.
The accurate playback of saturation data: first, carry out saturation analysis with every meter not less than 3 data lacunarity analysis data
Data combine playback；Then, the fine setting of core Location is carried out by micro resistor data, it is ensured that playback error is little
In 0.1m (Fig. 3).In Fig. 3, first is respectively CAL, Natural Gamma-ray Logging Curves and nutural potential logging curve,
Second is depth track, and the 3rd road is three Resistivity log, and the 4th road is tri-porosity logging curve, and the 5th road is nuclear-magnetism
Resonance log and core analysis porosity (scatterplot), the 6th road is the oil saturation (scatterplot) analyzed.
Iterative method determines saturation computation T2 T2Threshold value: by following equation be iterated calculate mean square error:
In formula: AT2 (j) is jth iteration T2The mean square calculating error of threshold value；
N is the number of oil saturation experimental data；
SOi is the saturation measurement data of i-th sampling point；
SSOji is jth iteration T2The i-th of value calculates saturation.
Calculate the T that mean square error is minimum2Value is the T determined2Saturation computation threshold value A T2.This mouthful of well sealed coring well section of Fig. 4 should
The mean square error calculated by different AT2 values, AT2 corresponding during mean square error minimum is 6ms, complete with core experiment result
Unanimously.
4) continuous T that the AT2 determined and nuclear magnetic resonance log obtain, is applied2Wave spectrum presses every mouthful of well of the following equation Continuous plus whole district
Saturation:
Wherein, So is the oil saturation of described fine and close oil reservoir, φiIt is hole corresponding to i-th millisecond of NMR relaxation time
Relative volume, AT2 is described T2 threshold value, and ATS is the nuclear magnetic resonance, NMR transverse relaxation time to count of effecive porosity,
ms；ATD is that the nuclear magnetic resonance, NMR transverse relaxation of effecive porosity terminates time, ms.
Fig. 5 is the continuous T that a bite well applies that the AT2 determined and nuclear magnetic resonance log obtain2Wave spectrum calculates the example of saturation.Figure
In, the 6th road is the saturation and sealing core drilling saturation (scatterplot) comparison diagram calculated, and in terms of comparing result, two kinds of methods obtain
The saturation obtained has preferable concordance, and precision is higher.In Fig. 5, first is respectively CAL, natural gamma is surveyed
Well curve and nutural potential logging curve, second is depth track, and the 3rd road is three Resistivity log, and the 4th road is three holes
Porosity log, the 5th road is that nuclear magnetic resonance, NMR surveys T2Spectrum, the 6th road is the oil saturation and the oil saturation of analysis calculated
(scatterplot).
Local area 3 mouthfuls has sealing core drilling saturation analysis data, totally 59 pieces of samples.Calculate saturation and analysis saturation error is divided
Analysis figure is shown in Fig. 6.The saturation contrast that 3 mouthfuls of well sealing core drilling saturations and the present invention calculate, calculate error be mainly distributed on 10% it
In, mean absolute error is 2.3%, and computational accuracy is high, is fully able to meet the required precision that reserves calculate.
As can be seen from the above embodiments, the example that the present invention is above-mentioned achieves following technique effect: the method is surveyed with nuclear magnetic resonance, NMR
Based on well data, being determined by fine and close oil reservoir in the case of hydrophilic or neutral wetting, its oil-containing pore throat radius lower limit is corresponding
Nuclear magnetic resonance, NMR transverse relaxation minimum time threshold value, thus define the T2 T of a kind of brand-new application nuclear magnetic resonance, NMR2Spectrum
The method directly calculating fine and close oil reservoir oil saturation.
These are only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art,
The present invention can have various modifications and variations.All within the spirit and principles in the present invention, any amendment of being made, equivalent,
Improve, should be included within the scope of the present invention.
Claims (7)
1. the method for the oil saturation calculating fine and close oil reservoir, it is characterised in that comprise the following steps:
Application nuclear magnetic resonance logging instrument gathers data, and processes acquisition T2 T through data2Spectrum；
Determine that hydrophilic or neutral wetting described fine and close oil reservoir corresponds to the T2 threshold of oil-containing pore throat radius lower limit
Value；
According to described T2 T2Described in spectrum and described T2 threshold value Continuous plus, fine and close oil reservoir contains
Oil saturation.
Method the most according to claim 1, it is characterised in that obtain described T2 T2The method of spectrum includes: profit
It is spaced the described log data gathering fine and close oil reservoir by the minimum echo of nuclear magnetic resonance analyser, and utilizes described log data to obtain
Take the described T2 T of proton2Spectrum.
Method the most according to claim 1, it is characterised in that when selecting one of following two method to determine described transverse relaxation
Between threshold value, first method is hydrogen-free desaturation and the nuclear magnetic resonance, NMR translocation method of core sample, and second method is airtight taking
The heart analyzes oil saturation data and nuclear magnetic resonance log T2Threshold value iterative method.
4. according to right, require method described in 3, it is characterised in that the hydrogen-free desaturation of described core sample joins with nuclear magnetic resonance, NMR
The step of survey method includes:
Water conservation hydrogen-free rock core is processed；
Nuclear Magnetic Resonance Measurement during water conservation hydrogen-free desaturation；
The identification of water spectrum and T2The determination of threshold value.
Method the most according to claim 3, it is characterised in that described sealing core drilling is analyzed oil saturation data and surveyed with nuclear-magnetism
Well T2The step of threshold value iterative method includes:
Every meter is utilized to combine playback not less than what 3 data lacunarity analysis data carried out saturation analysis data；
Utilize micro resistor data that the rock core of described fine and close oil reservoir is carried out playback fine setting；
Use iterative method to determine described T2 threshold value, calculate T during mean square error minimum2Value is saturation computation
Threshold value A T2, and its computing formula is as follows:
Wherein, WT2 (j) is jth iteration T2The mean square calculating error of threshold value, n is oil saturation experimental data
Number, SOiFor the saturation measurement data of i-th sampling point, SSOjiFor jth iteration T2The i-th of threshold value calculates saturated
Degree, its computing formula is as follows:
Wherein AT (j) is that oil saturation starts at nuclear magnetic resonance, NMR T2,
ms；ATS is the nuclear magnetic resonance, NMR transverse relaxation time to count of effecive porosity；ATD is that the nuclear magnetic resonance, NMR of effecive porosity is horizontal
The time is terminated to relaxation.
Method the most according to claim 3, it is characterised in that described first method and described second method use respectively；
Or described first method and described second method apply simultaneously.
Method the most according to claim 1, it is characterised in that according to described T2 T2Spectrum and described transverse relaxation
The formula of the oil saturation that time threshold calculates described fine and close oil reservoir is:
Wherein, So is the oil saturation of described fine and close oil reservoir, φiIt is that i-th millisecond of NMR relaxation time is corresponding
Hole relative volume, AT2For described T2 threshold value, ATS is that the nuclear magnetic resonance, NMR transverse relaxation of effecive porosity rises
Evaluation time；ATD is that the nuclear magnetic resonance, NMR transverse relaxation of effecive porosity terminates the time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510202590.0A CN106154343B (en) | 2015-04-24 | 2015-04-24 | Method for calculating oil saturation of tight oil reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510202590.0A CN106154343B (en) | 2015-04-24 | 2015-04-24 | Method for calculating oil saturation of tight oil reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106154343A true CN106154343A (en) | 2016-11-23 |
| CN106154343B CN106154343B (en) | 2018-07-13 |
Family
ID=57347356
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510202590.0A Active CN106154343B (en) | 2015-04-24 | 2015-04-24 | Method for calculating oil saturation of tight oil reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106154343B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108572129A (en) * | 2017-03-13 | 2018-09-25 | 中国石油化工股份有限公司 | Fine and close oil effective reservoir hole threshold value confining method and system |
| CN108896598A (en) * | 2018-04-17 | 2018-11-27 | 中国石油天然气股份有限公司 | Method and device for determining content of compact oil |
| CN109779619A (en) * | 2019-01-21 | 2019-05-21 | 中国石油天然气集团有限公司 | A method of it is composed by nuclear-magnetism T2 and calculates stratum water saturation |
| CN110487835A (en) * | 2019-09-18 | 2019-11-22 | 西南石油大学 | A kind of new method calculating fine and close oil-gas reservoir reservoir saturation exponent |
| CN111007230A (en) * | 2019-11-21 | 2020-04-14 | 中国石油天然气股份有限公司 | Method for quantitatively evaluating oil content of low-porosity compact oil reservoir of continental-phase lake basin |
| CN114414609A (en) * | 2022-01-13 | 2022-04-29 | 东北石油大学 | Experiment method for calculating influence of invasion liquid on shale oil momentum based on nuclear magnetic T2 spectrum |
| CN114660668A (en) * | 2022-02-18 | 2022-06-24 | 中国石油大学(华东) | Shale porosity overbalance correction method based on crack action |
| CN115420766A (en) * | 2022-09-06 | 2022-12-02 | 东北石油大学 | Based on T 2 Spectral analysis method for determining distribution of residual oil in sandstone water flooding oil extraction |
| CN117761102A (en) * | 2023-12-19 | 2024-03-26 | 西南石油大学 | Device and method for measuring evaporation degree of humidifying natural gas to formation water |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101892837A (en) * | 2010-04-29 | 2010-11-24 | 中国石油天然气股份有限公司 | Formation factor determining method and oil saturation determining method |
| US20130193964A1 (en) * | 2012-01-31 | 2013-08-01 | Numatex, Inc. | Method for detecting oil and gas from the surface by nuclear magnetic resonance imaging |
| CN104453876A (en) * | 2014-11-03 | 2015-03-25 | 中国石油天然气股份有限公司 | Method and device for predicting oil and gas yield of horizontal well of compact oil and gas reservoir |
| US20150094960A1 (en) * | 2013-10-02 | 2015-04-02 | Schlumberger Technology Corporation | Quality metrics for tight oil reservoirs |
-
2015
- 2015-04-24 CN CN201510202590.0A patent/CN106154343B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101892837A (en) * | 2010-04-29 | 2010-11-24 | 中国石油天然气股份有限公司 | Formation factor determining method and oil saturation determining method |
| US20130193964A1 (en) * | 2012-01-31 | 2013-08-01 | Numatex, Inc. | Method for detecting oil and gas from the surface by nuclear magnetic resonance imaging |
| US20150094960A1 (en) * | 2013-10-02 | 2015-04-02 | Schlumberger Technology Corporation | Quality metrics for tight oil reservoirs |
| CN104453876A (en) * | 2014-11-03 | 2015-03-25 | 中国石油天然气股份有限公司 | Method and device for predicting oil and gas yield of horizontal well of compact oil and gas reservoir |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108572129B (en) * | 2017-03-13 | 2021-12-17 | 中国石油化工股份有限公司 | Method and system for defining pore threshold of compact oil effective reservoir |
| CN108572129A (en) * | 2017-03-13 | 2018-09-25 | 中国石油化工股份有限公司 | Fine and close oil effective reservoir hole threshold value confining method and system |
| CN108896598A (en) * | 2018-04-17 | 2018-11-27 | 中国石油天然气股份有限公司 | Method and device for determining content of compact oil |
| CN108896598B (en) * | 2018-04-17 | 2021-05-28 | 中国石油天然气股份有限公司 | Method and device for determining content of compact oil |
| CN109779619A (en) * | 2019-01-21 | 2019-05-21 | 中国石油天然气集团有限公司 | A method of it is composed by nuclear-magnetism T2 and calculates stratum water saturation |
| CN109779619B (en) * | 2019-01-21 | 2023-02-28 | 中国石油天然气集团有限公司 | Method for calculating stratum water saturation through nuclear magnetic T2 spectrum |
| CN110487835A (en) * | 2019-09-18 | 2019-11-22 | 西南石油大学 | A kind of new method calculating fine and close oil-gas reservoir reservoir saturation exponent |
| CN111007230B (en) * | 2019-11-21 | 2022-03-29 | 中国石油天然气股份有限公司 | Method for quantitatively evaluating oil content of low-porosity compact oil reservoir of continental-phase lake basin |
| CN111007230A (en) * | 2019-11-21 | 2020-04-14 | 中国石油天然气股份有限公司 | Method for quantitatively evaluating oil content of low-porosity compact oil reservoir of continental-phase lake basin |
| CN114414609A (en) * | 2022-01-13 | 2022-04-29 | 东北石油大学 | Experiment method for calculating influence of invasion liquid on shale oil momentum based on nuclear magnetic T2 spectrum |
| CN114414609B (en) * | 2022-01-13 | 2022-11-01 | 东北石油大学 | Experiment method for calculating influence of invaded liquid on shale oil momentum based on nuclear magnetic T2 spectrum |
| CN114660668A (en) * | 2022-02-18 | 2022-06-24 | 中国石油大学(华东) | Shale porosity overbalance correction method based on crack action |
| CN115420766A (en) * | 2022-09-06 | 2022-12-02 | 东北石油大学 | Based on T 2 Spectral analysis method for determining distribution of residual oil in sandstone water flooding oil extraction |
| CN115420766B (en) * | 2022-09-06 | 2023-05-19 | 东北石油大学 | Based on T 2 Method for measuring distribution of residual oil in sandstone water injection oil extraction by spectrum analysis |
| CN117761102A (en) * | 2023-12-19 | 2024-03-26 | 西南石油大学 | Device and method for measuring evaporation degree of humidifying natural gas to formation water |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106154343B (en) | 2018-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106154343A (en) | Method for calculating oil saturation of tight oil reservoir | |
| Xiao et al. | A method to determine nuclear magnetic resonance (NMR) T2cutoff based on normal distribution simulation in tight sandstone reservoirs | |
| CN101487390B (en) | Archie mode method for confirming initial oil saturation of oil layer | |
| CN105545301B (en) | A kind of complicated reservoirs flow net model method under different mud system | |
| CN106468172A (en) | A kind of Oil in Super-low Permeability sandstone oil reservoir low-resistance reservoir log interpretation method | |
| CN104912550A (en) | Method for quantitatively calculating reservoir fluid producing profile by nuclear magnetic resonance well logging information | |
| Xiao et al. | A new technique for synthetizing capillary pressure (Pc) curves using NMR logs in tight gas sandstone reservoirs | |
| CN112255688B (en) | Method for inverting formation pressure by three-dimensional earthquake based on rock physics theory | |
| US20100185424A1 (en) | Method, Program and Computer System for Conciliating Hydrocarbon Reservoir Model Data | |
| CN104453874A (en) | Glutenite reservoir oil saturation calculation method based on nuclear magnetic resonance | |
| CN105182424A (en) | Method and device for quantitatively predicting reservoir porosity based on patch saturation model | |
| CN105093313A (en) | Predicting method and apparatus for production capacity of single well in Karst oil-gas reservoir | |
| CN104514552A (en) | Method for identification and abundance prediction of coalbed methane reservoirs | |
| CN109117505A (en) | A kind of pore structure reservoir flushed zone water saturation calculation method based on dielectric experiment | |
| CN106223942A (en) | A kind of Conglomerate Reservoir shale content computational methods based on Well logging curve reconstruction | |
| CN104714252A (en) | Method for analyzing fluid factor sensibility | |
| CN103643949A (en) | Quantitatively forecasting method and device for oil-gas possibility of reservoirs | |
| CN100523431C (en) | Method for simulating oil-water two-phase cable formation testing | |
| Liang et al. | Tight gas sandstone reservoirs evaluation from nuclear magnetic resonance (NMR) logs: case studies | |
| CN105931125A (en) | Method for predicting yield of compact oil staged multi-cluster volume fracturing horizontal well | |
| CN105240006A (en) | Oil and water layer recognition method suitable for volcanic reservoir | |
| Feng et al. | Accurate determination of water saturation in tight sandstone gas reservoirs based on optimized Gaussian process regression | |
| Al-Sulami et al. | The unconventional shale reservoirs of jafurah basin: An integrated petrophysical evaluation using cores and advanced well logs | |
| Sabea et al. | Geological model of the Khabour Reservoir for studying the gas condensate blockage effect on gas production, Akkas Gas Field, Western Iraq | |
| CN104484573B (en) | Method for determining formation rigidity coefficient |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |