The Injection of a Hot Liquid into a Porous Medium

1969 ◽  
Vol 9 (01) ◽  
pp. 100-114 ◽  
Author(s):  
J.E. Cheppelear ◽  
C.W. Volek

Abstract The injection of a hot liquid into initially cool porous media, saturated with the same liquid and porous media, saturated with the same liquid and surrounded by two impermeable but beat conducting media (cap and base rock), bas been studied both experimentally and theoretically. The temperature dependence of the viscosity was included in the theoretical model, but it was assumed that the specific heats and densities of the various materials were independent of the temperature. Solutions to the theoretical model were approximated by numerical methods. Both theoretical and experimental results indicate that centerline temperatures are significantly higher than boundary temperatures. Comparison of experimental and theoretical results with a cold/hot viscosity ratio of 19:1 were in reasonable agreement. Theoretical calculations show, that the effect of the temperature dependence of viscosity was very significant at ratios of 100:1 to 1000:1, which are typical of those that occur when injecting hot water to flood heavy oil reservoirs. Introduction We consider the problem of prediction of fluid flow and temperature distribution in an initially cold-fluid-filled reservoir on the injection of the same hot liquid by the use of mathematical and physical models. The results reported are for a physical models. The results reported are for a two-dimensional rectangular section of the reservoir, as shown in Fig. 1. The injection and withdrawal faces are assumed to be equipotentials for fluid flow. The ultimate purpose of such models would be to predict hot-water injection performance. However, we note that in the work presented here, one of the most significant aspects of the problem - the instabilities resulting from two-phase, water and oil flow - is not included. We will not give a historical review, but refer instead to the paper of Spillette and Nielsen, which contains a rather complete bibliography and critical discussion. The physical problem is exactly the same as Spillette and Nielsen, except for certain simplifications in our assumptions. The mathematical details are somewhat different, and we will present the details of our method here. The mass flow equation, which is elliptic in character, is handled by successive overrelaxation. The heat flow difficulty arose in the over-all heat balance due to by a straightforward explicit approximation. Some difficulty arose in the over-all heat balance due to small errors in the solution of the mass flow equation, and we feel that a different formulation of the heat flow equation would be desirable for future work. In addition to the mathematical solutions for the temperature distributions in a porous medium due to the injection of hot liquids, experimental data are presented to check the validity of these solutions. presented to check the validity of these solutions. A schematic diagram of the model is shown in Fig. 2. Certain qualitative physical conclusions were obtained from our numerical and experimental results These are:Assuming high (infinite) conductivity normal to the bedding plane in the reservoir is a poor approximation, and may lead to overestimates of the total heat losses (to cap and base rock) of as much as 50 percent.More heat is retained in the reservoir (per unit of heat injected) for higher viscosity changes.Any particular temperature isotherm moves more rapidly along the centerline for higher viscosity changes. Consequently, the approximation of temperature independent viscosity is not suitable for obtaining quantitatively correct results. MATHEMATICAL MODEL Our mathematical model is that of a reservoir of thickness 2h. The problem is idealized from that of a linear hot-water drive, and we imagine that the input face is sufficiently far away from the injection wells that the stream lines enter it normally. Similarly the production wells are far enough from the outflow face that the flow lines leave it normally. SPEJ P. 100

1999 ◽  
Author(s):  
Mario F. Letelier ◽  
César E. Rosas

Abstract A theoretical study of the fully developed fluid flow through a confined porous medium is presented. The fluid is described by the Bingham plastic model for small values of the yield number. The analysis allows for many admissible shapes of the wall contour. The velocity field is computed for several combination of relevant parameters, i.e., the yield number, Darcy resistance coefficient and the boundary perturbation parameter. The wall effect is especially highlighted and the characteristics of the central plug region as well. Plots of isovel curves and velocity profiles are included for a variety of flow and geometry parameters.


Author(s):  
David Jon Furbish

The concept of conservation of mass holds a fundamental role in most problems in fluid physics. For a given problem this concept is cast in the form of an equation of continuity. Such an equation describes a condition—conservation of mass—that must be satisfied in any formal analysis of a problem. Thus an equation of continuity often is one of several complementary equations that are solved simultaneously to arrive at a solution to a flow problem, for example, the flow velocity as a function of coordinate position in a flow field. (Typically these complementary equations, as we will see in later chapters, involve conservation of momentum or energy, or both.) Although we did not explicitly use this idea in analyzing the one-dimensional flow problems at the end of Chapter 3, it turns out that continuity was implicitly satisfied in setting up each problem. We will return to these problems to illustrate this point. We will develop equations of continuity for three general cases: purely fluid flow, saturated single-phase flow in porous media, and unsaturated flow in porous media. The most general of the three equations is that for unsaturated flow, where pores are partially filled with the fluid phase of interest, such that the degree of saturation with respect to that phase is less than one. We will then show that this equation reduces, in the special case in which the degree of saturation equals one, to a simpler form appropriate for saturated single-phase flow. Then, this equation for saturated flow could be reduced further, in the special case in which the porosity equals one, to a form appropriate for purely fluid flow. For pedagogical reasons, however, we shall reverse this order and consider purely fluid flow first. In addition we will consider conservation of a solid or gas dissolved in a liquid, and take this opportunity to introduce Fick’s law for molecular diffusion. For simplicity we will consider only species that do not react chemically with the liquid, nor with the solid phases of a porous medium. Most of the derivations below are based on the idea of a small control volume of specified dimensions embedded within a fluid or porous medium.


1974 ◽  
Vol 14 (03) ◽  
pp. 271-278 ◽  
Author(s):  
Milos Kojic ◽  
J.B. Cheatham

Introduction A number of problems occur in the fields of drilling and rock mechanics for which consideration must be given to the interaction of fluid flow and rock deformation. Such problems include those of borehole stability, chip removal from under a drill bit, drilling in the presence of a fluid pressure gradient between the drilling fluid and formation fluid, and drilling by use of hydraulic jets. We have recently developed a general theory of the influence of fluid pressure gradients and gravity on the plasticity of porous media. The solution of the problem considered here serves as an example of the application of that theory. The illustrative problem is to determine the load required on a flat problem is to determine the load required on a flat axially symmetric punch for incipient plasticity of the porous medium under the punch when fluid flows through the bottom face of the punch. The rock is assumed to behave as a Coulomb plastic material under the influence of body forces plastic material under the influence of body forces due to fluid pressure gradients and gravity. Numerical methods that have been used by Cox et al. for analyzing axially symmetric plastic deformation in soils with gravity force are applied to the problem considered here. Involved is an iterative process for determining the slip lines. The fluid flow field ‘used for calculating the fluid pressure gradient is based upon the work by Ham pressure gradient is based upon the work by Ham in his study of the potential distribution ahead of the bit in rotary drilling. The effective stresses in the porous rock and the punch force for incipient plasticity are computed in terms of the fluid plasticity are computed in terms of the fluid pressure and the cohesive strength and internal pressure and the cohesive strength and internal friction of the rock. PLASTICITY OF POROUS MEDIA PLASTICITY OF POROUS MEDIA A recently developed general theory of plasticity of porous media under the influence of fluid flow is summarized in this section. The equation of motion for the porous solid for the case of incipient plastic deformation reduces to the following equilibrium equation:(1) where Ts is the partial stress tensor of the solid; Fs is the body force acting on the solid per unit volume of the solid material; P is the interaction force between the solid and the fluid; and is the porosity, which is defined as the ratio of the pore porosity, which is defined as the ratio of the pore volume to the total volume of the solid-fluid mixture. The partial stress tensor Ts can be considered as the effective stress tensor that is used in sod mechanics. With the acceptance of the effective stress principle defined in Ref. 5, the yield function, f, in the following form is satisfied for plastic deformation of the porous medium. plastic deformation of the porous medium.(2) where EP is the plastic strain tensor and K and the work-hardening parameter. From the equation of motion for the fluid, the interaction force P can be expressed in the form(3) where is the inertial force of the fluid per unit volume of the mixture and F is the body force acting on the fluid per unit volume of fluid. For the case of incipient plastic deformation the solid can be considered static (velocities of the solid particles are zero), and the problem of determining particles are zero), and the problem of determining the fluid flow field is the one usually analyzed in petroleum engineering. petroleum engineering. Consider a flow of be fluid such that the inertial forces of the fluid can be neglected and assume that Darcy's law is applicable. SPEJ P. 271


1974 ◽  
Vol 14 (03) ◽  
pp. 263-270 ◽  
Author(s):  
Milos Kojic ◽  
J.B. Cheatham

Abstract Plastic deformation of a porous medium containing moving fluid is analyzed as a motion of a solid-fluid mixture. The fluid is considered to be Newtonian, and the porous material consists of interconnected pore spaces and of solid particles that can deform pore spaces and of solid particles that can deform elastically. The effective stress principle and a general form of the yield function-including work-hardening characteristics-and general stress-strain relations are applied to describe the plastic deformation of the solid. The system of plastic deformation of the solid. The system of governing equations with the number of unknowns being equal to the number of equations is formed. A possible method of solution of a general problem is described. Some simplification such as problem is described. Some simplification such as the assumptions of quasi-static plastic deformation and incipient plastic deformation with the application of Darcy's law for the fluid flow are discussed. To illustrate an application of the theory, the problem of incipient plane plastic deformation of a Coulomb material is presented. Introduction The motion of fluid through a porous medium and the deformation of a porous medium containing fluid have been the subjects of many investigations. For problems concerning fluid flow through porous media in petroleum and civil engineering literature, the porous material is usually considered undeformable and Darcy's law is taken as the governing relation between the velocity and the pressure of the fluid. pressure of the fluid. Most of the effort concerning fluidization of porous media has been experimental; here the task porous media has been experimental; here the task is to find the critical pressure gradient or the critical velocity of the fluid that will cause fluidization. Only the one-dimensional equilibrium equation, which relates Ne pressure gradient of the fluid and densities of solid and fluid, has been analyzed in most fluidization studies. Recently, a more general theoretical approach has been taken and equations of motion of fluid and solid have been established. Some of the results of this theory are used in the present study. Previous investigations of the deformation of porous media containing fluid have been both porous media containing fluid have been both empirical and theoretical. In the domain of elastic deformation much of the published material has dealt with experimental work aimed at finding the relation between a change in fluid pressure and stresses and deformation of the solid phase. A general theory of elasticity of porous media containing moving fluid was established by Biot. However, that theory is approximate since Darcy's law is considered as a governing relation for the fluid, and the change of permeability with the deformation of the solid is neglected. A simplification of this theory was presented by Lubinski. Experimental work has been carried out in the domain of plastic deformation of porous media containing fluid. The effective stress principle has been established as a result of experiments using saturated sand and porous rocks with various pore pressures (fluid is static in these experiments. pressures (fluid is static in these experiments. This principle, which is considered as a fundamental principle in soil mechanics, states that the pore principle in soil mechanics, states that the pore pressure does not affect the yield criterion of the pressure does not affect the yield criterion of the solid. In other words, the yield condition of the solid depends only on stresses transmitted among the solid particles. The influence of fluid flow on plasticity of porous media was indicated by Lambe and Whitman porous media was indicated by Lambe and Whitman in the analysis of stability of an infinite slope of a soil. In the equilibrium equation of a so-called "free body" a term equal to the negative pressure gradient is added. There is no general theory for plasticity of porous media containing moving fluid. plasticity of porous media containing moving fluid. GENERAL THEORY Consider the motion of a solid-fluid mixture and suppose that the motion of the solid is a plastic deformation. Then the problem reduces to the following: define the motion of a solid-fluid mixture so that the yield criterion of the solid is satisfied. The mechanical model can be described as follows. 1. The system comprises one fluid and one should constituent. SPEJ P. 263


2011 ◽  
Vol 8 (1) ◽  
pp. 25-38
Author(s):  
A.T. Akhmetov ◽  
S.V. Lukin ◽  
D.M. Balapanov ◽  
S.F. Urmancheev ◽  
N.M. Gumerov ◽  
...  

There are the results of experimental and theoretical studies on the propagation of weak shock waves in the wet sand at different water saturation. There are mathematical model and numerical analysis of propagation of pressure pulses in porous media, taking into account capillary forces. Non-monotonic dependence of the amplitude of the wave resulting in a wet porous medium vs. the degree of water saturation is installed. The evolution of the fast, slow and filtration waves depends on the saturation of the system with water is analyzed. The influence of capillary forces on the propagation of longitudinal waves is evaluated.


2020 ◽  
Author(s):  
Yuri Perepechko ◽  
Konstantin Sorokin ◽  
Georgiy Vasilyev

<p>The aim of the research is to construct a mathematical model of the formation of a fracture system in magma intrusion in the permeable zones of the lithosphere and on this basis to study the formation of magmatic channels in the lithosphere and crust. The lithosphere substrate is modeled by a saturated porous medium in which the processes of small-scale destruction in the mantle magma intrusion lead to the formation of faults and, consequently, to a magmatic channel. Destruction and occurrence of micro-fracture fields can be associated with both magma flow and external seismic effect leading to the rock breaking. The process of small-scale destruction is described within the framework of the dynamics of the elastoplastic fracture-porous medium and causes variations in the rheological properties of the lithosphere substrate. A feature of this process is the destruction substrate in the compression zone represented by a narrow area with a sharply changing concentration of micro-fractures. The micro-fracture accumulation provides the conversion of the broken area into a macro-fissure. The elastoplastic porous matrix in the destruction zone contains both broken and intact substrate, the relative content of which is determined by relaxation of deformations, the speed of which depends on stress and yield stress point according to the power law. The obtained mathematical model provides investigation of currents in fractured-porous media and their effect on the small-scale destruction. Based on the TVD-Runge Kutta method numerical simulation of the compressible fluid infiltration into the fracture-porous permeable channel has shown that stresses in the compression domain can reach stress limits of breaking and result in fracture formation. Change in relaxation time does not result in a marked change in stress fields. The concentration of maximum stresses is observed in the channel center leading to an increase in its fracture porosity. The computational results show the appearance of high stress values in the compression domain in the process of a liquid phase injection, for instance, magma, into a low-permeable fracture-porous layer. The introduction of the destruction criterion will help to associate the occurrence of such regions to the local breaking of the porous matrix. Thus, the proposed micro-fracture generation mechanism can be used to describe the formation of fracture or channels in micro-fracture porous media. Work is done on state assignment of IGM SB RAS with partial support from the Russian Foundation for Basic Research, grants No. 16-29-15131, 19-05-00788.</p>


Sign in / Sign up

Export Citation Format

Share Document