Pore-scale micro-computed-tomography imaging: nonwetting-phase cluster-size distribution during drainage and imbibition

Study on the Occurrence Characteristics of the Remaining Oil in Sandstone Reservoirs with Different Permeability after Polymer Flooding

After polymer flooding, the heterogeneity between different layers intensifies, forming intricate seepage channels and fluid diversions, which results in decreased circulation efficiency and lower recovery rates, leaving a significant amount of residual oil trapped within the reservoir. Understanding the characteristics of residual oil occurrence is crucial for enhancing oil recovery post-polymer flooding. This study focused on sandstone reservoirs with varying permeability in the Saertu block of the Daqing oilfield. Using cryosectioning and laser scanning confocal microscopy, the occurrence characteristics of the residual oil in these sandstone reservoirs post-polymer flooding were investigated. Additionally, micro-CT and scanning electron microscopy were employed to analyze the impact of the pore structure on the distribution characteristics of the residual oil. The results indicate that laser scanning confocal images reveal that post-polymer flooding, the residual oil in high- and low-permeability sandstone reservoirs predominantly exists in a bound state (average > 47%), mostly as particle-adsorbed oil. In contrast, the residual oil in medium-permeability reservoirs is primarily in a free state (average > 49%), mostly as intergranular-adsorbed oil. In high-permeability sandstone reservoirs, heavy oil components are mainly in a particle-adsorbed form; in medium-permeability sandstone reservoirs, residual oil predominantly consists of heavy components, with most light components occurring in a clustered form; in low-permeability sandstone reservoirs, clustered residual oil exists in a balanced coexistence of light and heavy components, while the heavy components primarily exist in a particle-adsorbed form. Post-polymer flooding, the large pore-throat structure in high-permeability sandstone reservoirs results in effective displacement and less free residual oil; medium-permeability sandstone reservoirs, with medium-large pores and throats, have preferential channels and fine particles blocking the throats, leading to some unswept pores and more free residual oil; low-permeability sandstone reservoirs, with small pores and throats, exhibit weak displacement forces and poor mobility, resulting in more bound residual oil. The distribution and content of clay particles and clay minerals, along with the complex microscopic pore structure, are the main factors causing the differences in the residual oil occurrence states in sandstones with varying permeability.

Study on Oil Recovery Mechanism of Polymer-Surfactant Flooding Using X-ray Microtomography and Integral Geometry

Understanding pore-scale morphology and distribution of remaining oil in pore space are of great importance to carry out in-depth tapping of oil potential. Taking two water-wet cores from a typical clastic reservoir in China as an example, X-ray CT imaging is conducted at different experimental stages of water flooding and polymer-surfactant (P-S) flooding by using a high-resolution X-ray microtomography. Based on X-ray micro-CT image processing, 3D visualization of rock microstructure and fluid distribution at the pore scale is achieved. The integral geometry newly developed is further introduced to characterize pore-scale morphology and distribution of remaining oil in pore space. The underlying mechanism of oil recovery by P-S flooding is further explored. The results show that the average diameter of oil droplets gradually decreases, and the topological connectivity becomes worse after water flooding and P-S flooding. Due to the synergistic effect of “1 + 1 > 2” between the strong sweep efficiency of surfactant and the enlarged swept volume of the polymer, oil droplets with a diameter larger than 124.58 μm can be gradually stripped out by the polymer-surfactant system, causing a more scattered distribution of oil droplets in pore spaces of the cores. The network-like oil clusters are still dominant when water flooding is continued to 98% of water cut, but the dominant pore-scale oil morphology has evolved from network-like to porous-type and isolated-type after P-S flooding, which can provide strong support for further oil recovery in the later stage of chemical flooding.

Impact of grain roughness on residual nonwetting phase cluster size distribution in packed columns of uniform spheres

We imaged the pore-scale distribution of air and water within packed columns of glass spheres of different textures using x-ray microcomputed tomography after primary drainage and after secondary imbibition. Postimbibition residual air saturation increases with roughness size. Clusters larger than a critical size of about 15 to 40 pores are distributed according to a power law, with exponents ranging from τ=2.29±0.04 to 3.00±0.13 and displaying a weak negative correlation with roughness size. The largest cluster constitutes 7 to 20% of the total residual gas saturation, with no clear correlation with roughness size. These results imply that activities that enhance grain roughness by, e.g., creating acidic conditions in the subsurface, will promote capillary trapping of nonwetting phases under capillary-dominated conditions. Enhanced trapping, in turn, may be desirable in some engineering applications such as geological CO_{2} storage, but detrimental to others such as groundwater remediation and hydrocarbon recovery.

Linking continuum-scale state of wetting to pore-scale contact angles in porous media

HYPOTHESIS

Wetting phenomena play a key role in flows through porous media. Relative permeability and capillary pressure-saturation functions show a high sensitivity to wettability, which has different definitions at the continuum- and pore-scale. We hypothesize that the wetting state of a porous medium can be described in terms of topological arguments that constrain the morphological state of immiscible fluids, which provides a direct link between the continuum-scale metrics of wettability and pore-scale contact angles.

EXPERIMENTS

We perform primary drainage and imbibition experiments on Bentheimer sandstone using air and brine. Topological properties, such as Euler characteristic and interfacial curvature are measured utilizing X-ray micro-computed tomography at irreducible air saturation. We also present measurements for the United States Bureau of Mines (USBM) index, capillary pressure and pore-scale contact angles. Additional studies are performed using two-phase Lattice Boltzmann simulations to test a wider range of wetting conditions.

FINDINGS

We demonstrate that contact angle distributions for a porous multiphase system can be predicted within a few percent difference of directly measured pore-scale contact angles using the presented method. This provides a general framework on how continuum-scale data can be used to describe the geometrical state of fluids within porous media.

Velocity distributions in trapped and mobilized non-wetting phase ganglia in porous media

Understanding the mobilisation of trapped globules of non-wetting phase during two-phase flow has been the aim of numerous studies. However, the driving forces for the mobilisation of the trapped phases are still not well understood. Also, there is little information about what happens within a globule before, at the onset and during mobilization. In this work, we used micro-particle tracking velocimetry in a micro-fluidic model in order to visualise the velocity distributions inside the trapped phase globules prior and during mobilisation. Therefore, time-averaged and instantaneous velocity vectors have been determined using fluorescent microscopy. As a porous medium, we used a polydimethylsiloxane (PDMS) micro-model with a well-defined pore structure, where drainage and imbibition experiments were conducted. Three different geometries of trapped non-wetting globules, namely droplets, blobs and ganglia were investigated. We observed internal circulations inside the trapped phase globules, leading to the formation of vortices. The direction of circulating flow within a globule is dictated by the drag force exerted on it by the flowing wetting phase. This is illustrated by calculating and analyzing the drag force (per unit area) along fluid-fluid interfaces. In the case of droplets and blobs, only one vortex is formed. The flow field within a ganglion is much more complex and more vortices can be formed. The circulation velocities are largest at the fluid-fluid interfaces, along which the wetting phase flows and decreases towards the middle of the globule. The circulation velocities increased proportionally with the increase of wetting phase average velocity (or capillary number). The vortices remain stable as long as the globules are trapped, start to change at the onset of mobilization and disappear during the movement of globules. They reappear when the globules get stranded. Droplets are less prone to mobilization; blobs get mobilised in whole; while ganglia may get ruptured and get mobilised only partially.

High pressure-elevated temperature x-ray micro-computed tomography for subsurface applications

Physical, chemical and mechanical pore-scale (i.e. micrometer-scale) mechanisms in rock are of key importance in many, if not all, subsurface processes. These processes are highly relevant in various applications, e.g. hydrocarbon recovery, CO₂ geo-sequestration, geophysical exploration, water production, geothermal energy production, or the prediction of the location of valuable hydrothermal deposits. Typical examples are multi-phase flow (e.g. oil and water) displacements driven by buoyancy, viscous or capillary forces, mineral-fluid interactions (e.g. mineral dissolution and/or precipitation over geological times), geo-mechanical rock behaviour (e.g. rock compaction during diagenesis) or fines migration during water production, which can dramatically reduce reservoir permeability (and thus reservoir performance). All above examples are 3D processes, and 2D experiments (as traditionally done for micro-scale investigations) will thus only provide qualitative information; for instance the percolation threshold is much lower in 3D than in 2D. However, with the advent of x-ray micro-computed tomography (μCT) - which is now routinely used - this limitation has been overcome, and such pore-scale processes can be observed in 3D at micrometer-scale. A serious complication is, however, the fact that in the subsurface high pressures and elevated temperatures (HPET) prevail, due to the hydrostatic and geothermal gradients imposed upon it. Such HPET-reservoir conditions significantly change the above mentioned physical and chemical processes, e.g. gas density is much higher at high pressure, which strongly affects buoyancy and wettability and thus gas distributions in the subsurface; or chemical reactions are significantly accelerated at increased temperature, strongly affecting fluid-rock interactions and thus diagenesis and deposition of valuable minerals. It is thus necessary to apply HPET conditions to the aforementioned μCT experiments, to be able to mimic subsurface conditions in a realistic way, and thus to obtain reliable results, which are vital input parameters required for building accurate larger-scale reservoir models which can predict the overall reservoir-scale (hectometer-scale) processes (e.g. oil production or diagenesis of a formation). We thus describe here the basic workflow of such HPET-μCT experiments, equipment requirements and apparatus design; and review the literature where such HPET-μCT experiments were used and which phenomena were investigated (these include: CO₂ geo-sequestration, oil recovery, gas hydrate formation, hydrothermal deposition/reactive flow). One aim of this paper is to give a guideline to users how to set-up a HPET-μCT experiment, and to provide a quick overview in terms of what is possible and what not, at least up to date. As a conclusion, HPET-μCT is a valuable tool when it comes to the investigation of subsurface micrometer-scaled processes, and we expect a rapidly expanding usage of HPET-μCT in subsurface engineering and the subsurface sciences.

X-ray Microtomography of Intermittency in Multiphase Flow at Steady State Using a Differential Imaging Method

We imaged the steady state flow of brine and decane in Bentheimer sandstone. We devised an experimental method based on differential imaging to examine how flow rate impacts impact the pore-scale distribution of fluids during coinjection. This allows us to elucidate flow regimes (connected, or breakup of the nonwetting phase pathways) for a range of fractional flows at two capillary numbers, Ca, namely 3.0 × 10-7 and 7.5 × 10-6. At the lower Ca, for a fixed fractional flow, the two phases appear to flow in connected unchanging subnetworks of the pore space, consistent with conventional theory. At the higher Ca, we observed that a significant fraction of the pore space contained sometimes oil and sometimes brine during the 1 h scan: this intermittent occupancy, which was interpreted as regions of the pore space that contained both fluid phases for some time, is necessary to explain the flow and dynamic connectivity of the oil phase; pathways of always oil-filled portions of the void space did not span the core. This phase was segmented from the differential image between the 30 wt % KI brine image and the scans taken at each fractional flow. Using the grey scale histogram distribution of the raw images, the oil proportion in the intermittent phase was calculated. The pressure drops at each fractional flow at low and high flow rates were measured by high-precision differential pressure sensors. The relative permeabilities and fractional flow obtained by our experiment at the mm-scale compare well with data from the literature on cm-scale samples.

Characterization of reactive flow-induced evolution of carbonate rocks using digital core analysis - part 2: Calculation of the evolution of percolation and transport properties

Percolation of reactive fluids in carbonate rocks affects the rock microstructure and hence changes the rock macroscopic properties. In Part 1 paper, we examined the voxel-wise evolution of microstructure of the rock in terms of mineral dissolution/detachment, mineral deposition, and unchanged regions. In the present work, we investigate the relationships between changes in two characteristic transport properties, i.e. permeability and electrical conductivity and two critical parameters of the pore phase, i.e. the fraction of the pore space connecting the inlet and outlet faces of the core sample and the critical pore-throat diameter. We calculate the aforementioned properties on the images of the sample, wherein a homogeneous modification of pore structure occurred in order to ensure the representativeness of the calculated transport properties at the core scale. From images, the evolution of pore connectivity and the potential role of micropores on the connectivity are quantified. It is found that the changing permeability and electrical conductivity distributions along the core length are generally in good agreement with the longitudinal evolution of macro-connected macroporosity and the critical pore-throat diameter. We incorporate microporosity into critical length and permeability calculations and show how microporosity locally plays a role in permeability. It is shown that the Katz-Thompson model reasonably predicts the post-alteration permeability in terms of pre-alteration simulated parameters. This suggests that the evolution of permeability and electrical conductivity of the studied complex carbonate core are controlled by the changes in the macro-connected macroporosity as well as the smallest pore-throats between the connected macropores.

Influence of phase connectivity on the relationship among capillary pressure, fluid saturation, and interfacial area in two-fluid-phase porous medium systems

Multiphase flows in porous medium systems are typically modeled at the macroscale by applying the principles of continuum mechanics to develop models that describe the behavior of averaged quantities, such as fluid pressure and saturation. These models require closure relations to produce solvable forms. One of these required closure relations is an expression relating the capillary pressure to fluid saturation and, in some cases, other topological invariants such as interfacial area and the Euler characteristic (or average Gaussian curvature). The forms that are used in traditional models, which typically consider only the relationship between capillary pressure and saturation, are hysteretic. An unresolved question is whether the inclusion of additional morphological and topological measures can lead to a nonhysteretic closure relation. Relying on the lattice Boltzmann (LB) method, we develop an approach to investigate equilibrium states for a two-fluid-phase porous medium system, which includes disconnected nonwetting phase features. A set of simulations are performed within a random close pack of 1964 spheres to produce a total of 42 908 distinct equilibrium configurations. This information is evaluated using generalized additive models to quantitatively assess the degree to which functional relationships can explain the behavior of the equilibrium data. The variance of various model estimates is computed, and we conclude that, except for the limiting behavior close to a single fluid regime, capillary pressure can be expressed as a deterministic and nonhysteretic function of fluid saturation, interfacial area between the fluid phases, and the Euler characteristic. To our knowledge, this work is unique in the methods employed, the size of the data set, the resolution in space and time, the true equilibrium nature of the data, the parametrizations investigated, and the broad set of functions examined. The conclusion of essentially nonhysteretic behavior provides support for an evolving class of two-fluid-phase flow in porous medium systems models.

CO2/water two-phase flow in a two-dimensional micromodel of heterogeneous pores and throats

Small capillary fingerings can be collapsed by capillary pressure in the beginning of displacement. Larger fingerings tend to form with developing displacement and cannot be overcome anymore.

Capillary saturation and desaturation

Capillary desaturation experiments produce disconnected (trapped) ganglia of mesoscopic sizes intermediate between pore size and system size. Experimental evidence for interactions between these mesoscale clusters during desaturation is analyzed and discussed within the established microscopic and macroscopic laws of Newton, Young-Laplace, and Darcy. A theoretical expression for capillary number correlations is introduced that seems to have remained unnoticed. It expresses capillary desaturation curves in terms of stationary capillary pressures and relative permeabilities. The theoretical expression shows that the plateau saturation in capillary desaturation curves may in general differ from the residual nonwetting saturation defined through the saturation limit of the main hysteresis loop. Hysteresis effects as well as the difference between wetting and nonwetting fluids are introduced into the analysis of capillary desaturation experiments. The article examines experiments with different desaturation protocols and discusses the existence of a mesoscopic length scale intermediate between pore scale and sample scale. The theoretical expression is derived entirely within the existing traditional theory of two-phase flow in porous media and compared to a recent experiment.