The Onset of Instability During Two-Phase Immiscible Displacement in Porous Media

1981 ◽  
Vol 21 (02) ◽  
pp. 249-258 ◽  
Author(s):  
Ekwere J. Peters ◽  
Donald L. Flock

Abstract This paper presents a dimensionless number and its critical value for predicting the onset of instability during immiscible displacement in porous media. The critical dimensionless number obtained from a stability theory for a cylindrical system successfully predicted the onset of instability in laboratory floods. Therefore, this number can be used to classify the stability of two-phase incompressible displacements in homogeneous porous media. Introduction When a fluid displaces a more viscous fluid, the displacement front may become unstable, resulting in viscous fingering. This phenomenon raises both practical and theoretical concerns. Apart from further reducing the displacement efficiency of an already inefficient displacement arrangement, instability may invalidate the usual method of simulating immiscible displacement performance based on relative permeability and capillary pressure concepts. Also, it introduces an additional scaling requirement for using model tests to forecast prototype displacement results. Therefore, it would be most beneficial to predict the onset of instability, so as to avoid viscous fingering, or, where it is unavoidable, to be able to recognize it as a factor in the displacement.The onset of instability call be predicted by a stability analysis of the displacement. The objective of such an analysis is to determine the conditions under which small disturbances or perturbations of the displacement front will grow to become viscous fingers. Ideally, the analysis should give a universal dimensionless scaling group together with its critical value above which instability will occur. The stability classification then would entail no more than the calculation of one dimensionless number in a manner analogous to the calculation of a Reynolds number to distinguish between laminar and turbulent flow.Several stability studies of immiscible displacement have been reported in the literature. Collectively, they show that these variables are pertinent to the stability problem:mobility (or viscosity) ratio,displacement velocity, system geometry and dimensions,capillary and gravitational forces, andsystem permeability and wettability. However, none of the previous studies have combined these variables into one dimensionless number that can be used to quantify the stability classification.The objective of this study was to obtain, by means of a stability analysis, a universal dimensionless scaling group and its critical value for predicting the onset of instability during immiscible displacement in porous media. This paper shows how the stability theory of Chuoke et al. was extended to achieve this objective and presents the results of laboratory floods that confirm the predicted onset of instability in cylindrical cores. Theory The pertinent dimensionless number for predicting the onset of instability was obtained by extending the stability theory of Chuoke et al. Their theory was based on a piston-like unperturbed displacement model in which the oil and water zones are separated by a planar interface. Details of the theory and our extension of it are presented in the following sections. SPEJ P. 249^

1977 ◽  
Vol 17 (01) ◽  
pp. 79-91 ◽  
Author(s):  
D.W. Peaceman

Abstract The usual linearized stability analysis of the finite-difference solution for two-phase flow in porous media is not delicate enough to distinguish porous media is not delicate enough to distinguish between the stability of equations using semi-implicit mobility and those using completely implicit mobility. A nonlinear stability analysis is developed and applied to finite-difference equations using an upstream mobility that is explicit, completely implicit, or semi-implicit. The nonlinear analysis yields a sufficient (though not necessary) condition for stability. The results for explicit and completely implicit mobilities agree with those obtained by the standard linearized analysis; in particular, use of completely implicit mobility particular, use of completely implicit mobility results in unconditional stability. For semi-implicit mobility, the analysis shows a mild restriction that generally will not be violated in practical reservoir simulations. Some numerical results that support the theoretical conclusions are presented. Introduction Early finite-difference, Multiphase reservoir simulators using explicit mobility were found to require exceedingly small time steps to solve certain types of problems, particularly coning and gas percolation. Both these problems are characterized percolation. Both these problems are characterized by regions of high flow velocity. Coats developed an ad hoc technique for dealing with gas percolation, but a more general and highly successful approach for dealing with high-velocity problems has been the use of implicit mobility. Blair and Weinaug developed a simulator using completely implicit mobility that greatly relaxed the time-step restriction. Their simulator involved iterative solution of nonlinear difference equations, which considerably increased the computational work per time step. Three more recent papers introduced the use of semi-implicit mobility, which proved to be greatly superior to the fully implicit method with respect to computational effort, ease of use, and maximum permissible time-step size. As a result, semi-implicit mobility has achieved wide use throughout the industry. However, this success has been pragmatic, with little or no theoretical work to justify its use. In this paper, we attempt to place the use of semi-implicit mobility on a sounder theoretical foundation by examining the stability of semi-implicit difference equations. The usual linearized stability analysis is not delicate enough to distinguish between the semi-implicit and completely implicit difference equation. A nonlinear stability analysis is developed that permits the detection of some differences between the stability of difference equations using implicit mobility and those using semi-implicit mobility. DIFFERENTIAL EQUATIONS The ideas to be developed may be adequately presented using the following simplified system: presented using the following simplified system: horizontal, one-dimensional, two-phase, incompressible flow in homogeneous porous media, with zero capillary pressure. A variable cross-section is included so that a variable flow velocity may be considered. The basic differential equations are (1) (2) The total volumetric flow rate is given by (3) Addition of Eqs. 1 and 2 yields =O SPEJ P. 79


2013 ◽  
Vol 437 ◽  
pp. 716-721 ◽  
Author(s):  
Yue Ming Li ◽  
Ying Hao Zhang ◽  
Guo Cheng Zhang ◽  
Zhong Hui Hu

This paper addresses the stability analysis on S Plane Control in terms of both position and velocity control. Employing Lyapunov stability theory and T-passivity theory, this paper proves the stability of the position controller based on S Plane Control, and on this ground, the stability analysis of the velocity controller based on S Plane Control is done. Finally, the S Plane Control results obtained from the sea trials are given.


2004 ◽  
Author(s):  
A. Khazaei ◽  
M. Rastgaar Aagaah ◽  
M. Mahinfalah ◽  
N. Mahmoudian ◽  
G. Nakhaie Jazar

This paper presents the stability theory and dynamic behavior of a micro-mechanical parametric-effect resonator. The device is a MEMS time-varying capacitor. The nonlinear dynamics of the MEMS are investigated analytically, and numerically. Applying perturbation methods, and deriving an analytical equation to describe the frequency response of the system enables the designer to study the effect of changes in the system parameters that can be used for design and optimization of the system.


2008 ◽  
Vol 609 ◽  
pp. 305-317 ◽  
Author(s):  
FRANÇOIS BLANCHETTE ◽  
THOMAS PEACOCK ◽  
RÉMI COUSIN

We present the results of a combined theoretical and experimental study of the stability of a uniformly stratified fluid bounded by a sidewall moving vertically with constant velocity. This arrangement is perhaps the simplest in which boundary effects can drive instability and, potentially, layering in a stratified fluid. Our investigations reveal that for a given stratification and diffusivity of the stratifying agent, the sidewall boundary-layer flow becomes linearly unstable when the wall velocity exceeds a critical value. The onset of instability is clearly observed in the experiments, and there is good quantitative agreement with some predictions of the linear stability analysis.


1960 ◽  
Vol 38 (2) ◽  
pp. 153-162 ◽  
Author(s):  
A. E. Scheidegger

It is well known that, during the displacement of a fluid contained in a porous medium by another less viscous one, the displacement front may become unstable: Fingers occur which proceed rapidly through the system.The Muskat–Aronofsky model of displacement in porous media, in which it is assumed that a sharp front exists with maximum saturation by the respective fluid being present on either side of the front, is analyzed in the light of the phenomenon of fingering. It is shown that the Muskat–Aronofsky model, in fact, demands that fingering occurs for mobility ratios (displaced/displacing fluid) smaller than one. This model should, therefore, not be used for the calculation of the steady progress of a front for such mobility ratios. The Muskat–Aronofsky model also yields some conditions regarding the geometry of fingers; the latter are deduced. It does not, however, describe the fingering process completely. In this connection, one would have to take recourse to the statistical geometry of porous media. This will be done in a separate paper.


1979 ◽  
Vol 46 (2) ◽  
pp. 259-262 ◽  
Author(s):  
P. C. Parks

In this paper the stability theory of Poincare´ and others, expounded in the last few pages of Sir Horace Lamb’s textbook Hydrodynamics, is extended by a stability analysis using Liapunov functions to cover the full nonlinear equations of motion of spinning liquid-filled spheroids. The linearized criterion is confirmed, that is for stability c/a < 1 or c/a > 3 where the semiaxes of the spheroid are a, a, c, the c-axis being also the spin axis. The unstable motion for 1 < c/a < 3 is examined, and a conjecture that viscosity will destabilize spinning spheroids with c/a > 3 is proposed.


1997 ◽  
Vol 1 (1) ◽  
pp. 1-8 ◽  
Author(s):  
Hermann Haken

This article extends the slaving principle of synergetics to processes with discrete time steps. Starting point is a set of nonlinear difference equations which contain multiplicative noise and which refer to multidimensional state vectors. The system depends on a control parameter. When its value is changed beyond a critical value, an instability of the solution occurs. The stability analysis allows us to divide the system into stable and unstable modes. The original equations can be transformed to a set of difference equations for the unstable and stable modes. The extension of the slaving principle to the time-discrete case then states that all the stable modes can be explicitly expressed by the unstable modes or so-called order-parameters.


Author(s):  
Avinash Vaidheeswaran ◽  
William D. Fullmer ◽  
Krishna Chetty ◽  
Raul G. Marino ◽  
Martin Lopez de Bertodano

The one-dimensional fixed-flux two-fluid model (TFM) is used to analyze the stability of the wavy interface in a slightly inclined pipe geometry. The model is reduced from the complete 1-D TFM, assuming a constant total volumetric flux, which resembles the equations of shallow water theory (SWT). From the point of view of two-phase flow physics, the Kelvin-Helmholtz instability, resulting from the relative motion between the phases, is still preserved after the simplification. Hence, the numerical fixed-flux TFM proves to be an effective tool to analyze local features of two-phase flow, in particular the chaotic behavior of the interface. Experiments on smooth- and wavy-stratified flows with water and gasoline were performed to understand the interface dynamics. The mathematical behavior concerning the well-posedness and stability of the fixed-flux TFM is first addressed using linear stability theory. The findings from the linear stability analysis are also important in developing the eigenvalue based donoring flux-limiter scheme used in the numerical simulations. The stability analysis is extended past the linear theory using nonlinear simulations to estimate the Largest Lyapunov Exponent which confirms the non-linear boundedness of the fixed-flux TFM. Furthermore, the numerical model is shown to be convergent using the power spectra in Fourier space. The nonlinear results are validated with the experimental data. The chaotic behavior of the interface from the numerical predictions is similar to the results from the experiments.


2016 ◽  
Vol 94 (1) ◽  
pp. 378-381 ◽  
Author(s):  
A. V. Akhmetzyanov ◽  
A. G. Kushner ◽  
V. V. Lychagin

Sign in / Sign up

Export Citation Format

Share Document