Search Documents

By G. L. Chierici

Starting from Muskat's theory of water and gas coning, maximum permissible oil production rates without water and/or free-gas production have been determined, in a broad range of reservoir and well parameters, using the potentiometric model technique. The main assumptions made are as follows: (1) the reservoir rock is homogeneeous (either isotropic or anisotropic); (2) the volume of the aquifer underlying the oil zone is very small, so that it does not contribute to reservoir energy; and (3) the gas cap expands at a very low rate, so that it can be assumed to be in quasi-static conditions. The results obtained are presented in the form of diagrams which can be used for solving two types of problems: (1) given the reservoir and fluids characteristics, as well as the position and length of the perforated interval, determine the maximum oil production rate without water and/or free-gas production; and (2) given the reservoir and fluids characteristics only, determine the position and length of the perforated interval which optimize the maximum permissible oil production rate, without water and/or free-gas production. INTRODUCTION In oil reservoirs where the oil-bearing formation is underlain by an aquifer which does not participate in the production mechanism, water-coning is a limiting factor to the flow rates of producing wells. Production rates are usually kept to a value that will prevent the water from entering the wells. The entry of water into a well lowers its productivity by increasing the weigbt of the fluid column; moreover, the separation of water from the effluent, at the surface, may constitute a very difficult problem in cases of heavy viscous oils. A similar situation is encountered in oil reservoirs with a gas cap overlying the oil-saturated zone; here a downward gas cone is induced by the flow of oil towards the producing wells. Production rates must be low enough to prevent the gas from being produced; producing gas from the gas cap would be a waste of energy. Of course, water-coning and gas-coning phenomena can occur at the same time in the same reservoir if the oil-producing formation is both overlain by a gas zone and underlain by a water zone. Due to its relevant practical importance, the mechanism of coning was studied by many people.2,3,5-8 Defining the conditions for getting the maximum water-free and/or gas-free oil production rate is a difficult problem, often encountered under one of the following aspects: 1. Predict the maximum flow rate that can be assigned to a completed well without the simultaneous production of water and/or free-gas. 2. Define the optimum length and position of the interval to be perforated in a well, in order to obtain the maximum water and gas-free production rate. A systematic study of these problems was made by means of the electrical analog technique. The results of this study are presented here, under the form of a set of curves providing solutions for the above stated problems. These curves are valid only for homogeneous forrnations, either isotropic or anisotropic. Should the formation be non-homogeneous (by horizontal or vertical variation of permeability, shale diaphragms, fractures, etc.), a specific potentiometric study would be required for each specific case. Especially when shale diaphragms of some radial extension are present, the critical rates observed are much larger than would be expected from the diagrams. STATEMENT OF THE PROBLEM In the present study the aquifer is supposed to be of such limited volume that it does not contribute to the energy of the reservoir. Moreover, the gas cap is supposed to expand at such a low rate that the potential gradient in the gas cap is negligible. Under static conditions water-oil and gas-oil interfaces (T1 and T2) are both horizontal. When the reservoir production starts, below each well these interfaces take a cone-like shape (Fig. 1) having as an axis the axis of the well. This shape results from the equilibrium between potential gradients in the oil zone and gravitational forces due to density differences between oil and water and between oil and gas. Assuming the oil-bearing formation to be homogeneous and the oil to be incompressible, the analysis of the problem (see Appendix) shows that the oil-water and gas-oil interfaces are stable only if the oil production rate of the well is not higher than the following values.

Jan 1, 1965

By R. W. Heins, N. Street

An investigation was made of the Rehbinder-Kuznetsov pendulum as an instrument to aid in the study of rock drillability and of the effect of surface active agents on rock drillability. The "critical weight", i.e., the pendulum weight at greatest sensitivity, was determined for several solids and it was shown that a change in surface environment changed the maximum value of the apparent hardness but not the critical weight, and the critical weight increased with increase of hardness of the solid used. Using silver halides as the solid material, and solutions of either silver nitrate or the appropriate alkali halide, it was shown that there is a hardness maximum at the zero point of charge of the particular silver halide. INTRODUCTION The development of relative hardness measurement by the method of damped oscillations, or what has been called the pendulum sclerometer, is due largely to the Russian workers P. A. Rehbinder and V. D. Kuznetsov and their co-workers.1,2 The apparent simplicity of the method has led to its use in a number of fields, particularly those where relative hardness has to be determined on a brittle solid. As will be shown, the simplicity of the method belies its complexities. THEORY OF THE PENDULUM SCLEROMETER The operation of the pendulum is quite simple. Make a sharp point such as a diamond, hardened steel or tungsten carbide, as the fulcrum of a pendulum, with the center of gravity of the system below the point of support. The attenuation of this pendulum is a function only of air friction, provided no fragmentation, breakage or other work is performed by the point. However, when the solid surface breaks or disintegrates under the point, and the latter penetrates the solid, part of the energy of oscillation is absorbed in the work of creating new surface, and the rate of attenuation increases. In the case of brittle solids, the point of the pendulum produces fragmentary particles, whereas the point produces plastic deformation in a plastic solid; but in both cases damping occurs. This pendulum damping is a function of the hardness or strength of the solid. Whereas a pendulum oscillating on a hard material attenuates primarily by air friction, in the case of a soft material, the points penetrate the surface of the solid, producing new surface, and the pendulum damps much faster. Thus the attenuation rate is a function of the hardness or strength of a material in both plastic and brittle substances. Russian workers initially used a single support for the pendulum, but this created many problems because of oscillation of the pendulum in more than one plane. A revised version used two points, or sometimes a prism as support, successfully confining oscillation to one plane. The version using two points was constructed and used throughout this investigation, making the recording of the attenuation as a function of time relatively easy. That both points should be in contact with the solid during the damping is implied. Rehbinder (c.f. Kuznetsov2) proposed that the hardness H as measured by this method could be expressed by a quantity which is inversely proportional to the relative initial rate of attenuation of amplitude. where H = hardness parameter (arbitrary units) a, = initial amplitude da/dt = rate of attenuation of amplitude with time. In the simplest case of damping with a constant logarithmic decrement, H can be regarded as a constant for the duration of the whole damped oscillation process. That is

Jan 1, 1966

By Caudle B. H., Cooper H. E. Jr., O. K. Kimbler

The results of studies of areal sweepout behavior in nine-spot injection patterns are presented in this paper. The ratio of the production rates of the two types of producing wells in a nine-spot pattern is considered in its effect on breakthrough. A model study was used to develop a method to predict the production histories of both side and corner wells and to develop the production history of the pattern or the field. INTRODUCTION The importance of areal sweepout behavior in injection operations has been recognized for a number of years. We now know that the fraction of a reservoir which will be reached by an injected fluid is a function of the geometric arrangement of the production and injection wells and of the mobility ratio of the fluids. Of these two, the mobility ratio (which is a measure of the relative ease with which the oil and the injected fluid can move through the reservoir) appears to have the greatest effect on recovery. The effect of mobility ratio on production history has been reported for most of the dispersed injection patterns (i. e., five-spot, staggered and straight-line drives).4-6 ere has been very little published on the nine-spot injection pattern. This may be because the two types of producing wells in the nine-spot pattern can—and sometimes should —be produced at different rates. This adds a second variable (the fist being mobility ratio) to be covered in a study of the nine-spot flooding network. Muskat1 has reported the sweep pattern efficiency at breakthrough for a mobility ratio of one at various ratios of producing rates for the two types of wells. For these conditions he found maximum oil production at breakthrough to occur when the corner well (the producing well farthest from the injection well) is produced 10.7 times as fast as the side well. Other investigatorsa2,3 have studied a few of the many possible operating conditions possible with this injection pattern. This paper presents the results of studies which are strictly applicable to a nine-spot injection pattern surrounded by other nine-spots in a homogeneous reservoir and which has no flowable gas present (no oil bank is formed). Although the full ranges of the variables mentioned above have not been covered, the results reported here should allow reasonable estimates of field operations to be made. EXPERIMENTAL The well spacing geometry of the nine-spot pattern is shown in Fig. 1. Here we see that each injection well is surrounded by eight production wells arranged in a square. Four of these producing wells are located at the corners of the square; the other four are in the midpoints of the four sides. These two types of producing wells are called "comer wells" and "side wells" respectively. It is the ratio of the producing rates of these two types of wells that affects the breakthrough recovery, as calculated by Muskat.' This ratio is defined as the total producing rate of a corner well divided by the total producing rate of a side well. The model used in this study represented one-eighth of a nine-spot pattern, as shown in Fig. 1. It contained the equivalent of one-fourth of a side well, one-eighth of a comer well, and one-eighth of an injection well. Because of these proportions, the actual producing rate ratios in the model were one-half of the ratios of producing rates as

Jan 1, 1965

By C. D. Stahl, R. Sandrea

This investigation was made to evaluate the comparative effects of the viscosity and the phase behavior of the buffer fluid in the composite solvent displacement of Bradford crude from waterflooded sandstone cores. Buffer slugs of propane, naphtha and other refined hydrocarbons exbihiting increasingly favorable viscosity ratios and decreasing solubility relationships, respectively, with Bradford crude, were used in long Berea sandstone cores. The secondary slug was isopropyl alcohol. The results indicate that higher oil recoveries are obtained for increasingly favorable phase relationships even when these are accompanied by unfavorable viscosity ratios within the range studied. Moreover, when propane is used as a buffer slug with an adverse viscosity ratio of 36, it gives higher oil recoveries than at a similar size slug of an any1 alcohol having a viscosity ratio of 1.09. The investigation was extended to the study of the effect of flooding rate on oil recovery. Residual Bradford crude was displaced from a 6-ft Berea sandstone core at rates varying from 0.3 to 30 ft/ day. The results show that as flooding rates were increased above or decreased below a minimum range of 1 to 2 ft /day, displacement efficiency increased considerably. INTRODUCTION In recent years, the alcohol slug process as a means of tertiary recovery has been the subject of many investigations. Gatlin and Slobod showed the advantages of the combination solvent slugs over the previous single-slug process.1 The combination slugs require less total slug material to obtain similar oil recoveries. The basic combination process involves injecting a primary or buffer solvent which is preferentially oleophilic, followed in turn by a water-driven slug of isopropyl alcohol. This procedure ensures that the buffer slug, which is miscible only with the reservoir oil, displaces and replaces the oil and is in turn displaced, together with the water, by the isopropyl alcohol slug miscible with both fluids. Likewise, for those reservoir oils that show some degree of phase incompatibility with isopropyl alcohol, as does Bradford crude, a proper buffer solvent can artificially make the alcohol process technically feasible. From an economic standpoint, propane or LPG would be the most appropriate primary solvent to be used in the field; however, Taber, et a12 have indicated that the alcohol process becomes highly inefficient in the displacement of light hydrocarbons because of the entrapment of these hydrocarbons during the stabilized bank formation. Moreover, since materials such as propane have such low viscosities, adverse viscosity ratios with most reservoir oils would be penerated, thereby possibly further reducing the efficiency of the buffer solvent. This paper presents the results of an experimental study of the comparative effects of viscosity ratio and phase relationships on the displacement efficiency of Bradford crude by means of dual solvents. The work was extended to include the effect of flooding rate on oil recovery. Rates within the range of field and laboratory values were used. EXPERIMENTAL PROCEDURE APPARATUS The porous media used in this study consisted of two 10-ft and one 6-ft Berea sandstone cores, 3-in. in diameter. These were cast in galvanized steel pipe of 31/2-in. ID and the annulus was filled with Armstrong cement. The 10-ft cores were assembled by connecting two 5-ft units, which were joined together by 8-in. steel flanges fitted to the pipe. Intimate contact between the two sections was obtained by inserting a thin disk of the same sandstone between the core faces. The outlet and inlet ends of the cores were covered with threaded pipe caps. The injected fluid was gravity fed to two variable-rate positive-displacement pumps. The fluids were filtered before entering the core and the effluent was collected in an automatic fraction collector and

Jan 1, 1966

By B. D. Smith, K. H. Coats

Experiments in which calcium chloride displaced sodium chloride from four cores showed the extent of asymmetry in the resulting effluent concentration profiles. These results provided a check on how validly the mixing process is modeled by a differential (i.e., not finite-stage) capacitance mathematical model. The effluent concentration profile from two consolidated cores exhibited considerable asymmetry, while two unconsolidated cores yielded nearly symmetrical profiles. All runs resulted in breakthrough of the 50 per cent concentration significantly before one pore volume was injected. In addition, velocity appreciably affected the effluent concentration profile from a Torpedo sandstone core. The differential capacitance model matched the data significantly better than a simple diffusion model. The capacitance model allows determination of the amount of dead-end pore space in a porous matrix and the effect of velocity on the rate of diffusion into this space. An experimental program yielding insight into the physical validity of the capacitance effect is described. INTRODUCTION Axial dispersion — the mixing accompanying the flow of miscible fluids through porous media — has been the subject of many relativeiy recent studiesl-13 and a comprehensive review of the topic has been given by Perkins and Johnston.14 This dispersion is of practical interest in studies of the miscible displacement process, fixed-bed chemical reactors, and the adsorption of solutes from a flowing stream onto the surface of a porous medium. In the latter case, the effect of dispersion must be considered when adsorption parameters are determined from the nature of concentration profiles. In general, early studies of dispersion assumed applicability of a simple diffusion equation and were concerned with correlation of the experimentallv determined "effective" diffusion coefficient with system properties over a large range of the latter. Recent have been concerned with the deviations between the asymmetrical effluent concentration profiles observed and the symmetrical ones predicted by the diffusion model. In the present study, effluent concentration profiles were obtained from consolidated and unconsolidated cores. These profiles were compared with those predicted by a differential (i.e., not finite-stage) capacitance model. Solutions to the simple diffusion model, for three sets of boundary conditions, were compared with one another and with the experimental profiles. SUMMARY OF PREVIOUS WORK The reader is referred to Perkins and johnston14 for an extensive review of studies of dispersion in porous media. Many investigators have employed the simple diffusion model characterized by Eq. 1 below: The dispersion coefficient D for unconsolidated systems is correlated by u is interstitial velocity and dp is particle diameter. Since heterogeneity of the sand pack affects the mixing, this equation is also expressed as where a is proportional to the degree of heterogeneity and is about 3.5 for random packs of unconsolidated sand.l4 Eq. 3 holds for 2 For a Do homogeneous (regular) type of packing, should be 1 or less.6 Aris and Amundson6 and Carberry and Bretton consider o to be the number of particle lengths per mixing cell in the finite-stage model. Data from consolidated cores indicate dp to be about 0.36 cm for outcrop rocks, Torpedo sandstone having a reported value of 0.17 cm.15 For v /Do less than about 2, the effect of molecular diffusion must be considered and the total dispersion coeffi-

Jan 1, 1965

By J. E. Warren

In their analysis of Eqs. 29 and 30, the authors correctly deduce the following behavior limits for the "differential capacitance model". 1. When the rate group a is sufficiently large, the model reduces to that implied by Eq. 1, i.e., the simple convection-diffusion model. 2. When the rate group a is sufficiently small, the model reduces to that implied by Eq. 1 with 1 replaced by 1, i.e., the simple convection-diffusion model with a total pore volume equal to that occupied by the mobile fluid in the original model. Unfortunately, the authors' solution, Eq. C-1, to the transformed versions of Eqs. 29 and 30 does not exhibit the same limiting behavior. For example, as a approaches zero, Eq. C-1 reduces to the solution to Eq. 1 or, equivalently, the solution to Eqs. 29 and 30 for sufficiently large values of a. Therefore, it can be concluded that Eq. C-1 is not a general solution to Eqs. 29 and 30; this conclusion can be easily verified by substitution. I suggest that the correct form of Eq. C-1 is the following: It is generally realized that accounting for dispersion effects is an extremely tricky business; therefore no one should criticize the pragmatic approach of Coats and Smith who proceeded by accessible theory, made use of some available phenomenology, and who then took the trouble to check the consequences, against their own carefully undertaken experimental results. Here I should like to defend three assertions: (1) Coats and Smith copy the error of previous authors in the formulation of the so-called simple diffusion model (which error is then retained in their differential capacitance model); (2) improving the fit between data and the predictions of an arbitrary model by increasing the number of parameters is a triviality unless the parameters are disposable by being independently determined; and (3) existing experimental procedures lack sufficient resolution to clarify some current theoretical uncertainties. Taking first the view that all theoretical approaches (whether constructed by statistical, finite difference, or continuum arguments) tend to reduce to the same model form as fixed by the number and nature of the disposable parameters, one notices that the "finite-stage" model of Deans1 is equivalent to the Coats and Smith so-called differential capacitance model for the pertinent case where a

Jan 1, 1965

By R. J. Blackwell, A. L. Pozzi

Scaled laboratory-model studies provide a powerful method for evaluation of a proposed oil-recovery process. In recent years, models have been used extensively to evaluate processes in which solvents displace oil, both for general cases and for specific reservoir conditions. Since the performance of a miscible flood in a horizontal reservoir can be significantly affected by transverse mixing between solvent and oil, this displacement mechanism must be accurately simulated in the scaled model studies. Unfortunately, precise scaling of transverse dispersion coupled with the requirement of geometric similarity requires impractically large laboratory models and long times for experiments. If scaling requirements for miscible displacements could be relaxed while accurate simulation of essential displacement mechanisms is maintained, the utility of model studies would be greatly enhanced. The purpose of the work reported herein was to evaluate the relative importance of various mechanisms affecting miscible displacement and to ascer tain whether the essential features of the displacement process can be simulated even though some scaling groups are not satisfied. These studies were performed with completely miscible systems in linear, horizontal models packed with unconsolidated media. From the experimental results, a set of relaxed scaling criteria was formulated which allows the requirements of geometric similarity and equality of the ratio of viscous to gravity forces to be omitted for specified conditions. The relaxed criteria are valid whether transverse mixing is by molecular diffusion or by convective dispersion. Correlations which permit prediction of vertical sweep efficiencies in linear, horizontal reservoirs were developed from the experimental data when transverse mixing is by molecular diffusion. These same correlations may be used when transverse mixing is by convective dispersion if an empirically defined, effective, transverse dispersion coefficient is used in the description of the mixing process. The effective transverse dispersion coefficient correlation essentially duplicates the dispersion coefficient correlation for equal-viscosity, equal-density fluid systems. Experimental values for the effective transverse dispersion coefficient can be measured readily. INTRODUCTION One of the most effective methods for evaluation of miscible-displacement oil-recovery processes is that of displacements in laboratory models scaled to simulate reservoir conditions. For these laboratory studies to be meaningful, however, the essential displacement mechanisms affecting reservoir performance must be accurately simulated. Since the performance of a miscible flood in horizontal reservoirs, or in dipping reservoirs at high rates, can be significantly affected by transverse mixing of solvent and oil, this mechanism must be considered in the design of laboratory experiments. Unfortunately, precise scaling of transverse dispersion coupled with the requirement of geometric similarity requires impractically large laboratory models and long experiment times. This difficulty seriously limits the utility of laboratory model studies. Craig,et al,1 demonstrated that geometric similarity is not required when mixing is unimportant. Their experimental data indicate, for the cases studied, that the displacement is sufficiently characterized by scaling the ratio of viscous-to-gravitational forces. This work suggests that relaxation of the requirement of geometric similarity and, possibly, other criteria might also be permissible when mixing is important, provided suitable groups describing the mixing process are scaled. The purpose of the work reported here was to evaluate the relative importance of various mechanisms affecting miscible displacement and to ascertain whether the essential features of the displacement process can be simulated even though some scaling groups are not satisfied. Consequently, laboratory investigations were made of the dependence of recovery on dispersion and other variables,

By W. B. Gogarty, H. R. Bailey

Methods are presented for determining molecular diffusion coefficients by using data from, capillary flow experiments. These methods are based on a numerical solution (presented in a previous paper) of the partial differential equation describing the combined mechanisms of flow and diffusion. Results from this numerical solution are given and compared with the approximate analytical solution of G. I. Taylor. The numerical solution is valid over a much larger time range. These methods are applied to experimental results for the fluid pairs water-potassium permanganate solution and amyl acetate-orthoxylene. Both of these fluid pairs have approximately equal densities and viscosities. Graphical and numerical techniques are presented for deter mining diffusion coefficients from the flow data. Values obtained by these techniques are compared with values obtained by other methods. INTRODUCTION The molecular diffusion coefficient is known to be a variable in determining the amount of mixing in a miscible displacement process. The effect of molecular diffusion on dispersion in longitudinal flow through porous media -has been examined by different investigators. These investigators concluded that at low velocities of flow, the amount of dispersion is approximate1y proportional to the molecular diffusion coefficient. The influence of diffusion on fingering, channeling, and overriding has been mentioned by other investigators.4 Recent studies5,6 have been made on the effects of molecular diffusion in connection with the problem of gravity segregation. Many different methods have been developed for the experimental determination of molecular diffusion coefficients? These methods differ mainly according to boundary conditions selected and analytical procedures used. Nevertheless, all of these methods have the condition in common that the bulk fluids in which diffusion is occurring are stationary with respect to each other. In connection with a series of papers on mixing in capillary flow,8-10 Taylor suggested the use of a flow method for determining molecular diffusion coefficients. Additional studies 11-13 have been conducted on miscible displacements in capillary tubes, but the data from these studies were not used for the specific purpose of determining diffusion coefficients. The flow method proposed by Taylor results in a single value of the diffusion coefficient for the fluid pair used in the displacement experiments. This single value represents the true value for the fluid pair when the diffusion coefficient is independent of concentration. If the diffusion coefficient is a function of concentration, the single value obtained by the flow method gives an average value for the coefficient of the fluid pair. These average values are based on diffusion taking place over the entire range of concentration, i.e., from 0 per cent of one fluid to 100 per cent of that same fluid. In field applications of the miscible displacement process, gradients occur over the same range of concentration as are found in the displacements in capillary tubes. Molecular diffusion coefficients obtained from the capillary flow method should, therefore, be especially relevant to field operations. This investigation was undertaken to evaluate the feasibility of obtaining molecular diffusion coefficients from capillary -flow experiments. In making this evaluation, diffusion coefficients were first determined for two systems from data obtained in capillary flow experiments. These values of the diffusion coefficient were then compared to values obtained by other methods. MIXING IN CAPILLARY FLOW—THEORETICAL The theoretical basis for determining molecular diffusion coefficients from capillary flow experiments is the partial differential equation relating the mechanisms of flow and diffusion. For the case of laminar flow and with a fluid pair having equal densities and viscosities, the equation is given by

By C. R. Johnson, R. A. Greenkorn, L. K. Shallenbarger

This paper describes a study, based on core data, of the directional permeability of a sandstone reservoir. Directional air permeabilities are explained and correlated with lithology by the tensor theory of permeability, which is extended to the more general case of heterogeneous anisotzopic porous media. In this case, the permeability tensor is made up of two components: (I) anisotropy (variation around a point) which correlates with bedding, and (21 heterogeneity (variation from point-to-point) which correlates with grain size. INTRODUCTION Previous studies of directional permeabilityl.5-8 have concerned themselves only with anisotropy at a point, rather than with both point-to-point differences and local anisotropies. In this paper, we present the existing tensor theory of permeability for anisotropic porous media and then extend it to the more general case of heterogeneous anisotropic porous media. The laboratory measurements of directional permeability are explained in view of this extended theory and the data are discussed in terms of the lithologic factors that correlate with it. The results show what the heterogeneity and anisotropy of the reservoir element are, and that the directional permeability correlates with lithology. Data used in this paper are measurements of air permeability in eight directions, spaced at 45° intervals, on 142 2-in. vertical plugs from 30 cores. The core data are meaningful in terms of directional permeability because the cores were oriented to within 45O during drilling and coring.3 Although we initially thought that only 60 per cent of the core material was reliable, subsequent study showed that almost all of it was reliably oriented. After determining the air permeabilities, the data were reduced to three independent variables for each plug: the major and minor permeability axes, and the direction of the major axis. These were obtained by converting the permeability data to the reciprocal square root of permeability and fitting the transformed data with ellipses according to the tensor theory of permeability.8 The point-to-point areal variation of the minimum permeability axis is related to grain size. The direction of the permeability axes, where the major axis exceeds the minor axis by at least 5 per cent (measured variation is about 4 per cent), correlates with the bedding. The permeability tensor used in this study must be considered as the sum of a scalar and a tensor, with the scalar being the minor axis permeability as a function of position, and the tensor the directional effect, which is additive permeability over the minor axis. In this case, the point-to-point variation or heterogeneity (minor axis) is substantially larger than the variation at a point or anisotropy. It may be important that this separation of "directional permeabilities" be recognized when considering migration of fluids due to permeability variation. Local migration may be due to anisotropy and point-to-point migration may be due to heterogeneity, but the direction and magnitude of these may not be the same. Furthermore, in truly heterogeneous systems, one would expect that anisotropy would be the smaller of the two effects. THEORY TENSOR THEORY OF PERMEABILITY FOR ANISOTROPIC POROUS MEDIA8 Darcy's law for flow in porous media in its usual form is where qis the flow rate vector, k is the permeability, is the viscosity, is the vector differential operator, and P is the Pressure. Consider flow, q, through a unit area in the direction of the unit vector n from elementary flow tubes parallel

Jan 1, 1965

By R. W. Parsons

This discussion concerns fitting an ellipse to angle, where the permeability k is measured for a series of plugs cut from varying angles in the horizontal plane. The conclusions and procedures given by Greenkorn, et al. are in no way affected by these comments, since the degree of anisotropy (ratio of maximum to minimum permeabilities) they found is generally quite close to one (see their Table A-1). However, since large anisotropies may arise, I feel that the comments should be made. Both Marcusl and Scheidegger2 have shown that there are two limiting expressions for the permeability of an unbounded anisotropic media. If the permeability is measured in the direction of the pressure gradient, the result is equivalent to the author's Eq. 11. This is the equation used to find kl1 and k22 from the experimental data. If, however, the permeability is measured in the direction of flow, then the resulting equation for the measured permeability is From these equations it follows that plots on angle will be ellipses. The principal axes will correspond with the major and minor permeability axes. The question arises as to what permeability is actually being measured on an anisotropic sample in the laboratory. An overall pressure gradient and a gross flow rate are determined and a permeability calculated by Darcy's law. The details of the macroscopic flow within the sample will be determined by the boundary conditions of the sample and the orientation and magnitude of the major and minor permeability axes. The situation is shown in Fig. 1 for the two-dimensional case. Using the method of conformal mapping, Marcusl has solved the problem of determining the overall apparent (or measured) permeability for this experiment for any condition of W-L ratios and orientation and magnitude of the principle permeability axes. For W-L ratios approaching zero, the streamlines are all forced by the boundary conditions to be in the direction of the measured permeability which is equal to kg. For W-L ratios approaching infinity, the effects of the no-flow bounds are negligible, and the pressure gradient and measured permeatility direction are parall;l. The streamlines are

Jan 1, 1965

By M. A. Torcaso, P. Raimondi

The distribution of the oil phase in Berea sandstone resulting from increasing and decreasing the water saturation by imbibition was investigated Three types of distribution were recognized: trapped, normal and lagging. The amount of oil in each of these distributions was determined as a function of saturation by carrying out a miscible displacement in the oil phase under steady-state conditions of saturation. These conditions were maintained by flowing water and oil simultaneously in given ratios and by using a displacing solvent having essentially the same density and viscosity as the oil. A correlation shows the amount of trapped oil at any saturation to be directly proportional to the conventional residual oil saturation Sor The factor of proportionality is related to the fractional permeability to the water phase. Part of the oil which was not trapped was displaced in a piston-like manner (normal part) and part was eluted gradually (lagging part). The observed phenomena are more than of mere academic importance. Oil which is trapped may well provide the fuel essential for forward combustion and thus be beneficial. On the contrary, in tertiary recovery operations, it is this trapped oil which seems to make current techniques uneconomic. INTRODUCTION A typical oilfield may initially contain connate water and oil. After a period of primary production water often enters the field either from surrounding aquifers or from surface injection. During primary production evolution and establishment of a free gas saturation usually occurs. The effect and importance of this third phase is fully recognized. However, this investigation is limited to a two-phase system, one wetting phase (water) and one non-wetting phase (oil). The increase in water content of a water-wet system is termed imbibition. In a relative permeability-saturation diagram such as the one shown in Fig. 1, the initial conditions of the field would be represented by a point below a water saturation of about 35 per cent, i.e., where the imbibition and the drainage curves to the non-wetting phase nearly coincide. When water enters the field the relative permeability to oil decreases along the imbibition curve. At watered-out conditions the relative permeability to the oil becomes zero. At this point a considerable amount of oil. called residual oil, (about 35 per cent in Fig. 1) remains unrecovered. Any attempt to produce this oil will require that its saturation be increased. In Fig. 1 this would mean retracing the imbibition curve upwards. In addition, processes like alcohol and fire flooding, which can be employed at any stage of production, involve the complete displacement of connate water and an increase, or imbibition. of water saturation ahead of the displacing front. Thus, in several types of oil production it is the imbibition-relative permeability curve which rules the flow behavior. For this reason a knowledge of the distribution of the non-wetting phase, as obtained through imbibition, whether "coming down" or "going up" on the imbibition curve, is important.

Jan 1, 1965

By A. S. Michaels, M. C. Porter, A. Stancell

Previous laboratory studies have demonstrated that the injection of small quantities of reverse wetting agents during water displacement can increase oil recovery from unconsolidated porous media.9 In the present investigation, an attempt has been made to determine more fully the effects of reverse wetting treatments and to clarify the mechanism by which increased oil recovery is effected. Water-oil displacements were performed in beds of 140-200 mesh silica sand. Hexylamine slugs (injected after 0.25 pore volume of water throughput), when adequate in size and concentration, were effective in promoting additional oil recovery. Their effectiveness increased with the quantity of amine injected. However, slugs of sufficient size and concentration to stimulate oil production at water flow rates of 34 ft/day did not do so at 4 ft/day. Visual studies in a glass grid micromodel have shown that the stimulation of oil production, via ayueous hexylamine, is a result of transient changes in the oil wettability of the pore walls. if the amine slug is of sufficient size and concentration to induce significant changes in the adhesion-tension, large continuous oil masses will be formed. if the superficial water velocity is high enough to result in rapid desorption of the amine, a favorable "wettobility gradient" may be established across the masses; under such conditions, high oil mobility is observed, and increased oil recovery results. INTRODUCTION It is generally agreed that the efficiency of oil displacement by water in porous media is limited in part by capillary forces which cause the retention of isolated masses of oil — resulting in the so-called "irreducible minimum oil saturation". Recent estimates10 indicate that there are about 220 billion bbl of petroleum in United States reservoirs which are not economically recoverable with present techniques (such as water flooding). This amounts to almost five times the known recoverable reserves. It has been recognized for some time that a suitable alteration in the water-oil interfacial tension and/or the contact angle, as measured between the water-oil interface and the solid surface, should result in better displacement efficiency. Surface active agents can be used as interfacial tension depressants to accomplish this objective,1,2,7 but unfortunately, the additional oil recovery is seldom commensurate with the treatment cost. In contrast to interfacial tension depressants, the effect of contact angle alterations on water-oil displacements has received relatively little attention in the literature. It is known that the wettability affects the displacement process.3,4,11 Displacements in water-wet systems generally result in lower residual oil saturations than those in oil-wet systems. The effect of "'transient" wettability alterations concurrent with the displacement process have been investigated by Wagner, Leach and coworkers,4,11 wherein it has been demonstrated that the establishment of water-wet conditions during water flooding of oil-wet, oil-saturated porous media is accompanied by significant increase in oil displacement efficiency. Michaels and Timminss studied the effects of transient contact angle alterations resulting from chromatographic transport of reverse wetting agents through unconsolidated sand. It was demonstrated that chromatographic transport of short-chain (C4 through C8) primary aliphatic amines can improve oil recovery and that the recovery increases with the quantity of amine injected (i.e., with either the amine concentration or the volume of the slug injected). Circumstantial evidence indicated that the increased displacement efficiency resulted primarily from transient changes in wettability of the porous medium. In the Present investigation. additional information has been obtained on the effects of reverse wetting treatments and the mechanism by which increased oil recovery is accomplished. Hexylamine slugs of varying size and concentration were chromato-

Jan 1, 1965

By B. L. Landrum, P. B. Crawford

Theoretical and potentiometric model studies have been made of the effect of non-uniform lateral permeabilities on pattern sweep efficiency and production capacity in waterflood and gas-cycling programs. It is shown that a diflerence in directional permeability by a factor of three may result in a sweep efficiency of only 43 per cent for a five-spot pattern or a sweep of either 79 or 38 per cent for a direct line-drive square pattern, depending on the direction of the line-drive flood. Changes in the pattern conductivity varied from about 0.8 to 1.34 over this same permeability variation, depending on the pattern used. It is suggested that measurements be made to determine the possible magnitude and extent of the directional permeability phenomenon early in the field developinent and certainly prior to the initiation of any fluid-injection progmm. INTRODUCTION Irregularities in reservoir sand properties long have been a major difficulty to anyone attempting to explicitly describe the field characteristics of oil production. In particular, it is well known that vertical and lateral permeabilities often differ appreciably; however, the existence of large regions with lateral permeabflity variation is not widely recognized. A number of years ago, extensive studies were conducted by the Secondary Recovery Research Labcratory of the Pennsylvania Grade Crude Oil Association, primarily on the Bradford field. Johnson and Hughes' reported a permeability trend in the northeast-southwest direction. They indicated that flow in the preferred direction may be 25 to 30 per cent greater in that direction than in the northwest-southeast direction. They also reported that similar effects may be found in other nearby fields. The origin of the permeability variation has been discussed by Griffith Hutchinson' described the results of labora- tory tests on limestone cores, pointing out that per-ferential directional permeabilities were significant in one-half of 10 formations studied and that the average permeability ratio was 16:1. Barfield, Jordan and Moore analyzed the fluid-injection displacement performance in the Spraberry and reported that the effective ratio of permeability parallel to the major axis to the permeability perpendicular to this axis was 144:l and that it inclined along an angle of about 30". To prevent the reduction of sweep efficiency which would normally result from flooding or cycling in such reservoirs, well spacing could be varied from the customary square to a diamond-shaped pattern. Johnson and Hughes' presented a method for determining the desired shape. Another method will be discussed in detail in the next section. It is the purpose of this paper to investigate quantitatively the effect of directional permeability upon the sweep efficiency and production capacity of the five-spot and line-drive well patterns. Information currently in the literature will be applied to this particular problem and augmented by a potentiometric model study. It is stressed that the effects herein discussed may or may not be of great importance because field data on this phenomenon are far too scant to form such an opinion. It is suggested that, if conditions similar to those found by Johnson and Hughes frequently exist (as may well be the case), more consideration cdrtainly should be given to this problem.

By G. C. Bernard, L. W. Bernard, L. W. Holm, W. L. Jacobs

The effect of loam on the permeability of porous media to water was studied as a {unction of foaming agent concentration, specific permeability, pressure gradient, length of a porous medium and its oil saturation. At a given fluid saturation in a porous medium, the permeability to water was found to be the same whether loam was present or not. Foam decreases the permeability to water by developing a higher trapped gas saturation than that obtained by water flooding without loam present, increasing the concentration of foaming agent increased the trapped gas saturation and thereby decreased the permeability to water. The presence of oil reduced the capability of most foaming agents to decrease the permeability of a porous medium to water. A few surfactants were found to be effective loaming agents even in the presence of oil. These results are similar to those reported in a previous paper on the effect of loam on the permeability of porous media to gas. The effect of loam was found to persist in long porous media at moderately high reservoir temperatures and during the passage of many pore volumes of surfactant-free water. INTRODUCTION This paper describes part of a study on a novel approach in the use of surfactants for oil recovery; the use of foam rather than water to displace oil. Previously it was found that foam can displace oil which normally is not displaced by water.4 The foam is formed by successively injecting a suitable surfactant solution and gas into a porous medium. Foam appears to have at least two uses in the field: (1) it shows promise as a superior oil recovery agent, and (2) it shows promise as a selective permeability reducing agent. Foam may be very useful in water floods, or in other oil recovery processes, where highly permeable streaks or unfavorable mobility ratios are a problem. A previous paper reported the effect of foam on the permeability of porous media to gas.1 In the present study the effect of foam on the permeability of porous media to water is reported. The specific objectives of the study were to determine: (1) the effect of foam on the permeability to water in porous media of various specific permeabilities, (2) the effect of foam on the permeability to water in the presence of oil, (3) the effect of foam and crude oil on the trapped gas saturation, (4) the effect of foam on permeability to water at trapped gas saturation, (5) the effect of pressure gradient on the permeability to water under foaming conditions, (6) the persistence of foam during the passage of surfactant-free water through the porous medium, and (7) the effect of various foaming agents, length of the porous medium and temperature on the permeability reduction caused by foam. EXPERIMENTAL PROCEDURES EQUIPMENT AND MATERIALS The experimental apparatus consisted of consolidated and unco n soli dated porous media, wet test meters and constant delivery pumps. The porous media consisted of consolidated sandstone cores (6 to 36 in. long), and unconsolidated sand packs (3 to 30 ft long). The consolidated cores had permeabilities of 32 and 1,000 md and porosities of about 20 per cent. The sand packs had permeabilities of 3,500 to 211,000 md and porosities of about 40 per cent. (Throughout this report a term such as 100 md sand" is used. This term means that the porous medium had a dry, nitrogen permeability of 100 md.) Fluids used in the experiments were distilled water, 1 per cent NaCl solution, aqueous solutions of foaming agents, nitrogen gas, air and crude oil.

Jan 1, 1966

By L. W. Holm, G. G. Bernard

Laboratory experiments were conducted to determine the effect of foam on gas flow in porous media. Previous studies have indicated that foam may be applicable as a restrictive agent in influencing underground gas flow. Foam was found to be exceedingly effective in reducing the permeability of porous media to gas. Consolidated and unconsolidated sands with specific permeabilities of 100 to 146,000 md had, in the presence of foam, gas permeabilities that were less than 1 per cent of the specific permeability; in many cases the gas permeability was practically zero. Foam reduced the gas permeability of loose sand to a much greater degree than that of a tight sand. For example, the permeability of a 125,000-md sand was reduced to 3 md while the permeability of a 4,000-md sand was reduced to 7 md. This effect should cause, to some degree, selective plugging of high permeability channels in various oil displacement processes. The presence of oil in a porous medium decreased the effectiveness of foam in reducing gas permeability; apparently oil acts as a foam depressant. However, it was found that certain foaming agents were very effective in reducing permeability even. in the presence of oil. Also, continuous injection of other foaming agents increased their effectiveness, when oil was present. The effect of foam on permeability of porous media to gas was studied as a function of foaming agent concentration and injection rate, absolute permeability, total pressure, pressure gadient, length of porous system, brine concentration and time. INTRODUCTION The use of surface active agents in flood water to increase the recovery of oil has been studied in the laboratory and in the field for decades, with rather limited success. In recent years a new approach to the problem was proposed. Instead of using an aqueous solution of surfactant as an oil recovery agent, Bond and Holbrook2 proposed that the oil recovery agent be a mixture of surfactant solution and gas. In their method, a water-soluble surface active agent with foam-producing characteristics is injected into an underground formation as an aqueous slug. This slug is followed by gas to produce a foam within the rock. Foams are defined as "agglomerations of gas bubbles separated from each other by thin liquid films'! A foam is fundamentally an unstable system. The foam process for oil recovery has since been studied by other investigators. Fried 1 has shown that foam can displace from porous structures oil that normally is unrecoverable by conventional water or gas drives. This superior oil displacing action is believed to be the result of several factors: (1) foam introduces into the reservoir many resilient interfaces of various sizes and curvatures, which increase the probability that a proper combination of forces for oil displacement will be created; (2) foam has appreciable viscosity which improves mobility ratio and contact efficiency; (3) foam accentuates the trapped gas effect because high gas saturations are possible without producing high gas/oil ratios. Deming4 studied the effect of various foam properties on the displacement of liquid. He found that (1) high foaming ability favors high displacement efficiency; (2) high foam stability is not necessary for high displacement efficiency; (3) displacement efficiency decreases with increase in plasticity of foam and (4) displacement efficiency is unaffected by surface tension of foaming agent solution. One of the important aspects of this oil recovery process is the effect of foam on gas permeability. The simultaneous flow of gas and liquid in porous media has been studied by numerous investigators and a large amount of literature exists on the subject. In this previous work, gas and liquid have generally been considered as essentially independent phases, whose flow characteristics are related through the saturation parameter. Foam, however, is a material with properties that are considerably different from those of its components; for example, the viscosity of a foam is greater than either of its components. Gas flow in a porous medium in the presence of foam should be considerably different

Jan 1, 1965

By Oren C. Baptist, Carlon S. Land

Laboratory research has been conducted to evaluute the effect of clay hydration on the permeability to gas of water-sensitive reservoir sands. Samples of a .sandstone containing trace amounts of montmorillonite and a sample of montmorillonite were .studied in the laboratory to detertnine whether swelling or dispersion was the cause of permeability reduction in these samples. Heliuin, containing various amounts of water vapor, was used to hydrate the clay minerals and to determine the gas permeability at various stages of clay hydration. The amount of water adsorbed by the samples using this method is small. The nonwetting-phase permeability at higher water saturations war investigated by saturating the with water and measuring the permeability to humid helium while decreasing the water saturation, Relative-permeability curves obtained from results of these procedures were used to estimate the effect of the swelling of trace amounts of mont/tlorillonite on the permeability of the .samples. Most of the damage to the permeability when reservoir sands containing trace amounts of montmorillonite are exposed to fresh water is due to dispersion and movement of clays. Blockage of pores by the increased volume of expanded montmorillonite is believed to result in permeability damage that is small in comparison to the observed damage to the samples tested. INTRODUCTION Studies have shown that permeability is severely damaged when sands containing only small amounts of montmorillonite are contacted by fresh water.15 When samples of sands containing large amounts of montmorillonite are placed in fresh water in the laboratory, these samples may completely disintegrate, forming an unconsolidated mass of larger volume than that occupied by the dry sample." In this case, it is apparent that the swelling of montmoril-lonite has destroyed the pore structure of the sand. If only a trace of montmorillonite is present in a sand. samples may remain intact when saturated with water, although the permeability to water is a small fraction of the gas permeability of the dry sample. Many workers in the field of water sensitivity have attributed this reduction in permeability to the blocking of pores and reduction of pore size by the increased volume occupied by expanded mont- niorillonite. if the sand contains a detectable amount of montmorill'onite or mixed-layer clay containing rnontmorillonite. Logically3 the smaller amount of montmorillonite present in a sand, the smaller should he the effect of montnlorillonite swelling on permeability; however, the quantity of montmorillonite sufficient to cause severe damage by swelling is not known. Although hundreds of samples have been tested in our laboratory, no correlation has been established between the amount of montmorillonite in samples and the permeability reduction caused by fresh water. To many petroleum engineers, the phrase "clay swelling" is synonymous with "water sensitivity", or "permeability reduction" implying that any formation damage due to the hydration of clays is caused by swelling. Although all clays adsorb water on their surfaces, montmorillonite is the only clay mineral commonly found in reservoir rocks which adsorbs water between intercrystalline layers, resulting in expansion of the clay particle. As montmoril-lonite swells, the first few layers of water adsorbed between platelets are strongly held and well oriented, and the montmorillonite retains its crystalline structure, although expanded. As swelling of sodium montmorillonite continues, the platelets become farther apart and the forces orienting the platelets in the crystalline structure become weaker, resulting in a less orderly orientation of platelets. In an abundance of water, small groups of platelets may become detached from the original monl-rnorillonite particle and may be dispersed throughout the water phase. Because of its swelling properties, sodium montmorillonite is very easily dispersed in water. Particles of other clay minerals. such as illite and kaolinite may also be dispersed in water. causing water sensitivity of sands not containing montmorillonite. The presence of an immobile layer of water adsorbed on the surface of clays has been considered a possible cause of the low permeability to water of dirty sands. Grim states that the thickness of the layer of immobile water held by sodium montrnorillonite is three nlolecular layers or 7.5 A (angstroms), with some orientation of water extending to 100 A. Assuming a very thick, immobile water layer adsorbed on the surface of a pore represented by a capillary tube, the maximum effect of the water layer on permeability can be calculated. Using a pore radius of 10 ' cm and an immobile water layer of 50 A. the calculation shows the permeability to be reduced only 2 per cent. Similar calculations can be used to show that the effect of electro-osmotic counterflow is of the same order of magnitude as that of bound water. The reduction of the permeability to water by either an immobile water layer

Jan 1, 1966

By E. P. Miesch, P. B. Crawford

A study was made of the effect of permeable and impermeable lenses in a reservoir on the production capacity of a well. Both steady-state and unsteady-state data were obtained. An electrical resistivity model was used to obtain the steady-state data and thermal models. were constructed to obtain the unsteady-state data. The productivity of a well is affected very greatly only when the lenses are close to the well. The effect of circular lenses on the productivity ratio can be correlated with the distance from the center of the lens to the center of the well divided by the radius of the lens. When this dimensionless distance is equal to six or greater, the effect of the lenses on production capacity will be negligible. The pseudo steady-state productivity of a heterogeneous reservoir can be predicted using steady-state data. INTRODUCTION Many analytical solutions of reservoir behavior assume that reservoir rock is uniform and homogeneous. Although this assumption is used, all of the data from core analyses and well logging indicate that the reservoirs are heterogeneous. Very little work has been done on the performance of heterogeneous reservoirs. The work of Landrum,l et al. showed that transient phenomena in oil reservoirs could be studied with thermal models. Pickering and Cotman2 used thermal models to study flow in stratified reservoirs and investigated the effect of inhomo-geneities in oil reservoirs on transient flow performance. Loucks3 made a mathematical study of the pressure build-up in a system composed of two concentric regions of different permeability. Root, Silberberg and Pirson4 studied the effect of the growth of the flooded region on water influx predictions using a thermal model consisting of three concentric cylindrical regions of different thermal properties which simulated the aquifer, the flooded region and the unflooded portion of the original hydrocarbon region. Tomme, et a1.5 made a mathematical study of vertical fractures. The object of this investigation was to study the effect of highly permeable and impermeable lenses in the vicinity of the wellbore on the pressure depletion history of the well. Steady-state data were obtained for both conductive and nonconductive lenses that completely penetrated the formation. The lenses were symmetrically located at various distances from the wellbore. The unsteady-state data were obtained on seven thermal models. EXPERIMENTAL EQUIPMENT AND PROCEDURE STEADY-STATE DATA The steady-state data were obtained from an electrical resistivity model 30 in. in diameter and approximately 11/2 in. deep. The outside of the model was lined with a 30-in. diameter copper strip, which served as the outer boundary of the reservoir. The bottom was covered with a sheet of plexiglass so that it would be nonconductive. The model was filled with a slightly saline solution. The well size was varied from an 0.064-in. diameter copper wire to a 10-in. diameter copper cylinder. Readings were taken with an impedance bridge using AC current to prevent polarization at the contacts. Copper and wax lenses were used to represent infinitely conductive and nonconductive lenses, respectively. The resistance was first measured for each well diameter with no lenses in the reservoir. Then the conductive and nonconductive lenses were spaced symmetrically at various distances from the well and the resistance read from each lens location. The diameters of the conductive lenses were 3, 1.022 and 0.624 in., and those of the nonconductive lenses were 3, 2.25 and 1.563 in. The ratio of each resistance reading to the resistance reading with no lenses in the reservoir

Jan 1, 1965

By D. M. Waldorf

Steam permeability measurements have been made in the laboratory on several samples of natural reservoir materials. The steam temperatures and pressures were selected to simulate conditions which might exist in a reservoir during the injection of steam. For each sample tested, the experimental permeability to superheated steam was comparable to that measured with air and no evidence of plugging was detected. Some samples were exposed to water at various temperatures and plugging was found to occur in materials which contained significant quantities of monmorillonite clay. Temperature had little effect on the degree of plug-ning between 75 and 325 F. The measured pemeabilities tended to increase slightly with temperature, but the changes were small compared with the initial loss of per~neability on wetting. Sequential pemzeability measurements were made on two samples using air, water, steam, water and air, in that order. Both samples were water-sensitive and plugged extensively after the initial injection of water. Upon exposure to superheated steatm the samples dehydrated and their permenbilities to superheated steam were comparable to those initially measured with air. The remaining measuretnetzts with water and air confirmed that the water plugging was reversible and that the samples were not seriorrsly damaged during the tests. INTRODUCTION The swelling of water-sensitive clays during water floods has long been recognized as a potential source of reservoir damage. The recent extensive application of steam injection and stimulation has compounded this problem since both hot water and steam (as well as fresh water at reservoir temperatures) are, at sume time, in contact with the producing zone adjacent to the bore of a steam injection well. The purpose of this paper is to present data which compare the sensitivity of some natural sedimentary rock samples to water at various temperatures, and to super-heated steam. Some properties of montmorillonite clay are briefly reviewed, and comparisons are drawn between empirical data and the predicted behavior of the montmorillonite known to be present in the samples. PROPERTIES OF MONTMORILLONIT E CLAY Water initially adsorbs on dry N a -montmorillonite clay in discrete layers in the interlaminar space between clal platelets. The platelet spacing, which is 9.6 A (angstroms) for a dehydrated clay, has been observed to expand in discrete steps to 12.4, 15.5, 18.4 and 21.4 A spacings, indicating the formation of four discrete layers of regularly oriented water molecules.' The first two layers are easily formed by hydrating a dry sample to equilibrium in an atmosphere with carefully controlled humidity. The formation of the higher layers is more difficult. The usual X-ray diffraction patterns of the more highly hydrated samples indicate a gradual increase in the average spacing betwcen 15.5 and 19.2 A, followed by a discontinuous expansion to 31 A when the weight ratio of water to dry clay is between 0.5 and 1.2.' Platelet expansion above 31 A proceeds monotonically as the moisture is increased and no regular arrangement of the platelets ib observed. Water-sensitivity in sedimentary rocks is usually associated with Na-montmorillonite clay when it is in the noncrystal-line state. Mering3 found that the average lattice spacing of sodium montmorillonite hydrated at 68 F and 70 per cent relative humidity was 15.5 A, and that the spacing, at 92 per cent humidity was 16.5 A. The water adsorbed at the higher humidity has the same free energy as liquid water at 65.6 F. Kolaian and Low' used a tensiometer to measure the thermodynamic properties of water in diffuse suspensions of montmorillonite clays relative to pure water. They observed that water in suspensions as dilute as 6 per cent clay became partially oriented when left undisturbed. The bonding associated with this orientation was not extensive because the free energy difference between the water in suspension and pure water was only a few millicalories per mole. They also found that the measured free energy difference decreased rapidly with temperature and became negligible above 100 F. This evidence indicates that montmorillonites contained in sedimentary rocks would dehydrate to a crystalline structure when exposed to superheated steam, and that the rock permeability measured with steam would be equivalent to that measured with air. The effect of elevated temperatures on the swelline of montmorillonite clays in aqueous suspensions has not been investigated. The Gouy-Chapman diffuse-ion-layer theory has been used to predict the swelling pressure of clay suspensions in dilute salt solutions at room temperature with reasonable success. theory also correctly predicts the direction of the thermal response of Na-mont-morillonite swelling pressures in dilute salt suspensions, 9 Over the temperature range of 33 to 68 F, an increase in

Jan 1, 1966

By R. A. Thomas, R. L. Solbod

The importance of transverse diffusion on the finger development in a miscible-phase displacement at an adverse mobility ratio of tbree was studied in a porous plate 1/4-in. thick, 3-in. wide and 18-in. long. Fast displacement rates (29 ft/D) and slow rates (1.6 ft/D) were used to determine the effect of residence time on the geometry of the fingers. The shape of the fingers was observed directly by use of the X-ray technique. At fast rates numerous narrow fingers were observed, but at slow rates a single somewhat bulging finger was produced. The amount of material moved transversely by diffusion across the plate was sufficient to modify the finger geometry in the slow-rate run because of the long residence time. These results are in contradiction to some of the postulates in the literature. The composition of the effluent stream, however, was not affected by the flow rate. This result is not inconsistent with the observed change in the shape of the finger in a short model, but it seems likely that a short model does not offer adequate and proper scaling of the reservoir. The model used was probably a valid one for studying the effect of transverse diffusion on the finger geometry, but a longer model would be needed for proper scaling of the effect of the change in the finger shape on the efficiency of displacement as measured by the composition of the effluent stream. INTRODUCTION Fingering can be defined as the uneven advance of the injected phase as it moves into a porous medium displacing the resident phase from the pore spaces of the rock. The use of this term is usually restticted to the situation in which the displacing phase is less viscous or more mobile than the fluid being displaced. Under these conditions, not only are fingers formed, but the length and width of the fingers grow with distance traveled in the porous medium. This subject has become one of great interest to the oil industry because of the present trend toward the use of various forms of miscible-phase displacement to increase oil recovery. Since in nearly all of the known modifications of the miscible-phase displacements an unfavorable mobility ratio exists (the displacing phase has a lower viscosity than that of the crude oil), the conditions are proper for fingering to develop. An appreciable amount of fingering appears to be a severe handicap to these processes for it increases the volume of agent required for the process to be a success, and such an increase puts a severe strain on the economics of the proposed processes. In some cases, such as for a mobility ratio of 200 unfavorable, it has already been demonstrated that the proposed process would not be economic if the fingering in the field were to be of the same magnitude as that observed in the laboratory. A number of aspects of fingering have been studied and reported in the literature. While the phenomenon of fingering cannot be regarded as a completely understood subject, considerable information exists on the effect of the path length and the mobility ratio on the growth of fingers. Less-complete data are available on the effect of the diameter of the flow path on the character and amount of the fingering, and even less agreement in results exists on the effect of rate of flow on the nature of fingering. This paper deals with one aspect of this latter subject. OBJECTIVE The objective of this study was narrowed down to one rather specific feature of the behavior of fingers in miscible-phase displacement in porous media. The variable studied was the effect of rate of flow on the nature and the development of fingers. It should be made clear at this point that, while rate was the apparent variable, the real variable was residence time; that is, at low rates the fluids are present at a given spot in the porous medium for a longer time interval than at fast rates. The purpose of the study, therefore, was to determine the changes which occur in the fingering and the possible benefits which might accrue from a longer residence time during that period when fingers are

By J. R. Kyte, R. L. Dalton, V. O. Naumann, R. L. Ridings, R. W. Greene

This paper presents a one-dimensional numerical method for calculating dissolved-gas-drive behavior for a single well or a linear reservoir. Pressure and saturation distributions can be calculated, including the effects of capillary and gravity forces. The analysis has been programmed for solution on the IBM 7090 digital computer. A comparison of experimental and calculated behavior for two linear models demonstrates that the method is suitable for calculating dissolved-gas-drive behavior under reservoir conditions of flow. Dissolved-gas-drive tests were run, at different production rates and with oils of different viscosities, on short (6.ft) and long (173ft) flow models. For tests approximating reservoir pressure gradients and flow rates, there was excellent agreement between experimental and calculated results. A series of calculations of well behavior for thin, horizontal, homogeneous, dissolved-gas-drive reservoirs shows that ultimate recovery is independent of production rate or well spacing in such cases. These calculations also show that Muskat's method and similar material-balance analyses are not applicable under some conditions particularly for high pressure drawdowns. INTRODUCTION The problem of calculating dissolved-gas-drive reservoir performance has received a great deal of attention throughout the oil-producing indusfry for many years. Really significant advances, however, have been made only in recent years. The rapid evolution of high-speed computing equipment and techniques has opened up new approaches to the problem and led to realistic, quantitative analyses free of many of the assumptions and limitations of earlier methods. The one-dimensional numerical method described in this paper includes the effects of both gravity and capillary forces. A number of test calculations, based on different rates, spacings, and rock and fluid properties, have shown that the method is adequate for systems that can be reasonably represented by a one-dimensional model (linear or radial) and in which fluid segregation perpendicular to the direction of flow is not important. Along with the development of mathematical methods, an experimental program was conducted and the data were used to test the validity of the theoretical basis of the mathematical method. Linear models were used and all experiments were run using 1,000-psig bubble-point oils. The use of long, high-pressure models enabled us to more closely approach reservoir gradients and flow rates. MATHEMATICAL FORMULATION DIFFERENTIAL EQUATIONS The differential equations describing dissolved-gas-drive behavior in one dimension dare summarized here. They include Darcy's law for flow of oil and gas (Eqs. 1 and 2) and the equations for conservation of mass for oil and gas (Eqs. 3 and 4). The oil and gas pressures are related by capillary pressure (Eq. 51.