CN109916934A - A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data - Google Patents
A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data Download PDFInfo
- Publication number
- CN109916934A CN109916934A CN201910265525.0A CN201910265525A CN109916934A CN 109916934 A CN109916934 A CN 109916934A CN 201910265525 A CN201910265525 A CN 201910265525A CN 109916934 A CN109916934 A CN 109916934A
- Authority
- CN
- China
- Prior art keywords
- porous media
- micro
- phase
- water
- oil
- 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.)
- Pending
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 title claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 42
- 238000006073 displacement reaction Methods 0.000 claims abstract description 12
- 239000012530 fluid Substances 0.000 claims abstract description 5
- 239000011435 rock Substances 0.000 claims description 40
- 238000002474 experimental method Methods 0.000 abstract description 11
- 238000012545 processing Methods 0.000 abstract description 5
- 238000011160 research Methods 0.000 abstract description 3
- 239000012071 phase Substances 0.000 description 39
- 238000000605 extraction Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 239000008346 aqueous phase Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 239000008367 deionised water Substances 0.000 description 1
- 229910021641 deionized water Inorganic materials 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Abstract
The present invention discloses a kind of method for calculating three-phase saturation in porous media based on micro- CT test image gradation data, in laboratory research, is vacuumized to porous media sample, full water, after full oil processing, utilizes CO2Fluid carries out three-phase displacement test.The porous media sample after each experimental procedure is scanned using micro- CT equipment, obtain micro- CT image of the porous media sample after testing each step, pass through the data processing to micro- CT test image gray value and combine related experiment test data, oil, gas and water three-phase saturation numerical value in porous media can be calculated.
Description
Technical field
The invention patent relates to using micro- CT test image gradation data calculate porous media in three saturation degrees method,
For CO2The displacement test process in porous media coexists in oil/water, utilizes micro- CT test image gray value and related experiment number
According to oil/gas in calculating porous media/water three-phase saturation distribution new method, affiliated micro- CT test image data processing field.
Background technique
Utilizing CO2Strengthen oil recovery factor laboratory research in, can accurately obtain CO2 coexisted in oil/water it is porous
Oil, the three-phase saturation parameter of gas and water are extremely important in medium displacement process.In laboratory research, to porous media sample
It vacuumized, full water, after full oil processing, utilize CO2Fluid carries out three-phase displacement test.Using micro- CT equipment to each experiment
Porous media sample after step is scanned, and is obtained micro- CT image of the porous media sample after testing each step, is passed through
Data processing to micro- CT test image gray value simultaneously combines related experiment test data, can calculate oil in porous media,
Gas and water three-phase saturation numerical value.
Summary of the invention
This method technical scheme applied to solve the technical problem is:
Using CO2During oil displacement experiment, it would be desirable to be vacuumized to rock sample, carry out CO after full water, full oil2Stream
Body displacement operation.CT scan is carried out to the rock core after each experimentation, obtains the gray value of rock core after each experimentation.Again
According to the certain physical characteristics of water yield, oil pump capacity and the rock core obtained in experimentation itself, so that it may calculate three-phase
Saturation degree.Due to
Utilize CTrock、CTgas、CTwaterAnd CTrockRespectively represent rock, gas phase, the gray scale of water phase and oily phase in rock core
It is worth, the calculating process of three-phase saturation distribution is as follows in rock core: (due to density of gas phase very little, it is assumed that CTgas=0).
Gray value CT at the end of vacuumizingdryAre as follows:
CTdry=(1- φ) CTrock+φCTgas (1)
It is made of in rock sample the gas component in non-porous rocks component and gap in dry sweep.Wherein φ is the hole of rock core
Porosity is known quantity.
The gray value CT that full water terminateswetAre as follows:
CTwet=(1- φ) CTrock+φCTwater (2)
After rock core satisfies water, the gray value in rock sample is made of the deionized water in non-gap rock and gap.
Simultaneous (1), (2) can obtain the gray value of water phase:
The gray value CT of CT test obtains after full oil oil, water two-phase coexistent rock coreoil+waterIt can be indicated by (4) formula
(wherein Sw、SoilRespectively represent water phase and oil-phase saturation):
CToil+water=CTdry+φSwCTwater+φSoCToil (4)
S in above formulaoil,SwIt can be obtained, be calculated public according to test data of experiment of the experiment front and back to rock core weight change
Formula is (m in formulawetFor rock core weight after full water, moil+waterFor rock core weight after full oil, ρwaterFor aqueous phase densities, ρoilFor oily phase
Density, V are the pore volume of rock core):
And water phase saturation can then be calculated by (6) formula:
Sw=1-So (6)
Simultaneous (1) (2) (3) (4) (5) (6) can be calculated the gray value CT of fuel-displaced phase by (7) formulaoil:
Carry out CO2Fluid displacement experiment, core sample after experiment be water, gas, oily three-phase coexistence porous media,
Core sample is scanned using micro- CT equipment at this time, obtains its gray value CToil+water+gas.Its value can be as shown in (8) formula
Relationship indicates (assuming that CTgas=0):
CToil+water+gas=CTdry+φSwCTwater+φSoCToil (8)
As shown in (9) formula, the oil-phase saturation So in (8) formula can be by total oil content amount (m in rock coretotoil) subtract displacement reality
Test the oily phase quality (m of middle extractionprodoil) obtain (ρ in formulaoilFor oily phase density):
(9) formula substitution (8) formula be can be obtained into the water phase saturation S in core samplew:
In turn, the gas phase saturation S in core sample can be acquired by formula (11)g:
Sg=1.0-Sw-Soil (11)
Based on above method, i.e., it is full the three-phase in porous media to be found out using the rock core gray values that micro- CT test obtains
With degree evidence.
The physical significance and unit list of each symbol are as follows:
V --- the total volume of rock sample void, ml
The porosity of φ --- rock sample, %
CTrock--- the gray value of rock sample skeleton
CTdry--- the gray value of dry core sample
CTwet--- the gray value of full water rock sample
CToil+water--- the gray value of full oil rock sample
CToil+water+gas--- the gray value of rock sample when gas displacement is completed
CTgas--- the gray value of gas phase
CTwater--- the gray value of water phase
CToil--- the gray value of oily phase
Sw--- water phase saturation, %
Sg--- the saturation degree of gas phase, %
So--- the saturation degree of oily phase, %
mwet--- after the completion of full water, the quality of rock sample, g
moil+water--- after the completion of full oil, the quality of rock sample, g
ρwater--- the density of water, g/ml
ρoil--- oily density, g/ml
mtotoil--- total oily phase quality, g in rock core after full oil
mprodoil--- gas drives the oily phase quality of extraction, g
Specific embodiment
In CO2In fluid three-phase displacement process, we to core sample, sequentially take vacuumize, full water, full oil, CO2It drives
For the step of simulate the mining method for truly descending oil reservoir rock, be scanned after each step terminates primary, and survey
Obtain the gray value data of each process rock core.
The central idea of this calculation method is: being directed to the same core sample, obtains gas respectively in the different phase of experiment
The gray value data of phase, water phase and oily phase, in this way can be to avoid due to sample difference bring gray value test error, Jin Erke
Accurately to obtain the saturation data of water during rock core displacement, gas, oily three-phase.
Embodiments above is merely to illustrate the technical solution for the patent that is not intended to limit the present invention.The invention patent essence is not departed from
Any modification or partial replacement of mind, should all cover in the scope of the claims of the invention patent.
Claims (2)
1. a kind of method for calculating three-phase saturation in porous media based on micro- CT test image gradation data described in, condition are
Rock sample is vacuumized, carries out CO after full water, full oil2Fluid displacement operation；CT is carried out to the rock core after each experimentation
Scanning, obtains the gray value of rock core after each experimentation；Further according to obtained in experimentation water yield, oil pump capacity, and
The certain physical characteristics of rock core itself, so that it may calculate three-phase saturation.
A kind of method that three-phase saturation in porous media is calculated based on micro- CT test image gradation data described in 2., due to gas
Phase density very little, it is assumed that CTgas=0.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910265525.0A CN109916934A (en) | 2019-04-03 | 2019-04-03 | A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910265525.0A CN109916934A (en) | 2019-04-03 | 2019-04-03 | A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109916934A true CN109916934A (en) | 2019-06-21 |
Family
ID=66968393
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910265525.0A Pending CN109916934A (en) | 2019-04-03 | 2019-04-03 | A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109916934A (en) |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5164590A (en) * | 1990-01-26 | 1992-11-17 | Mobil Oil Corporation | Method for evaluating core samples from x-ray energy attenuation measurements |
| US5430291A (en) * | 1992-05-01 | 1995-07-04 | Texaco Inc. | X-ray CT measurement of fracture widths and fracture porosity in reservoir core material |
| CN102183532A (en) * | 2011-03-03 | 2011-09-14 | 中国石油天然气股份有限公司 | Method for improving computed tomography (CT) measured fluid saturation precision |
| CN102262041A (en) * | 2011-04-20 | 2011-11-30 | 中国石油天然气股份有限公司 | Multispectral pore structural analysis-based saturation determining method |
| CN102809529A (en) * | 2012-08-03 | 2012-12-05 | 中国石油天然气股份有限公司 | Three-phase relative permeability testing method based on CT (computed tomography) scanning |
| CN103498669A (en) * | 2013-09-04 | 2014-01-08 | 中国石油天然气股份有限公司 | Quantitative determination method of interbedded cross flows of heterogeneous rock core models |
| CN104374682A (en) * | 2014-11-12 | 2015-02-25 | 中国石油天然气股份有限公司 | Rock core CT (Computed Tomography) scanning analysis method and device |
| CN104914017A (en) * | 2015-04-27 | 2015-09-16 | 大连理工大学 | Device and method using CT (computed tomography) technology to detect CO2 dispersion in porous media |
| CN105241901A (en) * | 2015-10-16 | 2016-01-13 | 中国石油天然气股份有限公司 | Method for determining three phase fluid saturation |
| CN106596596A (en) * | 2016-12-13 | 2017-04-26 | 中国海洋石油总公司 | Core X-ray scanning testing experimental method |
| CN107014839A (en) * | 2017-04-14 | 2017-08-04 | 中国石油天然气股份有限公司 | Bubble point pressure method of testing and device |
| CN107271460A (en) * | 2017-04-20 | 2017-10-20 | 东南大学 | A kind of porous material internal moisture saturation degree changes the quantitatively characterizing method of spatial distribution |
| CN107729590A (en) * | 2017-05-19 | 2018-02-23 | 青岛科技大学 | A kind of new three-phase saturation computational methods based on dual intensity CT numerical value |
| CN108613991A (en) * | 2016-12-13 | 2018-10-02 | 中国石油天然气股份有限公司 | The evaluation method and device of foam power of regeneration in heterogeneous core displacement process |
| CN108802073A (en) * | 2017-05-05 | 2018-11-13 | 中国石油化工股份有限公司 | A kind of rock electrical parameters acquisition methods and device based on digital cores |
| CN109376481A (en) * | 2018-08-16 | 2019-02-22 | 清能艾科（深圳）能源技术有限公司 | Calculation method, device and the computer equipment of digital cores phase percolation curve based on more GPU |
-
2019
- 2019-04-03 CN CN201910265525.0A patent/CN109916934A/en active Pending
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5164590A (en) * | 1990-01-26 | 1992-11-17 | Mobil Oil Corporation | Method for evaluating core samples from x-ray energy attenuation measurements |
| US5430291A (en) * | 1992-05-01 | 1995-07-04 | Texaco Inc. | X-ray CT measurement of fracture widths and fracture porosity in reservoir core material |
| CN102183532A (en) * | 2011-03-03 | 2011-09-14 | 中国石油天然气股份有限公司 | Method for improving computed tomography (CT) measured fluid saturation precision |
| CN102262041A (en) * | 2011-04-20 | 2011-11-30 | 中国石油天然气股份有限公司 | Multispectral pore structural analysis-based saturation determining method |
| CN102809529A (en) * | 2012-08-03 | 2012-12-05 | 中国石油天然气股份有限公司 | Three-phase relative permeability testing method based on CT (computed tomography) scanning |
| CN103498669A (en) * | 2013-09-04 | 2014-01-08 | 中国石油天然气股份有限公司 | Quantitative determination method of interbedded cross flows of heterogeneous rock core models |
| CN104374682A (en) * | 2014-11-12 | 2015-02-25 | 中国石油天然气股份有限公司 | Rock core CT (Computed Tomography) scanning analysis method and device |
| CN104914017A (en) * | 2015-04-27 | 2015-09-16 | 大连理工大学 | Device and method using CT (computed tomography) technology to detect CO2 dispersion in porous media |
| CN105241901A (en) * | 2015-10-16 | 2016-01-13 | 中国石油天然气股份有限公司 | Method for determining three phase fluid saturation |
| CN106596596A (en) * | 2016-12-13 | 2017-04-26 | 中国海洋石油总公司 | Core X-ray scanning testing experimental method |
| CN108613991A (en) * | 2016-12-13 | 2018-10-02 | 中国石油天然气股份有限公司 | The evaluation method and device of foam power of regeneration in heterogeneous core displacement process |
| CN107014839A (en) * | 2017-04-14 | 2017-08-04 | 中国石油天然气股份有限公司 | Bubble point pressure method of testing and device |
| CN107271460A (en) * | 2017-04-20 | 2017-10-20 | 东南大学 | A kind of porous material internal moisture saturation degree changes the quantitatively characterizing method of spatial distribution |
| CN108802073A (en) * | 2017-05-05 | 2018-11-13 | 中国石油化工股份有限公司 | A kind of rock electrical parameters acquisition methods and device based on digital cores |
| CN107729590A (en) * | 2017-05-19 | 2018-02-23 | 青岛科技大学 | A kind of new three-phase saturation computational methods based on dual intensity CT numerical value |
| CN109376481A (en) * | 2018-08-16 | 2019-02-22 | 清能艾科（深圳）能源技术有限公司 | Calculation method, device and the computer equipment of digital cores phase percolation curve based on more GPU |
Non-Patent Citations (3)
| Title |
|---|
| 宋晓东: "基于数字岩心的岩石声学特性微观数值模拟研究", 中国优秀硕士学位论文全文数据库 基础科学辑, no. 07, 15 July 2017 (2017-07-15), pages 9 - 18 * |
| 董云鹏等: "油气水三相相对渗透率物理模拟实验", 大庆石油地质与开发, vol. 37, no. 04, 31 August 2018 (2018-08-31), pages 56 - 57 * |
| 黄瑞瑶等: "用工业CT技术测量储层岩芯饱和度分布", 西南石油学院学报, vol. 12, no. 02, 31 May 1999 (1999-05-31), pages 25 - 26 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yassin et al. | Organic shale wettability and its relationship to other petrophysical properties: A Duvernay case study | |
| Iglauer et al. | X-ray tomography measurements of power-law cluster size distributions for the nonwetting phase in sandstones | |
| Liu et al. | Quantification of the microstructures of Bakken shale reservoirs using multi-fractal and lacunarity analysis | |
| Aljamaan et al. | Multiscale imaging of gas storage in shales | |
| Zhang et al. | Characterization of coal porosity and permeability evolution by demineralisation using image processing techniques: A micro-computed tomography study | |
| Suleimenova et al. | Acid demineralization with critical point drying: A method for kerogen isolation that preserves microstructure | |
| Ma et al. | Linking multi-scale 3D microstructure to potential enhanced natural gas recovery and subsurface CO 2 storage for Bowland shale, UK | |
| CN110231272A (en) | Tight sand aperture and nuclear magnetic resonance T2It is worth the determination method and system of transformational relation | |
| Hosseini et al. | The effect of heterogeneity on NMR derived capillary pressure curves, case study of Dariyan tight carbonate reservoir in the central Persian Gulf | |
| Thomson et al. | The influence of microporous cements on the pore network geometry of natural sedimentary rocks | |
| CN110470584B (en) | Method for evaluating comprehensive effect of imbibition and water lock | |
| Farokhpoor et al. | Permeability and relative permeability measurements for CO2‐brine system at reservoir conditions in low permeable sandstones in Svalbard | |
| CN105137033A (en) | Gypsum salt-carbonate reservoir rock supergene karst effect quantitative evaluation method and device | |
| Aljamaan et al. | Multiscale imaging of gas adsorption in shales | |
| Kim et al. | Data-driven estimation of three-phase saturation during gas hydrate depressurization using CT images | |
| Tang et al. | Generalizable framework of unpaired domain transfer and deep learning for the processing of real-time synchrotron-based x-ray microcomputed tomography images of complex structures | |
| Saputra et al. | Comprehensive Workflow for Laboratory to Field-Scale Numerical Simulation to Improve Oil Recovery in the Eagle Ford Shale by Selective Testing and Modelling of Surfactants for Wettability Alteration | |
| Drexler et al. | Wettability alteration by carbonated brine injection and its impact on pore-scale multiphase flow for carbon capture and storage and enhanced oil recovery in a carbonate reservoir | |
| CN109916934A (en) | A method of three-phase saturation in porous media is calculated based on micro- CT test image gradation data | |
| Azenkeng et al. | An image-based equation for estimating the prospective CO2 storage resource of organic-rich shale formations | |
| Wei et al. | Mechanism model for shale gas transport considering diffusion, adsorption/desorption and Darcy flow | |
| Li et al. | In-situ imaging of CO2 trapping and oil recovery in three-phase systems: Dependence on pore geometry and wettability | |
| Dalton et al. | Characterizing the evolution of trapped scCO2 curvature in Bentheimer and Nugget sandstone pore geometry | |
| Ansari et al. | The gas-in-place and CO2 storage capacity of shale reservoirs at subsurface conditions | |
| Bhattad et al. | X-ray micro-CT assisted drainage rock typing for characterization of flow behaviour of laminated sandstone reservoirs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |