Analysis of Chalk Petrophysical Properties by Means of Submicron-Scale Pore Imaging and Modeling
- Liviu Tomutsa (Lawrence Berkeley Laboratory) | Dmitriy Borisovich Silin (Lawrence Berkeley Laboratory) | Velimir Radmilovic (Lawrence Berkeley Laboratory)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- June 2007
- Document Type
- Journal Paper
- 285 - 293
- 2007. Society of Petroleum Engineers
- 4.3.1 Hydrates, 1.2.3 Rock properties, 5.1 Reservoir Characterisation, 4.3.4 Scale, 4.1.2 Separation and Treating, 5.8.1 Tight Gas, 4.1.5 Processing Equipment, 5.3.2 Multiphase Flow, 5.3.1 Flow in Porous Media
- 2 in the last 30 days
- 1,001 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
For many rocks of high economic interest such as chalk, diatomite, shale, tight gas sands, or coal, a submicron-scale resolution is needed to resolve the 3D pore structure, which controls the flow and trapping of fluids in the rocks. Such a resolution cannot be achieved with existing tomographic technologies. A new 3D imaging method based on serial sectioning, which uses the focused-ion-beam (FIB) technology, has been developed. FIB technology allows for the milling of layers as thin as 10 nm by using accelerated gallium (Ga+) ions to sputter atoms from the sample surface. After each milling step, as a new surface is exposed, a 2D image of this surface is generated, and the 2D images are stacked to reconstruct the 3D pore structure. Next, the maximum-inscribed-spheres (MIS) image-processing method computes the petrophysical properties by direct morphological analysis of the pore space. The computed capillary pressure curves agree well with laboratory data. Applied to the FIB data, this method generates the fluid distribution in the chalk pore space at various saturations.
Field-scale oil-recovery processes are the result of countless events happening in individual pores. To model multiphase flow in porous media at pore scale, the resolution of the 3D images must be adequate for the rock of interest. Chalk formations in the oil fields of Texas, the Middle East, the North Sea, and other areas hold significant oil reserves. The extremely small typical pore sizes in chalk impose very high requirements on imaging resolution. In the last decade, X-ray microtomography has been used extensively for direct visualization of the pore system and the fluids within sandstone (Jasti et al. 1993; Coles et al. 1998; Wildenschild et al. 2003; Seright et al. 2003). While this approach is fast and nondestructive, its applicability is limited mostly to micron resolutions, although recent developments are bringing the resolution to submicron range (Stampanoni et al. 2002). For chalk pore systems, which are characterized by submicron- to nanometer-length scales, 3D stochastic methods based on 2D scanning-electron-microscope (SEM) images of thin sections have been used to reconstruct the pore system (Talukdar et al. 2001).
The advent of FIB technology has it made possible to reconstruct submicron 3D pore systems for diatomite and chalk (Tomutsa and Radmilovic 2003) (Fig. 1). FIB technology is used in microelectronics to access individual components with nanoscale accuracy for design verification, failure analysis, and circuit modification (Orloff et al. 2002). FIB has been used in material sciences for sectional sample preparation for SEM and for 3D imaging of alloy components (Kubis et al. 2004). In earth sciences, the FIB also has been used for sample preparation for SEM and to access inner regions for performing microanalysis (Heaney et al. 2001).To access the pore structure at submicron scale, the FIB mills successive layers of the rock material as thin as 10 nm. As successive 2D surfaces are exposed, they are imaged with either the electron or the ion beam. After processing, the images are stacked to reconstruct the 3D pore structure. The geometry of the pore space of the obtained structure can be analyzed further to estimate petrophysical rock properties through computer simulations. To analyze the 3D chalk images obtained by the FIB method, we applied the MIS technique (Hazlett 1995; Silin et al. 2003, 2004; Silin and Patzek 2006). The MIS method analyzes the 3D pore-space image directly, without construction of pore networks. It bypasses the nontrivial task of extracting a simple but representative network of pore throats linking pore bodies from the 3D data (Lindquist 2002). Moreover, the pore-network extraction methods, which are based on relatively simple grain and pore shapes in sandstones (Øren and Bakke 2002), may not always be feasible for the complex pore structures of carbonates. Although a pore-network-based flow-modeling approach enjoyed a significant interest from the researchers and resulted in theoretically and practically sound conclusions (Øren et al. 1998; Xu et al. 1999; Patzek 2001; Blunt 2001), we believe that direct pore-space analysis deserves more attention. In addition, direct analysis of the pore space provides an opportunity to study alteration of the rock flow properties (e.g., those resulting from mechanical transformations or mineralization) (Jin et al. 2003).
|File Size||1 MB||Number of Pages||9|
Barenblatt G.I., Patzek, T.W., and Silin D.B. 2003.The MathematicalModel of Nonequilibrium Effects in Water-Oil Displacement.SPEJ 8 (4): 409-416. SPE-87329-PA. DOI: 10.2118/87329-PA.
Coles, M.E., Hazlett, R.D., Spanne, P., Soll, W.E., Muegge,E. L., and Jones, K.W. 1998. Pore Level Imaging of Fluid Transport UsingSynchrotron X-Ray Microtomography. J. Pet. Sci. & Eng.19(1): 55-63. DOI: 10.1016/S0920-4105(97)00035-1.
Heaney, P.J., Vicenzi, E.P., Gianuzzi, L.A., and Livi, J.T. 2001. Focused Ion BeamMilling: A Method of Site Specific Sample Extraction for Microanalysis of Earthand Planetary Materials. American Mineralogist 86:1094.
Image-J. 2007.http://rsb.info.nih.gov/ij/. Accessed on 11 January 2007.
Jasti, J.K., Jesion, G., and Feldkamp, L. 1993.MicroscopicImaging of Porous Media With X-Ray Computed Tomography.SPEFE 8 (3): 189-193. SPE-20495-PA. DOI: 10.2118/20495-PA.
Jin, G., Patzek, T.W., and Silin, D.B. 2003.Physics-BasedReconstruction of Sedimentary Rocks. PaperSPE 83587 presented at the SPE Western Regional/AAPG Pacific Section JointMeeting, Long Beach, California,19-24 May. DOI: 10.2118/83587-MS.
Jin, G., Patzek, T.W., and Silin, D.B. 2004.Direct Predictionof the Absolute Permeability of Unconsolidated and Consolidated ReservoirRock. Paper SPE 90084 presented at the SPE Annual TechnicalConference and Exhibition,Houston, 26-29 September. DOI:10.2118/90084-MS.
Lindquist, W.B. 2002. Quantitative Analysis of ThreeDimensional X-Ray Tomographic Images. In Development of X-ray TomographyIII, ed. U.Bouse; Proc., SPIE 4503, 103 SPIE; available atftp://ams.sunusb.edu/pub/papers/2001/susb01_06.pdf.
Magnani, F.S.,Philippi, P.C., Liang, Z.R., and Fernandes, C.P.2000. Modelling Two-Phase Equilibrium in Three-Dimensional PorousMicrostructures. Intl. J. of Multiphase Flow 26 (1): 99-123. DOI:10.1016/S0301-9322(99)00008-7.
Meijering, E.H., Niessen, W.W.J., and Viergever, M.A. 2001.Quantitative Evaluation of Convolution-Based Methods for Medical ImageInterpolation. Medical Image Analysis 5 (2):111-126. DOI: 10.1016/S1361-8415(00)00040-2.
Øren, P.-E., Bakke, S., and Arntzen, O.J. 1998.ExtendingPredictive Capabilities to Network Models.SPEJ 3 (4): 324-336.SPE-52052-PA. DOI: 10.2118/52052-PA.
Orloff, J., Swanson, L., and Utlaut, M. 2002. HighResolution Focused Ion Beams: FIB and Its Applications. Dordrecht, The Netherlands: Kluwer AcademicPublishers.
Patzek, T.W. 2001. Verification of a Complete PoreNetwork Simulator of Drainage and Imbibition. SPEJ 6 (2):144-156. SPE-71310-PA. DOI: 10.2118/71310-PA.
Seright, R.S.,Liang, J., Lindquist, W.B., and Dunsmuir, J.H.2003. Use of X-Ray Computed Microtomography to Understand Why Gels ReduceRelative Permeability to Water More Than That to Oil. J. Pet. Sci. &Eng. 39: 217-230. DOI: 10.1016/S0920-4105(03)00064-0.
Silin, D.B. and Patzek, T.W. 2003. Object-Oriented ClusterSearch for an Arbitrary Pore Network. Technical ReportLBNL-51599, Lawrence BerkeleyNational Laboratory, Berkeley, California (January 2003).
Silin, D.B. and Patzek, T.W. 2004. OnBarenblatt's Model of Spontaneous Countercurrent Imbibition. Transport inPorous Media 54 (3): 297-322. DOI: 10.1023/B:TIPM.0000003678.85526.b1.
Silin, D.B., Jin, G., and Patzek, T.W.2003. Robust Determination of thePore Space Morphology in Sedimentary Rocks. Paper SPE 84296 presented atthe SPE Annual Technical Conference and Exhibition, Denver, 5-8 October. DOI:10.2118/84296-MS.
Spinler, E.A. and Baldwin, B.A. 1997. Capillary Pressure ScanningCurves by Direct Measurement of Saturations. Paper SCA-9705 presented at theInternationalSymposium of the Society of Core Analysts, Calgary, 8-10September.
Stampanoni, M., Borchart, G., Abella, R., and Ruegsegger, P.2002. Bragg Magnifier: A Detector for Submicrometer X-Ray Computer Tomography.J. of Applied Physics 2 (12): 7630-7635. DOI: 10.1063/1.1520722.
Talukdar, M.S., Ioannidis, M., Howard, J., and Torsaeter, O.2001. Network Modeling as a Tool for Petrophysical Measurements in Chalk.Proc., 6th Nordic Symposium on Petrophysics, Trondheim, Norway, 15-16 May.
Tomutsa, L. and Radmilovic, V. 2003. Focused Ion BeamAssisted Three-Dimensional Rock Imaging at Submicron Scale. Paper SCA2003-47presented at the International Symposium of the Society of CoreAnalysts, Pau, France, 21-24September.
Willhite, G.P. 1986. Waterflooding. Textbook Series,SPE, Richardson, Texas,3.
Xu, B., Kamath, J., Yortsos, Y.C., andLee, S.H. 1999. Use ofPore-Network Models To Simulate Laboratory Corefloods in a HeterogeneousCarbonate Sample. SPEJ 4 (3): 178-186. SPE-57664-PA. DOI:10.2118/57664-PA.