CN102183532B - Method for improving computed tomography (CT) measured fluid saturation precision - Google Patents

Method for improving computed tomography (CT) measured fluid saturation precision Download PDF

Info

Publication number
CN102183532B
CN102183532B CN 201110051246 CN201110051246A CN102183532B CN 102183532 B CN102183532 B CN 102183532B CN 201110051246 CN201110051246 CN 201110051246 CN 201110051246 A CN201110051246 A CN 201110051246A CN 102183532 B CN102183532 B CN 102183532B
Authority
CN
China
Prior art keywords
core
value
water
oil
saturation
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
CN 201110051246
Other languages
Chinese (zh)
Other versions
CN102183532A (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.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas 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 China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN 201110051246 priority Critical patent/CN102183532B/en
Publication of CN102183532A publication Critical patent/CN102183532A/en
Application granted granted Critical
Publication of CN102183532B publication Critical patent/CN102183532B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention relates to a method for improving computed tomography (CT) measured fluid saturation precision. The method comprises the following steps of: a, performing CT scanning on a dry core sample inside and outside a core holder respectively to acquire CT values; b, performing CT scanning on a core of complete saturated brine inside and outside the core holder respectively to acquire CT values; c, making bound water by using the core in the holder, measuring the CT value of the core in a bound water state, and calculating the saturation of the bound water according to the water or liquid outlet volume after the bound water is made completely; d, performing CT scanning on brine for saturation, oil for experiments and peripheral air respectively to acquire CT values, wherein the CT value of the air is a constant value in any medium; e, calculating the CT value of the brine and the CT value of the oil in the core holder by using formulas III and V; and f, calculating the oily water saturation in the displacement process by using the formula V according to the corrected brine and oil CT values acquired in the step e, wherein the relative error of the saturation is within 1 percent.

Description

A kind of method that improves CT measurement fluid saturation precision
Technical field
The present invention relates to a kind of method of testing rock core fluid saturation, be specifically related to the method that a kind of CT accurately measures fluid saturation.
Background technology
The CT technology is as measuring technology conventional in the core analysis, is widely used in nonuniformity mensurations, core sample selection, crack quantitative test of core description, rock core, the online aspects such as measurement, flowing experiment research of saturation degree.And utilize CT that the saturation degree of fluid in the rock is measured, become that the CT technology is the most general in petroleum industry, one of the most effective experimental technique.
The ultimate principle of CT scan rock is: the X-ray beam that CT machine internal X-ray pipe produces from a plurality of directions along object certain selected tomography aspect shine, the X ray amount that sees through by mensuration, organize the absorption coefficient of constituent parts volume after the digitizing through calculating this aspect, these absorption coefficients can consist of different character matrixs; Carry out digital-to-analog conversion by high-speed computer in the machine, can show or take photos at screen, the image of reconstruction can also provide each pixel X ray attenuation coefficient, usually uses the CT value representation.
For a fault surface of rock, the CT value that image provides is the mean value of each volume element, that is:
CT dry=(1-Φ)CT grain+ΦCT air (1)
CT DryThe CT value of-dried rock fault surface
CT GrainThe CT value of-rock skeleton
CT AirThe CT value of-air
The factor of porosity of Φ-rock, %
Want the saturation degree of fluid in the computing rock, must be first with rock 100% saturated a kind of fluid.Such as elder generation's 100% saturation water, then obtain:
CT waterwet=(1-Φ)CT grain+ΦCT water (2)
CT WaterwetThe CT value of fault surface behind-rock 100% saturation water
CT GrainThe CT value of-rock skeleton
CT WaterThe CT value of-water
The factor of porosity of Φ-rock, %
And the factor of porosity of rock also can calculate by (1) (2) two formulas:
Φ = CT waterwet - CT dry CT water - CT air - - - ( 3 )
The CT value of a certain moment rock fault surface is in the displacement process:
CT x=(1-Φ)CT grain+Φ(S wCT water+S oCT oil) (4)
CT xThe CT value of-a certain moment rock fault surface
CT GrainThe CT value of-rock skeleton
CT WaterThe CT value of-water
CT OilThe CT value of-oil
The factor of porosity of Φ-rock, %
S w-water saturation, %
S o-oil saturation, %
According to S w+ S o=1 reaches (1) (2) (4) formula can calculate moisture, oil saturation:
S o = CT waterwet - CT x CT waterwet - CT dry · CT water - CT air CT water - CT oil
S w = 1 - CT waterwet - CT x CT waterwet - CT dry · CT water - CT air CT water - CT oil - - - ( 5 )
As can be seen from the above equation, calculate the CT value that the two-phase fluid saturation degree need to be known air, water and oil.But under normal conditions, rock core all is to test in core holding unit, and because ray hardened impact, it is different that water is directly exposed to airborne CT value with oil in the core holding unit neutralization; Oil and the CT value of water in core holding unit are difficult to again obtain, and all application are exposed to the CT value of airborne water and oil and calculate, and can bring certain error to the result of calculation of saturation degree.
Adopt another kind of method in some research work, be about to rock 100% saturation water post-drying, 100% saturated oil again obtains:
CT oilwet=(1-Φ)CT grain+ΦCT oil (6)
CT OilwetThe CT value of fault surface behind-rock 100% saturated oil
CT GrainThe CT value of-rock skeleton
CT OilThe CT value of-oil
The factor of porosity of Φ-rock, %
Again moisture, oil saturation are calculated:
S o = CT waterwet - CT x CT waterwet - CT oilwet
S w = CT x - CT oilwet CT waterwet - CT oilwet - - - ( 7 )
Avoid oil, water CT value to change the error of bringing with the method, also improved the degree of accuracy of measuring.But common rock core displacement test all needs to make irreducible water, and the method carries out needing twice also oven dry of saturated liquid before the displacement test so that corresponding experimentation is too loaded down with trivial details, and the core sample more intense to susceptibility also can cause certain injury.
Summary of the invention
The objective of the invention is to propose a kind of easy experiment test method that can improve CT measurement fluid saturation precision.
A kind of method that CT measures the fluid saturation precision that improves of the present invention, it comprises the steps:
A. after the core sample oven dry, carry out CT scan, obtain the CT value CT of this dried rock core fault surface Outside the dryAgain this rock core is packed in the core holding unit, and under the same condition of displacement process and under the same CT condition rock core is scanned, obtain the CT value CT of dried rock core fault surface Dry In
B. with behind this rock core 100% saturated brine, in the inside and outside CT scan of carrying out of core holding unit, experiment condition and step a are in full accord respectively, the outer CT value CT of core holding unit behind acquisition rock core 100% saturated brine Outside the waterwetWith the CT value CT in the core holding unit In the waterwet
C. rock core utilizes oily expelling water to make irreducible water in clamper, measures the CT value of core sample under the irreducible water state; And in endpiece metering water outlet, oil pump capacity, made behind the irreducible water according to water outlet, fluid volume irreducible water saturation S wCalculate;
D. carry out CT scan with oil and peripheral air and obtain CT value CT with salt solution, experiment saturated respectively Outside the water, CT Outside the oilAnd CT AirWherein, the CT value of air all is definite value in any medium;
E. utilize formula III that the CT value of salt solution in the core holding unit is calculated CT Water In
Figure BDA0000048721770000051
F. pass through the CT of dried core sample in the core holding unit In the dry, 100% saturated brine the CT of core sample In the waterwet, the CT of core sample under the irreducible water state In the x, the irreducible water saturation S that obtains of step c w, and the CT of salt solution in core holding unit that obtain of step e In the water, utilize formula V that the CT value of oil in the core holding unit is calculated CT In the oil
S w = 1 - CT waterwet - CT x CT waterwet - CT dry · CT water - CT air CT water - CT oil - - - ( V )
G. the corrected salt solution and the oily CT value that obtain according to step e, step f utilize formula V that the profit saturation degree in the displacement process is calculated.
Beneficial effect of the present invention is, final irreducible water saturation was easy to by the volumetric method accurate measurement when irreducible water was made in the method utilization, by this irreducible water saturation, and core sample scans the rock porosity value that obtains outside clamper, the CT value of water and oil is proofreaied and correct rear recycling formula convection cell saturation degree to be calculated, can greatly improve the degree of accuracy of measurement, so that the relative error of saturation degree is in 1%.On the other hand, the method only needs a saturated liquid, has reduced compared with the prior art operation steps.
Description of drawings
Fig. 1 is that oil-water is to permeability curve figure.
Embodiment
Embodiment 1 utilizes the CT scan technology accurately to measure profit two phase saturations and to make profit relative
Permeability curve
(1) sample: the I4D of Daqing oil field Portugal well Sandstone Cores sample 32#.
(2) instrument: GE LightSpeed 8 helical layer Medical CT scanners; The carbon fiber core holding unit; Quzix SP-5000 high-pressure metering pump.
(3) CT scan condition: scanning voltage: 100kV, sweep current: 100mA, scan mode: axial scan, scanning bed thickness: 1.25mm, each scan slice width of cloth number: 49 width of cloth.
(4) experimentation:
1. after core sample oven dry is weighed and measured factor of porosity and air permeability, be placed on core holding unit and scan outward, and the CT value of writing scan direction, position and each scan slice;
2. during core sample is packed core holding unit into, add confined pressure to 5MPa with distilled water, under the direction of scanning with step 1 rock core is scanned, and the CT value of writing scan position and each scan slice;
3. the NaI aqueous solution 24 hours of rock core vacuum saturated 5% is taken out and is claimed weight in wet base and buoyant weight, calculates volume of voids and degree of saturation.Degree of saturation is being carried out next step experiment again more than 99%;
Core sample respectively outside core holding unit and in the core holding unit (adding the 5MPa confined pressure) scan, direction of scanning, position respectively with step 1,2 in full accord, record the CT value of each scan slice;
5. rock core utilizes oily expelling water to make irreducible water in clamper, and in endpiece metering water outlet, oil pump capacity, according to water outlet, fluid volume irreducible water saturation is calculated after having made irreducible water;
6. the core sample (in clamper) of having made irreducible water is scanned, direction of scanning, position and step 2 are in full accord, record the CT value of each scan slice;
7. 5%NaI solution, the experiment of saturated usefulness are carried out CT scan with oil and peripheral air and obtain the CT value;
8. oil and 5%NaI solution are injected rock sample in the ratio of setting, after stablizing wait flowing, core sample (in clamper) is scanned, direction of scanning, position and step 2 are in full accord, record the CT value of each scan slice; Record simultaneously two ends pressure reduction;
9. change profit and inject ratio, repeating step 8 finishes experiment until last profit injects ratio;
10. utilize formula (III) and (V) each profit is injected profit two-phase saturation computation than lower each section, and average and be this profit and inject saturation degree than profit two-phase in the lower rock core; Utilize Darcy formula that the relative permeability of profit two-phase is calculated, and draw permeability-ratio curve as shown in Figure 1.
Experimental data and result are as shown in table 1.
Table 1
Reference examples 1 utilizes volumetric method, conventional CT method and the inventive method to measure profit two phase saturations
(1) sample: two in the I4D of Daqing oil field Portugal well Sandstone Cores sample.2# sample length 6.31cm wherein, diameter 2.541cm, factor of porosity 19.8%, air permeability 154mD; 34# sample length 5.717cm, diameter 2.538cm, factor of porosity 22.5%, air permeability 642mD.
(2) instrument: GE LightSpeed 8 helical layer Medical CT scanners; The carbon fiber core holding unit; Quzix SP-5000 high-pressure metering pump.
(3) CT scan condition: scanning voltage: 100kV, sweep current: 100mA, scan mode: axial scan, scanning bed thickness: 1.25mm.
(4) experimentation:
1. after the core sample oven dry, be placed on core holding unit and scan outward, and the CT value of writing scan direction, position and each scan slice;
2. during core sample is packed core holding unit into, add confined pressure to 5MPa with distilled water, under the direction of scanning with step 1 rock core is scanned, and the CT value of writing scan position and each scan slice;
3. the NaI aqueous solution 24 hours of rock core vacuum saturated 5% is taken out and is claimed weight in wet base and buoyant weight, calculates volume of voids and degree of saturation.Degree of saturation is being carried out next step experiment again more than 99%;
Core sample respectively outside core holding unit and in the core holding unit (adding the 5MPa confined pressure) scan, direction of scanning, position respectively with step 1,2 in full accord, record the CT value of each scan slice;
5. rock core utilizes oily expelling water to make irreducible water in clamper, and in endpiece metering water outlet, oil pump capacity, according to water outlet, fluid volume irreducible water saturation is calculated after having made irreducible water;
6. the core sample (in clamper) of having made irreducible water is scanned, direction of scanning, position and step 2 are in full accord, record the CT value of each scan slice;
7. 5%NaI solution, the experiment of saturated usefulness are carried out CT scan with oil and peripheral air and obtain the CT value;
8. with the flow velocity injection rock sample of 5%NaI solution with 0.5ml/min, 0.5min rear termination of pumping, the volume of oil and water scans core sample (in clamper) simultaneously in the recorded stream fluid, direction of scanning, position and step 2 are in full accord, record the CT value of each scan slice.
9. turn on pump again, every 0.5min stops pump No. one time, and repeating step 8 is until no longer fuel-displaced end experiment.
10. utilize formula (III) and (V) to the profit two-phase saturation computation of each section after the each time termination of pumping, and average and be the saturation degree of profit two-phase in this moment rock core; Calculate profit two phase saturations according to volumetric method simultaneously, computing formula is as follows:
S w = vt + S wi V p - Δ V w V p × 100 %
S w-water saturation, %
The v-injection rate, ml/min
The t-injection length, min
S Wi-irreducible water saturation, %
V p-volume of voids, ml
Δ V w-accumulative total water outlet volume, ml
Utilize volumetric method, conventional CT method (referring to the background technology part) and the inventive method to measure profit two phase saturation results as shown in table 2, take the volumetric method measurement result as standard, conventional CT method measuring error is 0.5~3.6%, and the measuring error of the inventive method only is 0~1.5%, as seen, method measuring accuracy of the present invention is higher.
Table 2
Figure BDA0000048721770000101

Claims (1)

1. one kind is improved the method that CT measures the fluid saturation precision, it is characterized in that, it comprises the steps:
A. after the core sample oven dry, carry out CT scan, obtain the CT value CT of dried rock core fault surface Outside the dryAgain this rock core is packed in the core holding unit, and under the same condition of displacement process and under the same CT condition rock core is scanned, obtain the CT value CT of dried rock core fault surface Dry In
B. with behind this rock core 100% saturated brine, in the inside and outside CT scan of carrying out of core holding unit, experiment condition and step a are in full accord respectively, the outer CT value CT of core holding unit behind acquisition rock core 100% saturated brine Outside the waterwetWith the CT value CT in the core holding unit In the waterwet
C. rock core utilizes oily expelling water to make irreducible water in clamper, measures the CT value CT of core sample under the irreducible water state In the xAnd in endpiece metering water outlet, oil pump capacity, made behind the irreducible water according to water outlet, fluid volume irreducible water saturation S wCalculate;
D. carry out CT scan with oil and peripheral air and obtain CT value CT with salt solution, experiment saturated respectively Outside the water, CT Outside the oilAnd CT AirWherein, the CT value of air all is definite value in any medium;
E. utilize formula III that the CT value of salt solution in the core holding unit is calculated CT WaterIn;
Figure FDA00002076152600011
F. pass through the CT of dried core sample in the core holding unit In the dry, 100% saturated brine the CT of core sample In the waterwet, the CT of core sample under the irreducible water state In the x, the irreducible water saturation S that obtains of step c w, and the CT of salt solution in core holding unit that obtain of step e In the water, utilize formula V that the CT value of oil in the core holding unit is calculated CT In the oil
G. the corrected salt solution and the oily CT value CT that obtain according to step e, step f Water InAnd CT In the oil, utilize formula V that the irreducible water saturation in the displacement process is calculated.
CN 201110051246 2011-03-03 2011-03-03 Method for improving computed tomography (CT) measured fluid saturation precision Active CN102183532B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110051246 CN102183532B (en) 2011-03-03 2011-03-03 Method for improving computed tomography (CT) measured fluid saturation precision

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110051246 CN102183532B (en) 2011-03-03 2011-03-03 Method for improving computed tomography (CT) measured fluid saturation precision

Publications (2)

Publication Number Publication Date
CN102183532A CN102183532A (en) 2011-09-14
CN102183532B true CN102183532B (en) 2013-02-13

Family

ID=44569753

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110051246 Active CN102183532B (en) 2011-03-03 2011-03-03 Method for improving computed tomography (CT) measured fluid saturation precision

Country Status (1)

Country Link
CN (1) CN102183532B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105241901A (en) * 2015-10-16 2016-01-13 中国石油天然气股份有限公司 Method for determining three phase fluid saturation

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102809529B (en) * 2012-08-03 2014-11-05 中国石油天然气股份有限公司 Three-phase relative permeability testing method based on CT (computed tomography) scanning
CN102809528B (en) * 2012-08-03 2015-02-25 中国石油天然气股份有限公司 Three-phase relative permeability testing system based on CT (computed tomography) scanning
CN103592319B (en) * 2013-10-23 2015-12-02 中国石油天然气股份有限公司 Be applicable to CT scan rock core displacement test and the heating constant-temperature system of energy simulating oil deposit condition
CN103592320B (en) * 2013-10-23 2015-11-18 中国石油天然气股份有限公司 Be applicable to CT scan rock core displacement test and the heated constant temperature method of energy simulating oil deposit condition
CN104316548A (en) * 2014-11-19 2015-01-28 中国石油大学(华东) Core holding unit based on Zeiss MCT-400 CT scanning and application thereof
CN105717274A (en) * 2016-01-19 2016-06-29 尹喜玲 Method for improving fluid saturation precision in CT (computed tomography) measurement
CN107014839B (en) * 2017-04-14 2019-11-08 中国石油天然气股份有限公司 Bubble point pressure test method and device
CN107255610B (en) * 2017-05-02 2020-05-08 中国石油天然气股份有限公司 Method and device for detecting fluid saturation of rock core
CN107729590A (en) * 2017-05-19 2018-02-23 青岛科技大学 A kind of new three-phase saturation computational methods based on dual intensity CT numerical value
CN107807078B (en) * 2017-09-13 2019-11-08 中国石油天然气股份有限公司 Rock core mobile water saturation online test method under formation condition
CN107918003B (en) * 2017-11-10 2020-11-06 中国石油天然气股份有限公司 Experimental system and method for measuring saturation front edge of core displacement process in real time
CN109916934A (en) * 2019-04-03 2019-06-21 青岛科技大学 A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data
CN111239006A (en) * 2020-02-03 2020-06-05 中国石油大学(北京) Method and system for detecting imbibition distance

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4868751A (en) * 1987-09-11 1989-09-19 Mobil Oil Corporation Method for determining relative permeability of a subterranean reservoir
US5297420A (en) * 1993-05-19 1994-03-29 Mobil Oil Corporation Apparatus and method for measuring relative permeability and capillary pressure of porous rock
US6178807B1 (en) * 1998-03-25 2001-01-30 Phillips Petroleum Company Method for laboratory measurement of capillary pressure in reservoir rock
CN101592027A (en) * 2009-06-18 2009-12-02 中国海洋石油总公司 A kind of combination flooding method according to correction of capillary number curve

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4868751A (en) * 1987-09-11 1989-09-19 Mobil Oil Corporation Method for determining relative permeability of a subterranean reservoir
US5297420A (en) * 1993-05-19 1994-03-29 Mobil Oil Corporation Apparatus and method for measuring relative permeability and capillary pressure of porous rock
US6178807B1 (en) * 1998-03-25 2001-01-30 Phillips Petroleum Company Method for laboratory measurement of capillary pressure in reservoir rock
CN101592027A (en) * 2009-06-18 2009-12-02 中国海洋石油总公司 A kind of combination flooding method according to correction of capillary number curve

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高建等.应用CT成像技术研究岩心水驱含油饱和度分布特征.《新疆石油地质》.2009,第30卷(第02期),269-271. *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105241901A (en) * 2015-10-16 2016-01-13 中国石油天然气股份有限公司 Method for determining three phase fluid saturation

Also Published As

Publication number Publication date
CN102183532A (en) 2011-09-14

Similar Documents

Publication Publication Date Title
CN102183532B (en) Method for improving computed tomography (CT) measured fluid saturation precision
Akin et al. Computed tomography in petroleum engineering research
CN103954544B (en) A kind of polymer control water increases experimental provision and the experimental technique of gas effect assessment
US5086643A (en) System and method for determining multi-phase relative permeability of a subterranean reservoir
US20080134760A1 (en) Method of characterizing the distribution of the absolute permeability of a heterogeneous sample
CN102809528B (en) Three-phase relative permeability testing system based on CT (computed tomography) scanning
CA2840931C (en) Method for estimating effective atomic number and bulk density of rock samples using dual energy x-ray computed tomographic imaging
CN102809529B (en) Three-phase relative permeability testing method based on CT (computed tomography) scanning
CN210264648U (en) Multifunctional rock core displacement device
Zhao et al. Visualization of CO 2 and oil immiscible and miscible flow processes in porous media using NMR micro-imaging
CN106153662A (en) The measuring method of rock core stress sensitivity
CN107709964A (en) For evaluating the nuclear magnetic resonance gas isothermal line technology of reservoir rock wetability
CN109100812B (en) Rock pore fractal dimension evaluation method and device based on nuclear magnetic resonance
CN108195732A (en) A kind of compact rock core imbibition experimental provision and imbibition weight testing method
CN202814845U (en) Computed tomography (CT) scan-based three-phase relative permeability test system
Rocchi et al. Experimental accuracy of the initial specific volume
RU2467316C1 (en) Method to detect spatial distribution and concentration of component in pore space of porous material
US11112373B1 (en) Systems and methods for slice selective nuclear magnetic resonance testing of fractured core plugs to determine in-situ pore volume
Jang et al. Novel apparatus to measure the low-permeability and porosity in tight gas reservoir
CN103900755A (en) Device and method for measuring minimum miscibility pressure of oil and gas through CT
CN201965077U (en) Rock core computed tomography (CT) scanning fixing device
CN112577979B (en) Quantitative analysis device and method for rock internal fluid saturation spatial distribution
Andersen Comparison of intercept methods for correction of steady-state relative permeability experiments for capillary end effects
CN104914017B (en) One kind utilizes CO in CT technology for detection porous medias2The method of disperse
Hu et al. Unsteady-State coreflooding monitored by positron emission tomography and X-ray computed tomography

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant