Search Documents

By C. W. Carpenter, P. T. Bail, J. E. Bobek

This paper presents the results of model studies made to demonstrate the validity of applying the "basic" waterflooding scaling relationships of Rapoport1 to water-wet reservoirs composed of communicating strata of different permeability. Results demonstrate that these relationships can be validly applied, and also show the combined effect of variables on recovery. This demonstration was needed to confirm that the sizable effects of communicating strata on waterflood behavior, observed in model studies, were actually representative of field behavior. INTRODUCTION Rapoport1 Collins2 and others have presented scaling laws for modeling water-oil displacements. The validity of these laws has been experimentally demonstrated for homogeneous systems only. Recently, however, several investigators have presented the results of multilayer model studies made to determine the effects of viscous, capillary and gravitational forces on oil recovery in systems composed of communicating layers of different permeability. Dyes and Braun,3 using miscible displacements to simulate water flooding, showed that vertical communication adversely affects recovery at unfavorable viscosity ratios in the absence of capillarity. Ogandzhanyants and Egorova, 4 Gaucher and Lindley,5 and Richardson and Perkins6 modeled water-oil displacements in layered systems using immiscible fluids with capillary forces present. These investigators found more efficient behavior than in the miscible system of Dyes and Braun. However, since the scaling laws had not been verified for heterogeneous systems prior to or during these tests, it was necessary to assume that the results were representative of field behavior. This report presents the results of model studies made at Jersey Production Research Co. to find out if the waterflooding scaling relationships are valid when applied to water-oil displacements in heterogeneous reservoirs with vertical communication. Data obtained in the study show the combined effects of viscous, capillary and gravitational forces on oil-recovery behavior. The studies were performed in strongly water-wet systems where strong imbibition forces were present. For systems of different wettability, the effects of capillary imbibition would be somewhat different than those presented. APPROACH The scaling laws studied were those developed by Rapoport1 and designated by him as "basic" scaling laws. These relationships require that the model be scaled to the geometry of the prototype, and that the model and prototype be operated at the same viscosity ratio and have identical relative permeability curves, contact angles and boundary conditions. It is also necessary that the groups** be identical in the two systems, although any of the parameters making up the groups may be different. Finally, the capillary pressure of one system must be directly proportional to that of the second. Studies were conducted in two four-layer models, representing quadrants of five-spots, which were geometrically similar and which had the same glass-bead packing. The fluids used assured the same wettability in both models. Thus, the contact angles and relative permeabilities of comparative sand layers in the two models were identical and the capillary pressures were directly proportional. However, one model had a characteristic length three times that of the other. To make the scaling groups equal and to check the scaling laws, it was necessary to adjust the rate, water-viscosity level, interfacial tension and water-oil density difference to compensate for the difference in model sizes.

By W. D. Von Gonten, R. L. Whiting

Regression analysis was used to correlate the physical properties of 478 sandstone and 90 carbonate core samples. Porosity, permeability, electrical formation resistivity factor, capillary pressure and the sonic velocity of the shear and compressional waves were measured. Prediction equations for porosity, permeability and electrical formation resistivity factor were found which should be useful in understanding the relationships between the physical properties of porous media in formation evaluation. INTRODUCTION The physical properties of porous media are important to the petroleum engineer and geologist. To evaluate fully the potential and behavior of a subsurface formation such as a petroleum reservoir, certain physical properties of the porous medium must be known. The problem of accurately determining physical properties of subsurface formations has not been solved because many determinations must be made by indirect measurements and because of the difficulties caused by complex pore structure and presence of clays in most naturally occurring porous media. Due to the difficulty in measuring some of the physical properties of porous media, it would be advantageous to be able to predict a certain physical property of a rock from other physical properties of the rock which could be measured more easily and more accurately. Since most potential reservoir rocks are heterogeneous, relationships between the physical properties are very complex and thus far no satisfactory correlations based on theory or laboratory models have been developed. It appears that empirical relationships obtained by measuring the physical properties of a large number of samples of naturally occurring porous media and applying regression analysis to develop equations for one physical property in terms of other rock properties is the best approach. REVIEW LITERATURE The three physical properties used as dependent variables for correlating purposes were porosity, permeability and formation factor. Since formation factor is more difficult to determine, a brief review of the literature is provided on this property. The first work on determining formation factors was published by Archie1 in 1942. He defined this property of a porous medium as F - Ro/Rw......(1) where Ro is the resistivity of the porous medium when completely saturated with a brine of resistivity Rw. Archie found the best correlation between formation factor and porosity was the following equation, where is the porosity fraction and m is the constant characteristic of the rock. The value of m was 1.3 for unconsolidated sand packs, and ranged from 1.8 to 2.0 for consolidated sandstones. In 1950 Patnode and Wyllie2 and De Witte3 observed that the formation factor as determined by Eq. 1 was valid only when the porous medium contained no conductive solids such as clay or shale. When conductive solids are present the formation factor is also dependent on the resistivity of the saturating fluid. Therefore, in samples containing conductive solids the formation factor decreased as the resistivity of the saturating fluid increased. Because of this. the measured formation factor was called the apparent formation factor and was designated Fa. Patnode and Wyllie proposed the following equation,

By Jean-Jacques Prompsy, J. S. Rinehart

Some time ago a study was initiated at the Colorado School of Mines in an effort to arrive at a better understanding of the stress fields developed within drill bits under dynamic loading and the influence that these stress fields could have on the failure of the bits. Particular emphasis has been given to the influence that bit shape has on the establishment of highly localized, potentially destructive stress inhomogeneities within the bit. The study has been divided into three phases involving three velocity regimes: impacts at very low velocities, 0 to 20 ft/sec, a velocity range in which dynamic effects are just beginning to be found; impacts at velocities ranging from 20 to several hundred ft/sec, a range commonly encountered in practical drilling operations; and impulsive loading through the detonation of explosives, a region in which the dynamic effects are greatly exaggerated and made more identifiable. The ultimate objective is to provide basic data which will enable the mitigation through judicious design of the frequency and severity of these transient concentrations of intense stress, thereby prolonging and increasing drill efficiency. This paper presents results of the third phase of this study — the development of a better understanding of the dynamics and mechanics of failures caused by transient stress-wave interactions. The impulsive loading or transient waves in these experiments were generated by detonating small explosive charges placed on the ends of Plexiglas cylindrical rods terminating in truncated cones. Under such intense loading, dynamic effects are greatly exaggerated and, hence, made more identifiable. The explosion generates a high-intensity transient stress disturbance which moves through the rod, reflects from surfaces which lie in its path, and establishes momentarily narrow regions of high concentrations of stress where failure of fracture may occur. The peregrinations of these waves, their interactions and their effects have already been described in the literature (see, for example, Refs. 1 and 2). The locations and extent of the regions quite apparently are strongly dependent on specimen shape, making it possible through judicious shaping of specimens to relocate the highly stressed regions, minimize their extent, or eliminate them entirely. These experiments were designed, first, to identify more clearly the vulnerable regions in geometries relevant to drill-bit design and, then, to modify these regions in a predictable way by changing specimen shape. Through this latter process, it is anticipated that better bit design might evolve. The specimens used in the experiments were cylindrical rods, 1 1/2 and 1 1/8 in. in diameter, ended by a truncated cone, with the cone angles ranging from 45° to 130°. Upper cone face diameters were 1/2, 3/4, or 1 in. A small plastic-covered electric blasting cap detonated on the axis of the specimen produced the spherical shock wave. The detonators used were Olin Mathieson No. 6, which induced into the rod a transient saw-toothed stress disturbance of about 3-microseconds duration, corresponding to a wave length of 0.3 in. in the Plexiglas. Five distinct systems of fracture were observed: (1) along the axis of the rod, a lower system of fractures, the cracks focusing around one point for specimens having cone angles of 80° and 90° and, as the cone angle was increased, spreading downward; (2) just above these fractures, an upper system of fractures composed of horizontal cracks observable in some specimens but not of any appreciable extent; (3) a third system of axial cracks still above the latter, observable in the small cone-angle specimens; (4) in the upper part of the specimen, a system of radial fractures extending from the blasting point; and (5) circular spalling due to reflection of the wave on the boundary of the cylinder. Only the first three systems, the axial fractures, are examined here. LOWER TENSILE FRACTURE SYSTEM The first, lower system of tensile fractures, shown in the photograph of Fig. 1, was observed in almost all of the specimens. It was due to the reflected tensile wave coming in from the boundary of the cylinder as shown in Fig. 2. The reflected wave was cylindrically convergent and focused

By L. C. Kern

This paper presents information on the use of a new propping agent (malleable aluminum particles) which has been used successfully for producing high-conductivity fractures. The conductivity of a packed fracture depends on the width of the fracture as well as the permeability of the propping material. It is significant that the flow capacity of a propped crack is often increased many fold by propping the crack with a sparse distribution of particles. However, the flow capacity of a sparsely propped crack can be limited by two phenomena: (I) crushing of the propping particles, and (2) embedment of the propping particles in the formation walls. Malleable aluminum particles exhibit interesting properties which tend to overcome these limitations. As high loads are applied to the particles they will deform slightly, rather than shattering. The deformation results in a greater bearing area against the formation wall. This reduces the stress on the particle and reduces particle penetration into the formation. The behavior of aluminum particles under load and the bearing capacities of serveral formation rocks have been measured in the laboratory. By properly using these data, as explained in this paper, the flow capacities of aluminum-propped cracks can be calculated for a variety of conditions. Aluminum particles have been used successfully in the field to prop fractures. Where high-conductivity fractures were needed, the aluminum-propped fractures have given greater well stimulation than could have been expected for conventional fracturing processes. INTRODUCTION The purpose of a hydraulic fracture treatment is to create a conductive channel from the wellbore out into the oil or gas formation. The word conductive is emphasized because we believe that the fracture conductivity (flow capacity) is one of the most important factors influencing fracturing results. Previous papers have pointed out the additional well stimulation that can sometimes be achieved by increasing fracture conductivity. Several years ago, we investigated in some detail the effect of fracture conductivity on productivity increases to be expeced from vertical fractures. Both electrical analog and theoretical studies of unsteady-state flow behavior were made. A discussion of the electrical analog equipment and some of the results of the study were given in a previous paper' by Dyes, Kemp and Caudle. Further details were given recently by McGuire and Sikora.5 Fig. 1 shows the remarkable influence of relative fracture flow capacity and length on well stimulation. Briefly, the study indicates the following. 1. High-permeability formations require high-conductivity cracks for best results. 2. The longer the crack, the higher the fracture conductivity required for best results. 3. If the conductivity of the crack is very large compared to the formation (i.e., flow is not being limited by crack conductivity), then the longer the fracture, the better will be the results. Judging from recently published articles,1,314 there is current interest in methods for obtaining high fracture conductivities. Hence, we would like to present in this paper a discussion of a new fracture propping process (use of aluminum particle) which we have used successfully to obtain highly conductive fractures. There are several factors which influence the conductivity of a propped crack. Conventional fracturing techniques frequently result in a crack completely packed with fracturing sand. For such conditions, the crack conductivity is the product of the width of the packed fracture and the permeability of the propping material. It is significant that a large increase in conductivity sometimes can be obtained by propping with particles which are distributed sparsely (less than a

By N. H. van Lingen

A laboratory stud], has been made to determine what factors affect the penetration rate of roller bits, diamond bits and drag bits in rock drilling with clay /water muds. The rather simple relations that exist when pressures in and around the borehole are equal become more conzplicated when under down-hole conditions the penetration rate is hampered by the existence of a pressure differential between the mud at the hole bottom and the pore liquid at cutting depth. Expressions that have been derived for both the penetration rate and the magnitude of the pressure differential in permeable rock together fully account for operating, rock, mud and bit variables. In impermeable rock a similar pressure differential is caused by the bit action itself. In all cases, the pressure difierential and the reduction in penetration rate increase with the effectiveness of the plastering at the bottom of the hole by mud particles. Where bits are employed whose action is largely that of crushing, however, the plastering may become even more effective owing to the addition of rock particles rubbed into the pores of the rock. With roller bits, a plastically hehaving layer may be formed, which causes a further reduction in penetration. In the case of bits whose action is chiefly scraping, moreover, penetration may be arrested by the bit's becoming balled-up. The various adverse effects are reduced by thorough scavenging of the hole bottom. This paper shows how with jet bits the efficiency of such scavenging may he improved by suitable choice of the position of the nozzles. INTRODUCTION The reduction in a bit's penetration rate with increasing depth of hole has been the subject of many investigations. Various investigators ',' have independently reached the conclusion that this reduction occurs not so much because the rock-breaking process becomes more difficult as it does because the lifting of the rock fragments is impaired by a fluid pressure differential holding the fragments down. When such a pressure differential can be avoided, as in air drilling and to a smaller extent in water drilling, penetration rates remain high. Although the range of applicability of these techniques is being extended, most oil wells still have to be drilled with mud. In this paper we present the results of laboratory experiments carried out to determine what factors govern the hold down effect encountered in mud drilling. In the course of this investigation, the well known importance of bottom scavenging came more and more into prominence. Not only may it reduce the magnitude of the hold-down forces, but also it appeared to be a means of preventing the accumulation of a plastic mass of cuttings and mud on the hole bottom and bit which may cause penetration to cease almost completely. So that the effects of fluid pressures may be more readily understood, the factors that govern penetration rate in the absence of these pressures will first be discussed. The effect of down-hole pressures will then be assessed for the case of permeable rock where fluid pressures are well defined, and it will be demonstrated how rock strengthening due to confining can be accounted for. Subsequently, fluid pressures in less-permeable rock will be examined; we will show that in impermeable rock it is the bit action that governs fluid-pressure distribution and hold-down. In the case of poor scavenging and/or high bit load. cutting cake manifests itself. The effect, like various others. appears to depend on the type of bit used. From the very beginning, therefore, we shall base our discussion on the three main types of bits used in rotary oil-well drilling— roller bits, diamond bits and drag bits. Since space is limited, it will be impossible to make more than a passing reference to many aspects of the problem. Moreover. certain relations that are not of primary irnportance to our argument will unavoidably have to be stated in a somewhat dogmatic fashion, and the experimental data on which these are based necessarily must be omitted. EQUIPMENT The drilling experiments were performed on three machines, of which one has been designed for the realistic simulation of down-hole pressure conditions. HIGH-PRESSURES MACHINE In this machine (diagramed in Fig. 1) the mud, pore and confining pressures can all be adjusted independently. The foil-covered rock sample contained in the inner pressure vessel is confined by oil that is pressurized by means of a hand pump. The pore space of the watersaturated sample is connected to an arrangement which keeps the pressure in the pore space constant regardless of the rate of filtrate flow through the sample.

By M. R. Tek, K. H. Coats, D. L. Katz

The performance of a gas reservoir su,bject to water drive is often affected by interference due to gas procluction or injection in neighboring reservoirs adjacent to a common aquifer. Field data available on production-pressure behavior of several reservoirs indicate that prediction of reservoir performence must include the effects due to mutual interference. A method developed to predict the performance of two or more interfering gas reservoirs, subject to specified production-iniection .rchedule.s, is presented. The method involves the comhinatiorz of the theory of unsteady-state water movernent with appropriate meterial balance consideration.s. The procedure which extends M. Mortada's treatment for interfering oil reservoirs makes use of his dimen-.sionless pressure drop quantities. prevented graphically in the literature. In order to evaluale the predictive accuracy and practical field applicability of the theory, the method has been applied to analyze the performance of two actual gas reservoirs in mutual interference. The calculations were performed on an IBA4 704 computer. The predicted performances were alro calcltlated with the interference effects ignored. The results show that considerably better agreement between calculated and observed behavior ic. obtained when the interference effects are included. INTRODUCTION Fluid communication between any number of oil or gas reservoirs adjacent to a common aquifer results in mutual interference effects on their production-pressure behavior. The theory, application and interpretation of pressure build-up behavior of oil wells, drawdown tests on oil or gas-producing formations and behavior of several storage reservoirs indicate ample experimental evidence of mutual interference. A method developed to compute the behavior of interfering oil reservoirs was published several years ago by Mortada.4 Studies conducted during the recent years on underground storage of natural gas clearly indicated the need for extending Mortada's work to the case of interfering gas reservoirs. The adoption of appropriate material balance equations in conjunction with the unsteady-state pressure interference theory yields pressure-explicit equations. Their solution, their application to a specific case history and an over-all evaluation of the developed method through the comparison of predicted and observed field performances were the main objectives of this study. DEVELOPMENT OF MATHEMATICAL RELATIONS The diffusivity Eq. 1 governs planar, radial, unsteady-state compressible liquid flow through a porous medium. ?p/?r2+1/r ?p/?r=µFc/k ?p/?t (1) The 1iterature1,2,3 contains derivations and solutions of

By M. R. J. Wyllie, Michael A. Torcaso

With the continued deep drilling of today, increasing numbers of high pressure and high temperature gas-condensate reservoirs are being discovered. Correspondingly, the ranges of properties of gas-condensate reservoir fluids are also being increased. The literature1*'." contains several excellent reviews of the property ranges in presently known condensate reservoirs. When one examines these property ranges, it becomes apparent that the need for reliable data is great if the engineer is to evaluate properly or appraise a gas-condensate reservoir. It is the purpose of this paper to present phase-behavior data that will point out some of the important variables that influence recovery from gas-condensate reservoirs produced by pressure depletion. All these variables may be resolved into two major categories: the hydroc:arbon composition of the reservoir fluid and the physical properties of the reservoir. But the variables which we will consider are those which are usually determined from well tests shortly after discovery and are, consequently, readily available to the engineer. They are as follows: (1) the producing separator gas-stock-tank liquid ratio, (2) stock-tank liquid gravity, and (3) reservoir temperature. Laboratory depletion studies have been performed on several reservoir fluids to illustrate the effect of changes in the cronsidered variables upon ultimate recoveries. The deviation factors, Z, of a series of reservoir fluids were determined experimentally to illustrate their order of magnitude and their effect upon the calculation of in-place reserves existing at high pressurcs and temperatures. The procedures followed in conjunction with the performance of these studies and their associated calculations are illustrated in Fig. 1. For the purpose of this paper, the depletion recoveries are expressed only in terms of pcr cent of initial stock-tank liquid, with the abandonment pressure assumed to be 500 psig. VARIATION OF RECOVERY WITH INITIAL PRODUCING SEPARATOR GAS-STOCK-TANK LIQUID RATIO The producing separator gas-stock-tank liquid ratio has long been one of the principal criteria in assumptions and generalized correlations' used for the estimation of ultimate recovery from gas-condensate reservoirs. Fig. 2 presents a plot of the cumulative recovery of stock-tank liquid at abandonment as a function of the initial producing separator gas-stock-tank liquid ratio. The graphical illustration represents the results of 44 laboratory depletion studies performed on actual gas-condensate fluids. These data are presented not only to illustrate the general trend of ultimate recovery increase with initial producing ratio increase but to emphasize the considerable variation in possible recovery at any given producing ratio. Table 1 presents the properties of two actual reservoir fluids chosen to illustrate the approximate extremities of producing ratios for gas-condensate systems. The fluids designated as A and B exhibited initial producing ratios of 3,400 and 114,000 scf/bbl, respectively. Fig. 3 presents a graphical comparison of the cumulative stock-tank liquid recoveries, calculated from the well-stream compositions and produced volumes for each fluid during depletion. It is to be noted that Fluid A exhibited a recovery of only 16.5 per cent when produced to an assumed abandonment pressure of 500 psig while the corresponding recovery for Fluid B was calculated to be 62.2 per cent. This difference in recovery, which exemplifies the general trend of recovery increase with initial producing ratio increase, is readily explained by consid-

By M. R. Tek, K. H. Coats, D. L. Katz

The nature and the limits of validity of Darcy's law US applied to the flow of natural gas through reservoirs has been considered in order to resolve some controversial aspects of the effect of turbulence on pressure drops. The equivalence between various concepts and viewpoints advanced in the past by several investigators to explain how and why a gas well does not necessarily perform according to Darcy's law is shown. Starting with generalized equations of flow of fluids through porous media, a partial differential equation has been derived which accurately represents the flow at all rates. This equation has been numerically solved using an IBM 704 digital computer. The results permit plots of unsteady radial pressure distribution curves from which specific isochronal backpressure curves may be constructed. These back-pressure curves show the effect of the ,8 factor on the slope of the back-pressure curve. The calculations further indicate that the drainage radius for a gas well in turbulent flow propagates at a rate dependent upon the rate of production at the wellbore. This is quite different from the case with liquid flow or natural-gas flow in laminar regime. Additionally, the effect of reservoir inhomogeneities and crossflow between layers of different permeability on the back-pressure performance of gas wells has been conridered. In light of the current numerical results the significance and linlitation of the rate of flow function Y proposed by Smith" has been discussed. INTRODUCTION The relationship between the pressure drop and flow rate in problems of fluid flow through porous media is known to be affected by the nature of flow through the porous matrix. It has been observed by many that, for a range of flow rates, the pressure drop remains proportional to the rate of flow. When some flow rate is reached, however, it is usually observed that the pressure drop gradually begins to increase more than proportionally to the flow rate. It is well known that this phenomenon was first observed by Osborne Reynolds in 1901 in experimenting with flow through pipes. In his classical experiments, Reynolds made visual observations on the condition of streamlines evidenced by injecting a dye into water flowing through glass tubes. In these experiments, the abrupt transition between steady, "streamline, laminar" flow and unsteady random turbulent flow was found to be a function of the dimensionless group (Dvp/u), now known as the Reynolds number. During these experiments, in addition to observations on the nature of flow regimes, the proportionality between flow rate and pressure drop in laminar flow was contrasted with the nonlinearity between these variables in turbulent flow. Fancher and Lewis' reported data on various consolidated and unconsolidated sands in 1933. Their conclusions were that ". . . the flow of fluids through these porous materials closely resembles that through pipes; that there is a condition of flow in porous systems which resembles viscous flow, another which corresponds to turbulent; that the change from one type to the other takes place at a definite and reproducible condition for each system". In 1947, Brownell and Katz published a method to predict the laminar and turbulent flow behavior from the particle size, bed porosity and the particle sphericity, employing the friction factor-Reynolds number charts for pipes. Several investigators have verified the work of Fancher, Lewis and Barnes and presented their data as friction factor-vs-Reynolds number plots.' The equation which would represent the pressure gradient over the whole range of velocity must have an added term over that represented by Darcy's law. Accordingly, the pressure gradient necessary to sustain flow at the velocity (v) through a porous medium may be represented by the following equation, suggested by Forscheimer." The nature and the range of validity of Darcy's law has been the subject of studies by many investigators over the past years. While everyone seemed to agree on the need for a quadratic correction term to Darcy's law to make it effective over the range of velocities, the concept of inception of turbulence and the use of the term "turbulent flow" remained controversial. Some fluid dynamicists define turbulence as a flow

By Ancel St. John

At the New York meeting of the Institute of Metals Division in February, 1927, Jeffries and Tarasov presented a paper on Tungsten and Thoria,' in which the experimental facts were interpreted in accordance with the prevailing point of view that the thorium is present as the oxide, thoria. In discussing that paper the writer of this article stated that his own investigations, which had been unavailable to Messrs. Jeffries and Tarasov, showed that the so-called "thoriated" tungsten filament failed to disclose any sign of the oxide when subjected to X-ray diffraction analysis, but did give plenty of evidence of the presence of metallic thorium and, in certain cases, of a material believed to be an alloy of tungsten in thorium. Since these results were completely at variance with the results reported in the paper under discussion and with other papcrs which had preceded it from time to time in the past, it was felt that further discussion should be withheld until the new evidence could be presented adequately. Arrangc-mcnts were therefore made to present that evidence at this meeting. Material Examined The material exalllined was high-ernission filament such as is used in radio tubes, supposed to contain about 0.75 per cent. thoria. The actual analysis was not stated. Scope of the Investigation The investigation was undertaken because of a feeling that the thorium present in filamentary material of this type exists as metallic thorium or as an alloy of thorium and tungsten and not as the oxide. Attention was therefore directed to determining: 1. The minimum amount of thoria detectible in a mechanical mixture of thoria powder and finely divided tungsten. 2. The amount of t.horia detectible in the filament. 3. The condition of any thorium present in other form than as thoria. 4. The amount of thorium dissolved in the tungsten.

By R. L. Perrine

This paper presents results obtained by consider ation of unstable miscible displacement within a porous medium as a quasi-turbulent displacement process, that is, a process characterized by fluctuations around a simpler, mean flow behavior. The instability usually results from the miscible displacing fluid having a viscosity lower than that of the displaced fluid. Such instability is evidenced in suitable experiments by the presence of viscous fingering. A unified theoretical relationship applicable to both stable and unstable displacement is derived. Parameters in the theory are obtained by empirical means. Although the theory describes an idealized conceptual model, agreement with experiment is quite reasonable up to moderately adverse viscosity ratios, such as 17:1. INTRODUCTION Many experimental studies of miscible displacement in porous media have been conducted with prediction of reservoir displacement efficiency as the ultimate objective. Most such studies have utilized lower displacing than displaced fluid viscosity, scaled to potential reservoir fluid pairs. Theoretical approaches have been largely limited to unit viscosity ratios, however, in spite of the necessity for an understanding of the mechanism of the unstable, adverse viscosity ratio displacement process. An obvious reason is the difficulty of describing in mathematical form the viscous fingering characteristic of these conditions. Observation of experiments conducted with dyed fluids in transparent systems suggests an analogy between unstable miscible displacement and turbulent flow in a pipe. In both cases there are fluctuations around a simpler, mean flow behavior. An important difference is that flow behavior of interest in the porous medium is entirely transient, contained within a transition zone between displacing and displaced fluids. Transient behavior complicates description in that coefficients in the equations become variables rather than constants. In the study reported here, the analogy with turbulent flow has been used in treating unstable miscible displacement as a quasi-turbulent displacement process. The purpose has been to derive, even if restricted to an idealized conceptual model, a unified theoretical relationship applicable to both stable and unstable displacement. A relationship meeting these specifications up to moderately adverse viscosity ratios, such as 17:1, is presented. One fluctuation parameter in the theory and dispersion coefficients are obtained by empirical means. The idealized theory is compared with experimental results. UNSTABLE MISCIBLE DISPLACEMENT AS QUASI-TURBULENT DARCY FLOW The miscible displacement of one fluid by another within a porous medium is considered in this study. Flow conditions are such that Darcy's law is applicable. It is further assumed that, by the stability criterion proposed by perrine,1 initial flow conditions lie well within the regime of instability. Thus substantial viscous fingering is assured. We wish to show how this flow regime can be represented as quasi-turbulent. That is, the Reynolds number for the established flow conditions is below that for inertial or turbulent flow within a porous medium, and lies within the regime for which Darcy's law is valid. Yet flow behavior can be described as the combination of some relatively simple average result, and characteristic fluctuations that are superimposed on the simpler behavior. Stated another way, flow behavior includes the movement of layers of fluids with different velocities past or over one another. Such descriptions are characteristic of turbulent pipe flow. The fluctuations in the present case are a direct consequence of the local viscous fingering which accompanies the unstable displacement process. Should displacement become stable, fluctuations would die out. It is of particular importance to note that fluctuations such as these can interact in a way that contributes to material transport by the basic flow. The basic transport equations required for engineering calculations must be modified to reflect this fact.

By F. W. Jessen, C. A. Johnson

The effect of treatment with ferrochrome lignosul-fonate on both sodium and calcium bentonites has been examined. In the early stages of treatment it appears that some base exchange of iron and chromium for sodium may occur with the sodium bentonite. Such exchange appears specific for exchange of sodium ions when the calcium ion concentration represents 40-5.5 per cent of the total exchangeable cations, and when additions of lignosulfonate are less than 1 lb per bbl. Adsorption appears entirely physical above 1 lb per bbl treatment level for all concentrations of calcium studied. Gel strength and plastic viscosity of both calcium and sodium bentonite systems were reduced with treatments of ferrochrome lignosulfonate up to 2 Ib per bbl, the greatest reduction occurring within a range of calcium concentration of 40-5.5 per cent of the total exchange capacity of the clay. INTRODUCTION The problem of control of viscosity and gel strength of drilling fluids by the action of thinners or dispersants has been studied by many investigators.1-7 In addition, much work has been done on formation of gel structure in clay suspensions. However, gel formation and structure still are incompletely understood, and the use of the various dispersing chemicals is based, to a large extent, on data almost entirely empirical in nature. In recent years the lignosulfonates have become increasingly important in the treatment and control of drilling fluids. More specifically, ferrochrome lignosulfonate has found rather wide acceptance, particularly in calcium systems.5-7 Lignosulfonates are refined lignin products made from spent sulfite liquor from which the fermentable sugars have been removed.15 The exact structure of lignin is not known; however, it is believed to be a phenolic propane type polymer which may exist as a branched chain or cross linked structure. The molecular weights of lignosulfonates probably range from 300 to 100,000 or more.15 A better understanding of the mechanism of adsorption of dispersants is essential before improvement and development of characteristics of drilling fluids and drilling fluid control is possible. This investigation was undertaken to determine the action of ferrochrome lignosulfonate on sodium and calcium montmOrillonite' LABORATORY INVESTIGATION A Stormer viscosimeter and Fann V-G meter were used to measure viscosity and gel strength. A Zeiss spectrophotometer was used for analysis of Fe and Cr and the flame attachment was employed in the analysis of Na, Ca and Mg. A General Electric X-ray diffraction unit utilizing Copper Ka radiation and a Ni-filter was utilized in the X-ray studies. PROCEDURE Two systems were studied. Predominately sodium clay, obtained from Wyoming, and a predominately calcium clay from Texas were used. These particular samples were chosen because of the extensive work of characterization completed by Mungan. Ionic base exchange capacity determined by flame photometry and total exchange capacity by the Kjehldahl method were obtained experimentally. A 7-per cent suspension of the sodium clay and a 20-per cent suspension of the calcium clay in distilled water were prepared by mixing with a high-speed mixer. The suspensions were allowed to stand for two weeks for complete hydration. Calcium chloride was added to portions of the sodium and calcium cla suspensions to obtain suspensions varying in Ca++ concentration as shown in Table I. A series of treatments was run with 500-ml portions of the original suspensions using a ferrochrome lignosulfonate. Treatment was carried out in the same way with each of the suspensions of varying Ca concentration. Viscosity and gel strength were measured and a filtrate obtained at each treatment level. Ten ml of filtrate were diluted to 50 ml, and aliquot portions of the 50 ml were taken for the Na, Ca, and Mg analysis by the flame photometer. The pH of these filtrates ranged from 7 to 7.3. Standard

By M. L. Millican

The objective of this study has been to develop a method of log interpretation in the Delaware sand whereby the effects of the shale contained within the sand can be recognized and accounted for in a quantitative solution for porosity and water saturation. Both theoretical and empirical considerations are discussed. A simplified interpretation chart is presented which provides a means of amending the conventional porosity resolution of the Sonic log to compensate for the presence of shale. The amount of shale within the sand and its degree of lamination may be inferred from a calibrated gamma ray log. The resulting porosity value may then be used in conjunction with a knowledge of the formation resistivity to determine water saturation. Of added interest, particularly in the Delaware sand, is an approximation of permeability which may be derived empirically from the simplified chart. The discussion is illustrated by field examples. INTRODUCTION The performance of the Sonic log in defining the porosity fraction in clean, consolidated formations is well known, particularly in sandstones that are well compacted and free of shale. An approach to Sonic log porosity determination in shaly sands has been recently developed.' This method involves correcting the Sonic log readings by means of an a factor, where a is the ratio of the pseudo-static SP to, the static SP, and is a function of the amount of shale present. This method is dependent upon a contrast between the salinities of the drilling fluid and the connate water in the formation, inasmuch as this contrast produces the SP effect. In wells drilled with salt mud, the contrast, and with it the means of obtaining the a correction, disappear. In the Delaware basin, the large majority of rotary wells drilled to the uppermost sands of the Bell Canyon group employ salt mud as the drilling fluid. These sands have been the objective of most of the drilling exploration in the basin to date. A comparison of porosity values in these sands, as determined by core analysis, to Sonic log readings has indicated that a correction for the shaliness of the sand is needed. In the absence of an SP curve to determine the correction factor, the porosity resolution of the Sonic log in these sands has been questionable.