Dimensionally Scaled Miscible Displacements in Heterogeneous Permeable Media

Understanding and predicting the performance of solvent drives and remediation of contaminated aquifers in heterogeneous reservoirs is of great importance to the petroleum and environmental industries. In this paper, a general method to scale flow through heterogeneous reservoirs is presented for a miscible displacement of oil by a solvent. Results show that scaling miscible displacements in a two-dimensional, heterogeneous, anisotropic vertical cross section requires the matching of 13 dimensionless scaling groups. These groups were derived using a general procedure of inspectional analysis. A detailed numerical sensitivity study was performed to reveal the relationship between the scaling groups and the fractional oil recovery of miscible displacements in heterogeneous reservoirs. This relationship was then mapped using an artificial neural network, which can be used as a quick prediction tool for the fractional oil recovery for any combinations of the scaling groups, thus eliminating the need for the expensive fine-mesh simulations. These results have potential applications in modeling miscible displacements and in the scaling of laboratory displacements to field conditions.     

A relationship between the fractal dimension and scaling groups of unstable miscible displacements

A detailed two-dimensional flow visualization study was performed to examine the dynamics of viscous fingering in miscible displacements. Detailed quantitative miscible displacement experiments using a microcomputer-based imaging workstation on a variety of oil recovery fluid systems were performed. The effect of two dimensionless scaling groups, namely gravity number and viscosity ratio, on the displacement behavior was investigated. Based on image analysis, the irregular fingering patterns of the flow visualization experiments were analyzed for fractal characteristics. Results indicate that the areal sweep efficiency of unstable miscible displacement follows a fractal scaling law with a fractal dimension and proportionality constant related to the gravity number and the viscosity ratio. The study shows that the fractal dimension decreases with decreasing gravity number and increasing viscosity ratio. This relationship was mapped by an artificial neural network model, which can be used to estimate the fractal dimension and the proportionality constant of miscible displacements as functions of the two scaling groups. These results have potential application in the mathematical modeling of unstable EOR displacements and in the scaling of laboratory displacements to field conditions. Received: 19 December 1999/Accepted: 26 January 2001     

Experimental displacement patterns in a 2 × 2 quadrant block with permeability and wettability heterogeneities—problems for numerical modelling

The effect of heterogeneities on miscible and immiscible flood displacements in 2D bead packs in quadrant form, 2 × 2 block heterogeneity, with either a permeability or a wettability contrast is the subject of this paper. The physical processes occurring during miscible and immiscible flow and displacement within permeability and wettability quadrant bead pack models have been studied experimentally. This geometry occurs in a number of situations relevant to hydrocarbon production: particularly faults where adjacent rocks have large permeability contrasts with rapid changes, in the laboratory with core butting, in reservoir simulation where grid blocks have different permeability and in reservoirs having near-wellbore damage problems. The model quadrants 1–4, had 1 and 4 and 2 and 3 with identical properties, either in permeability or wettability. Reported are complete unit mobility miscible displacements, then the effects of viscosity differences (mobility modifiers) and finally immiscible displacements on displacement patterns for initial linear injection. The experiments demonstrate that nodal flow occurs for both miscible and immiscible flow, but for immiscible flow there are boundary effects due to capillary pressure differences created by water saturation changes or wettability contrasts which can leave patches of isolated fluid within a quadrant. The displacement patterns for the different models and fluids change significantly with the viscosity and wettability changes, particularly for the immiscible displacements. This is due to the changing capillary pressure between the quadrant blocks as the water saturation change. These are difficult to address in numerical modelling but should be accounted for. Other effects include coupling of all physical processes governing the flow through the node and creations of microzones of trapped residual oil. Our displacement patterns can therefore be a valuable verification benchmark tool for numerical modelling and a calibration data source for those wishing to simulate the effects of capillary pressure under differing wettability conditions and for those investigating upscaling modelling procedures. However, the possible loss of physical reality when averaging must always be considered.     

Pore-Scale Imaging and Analysis of Wettability Order,Trapping and Displacement in Three-Phase Flow in Porous Media with Various Wettabilities

Three-phase flow in porous media is encountered in many applications including subsurface carbon dioxide storage, enhanced oil recovery, groundwater remediation and the design of microfluidic devices. However, the pore-scale physics that controls three-phase flow under capillary dominated conditions is still not fully understood. Recent advances in three-dimensional pore-scale imaging have provided new insights into three-phase flow. Based on these findings, this paper describes the key pore-scale processes that control flow and trapping in a three-phase system, namely wettability order, spreading and wetting layers, and double/multiple displacement events. We show that in a porous medium containing water, oil and gas, the behaviour is controlled by wettability, which can either be water-wet, weakly oil-wet or strongly oil-wet, and by gas–oil miscibility. We provide evidence that, for the same wettability state, the three-phase pore-scale events are different under near-miscible conditions—where the gas–oil interfacial tension is ≤?1 mN/m—compared to immiscible conditions. In a water-wet system, at immiscible conditions, water is the most-wetting phase residing in the corners of the pore space, gas is the most non-wetting phase occupying the centres, while oil is the intermediate-wet phase spreading in layers sandwiched between water and gas. This fluid configuration allows for double capillary trapping, which can result in more gas trapping than for two-phase flow. At near-miscible conditions, oil and gas appear to become neutrally wetting to each other, preventing oil from spreading in layers; instead, gas and oil compete to occupy the centre of the larger pores, while water remains connected in wetting layers in the corners. This allows for the rapid production of oil since it is no longer confined to movement in thin layers. In a weakly oil-wet system, at immiscible conditions, the wettability order is oil–water–gas, from most to least wetting, promoting capillary trapping of gas in the pore centres by oil and water during water-alternating-gas injection. This wettability order is altered under near-miscible conditions as gas becomes the intermediate-wet phase, spreading in layers between water in the centres and oil in the corners. This fluid configuration allows for a high oil recovery factor while restricting gas flow in the reservoir. Moreover, we show evidence of the predicted, but hitherto not reported, wettability order in strongly oil-wet systems at immiscible conditions, oil–gas–water, from most to least wetting. At these conditions, gas progresses through the pore space in disconnected clusters by double and multiple displacements; therefore, the injection of large amounts of water to disconnect the gas phase is unnecessary. We place the analysis in a practical context by discussing implications for carbon dioxide storage combined with enhanced oil recovery before suggesting topics for future work.

     

Phase Behaviour,Fluid Properties and Recovery Efficiency of Immiscible and Miscible Condensate Displacements by SCCO<Subscript>2</Subscript> Injection: Experimental Investigation

This paper presents a quantitative investigation of the interfacial tension dependent relative permeability (IFT-DRP) and displacement efficiency of supercritical CO₂ injection into gas-condensate reservoirs. A high-pressure high-temperature experimental laboratory was established to simulate reservoir conditions and to perform relative permeability measurements on sandstone cores at a constant reservoir temperature of 95°C and displacement velocity of 10 cm/h. This investigation covers immiscible displacements (1100 and 2100 psi), near-miscible displacement (3000 psi) and miscible displacements (4500 and 5900 psi). The coreflooding results demonstrated that displacement pressure is a key factor governing the attainment of optimum sweep efficiency. The ultimate condensate recovery increased by almost threefold when CO₂ was injected at near-miscible conditions (i.e., 23.40% ultimate recovery at 1100 psi compared to 69.70% at 3000 psi). Miscible flooding was found to give the optimum condensate recovery (9% extra ultimate recovery compared to near-miscible injection). Besides improving the ultimate recovery, miscible floods provided better mobility ratios and delayed gas breakthrough (0.62 PV BT at 5900 psi compared to 0.21 PV BT at 1100 psi). In addition to the elimination of IFT forces in miscible displacements, favourable ratios of fluid properties and phase behaviour relationships between the SCCO₂ and condensate were believed to be the driving force for the improved recovery as they provided a stabilising effect on the displacement front and stimulated swelling of the condensate volume. This paper incorporates the theoretical aspects of phase behaviour and fluid properties that largely affect the microscopic displacement efficiency and serves as a practical guideline for operators to aid their project designs and enhance their recovery capabilities.     

Microscopic Mechanisms of Oil Recovery By Near-Miscible Gas Injection

As gas flooding becomes a more viable means of enhanced oil recovery, it is important to identify and understand the pore-scale flow mechanisms, both for the development of improved gas flooding applications and for the predicting phase mobilisation under secondary and tertiary gas flooding. The purpose of this study was to visually investigate the pore-level mechanisms of oil recovery by near-miscible secondary and tertiary gas floods. High-pressure glass micromodels and model fluids representing a near-miscible fluid system were used for the flow experiments. A new pore-scale recovery mechanism was identified which significantly contributed to oil recovery through enhanced flow and cross-flow between the bypassed pores and the injected gas. This mechanism is strongly related to a very low gas/oil interfacial tension (IFT), perfect wetting conditions and simultaneous flow of gas and oil in the same pore, all of which occur as the gas/oil critical point is approached. The results of this study helps us to better understand the pore-scale mechanisms of oil recovery in very low-IFT (near-miscible) systems. In particular we show that in near-miscible gas floods, behind the main gas front, the recovery of the oil continues by cross-flow from the bypassed pores into the main flow stream and as a result almost all of the oil, which has been contacted by the gas, could be recovered. Our observations in high-pressure micromodel experiments have demonstrated that this mechanism can only occur in near-miscible processes (as opposed to immiscible and completely miscible processes), which makes oil displacement by near-miscible gas floods a very effective process.     

Unstable Displacement of Non-aqueous Phase Liquids with Surfactant and Polymer

In this paper, we study two-phase multicomponent displacement of two immiscible fluids in both homogeneous and heterogeneous porous media. In many applications such as enhanced oil recovery, fluid mixing and spreading can be detrimental to the efficacy of the process. Here, we show that when an initially immobile phase is being displaced by a finite-size slug of solvents (surfactant and polymer), viscous fingering significantly enhances mixing and spreading of solvents. These effects are similar to those caused by medium heterogeneity and lead to poor displacement efficiency. We first quantify the displacement efficiency subject to different mobility ratios, Peclet numbers, and levels of medium heterogeneity. We observe a non-monotonic behavior in displacement efficiency as a function of mobility ratio, indicating that although stable frontal interface is desirable, miscible viscous fingering on the rear interface will eventually disintegrate the solvents slugs and reduce the displacement efficiency. Then, we show that miscible viscous fingering developing on the rear interface of the chemical slug could be greatly suppressed when viscosity contrast is gradually decreased using exponential or linear functions, leading to 10% increase in displacement efficiency while using the same amount of chemicals. To elucidate this low displacement efficiency, we study the evolution of mixing, spreading, and interfacial length and show that while higher viscosity ratios are quite effective in mobilizing the initially immobile phase in 1D displacements, they are in fact detrimental in 2D unstable displacements since they enhance mixing and spreading of solvents.

     

Immiscible Displacement in Cross-Bedded Heterogeneous Porous Media

Direct insight into the mechanisms of flow and displacements within small-scale (cm) systems having permeability heterogeneities that are not parallel to the flow direction (cross-bedding and fault zones) have been carried out. In our experiments, we have used visual models with unconsolidated glass bead packs having carefully controlled permeability contrasts to observe the processes with coloured fluids and streamlines. The displacements were followed visually and by video recording for later analysis. The experiments show the significance that heterogeneities have on residual saturations and recovery, as well as the displacement patterns themselves. During a waterflood, high permeability regions can be by-passed due to capillary pressure differences, giving rise to high residual oil saturations in these regions. This study demonstrates the importance of incorporating reservoir heterogeneity into core displacement analysis, but of course the nature of the heterogeneity has to be known. In general, the effects created by the heterogeneities and their unknown boundaries hamper interpretation of flood experiments in heterogeneous real sandstone cores. Our experiments, therefore, offer clear visual information to provide a firmer understanding of the displacement processes during immiscible displacement, to present benchmark data for input to numerical simulators, and to validate the simulator through a comparison with our experimental results for these difficult flow problems.     

Influence of interfacial reaction and mobility ratio on the displacement of oil in a Hele-Shaw cell

A comparative experimental study was conducted to examine the effect of fluid-fluid interfacial reaction on immiscible displacements in a wide range of mobility ratios. The nonreactive fluid pair consisted of paraffin oil/decane mixtures and water. In order to formulate a reactive system, linoleic acid (10 mol/m³) was added to the oil and sodium hydroxide (25 mol/m³) was added to the water. The experiments were performed in a square Hele-Shaw cell, simulating a quarter of a reversed five-spot pattern. Displacement patterns produced by reactive and nonreactive fluid pair systems were fundamentally different from each other. The recovery in the reactive system was generally higher than in the non-reactive system. The recovery at breakthrough for both the reactive and the non-reactive systems increased with a decrease in the viscosity of oil. In the non-reactive system, the recovery at one hour after breakthrough increased with the decrease in oil viscosity and became constant at 63% when the oil viscosity reached 6.86 mPa · s. In the reactive system the recovery at one hour after breakthrough was nearly 100% regardless of oil phase viscosity. When the oil viscosity was lower than the viscosity of the displacing phase (favorable mobility ratio), the recovery percentage was high. However, the jagged character of the displacement front characteristic of the reactive system was still preserved.     

Enhancing Immiscible Fluid Displacement in Porous Media by Capillary Pressure Discontinuities

Fluid displacement in porous media plays an important role in many industrial applications, including biological filtration, carbon capture and storage, enhanced oil recovery, and fluid transport in fuel cells. The displacement front is unstable, which evolves from smooth into ramified patterns, when the mobility (ratio of permeability to viscosity) of the displacing fluid is larger than that of the displaced one; this phenomenon is called viscous fingering. Viscous fingering increases the residual saturation of the displaced fluid, considerably impairing the efficacy of fluid displacement. It is of practical importance to develop suitable methods to improve fluid displacement. This paper presents an experimental study on applying the discontinuity of capillary pressure to improve immiscible fluid displacement in drainage for which the displacing fluid (air) wets the porous media less preferentially than does the displaced fluid (silicone oil). The concept involves using a heterogeneous packing system, where the upstream region features large pores and small capillary pressure, and the downstream region features small pores and large capillary pressure. The increase in capillary pressure prevents fingering from directly crossing the media interface, thus enhancing the displacement. The experimental apparatus was a linear cell comprising porous media between two parallel plates, and glass beads of 0.6 and 0.125 mm diameter were packed to compose the heterogeneous porous media. The time history of the finger flow was recorded using a video camera. Pressure drops over the model from the inlet to the outlet were measured to compare viscous pressure drops with capillary pressures. The results show that the fluid displacement was increased by the capillary discontinuities. The optimal displacement was determined through linear regression by adjusting the relative length of the large- and small-pore region. The results may assist in the understanding of fingering flow across the boundaries of different grain-sized bands for the gas and oil reservoir management, such as setting the relative location of the injection and production wells. The findings may also serve as a reference for industrial applications such as placing the grain bands in an adequate series to improve the displacement efficacy in biological filtration.     

Characterizing the Role of Shale Geometry and Connate Water Saturation on Performance of Polymer Flooding in Heavy Oil Reservoirs: Experimental Observations and Numerical Simulations

Many heavy oil reservoirs contain discontinuous shales which act as barriers or baffles to flow. However, there is a lack of fundamental understanding about how the shale geometrical characteristics affect the reservoir performance, especially during polymer flooding of heavy oils. In this study, a series of polymer injection processes have been performed on five-spot glass micromodels with different shale geometrical characteristics that are initially saturated with the heavy oil. The available geological characteristics from one of the Iranian oilfields were considered for the construction of the flow patterns by using a controlled-laser technology. Oil recoveries as a function of pore volumes of injected fluid were determined from analysis of continuously recorded images during the experiments. We observed a clear bypassing of displacing fluid which results in premature breakthrough of injected fluid due to the shale streaks. Moreover, the results showed a decrease of oil recovery when shales’ orientation, length, spacing, distance of the shale from production well, and density of shales increased. In contrast, an increase of shale discontinuity or distance of the shale streak from the injection well increased oil recovery. The obtained experimental data have also been used for developing and validating a numerical model where good matching performance has been observed between our experimental observations and simulation results. Finally, the role of connate water saturation during polymer flooding in systems containing flow barriers has been illustrated using pore level visualizations. The microscopic observations confirmed that besides the effect of shale streaks as heterogeneity in porous medium, when connate water is present, the trapped water demonstrates another source of disturbance and causes additional perturbations to the displacement interface leading to more irregular fingering patterns especially behind the shale streaks and also causes a reduction of ultimate oil recovery. This study reveals the application of glass micromodel experiments for studying the effects of barriers on oil recovery and flow patterns during EOR processes and also may provide a set of benchmark data for recovery of oil by immiscible polymer flood around discontinuous shales.     

A Stochastic Approach to the Two-Phase Displacement Problem in Heterogeneous Porous Media

The problem of two-phase immiscible flow in heterogeneous porous media in the case of a horizontal displacement of some fluid by another, which is of practical importance in industrial oil recovery, is considered. Assuming that (a) the saturation jump on the displacement front is constant, (b) the log-permeability of the medium obeys Gaussian statistics, and (c) the case when the front is stable, the displacement front position and the saturation distribution are described analytically in terms of generalized functions. Note that in our analysis we do not assume that the front shape fluctuations are small, and in this respect our results may be regarded as exact. The assumption that the log-permeability fluctuations are small was only used in deriving the linear relation between the log-permeability of a porous medium and the total flow velocity (Nœtinger et al. in Fluid Dyn 41(5):830–842, 2006). By means of ensemble averaging, the mean saturation and saturation variance are found in the vicinity of the front. These characteristics are related to the variance of front displacements, which, in turn, can be calculated analytically. Next, a method for reconstructing the full solution for the saturation (rarefaction wave) is proposed. Such a full solution satisfies the mass conservation requirement. Finally, the theoretical predictions are compared with the results of numerical simulations carried out within the framework of Monte-Carlo method.     

Interpretation of pressure recovery curves recorded simultaneously in different sections of a horizontal well borehole

A new approach to carrying out and interpreting hydrodynamic investigations of horizontal wells is proposed. The problem of determining the flow parameters of an inhomogeneous oil reservoir is considered. As the initial information, the pressure recovery curves recorded simultaneously in different sections of the horizontal well borehole are used.     

An Experimental Study of Three-Phase Trapping in Sand Packs

The trapped saturations of oil and gas are measured as functions of initial oil and gas saturation in water-wet sand packs. Analogue fluids—water, octane and air—are used at ambient conditions. Starting with a sand-pack column which has been saturated with brine, oil (octane) is injected with the column horizontal until irreducible water saturation is reached. The column is then positioned vertically and air is allowed to enter from the top of the column, while oil is allowed to drain under gravity for varying lengths of time. At this point, the column may be sliced and the fluids analyzed by gas chromatography to obtain the initial saturations. Alternatively, brine is injected through the bottom of the vertical column to trap oil and gas, before slicing the columns and measuring the trapped or residual saturations by gas chromatography and mass balance. The experiments show that in three-phase flow, the total trapped saturations of oil and gas are considerably higher than the trapped saturations reported in the literature for two-phase systems. It is found that the residual saturation of oil and gas combined could be as high as 23 %, as opposed to a maximum two-phase residual of only 14 %. For very high initial gas saturations, the residual gas saturation, up to 17 %, was also higher than for two-phase displacement. These observations are explained in terms of the competition between piston-like displacement and snap-off. It is also observed that less oil is always trapped in three-phase flow than in two-phase displacement, and the difference depends on the amount of gas present. For low and intermediate initial gas saturations, the trapped gas saturation rises linearly with initial saturation, followed by a constant residual, as seen in two-phase displacements. However, at very high initial gas saturations, the residual saturation rises again.     

Investigating the Fracture Network Effects on Sweep Efficiency during WAG Injection Process

In this study, the main recovery mechanisms behind oil/water/gas interactions during the water-alternating-gas (WAG) injection process, in a network of matrix/fracture, were fundamentally investigated. A visual micromodel was utilized to provide insights into the potential applications of WAG process in fractured oil-wet media as well as the possibility of observing microscopic displacement behavior of fluids in the model. The model was made of an oil-wet facture/matrix network system, comprised of four matrix blocks surrounded with fractures. Different WAG injection scenarios, such as slug arrangements and the effects of fluid injection rates on oil recovery were studied. A new equation representing the capillary number, considering the fracture viscous force and matrix capillary force, was developed to make the experimental results more similar to a real field. In general, WAG tests performed in the fractured model showed a higher oil recovery factor compared with the results of gas and water injection tests at their optimum rates. The results showed that the presence of an oil film, in all cases, was the main reason for co-current drainage and double displacement of oil under applied driving forces. Furthermore, the formation of oil liquid bridges improved the recovery efficiency, which was greatly influenced by the size of fracture connecting the two matrix blocks; these connecting paths were more stable when there was initial water remaining in the media. Analyzing different recovery curves and microscopic view of the three phases in the transparent model showed that starting an injection mode with gas (followed by repeated small slugs of water and gas), could considerably improve oil recovery by pushing water into the matrix zone and increasing the total sweep efficiency.     

A Semi-Analytical Model for Two Phase Immiscible Flow in Porous Media Honouring Capillary Pressure

This article describes a semi-analytical model for two-phase immiscible flow in porous media. The model incorporates the effect of capillary pressure gradient on fluid displacement. It also includes a correction to the capillarity-free Buckley–Leverett saturation profile for the stabilized-zone around the displacement front and the end-effects near the core outlet. The model is valid for both drainage and imbibition oil–water displacements in porous media with different wettability conditions. A stepwise procedure is presented to derive relative permeabilities from coreflood displacements using the proposed semi-analytical model. The procedure can be utilized for both before and after breakthrough data and hence is capable to generate a continuous relative permeability curve unlike other analytical/semi-analytical approaches. The model predictions are compared with numerical simulations and laboratory experiments. The comparison shows that the model predictions for drainage process agree well with the numerical simulations for different capillary numbers, whereas there is mismatch between the relative permeability derived using the Johnson–Bossler–Naumann (JBN) method and the simulations. The coreflood experiments carried out on a Berea sandstone core suggest that the proposed model works better than the JBN method for a drainage process in strongly wet rocks. Both methods give similar results for imbibition processes.     

Effect of Well Orientation (Vertical vs. Horizontal) and Well Length on the Injection of CO<Subscript>2</Subscript> in Deep Saline Aquifers

Simulations of CO₂ injection into confined saline aquifers were conducted for both vertical and horizontal injection wells. The metrics used in quantifying the performances of different injection scenarios included changes in pressure near the injection well, mass of CO₂ dissolved into brine (solubility trapping), and storage efficiency, all evaluated with an assumed injection period of 50 years. Metrics were quantified as functions of well length, well orientation, CO₂ injection rate, and formation anisotropy (ratio of vertical to horizontal conductivity). When equal well lengths are compared, there is not a significant difference between the predicted performances of horizontal and vertical wells. However, the length of a horizontal well may exceed the length of a vertical well because the length of the horizontal well is not constrained to the vertical thickness of the geologic formation. Simulations show that, as the length of the horizontal well is allowed to increase, the geologic formation can receive a significantly higher CO₂ injection rate without exceeding a maximum allowable pressure. This result is observed in both isotropic and anisotropic formations, and suggests that horizontal wells may be advantageous under pressure-limited conditions. However, the use of horizontal wells does not significantly improve the storage efficiency, and under strongly anisotropic conditions, a vertical well provides higher storage efficiency than a horizontal well. We conclude that horizontal wells may be preferable if the goal is to sequester a large amount of CO₂ in a short period of time, but do not offer a significant advantage in terms of long-term capacity of a potential repository.     

A new simple third-order shear deformation theory of plates

An improved simple third-order shear deformation theory for the analysis of shear flexible plates is presented in this paper. This new plate theory is composed of three parts: the simple third-order kinematics of displacements reduced from the higher-order displacement field derived previously by the author; a system of 10th-order differential equilibrium equations in terms of the three generalized displacements of bending plates; five boundary conditions at each edge of plate boundaries. Although the resulting displacement field is the same as that proposed by Murthy, the variational consistent governing equations and the associated proper boundary conditions are derived and identified in this work for the first time in the literature. The applications and accuracy of the present shear deformation theory of plates are demonstrated by analytically solving the differential governing equations of a twisting plate, a bending beam and two bending plates to which the 3-D elasticity solutions are available, and excellent agreements are achieved even for the torsion of a plate with square cross-section as well the local effects of stresses at plate boundaries can be characterized accurately. These analytical solutions clearly show that the simple third-order shear deformation theory developed in this work indeed gives better results than the first-order shear deformation theories and other simple higher-order shear deformation theories, since the present third-order shear flexible theory is based on a more rigorous kinematics of displacements and consists of not only a system of variational consistent differential equations, but also a group of consistent boundary conditions associated with the differential equations. The present simple third-order shear deformation theory can easily be applied to the static and dynamic finite element analysis of laminated plates just like the applications of other popular shear flexible plate theories, and improved results could be obtained from the present simple third-order shear deformable theories of plates.     

High Resolution Imaging of Unstable,Forced Imbibition in Berea Sandstone

An experimental investigation is presented of immiscible, high-mobility ratio forced imbibition in a representative linear homogeneous sandstone. Water floods with mobility ratios (from 1 to 155) at various water injection rates were conducted. Fine-scale (order mm³) in situ water saturation history was collected via X-ray computed tomography (CT). Three-dimensional images were constructed of stable displacement and the initiation and growth of unstable water fingers. Interestingly, viscous fingers do not lead the displacement front by significant distances, counter to experience in miscible systems. In this homogeneous porous medium, both water (displacing phase) injection rate and oil (displaced phase) viscosity have an obvious effect on the stability of the water front. As the oil viscosity and displacement rate increase, the water front becomes less stable. In addition, the so-called shock mobility ratio, as computed from steady-state relative permeability, is found to be predictive regarding displacement front stability. When the shock mobility ratio is greater than 1, the displacement is always unstable. Steady-state relative permeability, however, is found to be a function of viscosity ratio for unstable displacements.