Hybrid Mixed Methods Applied to Miscible Displacements with Adverse Mobility Ratio

Abstract Areal sweep efficiency of oil displacement by enhanced-viscosity water exhibiting pseudoplastic behavior was measured in a Hele-Shaw model representing one-quarter of a five-spot pattern. The pseudoplasticity of polymer solutions and the velocity distribution in the five-spot pattern produced a condition under which the mobility ratio between the displacing and the displaced fluid could not be assigned a single value. Instead, the movement of the displacement front is governed by local mobility ratios which are also time dependent. The areal sweep at breakthrough with polymer solutions was poorer than the sweep obtained with Newtonian fluids of comparable viscosity. However, the areal sweep and 1 PV throughput was greatly improved as compared to flood water without polymer. It was also demonstrated that, even after the oil-cut had declined to a low value during a regular waterflood, switching to polymer flood efficiently swept out the oil remaining in the model. Introduction The behavior of fluid displacements in isotropic porous media for various patterns of injection and production wells has been extensively investigated. These investigations all concerned Newtonian fluids, i.e., the viscosity of each fluid was constant regardless of flow rate. The generally unfavorable influence on areal sweep efficiency of higher mobility of the displacing fluid as compared to the mobility of the displaced fluid has been established for both miscible and immiscible fluids. The principle was also established that a close correspondence exists between miscible and immiscible flood front behavior, although oil recovery in a waterflood at unfavorable mobility ratio may be less than that observed in a miscible displacement at the same mobility ratio. This is true even when oil recovery is expressed on the basis of movable oil. The reason is that oil saturation only slowly achieves its final value behind the waterflood front in accordance with the Buckley-Leverett simultaneous flow relations. It is convenient to use miscible displacements for laboratory simulation of waterflood frontal advance since the interfacial tension forces which are negligible in proportion to viscous forces on a reservoir scales are thus made nonoperative in the laboratory model. For miscible displacements, the Hele-Shaw type of model adequately represents a porous medium so long as the appropriate scaling rules are observed in its design and operation. During simulation of waterflood front behavior in the laboratory by using miscible displacements, the behavior of connate water may ordinarily be disregarded since it is usually indistinguishable from flood water in this process. However, when the flood water is deliberately thickened to improve the mobility ratio between water and oil, the effect on the sweep efficiency due to generation of a connate water bank during the process must be considered. In a uniform porous medium, such a bank is generated and efficiently displaced by injection of thickened water. The oil originally in-place at the start of the waterflood is then displaced by connate water followed by thickened water. If the flood water must be thickened to obtain a favorable mobility ratio, the mobility of the oil phase is appreciably less than that of the connate water. Hence, the oil phase is inefficiently displaced by the connate water bank, and a considerable proportion of the oil comes in contact with and is displaced by the thickened waterflood front. SPEJ P. 52ˆ

Download Full-text

Systematic Modelling of Unstable Displacement

10.2118/207493-ms ◽

2021 ◽

Author(s):

Arne Skauge ◽

Kenneth Stuart Sorbie ◽

Iselin Cecilie Salmo ◽

Tormod Skauge

Keyword(s):

Single Point ◽

Viscous Fingering ◽

Mobility Ratio ◽

Fluid Viscosity ◽

Fractional Flow ◽

Fine Grid ◽

Rock Heterogeneity ◽

Unstable Displacement ◽

The Impact ◽

Adverse Mobility Ratio

Abstract Modelling unstable displacement is a challenge which may lead to large errors in reservoir simulations. Field scale coarse grid simulations therefore need to be anchored to more reliable fine grid models which capture fluid displacement instabilities in a physically correct manner. In this paper, a recently developed approach for accurately modelling viscous fingering has been applied to various types of unstable displacement. The method involves estimation of dispersivity of the porous medium and length scale of the model to determine the required size of the simulation grid cell. Fractional flow theory is then applied to obtain the correct saturation of the injected phase in the unstable fingers formed due to the adverse mobility ratio. Unstable displacement experiments have been history matched using 2D-imaging of in-situ saturation as a calibration of our method, before carrying out sensitivity calculations on the effect of fluid viscosity, and rock heterogeneity. Our modelling approach allows us to carry out simulations using a conventional numerical simulator using elementary numerical methods (e.g. single-point upstreaming). The methods used to model instability (Sorbie et al, 2020) was originally developed for immiscible water/oil systems. The current paper now presents new results applying this approach to unstable gas displacements, where adverse viscosity ratios may be even higher than in water/oil systems. The displacement with injected gas is shown to be influenced by mass exchanges between the gas and oil as the alternating fluids (water and gas) are injected in WAG processes. Swelling of fingers delay the gas front and WAG processes divert the injected gas and improve sweep efficiency. We have also modelled water-oil displacement at adverse mobility and shown the benefit which is obtained by reducing the instability by adding polymers to viscosify the injected water. The impact of rock heterogeneity has different effect depending on buoyancy forces and the degree of crossflow into the high permeable zones. This paper extends our novel approach to modelling the fine scale distribution of the injected fluids in adverse mobility ratio displacements. This approach has now been applied to both, gas/oil and water/oil systems where viscous fingering is present, either at extremely adverse mobility ratios and/or for reservoirs where the permeability field is very heterogeneous.

Download Full-text

A Finite-Element Method for Reservoir Simulation

Society of Petroleum Engineers Journal ◽

10.2118/7413-pa ◽

1981 ◽

Vol 21 (01) ◽

pp. 115-128 ◽

Cited By ~ 29

Author(s):

Larry C. Young

Keyword(s):

Finite Element ◽

Finite Element Methods ◽

Reservoir Simulation ◽

Finite Differences ◽

Time Derivative ◽

Mobility Ratio ◽

Difference Approximation ◽

Central Difference ◽

Fine Grid ◽

Adverse Mobility Ratio

Abstract Several previous studies have applied finite-element methods to reservoir simulation problems. Accurate solutions have been demonstrated with these methods; however, competitiveness with finite difference has not been established for most nonlinear reservoir simulation problems. In this study a more efficient finite-element procedures is presented and tested. The method is Galerkin-based, and improved efficiency is obtained by combining Lagrange trial functions with Lobatto quadrature in a particular way. The simulation of tracer performance, ion exchange preflush performance, and adverse mobility ratio miscible displacements is considered. For the problems considered, the method is shown to yield accurate solutions with less computing expense than finite differences or previously proposed finite-element techniques. For the special case of linear trial functions, the method reduces to a five-point central difference approximation. In contrast to previously reported results, this approximation is found to simulate adverse mobility ratio displacements without grid orientation sensitivity, provided a sufficiently fine grid is used. Introduction In the past few years several studies have investigated the use of finite-element methods in reservoir simulation. These include single-phase two-component simulations in one1–3 and two4,5 spatial dimensions and two-phase immiscible calculations in both one6–8 and two9,10 dimensions. These studies have demonstrated that the method is capable of giving accurate solutions, particularly for small slug problems and adverse mobility ratio displacements. All these studies used what we term conventional Galerkin finite-element techniques,1 and, unfortunately, these methods have not proved to be cost competitive with finite differences for most nonlinear reservoir simulation problems. A reduction in computing requirements is, therefore, necessary to make finite-element methods truly useful for reservoir simulation. Relative to finite differences, the increased computing requirements of conventional Galerkin-based methods are due to the following.The approximation of time-derivative terms involves the same number of surrounding grid points as the approximation of flow terms; thus, implicit-pressure/explicit-saturation (IMPES) techniques are not possible (see Ref. 12, Chap. 7).The matrices which result from the approximation of flow terms are not nearly so sparse as in finite differences; thus, the solution of matrix problems requires more computation.The computational work required to generate matrix coefficients is considerably greater than with finite differences due to the number of numerical integrations which must be performed.

Download Full-text