CN109540765A - A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core - Google Patents

A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core Download PDF

Info

Publication number
CN109540765A
CN109540765A CN201910017275.9A CN201910017275A CN109540765A CN 109540765 A CN109540765 A CN 109540765A CN 201910017275 A CN201910017275 A CN 201910017275A CN 109540765 A CN109540765 A CN 109540765A
Authority
CN
China
Prior art keywords
image
rock core
micro
hole
model
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
Application number
CN201910017275.9A
Other languages
Chinese (zh)
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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN201910017275.9A priority Critical patent/CN109540765A/en
Publication of CN109540765A publication Critical patent/CN109540765A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N2015/0846Investigating permeability, pore-volume, or surface area of porous materials by use of radiation, e.g. transmitted or reflected light

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a kind of, and parameter prediction method is seeped in the hole based on the micro- CT image of rock core.This method is based on the micro- CT image of rock core, generate the surface mesh and volume mesh of rock core pore structure respectively by Mimics and ICEM software, to establish micron order hole unstructured grid model, and be conducted into Fluent software and carry out micro flow numerical simulation, it obtains and coincide preferable porosity, permeability and permeability saturation curve data with experimental result.This method is at low cost compared with traditional core laboratory experiment, and the test period is short, and can use grid model generated and carry out comparative simulation experiment.

Description

A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core
Technical field
The present invention relates to petroleum works fields more particularly to a kind of hole based on the micro- CT image of rock core to seep parameter prediction side Method.
Background technique
Porosity, permeability and the relative permeability of rock are to determine the key physical ginseng of oil reservoir development scheme and production capacity Number, and the pore structure characteristic of these parameters and rock is closely related.The present invention is based on the micro- CT images of rock, in conjunction with Mimics weight Structure software and ICEM grid dividing software propose and a kind of the micro- CT image of core three-dimensional are converted to unstructured grid model Modeling method, and the single-phase and water-oil phase microscopic seepage process in model is simulated using Fluent software, obtain hole The parameters such as degree, permeability and permeability saturation curve seep parameter prediction for rock pore and provide a kind of new method.
Summary of the invention
It is an object of the invention to combine the micro- CT image of rock core, Mimics and ICEM software building hole unstructured grid Model, and microcosmic single-phase and oil-water two-phase flow numerical simulation is carried out using Fluent software based on reconstruction model, to obtain rock The parameters such as heart porosity, permeability and permeability saturation curve.
In order to achieve the above object, the present invention is implemented as follows:
(1) based on the micro- CT image of rock core, hole is extracted using image processing techniques such as filtering, segmentations, is passed through Mimics and ICEM Software Create hole unstructured grid model.
(2) by statistics hole unstructured grid model volume, core porosity is acquired.
(3) based on the hole grid model of above-mentioned acquisition, carry out single-phase and grease in hole using Fluent software Two phase fluid flow numerical simulation, difference prediction model permeability and permeability saturation curve.
Compared with prior art, benefit of the invention is that:
1. can be used for the unstructured grid model of numerical value calculating based on the micro- CT picture construction of rock core, the model is preferably The rock labyrinth feature for presenting the micro- CT image of original rock core, meets requirement of the finite element software to mesh quality.
2. few based on the experiment capital investment of reconstruction model Numerical Simulation of Seepage compared with traditional rock core laboratory experiment, experiment Period is short, and can carry out the comparative test under different operating conditions on same model.
Detailed description of the invention
In order to illustrate more clearly of the technical solution of the method for the present invention, with reference to the accompanying drawings and detailed description to this Shen Please embodiment be described further.
Fig. 1 is the application method implementation flow chart.
Fig. 2 is the original micro- CT image of rock core used by this method.
Fig. 3 is hole unstructured grid model.
Fig. 4 is hole unstructured grid model degree of skewness distribution map.
Fig. 5 is the setting of numerical simulation boundary condition.
Fig. 6 is rock sample S2-1 single-phase flow analog result, and (a) is velocity field cloud atlas, (b) is pressure field cloud atlas.
Fig. 7 be rock sample S1 water-oil phase flow field simulation oil saturation cloud atlas, (a) 5 time step, (b) 22 time step, (c) 40 Time step, (d) 60 time step.
Fig. 8 is that the permeability saturation curve that numerical simulation acquires and experimental result compare.
Specific embodiment
In order to make the present invention realize technological means, reach purpose and model efficacy ease of explanation, below with reference to specific reality Applying details and attached drawing, the present invention is further illustrated.It should be noted that described embodiment is only the application part Embodiment is not whole embodiments.Based on the present processes, other researchers, technical staff of this field are not having The every other embodiment obtained under the premise of other innovative labors all should be the application protection scope.
Fig. 1 is the flow chart that parameter prediction embodiment of the method is seeped in a kind of hole based on the micro- CT image of rock core of the application, including Following steps.
S1: being illustrated in figure 2 the micro- CT scan image of the original rock core used in the present embodiment, and rock sample as shown in the figure is sandstone (S2-2), remaining rock sample further include: Berea sandstone (B1), sandstone (S1), sandstone (S2-1), synthesis silicon (SS) and carbonate rock (C2).The micro- CT image of rock core used by the embodiment of the present invention mainly passes through Southwest Petrol University's hydrocarbon reservoirs and Development Engineering The Zeiss Xradia MicroXCT-400 microscopic imaging systems of National Key Laboratory, which scan, to be obtained.
S2: importing Mimics software for the micro- CT image of rock core, generates gore grid in model surface.
S3: the veil lattice model of generation is imported into ICEM software, surface grids ingrowing is made to complete the generation of volume mesh.? In the partition process of volume mesh, using contraction-expansion algorithm, the production rate of volume mesh is improved, is rejected in (filling) model Sharp protrusion (slit), hole unstructured grid model generated are as shown in Figure 3.
S4: in ICEM software, by counting the volume of hole grid model, model porosity is acquired using formula (1-1).
Wherein, VpFor hole grid model volume, VtFor hole and skeleton total volume, L is grid model side length.
S5: boundary condition is applied to grid model in Fluent software.In single-phase flow simulation process, setting model is Inlet and outlet pressure boundary condition, by taking the direction z as an example (Fig. 5), model upper and lower surface is respectively set to pressure inlets and pressure export Boundary, remaining surface are set as impermeable boundary.Assuming that fluid is the laminar flow that property is kept constant, repaired using SIMPLEC pressure Positve term keeps the relaxation factor of default.This simulates the influence for not considering temperature, and temperature keeps 273K constant in simulation process. Then carry out single-phase flow numerical simulation, it is to be calculated reach convergence after, obtain rate of discharge, model can be acquired using formula (1-2) Permeability:
Wherein, QiFor mold exit flow, A is mold exit sectional area, and Δ p is model inlet and outlet pressure gradient.
By changing boundary condition, simulate fluid respectively along the flowing in x, y, z direction, calculate separately to obtain model along x, y, The permeability in the direction z, wherein kzThe permeability of representative model in the z-direction.Table 1 give unstructured grid model permeability with The comparison of experimental data.As known from Table 1, model Permeability Prediction value and experimental data have a different, this species diversity be by It is original big core permeability in experimental data and the permeability of model prediction is small size core permeability.Therefore, the deviation It is in acceptable range.
1 model permeability data of table and experimental data
In addition, it can be seen that pore-fluid pressure from the pressure cloud atlas and speed cloud atlas (Fig. 6) of the single-phase flow field simulation of rock sample S1 The main thoroughfare that the uneven distribution and fluid in the field of force flow through.
S6: based in Fluent software VOF model carry out water drive oil numerical simulation, the initial saturated oils of model, into Water is injected at mouthful to simulate the process of water drive oil.Apply the barometric gradient of 5MPa/m in model upper and lower surface, lap is not permeable Saturating boundary.Grease physical parameter, surface tension and the contact angle used in simulation process is shown in Table 2.By the way that different contacts is arranged The characterization of rock sample different wetting is realized at angle.
2 grease physical parameter of table
5th, 22,40, the model oleaginous saturation degree cloud atlas of 60 time steps it is as shown in Figure 7.It is every in simulation process by extracting The rate of discharge of the oil saturation of a time step and each phase fluid obtains the microcosmic water drive of model in conjunction with formula (1-3) and (1-4) The permeability saturation curve of oily process, as shown in Figure 8.
sii (1-4)
Wherein S0 be the displacement of primitive rock intra-ventricle as a result, S1, S2-1, S2-2 respectively correspond hygroscopicity, mixed wettability and Oil-wet rock sample.Original rock sample is identical with model S2-1 wetability, is mixed wettability.From permeability saturation curve it is found that The permeability saturation curve and primitive rock intra-ventricle displacement result S0 of model S2-1 coincide preferably, demonstrates the phase of the application proposition To the reasonability of permeability curve prediction technique.
Embodiment described above is only section Example of the invention, and basic principle for describing the present invention implements mesh And process, be not intended to limit use scope of the invention.Embodiment of above is done according to the technical essence of the invention Any modification, equivalent variations and modification, belong in the range of technical solution of the present invention.

Claims (7)

1. parameter prediction method is seeped in a kind of hole based on the micro- CT image of rock core, which is characterized in that this method comprises:
(1) unstructured grid model is rebuild based on the micro- CT image of rock core;
(2) core porosity is obtained by the volume of statistics pore model;
(3) using Fluent software as numerical value computing platform, the simulation of single-phase and water-oil phase microscopic seepage is realized respectively, obtains rock core Permeability and permeability saturation curve.
2. parameter prediction method is seeped in a kind of hole based on the micro- CT image of rock core described in accordance with the claim 1, it is characterised in that: institute State in step (1), based on the micro- CT image of rock core, the image such as scanning electron microscope of similar available rock micropore structure and Rock core casting body flake etc. can also be used as the basic data of reconstruction model.
3. parameter prediction method is seeped in a kind of hole based on the micro- CT image of rock core described in accordance with the claim 1, it is characterised in that: institute It states in step (1), by the image processing techniques such as filtering, dividing, extracts rock core skeleton and pore structure.
4. parameter prediction method is seeped in a kind of hole based on the micro- CT image of rock core described in accordance with the claim 1, it is characterised in that: institute It states in step (1), the surface grid model of pore structure is obtained by Mimics software, importing ICEM generation can be used for numerical value meter The hole unstructured grid model of calculation.
5. parameter prediction method is seeped in a kind of hole based on the micro- CT image of rock core described in accordance with the claim 1, it is characterised in that: institute It states in step (2), hole unstructured grid model volume is calculated by statistics, predicts core porosity.
6. parameter prediction method is seeped in a kind of hole based on the micro- CT image of rock core described in accordance with the claim 1, it is characterised in that: institute It states in step (3), based on Fluent software, simulates the single-phase flow in hole unstructured grid model, it is fixed by darcy Restrain solving model permeability.
7. parameter prediction method is seeped in a kind of hole based on the micro- CT image of rock core described in accordance with the claim 1, it is characterised in that: institute It states in step (3), based on the VOF model in Fluent software, simulates the grease two in hole unstructured grid model Xiang Liu passes through analog result solving model permeability saturation curve.
CN201910017275.9A 2019-01-08 2019-01-08 A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core Pending CN109540765A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910017275.9A CN109540765A (en) 2019-01-08 2019-01-08 A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910017275.9A CN109540765A (en) 2019-01-08 2019-01-08 A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core

Publications (1)

Publication Number Publication Date
CN109540765A true CN109540765A (en) 2019-03-29

Family

ID=65834605

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910017275.9A Pending CN109540765A (en) 2019-01-08 2019-01-08 A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core

Country Status (1)

Country Link
CN (1) CN109540765A (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109946328A (en) * 2019-04-10 2019-06-28 中国石油大学(北京) Fibrous filter Micro-Structure Analysis method and apparatus
CN110160932A (en) * 2019-06-03 2019-08-23 西南石油大学 A kind of oil-water relative permeability curve test device and test method
CN110609046A (en) * 2019-10-22 2019-12-24 西南石油大学 Capillary force prediction method based on micro-CT water flooding image
CN112394072A (en) * 2020-11-26 2021-02-23 西安石油大学 Micro-CT-based core broadband dielectric constant characterization method and device
CN113237809A (en) * 2021-04-16 2021-08-10 贵州电网有限责任公司 Composite insulator core rod porosity evaluation method
CN113295592A (en) * 2021-05-20 2021-08-24 中国石油大学(北京) Relative permeability determination system and method for low-permeability rock
CN113295594A (en) * 2021-05-20 2021-08-24 中国石油大学(北京) System and method for measuring relative permeability based on shunt model and CT scanning
CN113916916A (en) * 2021-09-29 2022-01-11 西南石油大学 Simulation method for three-dimensional seepage-particle flow coupling of shale digital core
CN114397235A (en) * 2022-01-17 2022-04-26 中国矿业大学 In-situ temperature-pressure two-phase flow analysis method for reconstructing digital core based on porosity
CN114428040A (en) * 2021-09-13 2022-05-03 中国石油化工股份有限公司 Quantitative characterization and parameter acquisition method for shale oil reservoir storage and seepage space
CN118150441A (en) * 2024-05-11 2024-06-07 中国地质大学(北京) Method for evaluating rock micro-pore structure based on fluid flow characteristics

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109946328A (en) * 2019-04-10 2019-06-28 中国石油大学(北京) Fibrous filter Micro-Structure Analysis method and apparatus
CN110160932A (en) * 2019-06-03 2019-08-23 西南石油大学 A kind of oil-water relative permeability curve test device and test method
CN110160932B (en) * 2019-06-03 2023-12-15 西南石油大学 Oil-water relative permeability curve testing device and testing method
CN110609046A (en) * 2019-10-22 2019-12-24 西南石油大学 Capillary force prediction method based on micro-CT water flooding image
CN110609046B (en) * 2019-10-22 2022-12-13 西南石油大学 Capillary force prediction method based on micro-CT water flooding image
CN112394072A (en) * 2020-11-26 2021-02-23 西安石油大学 Micro-CT-based core broadband dielectric constant characterization method and device
CN112394072B (en) * 2020-11-26 2021-10-22 西安石油大学 Micro-CT-based core broadband dielectric constant characterization method and device
CN113237809A (en) * 2021-04-16 2021-08-10 贵州电网有限责任公司 Composite insulator core rod porosity evaluation method
CN113295594B (en) * 2021-05-20 2022-07-12 中国石油大学(北京) System and method for measuring relative permeability based on shunt model and CT scanning
CN113295592B (en) * 2021-05-20 2022-06-07 中国石油大学(北京) Relative permeability determination system and method for low-permeability rock
CN113295594A (en) * 2021-05-20 2021-08-24 中国石油大学(北京) System and method for measuring relative permeability based on shunt model and CT scanning
CN113295592A (en) * 2021-05-20 2021-08-24 中国石油大学(北京) Relative permeability determination system and method for low-permeability rock
CN114428040A (en) * 2021-09-13 2022-05-03 中国石油化工股份有限公司 Quantitative characterization and parameter acquisition method for shale oil reservoir storage and seepage space
CN114428040B (en) * 2021-09-13 2022-07-22 中国石油化工股份有限公司 Quantitative characterization and parameter acquisition method for shale oil reservoir storage and seepage space
CN113916916A (en) * 2021-09-29 2022-01-11 西南石油大学 Simulation method for three-dimensional seepage-particle flow coupling of shale digital core
CN113916916B (en) * 2021-09-29 2024-03-19 西南石油大学 Simulation method for shale digital rock core three-dimensional seepage-particle flow coupling
CN114397235A (en) * 2022-01-17 2022-04-26 中国矿业大学 In-situ temperature-pressure two-phase flow analysis method for reconstructing digital core based on porosity
CN118150441A (en) * 2024-05-11 2024-06-07 中国地质大学(北京) Method for evaluating rock micro-pore structure based on fluid flow characteristics

Similar Documents

Publication Publication Date Title
CN109540765A (en) A kind of hole infiltration parameter prediction method based on the micro- CT image of rock core
Hálek et al. Groundwater hydraulics
Ramstad et al. Relative permeability calculations from two-phase flow simulations directly on digital images of porous rocks
Corapcioglu et al. Micromodel visualization and quantification of solute transport in porous media
Van Marcke et al. An improved pore network model for the computation of the saturated permeability of porous rock
Lu et al. Controlling capillary fingering using pore size gradients in disordered media
Corapcioglu et al. Glass bead micromodel study of solute transport
Chang et al. Scaling the impacts of pore-scale characteristics on unstable supercritical CO2-water drainage using a complete capillary number
Melnikov et al. Grain-scale modeling of arbitrary fluid saturation in random packings
CN109632604B (en) Method for coarsening relative permeability of polymer flooding from pore size to core size
CN104881522B (en) A kind of remaining oil characterizing method and device based on formulas calculating
Taylor et al. Sub-particle-scale investigation of seepage in sands
Shen et al. Spontaneous imbibition in asymmetric branch-like throat structures in unconventional reservoirs
Gao et al. Capillarity and phase-mobility of a hydrocarbon gas–liquid system
Mourzenko et al. Deposition in fractures
CN110987749B (en) Method for researching equivalent permeability coefficient microscopic scale of multiphase composite material
Li et al. Partial fractional differential model for gas-liquid spontaneous imbibition with special imbibition index: imbibition behavior and recovery analysis
Standnes Spontaneous imbibition of water into cylindrical cores with high aspect ratio: Numerical and experimental results
Zakirov et al. Simulation of two-phase fluid flow in the digital model of a pore space of sandstone at different surface tensions
Kanitz et al. Multiscale investigation of suffusion with coupled CFD-DEM-Influence of different drag force models
CN109387467B (en) Rapid testing method for oil-water comprehensive starting pressure gradient
Yao et al. Upscaling of carbonate rocks from micropore scale to core scale
Sarris et al. A novel methodology for synthesizing the wetting capillary pressure curve of rocks through sorptivity data
Zhang et al. Simulation of non-uniform sediment transport in a German Reservoir with the SSIIM Model and sensitivity analysis
Crestana 2-D and 3-D fingering in unsaturated soils investigated by fractal analysis, invasion percolation modeling and non-destructive image processing

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190329