Analysis of Chalk Petrophysical Properties by Means of Submicron-Scale Pore Imaging and Modeling

Abstract

Summary 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. Introduction 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).

Study on microheterogeneity and seepage law of reservoir after long-term water flooding

“…Moreover, after long-term flooding in reservoir of the Jidong Oilfield, the separation of pore and throat becomes better, and microheterogeneity of reservoir gradually becomes weaker. The rule has not been reported before (Elliott and Dover 1982;Sorbie et al 1989;Lymberopoulos and Payatakes 1992;Xu et al 1999;Fredrich et al 1995;Fredrich 1999;Delerue et al 1999;Lindquist and Venkatarangan 1999;Liang et al 2000;Blunt (2001); Radlinski et al 2002;Siddiqui and Khamees 2004;Coenen et al 2004;Arns et al 2004Arns et al , 2005Knackstedt et al 2006;Al-Kharusi and Blunt 2007;Tomutsa et al 2007;Hou et al 2007;Youssef et al 2008;Dong and Blunt 2009;Mustafa et al 2009;Xu et al 2010;Peng et al 2012;Saurabh et al 2013;Theis et al 2014). …”

“…Since the pore network model was proposed by Fatt (1956aFatt ( , b, c) in 1956, it has been a powerful tool to study microscopic reservoir structure and seepage characteristics. Elliott and Dover (1982), Sorbie et al (1989), Lymberopoulos and Payatakes (1992), Xu et al (1999), Fredrich et al (1995), Fredrich (1999, Delerue et al (1999), Lindquist and Venkatarangan (1999), Liang et al (2000), Blunt (2001), Radlinski et al (2002), Siddiqui and Khamees (2004), Coenen et al (2004), Knackstedt et al (2006), Al-Kharusi and Blunt (2007), Tomutsa et al (2007), Blunt (2009), Peng et al (2012), Saurabh et al (2013), Theis et al (2014), have researched the technologies and algorithms of constructing the digital core and 3D pore network model, and the 3D pore network model has been rapidly developed and widely applied. Among those technologies, the method of constructing 3D pore network model with digital core reconstructed by high-precision X-ray CT scanning is dominant, and it is future direction of researching microstructure of reservoir pore (Hou et al 2007;Arns et al 2004Arns et al , 2005Youssef et al 2008;Mustafa et al 2009).…”

Applications of digital core analysis and hydraulic flow units in petrophysical characterization

“…For high resolution requirements, FIB/SEM-first use FIB (focused ion beams) to make very fine slices through the samples (typically just a few µm across), that enable sequential SEM images to be obtained down to resolutions of 10 nm (Tomutsa et al, 2007;Lemmens et al, 2010).…”

A pore-scale model of two-phase flow in water-wet rock

“…With a due scaling, it produces a reasonable approximation of the capillary pressure curve [43]. The calculated fluid distribution yields a list of the voxel which are deemed to be occupied by water or oil.…”

“…However, a combination of flow simulations with Maximal Inscribed Spheres calculations provides an opportunity to compute relative permeability curves. MIS-calculations proved to be an efficient way of generating realistic capillary pressure curves [43,6]. It turns out that the MIS-computed relative permeability curves are in a good agreement with the coreflood measurements.…”