Hostname: page-component-7c8c6479df-8mjnm Total loading time: 0 Render date: 2024-03-26T13:59:39.667Z Has data issue: false hasContentIssue false

Transient effects in oilfield cementing flows: Qualitative behaviour

Published online by Cambridge University Press:  01 August 2007

M. A. MOYERS-GONZÁLEZ
Affiliation:
Département de mathématiques et de statistique, Université de Montréal, CP 6128 succ. Centre-Ville, Montréal, Quebec, CanadaH3C 3J7 email: moyers@dms.umontreal.ca
I. A. FRIGAARD
Affiliation:
Department of Mathematics and Department of Mechanical Engineering, University of British Columbia, 2054-6250 Applied Science Lane, Vancouver, British Columbia, CanadaV6T 1Z4 email: frigaard@mech.ubc.ca
O. SCHERZER
Affiliation:
Department of Computer Science, Universität Innsbruck, Technikerstraße 25, A-6020 Innsbruck, Austria email: Otmar.Scherzer@uibk.ac.at
T.-P. TSAI
Affiliation:
Department of Mathematics, University of British Columbia, 1984 Mathematics Road, Vancouver, British Columbia, CanadaV6T 1Z2 email: ttsai@math.ubc.ca

Abstract

We present an unsteady Hele–Shaw model of the fluid–fluid displacements that take place during primary cementing of an oil well, focusing on the case where one Herschel–Bulkley fluid displaces another along a long uniform section of the annulus. Such unsteady models consist of an advection equation for a fluid concentration field coupled to a third-order non-linear PDE (Partial differential equation) for the stream function, with a free boundary at the boundary of regions of stagnant fluid. These models, although complex, are necessary for the study of interfacial instability and the effects of flow pulsation, and remain considerably simpler and more efficient than computationally solving three-dimensional Navier–Stokes type models. Using methods from gradient flows, we demonstrate that our unsteady evolution equation for the stream function has a unique solution. The solution is continuous with respect to variations in the model physical data and will decay exponentially to a steady-state distribution if the data do not change with time. In the event that density differences between the fluids are small and that the fluids have a yield stress, then if the flow rate is decreased suddenly to zero, the stream function (hence velocity) decays to zero in a finite time. We verify these decay properties, using a numerical solution. We then use the numerical solution to study the effects of pulsating the flow rate on a typical displacement.

Type
Papers
Copyright
Copyright © Cambridge University Press 2007

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Andreu, F., Ballester, C., Caselles, V. & Mazón, J. M. (2001) The Dirichlet problem for the total variation flow. J. Funct. Analy. 180, 347403.CrossRefGoogle Scholar
[2]Barenblatt, G. I., Entov, V. M. & Ryzhik, V. M. (1990) Theory of Fluid Flows Through Natural Rocks. Kluwer Academic Publishers, New York.CrossRefGoogle Scholar
[3]Bittleston, S. H.Ferguson, J. & Frigaard, I. A. (2002) Mud removal and cement placement during primary cementing of an oil well; laminar non-Newtonian displacements in an eccentric annular Hele–Shaw cell. J.Eng.Math. 43, 229253.CrossRefGoogle Scholar
[4]Dacorogna, B. (1989) Direct Methods in the Calculus of Variations. AMS Vol. Springer-Verlag.CrossRefGoogle Scholar
[5]Dustrehoft, D., Wilson, G. & Newman, K. (2002) Field Study on the Use of Cement Pulsation to Control Gas Migration. Society of Petroleum Engineers paper number SPE 75689.Google Scholar
[6]Ekeland, I. & Témam, R. (1999) Convex Analysis and Variational Problems. SIAM.CrossRefGoogle Scholar
[7]Evans, L. C. (1998) Partial Differential Equations. AmericanMathematicalSociety.Google Scholar
[8]Fortin, M. & Glowinski, R. (1983) Augmented Lagrangian Methods. North-Holland.Google Scholar
[9]Galdi, P. (1991) An Introduction to the Mathematical Theory of the Navier Stokes Equations, Part 1. Springer-Verlag, 463 pp.Google Scholar
[10]Gilbarg, D. & Trudinger, N. S. (1983) Elliptic Partial Differential Equations of Second Order. Springer-Verlag, New York.Google Scholar
[11]Glowinski, R. (1983) Augmented Lagrangian Methods: Applications to the Numerical Solution of Boundary-Value Problems. North-Holland.Google Scholar
[12]He, J. W. & Glowinski, R. (2000) Steady Bingham fluid flow in cylindrical pipes: a time dependent approach to the iterative solution. Numer. Linear Algebra Appl. 7 (6), 381428.3.0.CO;2-W>CrossRefGoogle Scholar
[13]Moyers-González, M. A. (2006) Transient Effects in Oilfield Cementing flows. Ph.D. thesis, University of British Columbia.Google Scholar
[14]Moyers-González, M. A. & Frigaard, I. A. (2006) Kinematic instabilities in two-layer eccentric annular flows, part1: Newtonian fluids. J.Eng.Math., October, accepted for publication.Google Scholar
[15]Pelipenko, S. & Frigaard, I. A. (2004) On steady state displacements in primary cementing of an oil well. J.Eng.Math. 46 (1), 126.CrossRefGoogle Scholar
[16]Pelipenko, S. & Frigaard, I. A. (2004) Two-dimensional computational simulation of eccentric annular cementing displacements. IMA J. Appl. Math. 69, 557583.CrossRefGoogle Scholar
[17]Pelipenko, S. & Frigaard, I. A. (2004) Visco-plastic fluid displacements innear-vertical eccentric annuli: lubrication modelling. J.FluidMech. 520, 343377.CrossRefGoogle Scholar
[18]Tehrani, A., Bittleston, S. H. & Long, P. J. (1993) Flow instabilities during annular displacement of one non-Newtonian fluid by another. Exp. Fluids 14, 246256.CrossRefGoogle Scholar
[19]Zalesak, S. T. (2001) Fully multidimensional flux-corrected transport algorithms for fluids. J. Non-Newt. Fluid Mech. 100, 4976.Google Scholar