Search Documents

By R. L. Perrine

A stability theory is developed for miscible liquid-liquid displacement within a porous medium. In the usual case considered, a high-density high-viscosity "oil" is displaced downdip by a low-density low-viscosity "solvent'! Perturbation methods are used to find the conditions under which the spreading mechanism changes from the stable dispersion process to unstable viscous fingering. We find that instability is conditional, that there is a dependence on the shape of a disturbance leading to a"diameter" effect and that very difficult experimental scaling problems may result. A useful consequence is the definition of a minimum "slug size" for stable miscible displacement. This should make possible optimum use of the solvent process for oil recovery. The results may apply to many situations in which one fluid displaces another of somewhat different fluid properties within a porous medium. INTRODUCTION The process of miscible liquid-liquid displacement looked very favorable when it first received serious consideration as a means to increase petroleum recovery. Some early laboratory results were interpreted to mean that only a small "slug" of solvent was needed, perhaps 2 to 3 per cent of a hydrocarbon pore volume. This "slug" could recover the oil in an entire reservoir provided the fluids remained miscible.1,2 Only a slow growth of the "mixing zone", or gradual transition from oil to solvent, should occur. This would follow naturally from mixing by a dispersion mechanism such as that described by Scheidegger.3 And, indeed, laboratory cores have shown this kind of behavior provided solvent viscosity was at least as great as that of the oi1.2,4-6 The same kind of behavior has also been observed with solvent viscosity lower than that of the oil after the distance traversed became large. Not all laboratory results have been this favorable, however.7-9 Mixing zone growth may become proportional to the distance traveled rather than the square root of this quantity. Such behavior accompanies high rates in short systems of large diameter, provided the viscosity ratio is adverse. If this more rapid spreading were to persist in a reservoir, the solvent requirement for miscible "slug" displacement could exceed 30 per cent of a pore volume. There is considerable economic importance in the difference between 3 and 30 per cent solvent. A 1ogical conclusion is that two different kinds of flow behavior are possible, with each kind leading to a different result. One kind gives efficient displacement, the other does not. In the efficient case, the solvent bank spreads out only by the mechanism we have termed dispersion. Under these conditions, displacement is as near piston-like as possible. In the second kind of flow, found only with adverse viscosity ratios, viscous fingering occurs. That is, permeability variations cause a small finger of low-viscosity solvent-rich fluid to move ahead of its average position within the mixing zone. This creates a path of low resistance to flow, and an even greater amount of solvent-rich fluid follows. Thus, the process is autocatalytic. Once started, the fingering mechanism rapidly becomes dominant. The question to be answered is this. Under what conditions will each of these two different kinds of flow occur? In particular, what conditions denote the transition from dispersion to viscous fingering as the solvent spreading mechanism? The answer may tell us whether or not the desirable, near piston-like miscible displacement process is practical. An answer to this problem can be obtained from theory by the use of perturbation methods. The procedure is as follows. We first formulate a mathematical representation of the system. Then a small disturbance (or perturbation) in solvent concentration profile is introduced and observed to see what happens. Unstable flow is indicated when a small disturbance will grow larger. This will lead to eventual viscous fingering. If on the other hand the disturbance dies out, the displacement is stable. In this latter case, with dispersion as the spreading mechanism, near piston-like displacement is possible. Thus, this paper presents the development of a stability theory for miscible liquid-liquid displace-

By H. Rumfelt

The author demonstrates an approach for analyzing overcasting requirements for a stripping project. This approach to the problem employs indicated trends in the relationship of the weight of the machine (only electropowered draglines are considered) to its ability to do stripping work. Strip mining in the domestic coal industry is contributing a greater proportion of overall production than at any time in the past. U.S. Bureau of Mines statistics show that the percentage of total production mined in 1958 by stripping methods increased to more than 30 pct, up from 22 pct in 1951. The tendency to favor strip and open pit mining exists in other mining activities, but such trends are not as clearly outlined as in the coal industry. Significantly, the depth of overburden in coal properties considered suitable for stripping is also increasing. A depth of 50 to 60 ft of overburden was once considered to be about the maximum amount that could be handled economically. Today, a depth of 80 to 100 ft is frequently taken into consideration for simple overcasting operations (i.e., excavating overburden and dumping it in the spoil position where it remains indefinitely). Accompanying these changes, there have been continuing studies searching for procedures to appraise deeper stripping problems. Although the only known practicable way to accurately evaluate a proposed stripping venture is to make studies (including cost estimates) step by step, it is nevertheless desirable to have a quick method for making preliminary evaluations of potential strip mining operations. A recently devised digital computer method of providing a preliminary simple overcasting analysis for a stripping prospect relates the indicated trends in the relationship of the weight of the machine to its ability to do stripping work. In addition, this approach employs situations where the geometry of each cut-and-spoil section is assumed to take certain defined relationships for varying overburden depths. Slopes are necessarily considered to be stable, which means such practical factors as the mechanics of soils are neglected for the sake of convenience. Definition of Terms: Photographs A and B show a shovel-type and a dragline-type of simple overcast operation, plus terms used to denote certain parts of a pit. Pits that follow a straight line when projected on a horizontal plane are referred to as straightaway. It has been assumed that all section drawings accompanying this article pertain to straightaway cuts which, in turn, means that areas can be compared by relative volume. The volume of overburden is measured in virgin cubic yards. Material displaced from its virgin state to a spoil area normally occupies a larger volume than when in situ; the difference is termed swell, and it is expressed as a percentage of the original volume. For example, if the original cut volume is denoted V cu yd, and the spoil volume is 1.2 V cu yd, the swell is a positive 20 pct. MAXIMUM USEFULNESS FACTOR (MUF) CONCEPT Men experienced in overburden removal give emphasis to the weight of the machine in relation to its value as an excavator. From the writer's understanding, value has been only loosely defined and discussed generally and vaguely. The concept is considered important, however, since it is based upon the observations of many men. It provides a clue for the approach taken to solve deep stripping problems. In the procedure described below, value has been arbitrarily established as the product of the nominal dipper size of the shovel (or dragline) times a functional dumping reach. MUFs—Shovel: A shovel's capability to handle thick overburden is usually limited by its ability to dispose of the spoil. Thus the dumping reach, as well as its dumping height, is significant. The angle of repose of spoil material must also be considered. Although such slope angles vary among different mines, jobs, and materials, a slope of 1.25 to 1, frequently found in both planning and practice, is used by this author in the illustration of this analytical method. The relationship of the shovel's geometry to a stripping operation is shown in Fig. 1. To gain maximum advantage in constructing spoil piles, the shovel is placed so that its tracks adjacent to the spoil are as close as possible to the rib of the coal deposit. A vertical plane passing through the line

Jan 1, 1961

By W. S. Gillam

A need for new supplies of fresh water exists today and in many specific areas that need is urgent. One solution lies in saline water conversion, a problem complicated by cost factors. The principles involved in saline water conversion, the status of development, and the estimated costs (present and future) of several processes are presented. Among the methods discussed are distillation, electrodialysis, and freezing. In general, the costs presented are based on a standardized procedure for estimating conversion costs, permitting a valid comparison among the various processes. The need for new supplies of fresh water and the potential benefits to be derived from an abundant supply of converted water are recognized by practically everyone concerned with water problems. The water supply problem exists today; it is urgent in many specific areas in this country and also in the world, and it will become more acute in the future. One answer to the growing problem of adequate water supplies is the development of new sources. Very significant quantities of brackish underground and surface waters exist in certain areas and an inexhaustible supply of ocean water is available. Thus in many areas water resources can be extended through saline water conversion. Congress recognized the need for new sources of fresh water in 1952 and passed the Saline Water Act, Public Law 448, amended it in 1955, and in September 1958, enacted Public Law 85-883, calling for the construction of at least five demonstration plants. The program is administered by the Dept. of Interior through the Office of Saline Water, and its primary objective is to reduce the cost of converted water produced, whether it be by development of new processes or improvement of known processes. This is a most difficult problem and one that will require several years of prodigious effort. It is difficult—not because of any intricate or new chemistry, engineering, or physics involved—but because of the difficulty in converting water at low cost. Whatever the sources of the saline water, the salts which are held tenaciously in solution must be removed before the water becomes suitable for industrial or domestic uses. Saline water is a relatively simple system of salts dissolved in water. It has certain chemical and physical properties that determine the various methods by which the salts may be separated from the water. The system, although not complex, in most instances, has had countless years in which to reach equilibrium and is, therefore, comparatively stable. Because of its stability, separation of saline solutions requires relatively large quantities of energy. The unique properties of water depend on the fact that its molecules are chemically active. The chemical and physical properties of water are associated with the type of bonding involved in the water molecule. Chemical changes such as hydrolysis, or rusting of iron, involve the breaking of chemical bonds between the hydrogen and oxygen atoms. Physical changes, such as evaporation in a boiler, the melting of ice, or the viscous resistance to flow in a pipe, involve breaking of the hydrogen bonds. (The hydrogen nucleus is so small that it can attract two negative atoms.) Thus water molecules not only combine with molecules of other compounds but even with one another; e.g., each molecule may be bounded to four other molecules. Water molecules cling to the ions of dissolved salt to form water-encumbered hydrated ions and they cling to one another to form entangling networks through which hydrated ions can be propelled only by tearing the networks apart. That is one reason why considerable energy still needs to be expended in our simplest procedures for purifying water. If water molecules did not have this habit of clinging so tenaciously to other molecules, and to one another, it would be easy to push salt ions past the water molecule and get a separation. But the water would not then dissolve salt, so the problem would not exist.' Water when heated evaporates very slowly, relative to other liquids having simple molecules. Vaporization involves the separation of molecules from the liquid, and this means overcoming the attraction between molecules which is due to the hydrogen bonding. The heat of vaporization for water is high; consequently, the boiling point of water is also high. Water boils at 100" C; hydrogen sulfide (H2S) at -60" C; oxygen (02) at -183" C; nitrogen (N2) at -196" C; and methane (CH4) at-161°C, even though the latter has about the same molecular weight as water. Because of these peculiar properties of water, it exists as a liquid on earth instead of a gas such as hydrogen sulfide or nitrogen and oxygen.

Jan 1, 1961

By H. D. Outmans

Present first-order theory for frontal stability and viscous fingering of immiscible liquids is improved by including the nonlinear terms in the equations describing conditions at the interface of the liquids. This leads to the revision of several conclusions which were based on first-order theory. They concern the growth rate and changing shape of a sinusoidal disturbance, particularly in relation to the wave number of this disturbance. Nonlinear aspects of effective interfacial tension are discussed and it is shown that this tension is not simply a positive proportionality constant in a linear relation between pressure difference and curvature at the interface. Scaling requirements are determined from the dimensionless groups which govern fingering. Gravity and interfacial tension invalidate a previously formulated conclusion that the shape of a finger after a given displacement is independent of the displacement velocity. Also, similarity of fingering (and, hence, of sweep efficiency in case of an unstable front) requires geometrical similarity of the initial disturbance in the model and the reservoir with a scale factor which is the same as the one for scaling down the dimensions of the reservoir. INTRODUCTION The critical velocity at which the transition zone in the vertical displacement of immiscible liquids becomes unstable was calculated by ill. More recently, the same was done for nonvertical displacements and it was also shown that the stability is affected by an effective interfacial tension between the liquids. 2-4 Although in the oil reservoir a stable front is always desirable from a recovery point of view, the necessary velocity may have to be so low that the corresponding production rate is not economical. In that case, a knowledge of the fingering unstable front is necessary for predicting recovery at break- through. Studies in this direction have not gone beyond the very moment at which fingering first occurs.5,6 Conclusions about stability and fingering in these references are all based on linear theory. In this theory the nonlinear terms in the equations describing conditions at the front are neglected. The usual justification for this lack of mathematical rigor is that the nonlinear terms can be made small relative to the linear terms and, supposedly, small causes produce only small effects. It has been recognized, however, that this is not necessarily true in the nonlinear problems of hydrodynamics.7 The hydro-dynamic stability and fingering of the front between two liquids, accelerated in a direction normal to the front, for instance, is strongly influenced by the nonlinearity of the boundary conditions.8-10 Since these boundary conditions are not unlike those arising in the problem of stability and fingering of the front during a slow immiscible displacement in porous media, it was thought that there, too, the nonlinear terms should be taken into consideration. A summary of the results obtained in linear theory and a comparison with experimental data precede the nonlinear theory developed in this paper. This summary serves to introduce some quantities which will be used in the further develop ment. It should also be pointed out that the method of higher-order approximation by which the nonlinear stability problem is solved has application in other reservoir studies.11 SUMMARY OF LINEAR THEORY 1. A plane interface remains stable if its velocity is smaller than a critical velocity. vc = (p2-p1)gn/(l/A2- l/A1 .........(1.1) This equation is derived for uniform flow in a thin layer inclined to the horizontal plane. The displacement takes place in a direction normal to the intersection of this layer and the horizontal plane. The undisturbed (line) interface is normal to the velocity. Fluid 1, the upper fluid (gas), is displacing Fluid 2 (oil). (VcQS would be the injection rate in cc/

By AIME AIME

AT bottom of this page is a photograph recently taken by a student-John E. Stewart-of the most northerly situated college in the world, the Alaska Agricultural College and School of Mines. It is situated at College, near Fairbanks, Alaska, approximately on the 65th parallel of latitude, and therefore just south of the Arctic Circle. It is farther north than Hudson Bay, and in about the same latitude as southern Greenland.

Jan 1, 1932

By J. F. Abel

It can be concluded that rock mechanics instrumentation, geologic mapping, and operations research in combination will produce an accurate estimate of tunnel support requirements for establishing a steel-set load prediction model. This estimate can be used to increase safety and to reduce tunneling costs. Unnecessary supports can be eliminated; overloaded or insufficient supports can be replaced on a more scientific basis than previously possible. Only a limited number of measurements are required to develop a steel set load prediction model. The results of this study can be applied generally; but, every tunnel must be considered as a separate problem. The construction and geologic variables which were found to be significant in the prediction of rock loads in this tunnel were: steel section modulus, set blocking points, percentage of alteration, relative degree of faulting, relative water condition, distance to the nearest fault, and average joint spacing. INTRODUCTION A combined engineering geology and rock mechanics program was performed throughout the Straight Creek tunnel pilot bore in Colorado in order to develop a technique for estimation of steel tunnel support loads. The technique involves the construction of a statistical model from measured steel set loads, geologic and construction factors. The aim of the investigation was to produce a set of charts, based on these factors, which could be used by a trained man equipped with a hand lens and measuring tape, to specify a safe and economical steel size and set spacing. In order to be of practical value the measured parameters had to be rapidly measurable and capable of quantification. It was apparent from the beginning that quantitative values for many of the geologic factors would have to be arbitrary. It was impossible to ascertain in advance which of the readily determi-nable geologic and construction factors significantly affect steel set loading. As the investigation progressed the number of observations was systematically reduced in order to determine the minimum amount of data which would be necessary for a significant mathematical model. This paper presents an analysis of all the variables measured at each instrumented steel set. The statistical significance of the measured variables was determined. On the basis of this determination, a mathematical model was constructed and charts were prepared to select and space the steel tunnel supports at any given point in the tunnel. Once the effect of geologic and construction factors on set loading is sufficiently established, it is possible to exert close control over the selection of steel tunnel support size and spacing. The magnitude of this rock mechanics program can be appreciated when it is known that there were 44 instrumentation stations for which steel set loads were measured. For each of these sets there were three construction variables recorded, and nine geologic parameters determined (Table I). A rough approximation of the cost for both the collection and reduction of the data would be $50,000. This is in addition to the cost of geologic mapping.

Jan 1, 1967

By H. H. Evinger, M. Muskat

A methoD of calculating productivity factors for oil, gas, and water systems in the steady state is presented as an illustration of the quantitative application of the fundamental data on the flow properties of heterogeneous fluid systems. While the numerical results of the calculations cannot be applied directly to field conditions, because they are based on permea-bility data for unconsolidated sands and involve specific assumptions regarding the nature of the reservoir fluids, the methods of analysis and interpretation should havc significance in considering certain specific field situations when the necessary field and laboratory data become available. Moreover, it is felt that even under the idealized conditions to which the numerical calculations refer they should serve to give at least the orders of magnitude of the effects of gas-oil ratio and connate water on observed productivity factors. Within the indicated limitations, the computations imply that the foregoing factors can explain only a part of the apparent discrepancy between the homogeneous fluid productivity factors and those obtained recently in field measurements on the West Coast. IntroductioN Sufficient laboratory research has been carried out in the last several years on the fundamental features of the flow of heterogeneous fluids through porous media to call for consideration of what can be done with this type of information. It has, of course, already led to a, much clearer general understanding of the mechanics of oil production. However, examples of the quantitative application of the results of such laboratory studies are few in the literature. part of the explanation for this fact lies undoubtedly in the lack of sufficient laboratory data relating to practical fluids and sand systems such as would occur ill actual oil reservoirs. However, unless it can be shown that if such data were obtained quantitative applications could be made, serious doubt may be cast upon the desirability of carrying on further extensive experimental research work on fluid flow. The purpose of this paper is to present one approach to the problem of quantitative application of laboratory data to practical field problems. Because the published results on the basic experiments with three-phase fluid systems in porous media are limited entirely to uncon-solidated sands, the numerical applications made here will not pertain directly to any actual field situation. However, they should serve to illustrate what could be done with fundamental laboratory data if they were extended to cover the field of consolidated sands. From this point of view, it is hoped it will serve to encourage further laboratory research on this problem. The particular application to be treated here relates to the calculation of the productivity factors for oil-gas-water systems in the steady state. This specific problem seems especially appropriate for discussion in view of the recent publication by Johnston1 of the results of produc-

Jan 1, 1942

By M. Muskat, H. H. Evinger

A methoD of calculating productivity factors for oil, gas, and water systems in the steady state is presented as an illustration of the quantitative application of the fundamental data on the flow properties of heterogeneous fluid systems. While the numerical results of the calculations cannot be applied directly to field conditions, because they are based on permea-bility data for unconsolidated sands and involve specific assumptions regarding the nature of the reservoir fluids, the methods of analysis and interpretation should havc significance in considering certain specific field situations when the necessary field and laboratory data become available. Moreover, it is felt that even under the idealized conditions to which the numerical calculations refer they should serve to give at least the orders of magnitude of the effects of gas-oil ratio and connate water on observed productivity factors. Within the indicated limitations, the computations imply that the foregoing factors can explain only a part of the apparent discrepancy between the homogeneous fluid productivity factors and those obtained recently in field measurements on the West Coast. IntroductioN Sufficient laboratory research has been carried out in the last several years on the fundamental features of the flow of heterogeneous fluids through porous media to call for consideration of what can be done with this type of information. It has, of course, already led to a, much clearer general understanding of the mechanics of oil production. However, examples of the quantitative application of the results of such laboratory studies are few in the literature. part of the explanation for this fact lies undoubtedly in the lack of sufficient laboratory data relating to practical fluids and sand systems such as would occur ill actual oil reservoirs. However, unless it can be shown that if such data were obtained quantitative applications could be made, serious doubt may be cast upon the desirability of carrying on further extensive experimental research work on fluid flow. The purpose of this paper is to present one approach to the problem of quantitative application of laboratory data to practical field problems. Because the published results on the basic experiments with three-phase fluid systems in porous media are limited entirely to uncon-solidated sands, the numerical applications made here will not pertain directly to any actual field situation. However, they should serve to illustrate what could be done with fundamental laboratory data if they were extended to cover the field of consolidated sands. From this point of view, it is hoped it will serve to encourage further laboratory research on this problem. The particular application to be treated here relates to the calculation of the productivity factors for oil-gas-water systems in the steady state. This specific problem seems especially appropriate for discussion in view of the recent publication by Johnston1 of the results of produc-

Jan 1, 1942

By R. W. Parsons

wben a displacement front encounters an isolated permeability heterogeneity, it will be perturbed. The details of this perturbation will depend on the beterogeneity size and shape, the Permeability contrast, and the mobility ratio. This paper deals with the displacement front behavior for one, simple, mathematically tractable model: a circular permeability discontinuity and unit mobility ratio. The progress of an initially planar front is shown for four permeability contrasts: k2/k1 = 0, 5/2, 2, m. The front distortion is more pronounced than streamline and isopotential distortion and it approaches a constant shape downstream from the discontinuity. INTRODUCTION The impact of reservoir heterogeneities on oil displacing processes depends upon the variation of the heterogeneous property (frequency distribution), its disposition (spatial distribution), and the inherent stability of the displacement mechanism.l The interplay of these factors plus the impracticability, if not impossibility, of describing a natural rock system in any detail makes exact displacement front movements difficult to determine. It is useful, however, to be able to qualitatively visualize what heterogeneities will do to a displacement front. In this respect certain simple models that can be treated analytically are of value. The following analysis will show the displacement front behavior for one such model. MODEL DESCRIPTION The flow field is two-dimensional and infinite in extent. The fluid flow is single phase, steady state and unidirectional at x, y = 2 m, or before insertion of the discontinuity. A circular permea- bility discontinuity is placed at the origin. Porosity and thickness are the same within and exterior to the heterogeneity. This is exactly the same model that was used by Greenkorn et al. to check analytically Streamlines as determined by Heie-Shaw flow models and finite difference computer techniques.2 FRONT POSITION CALCULATIONS The results of displacement front calculations are presented in Figs. 1 through 4 for four different permeability contrasts: k2/kl = O, %, 2, m. The times corresponding to the front locations are the Same for the 0 and cases! and for the 1/2 and 2 cases. The determination of front positions at various times is based on known analytical expressions for velocity potentials ().2 and stream functions ().2 These expressions are in terms of the x and y

Jan 1, 1969

By P. E. Hildebrand, S. W. Mao, C. N. Crain

Hydrologic and economic criteria figure in many obvious ways in water resources development, but they are rarely linked quantitatively, and most of the applications are pertinent only to the case in point. In the preparation of a regional plan for the massive program of groundwater development in the Punjab, situations have arisen where the choice among two or more alternatives was made by analyzing the economic implications of critical hydrologic parameters. Two of these situations which are of more or less general interest are described herein. The first involves basic development policy, i.e., the question of whether groundwater development should be based on the demand of the economy, or the conservation concept of safe yield. The analysis demonstrated that mining of groundwater is clearly indicated for areas where alternative supplies are unavailable. The second problem involves integrating hydrologic data into an economic model to determine optimum well design. From an analysis of all pertinent factors, a simple nomograph was derived for determining well specifications on site during construction. The nomograph can be modified for use in any alluvial terrain. In the Northern Zone of the Indus Plains —the region known as the Punjab- Pakistan is proceeding with integrated development of water and land resources on a scale never before contemplated. The heart of the development plan is a massive program of groundwater development to supplement canal irrigation supplies and to provide drainage relief to the irrigated lands. By the close of 1966 about 2500 high-capacity wells were in operation and furnishing the water supply and drainage requirements for nearly two million irrigated acres. The development plan calls for the construction of a total of nearly 30,000 wells with an aggregate capacity of 120,000 cfs by 1980. This is part of an integrated program eventually affecting over 20 million acres of irrigated lands. [This development program is administered by the West Pakistan Water and Power Development Authority (WAPDA) for which the firm of Tipton and Kalm-bach, Inc. serves as regional planning consultant for the Northern Zone. The material presented herein is derived from planning studies which have been made in connection with the preparation of a regional development plan for the Northern zone.] Pakistan is proceeding with this large-scale development of land and water resources at a time when she is at the threshold of an almost calamitous population growth. Furthermore, she is attempting this development at a time when most of her usable land is already occupied and her food demands exceed current production. Here in the Northern Zone is one of the largest continuous areas of irrigated land in the world, easily irrigated, with a warm, dry climate and ample potential supplies of surface water and groundwater. Except for an abundance of land and labor, other primary resources are lacking or in very limited supply. The agriculture sector so dominates the economy that only by means of an effective development of the agricultural resources can other sectors grow and contribute to a viable economy. Because of inadequate irrigation, farming here has not effectively participated in the twentieth century revolution in agriculture, and only by means of massive development of potential supplies can Pakistan catch up. Capital for any kind of mass development must come from the outside. But the overriding restraint to development of the land and labor resources is the availability of water. ECONOMICS OF DEVELOPMENT The decision-making process in planning regional development becomes increasingly complex as the various sectors of the economy become more highly developed. Where, as in this case, one sector dominates the economy and its development can be considered the key to development of other sectors, the decision-making process is simplified and it is easier to identify alternatives and to trace the prob-

Jan 1, 1970

By G. W. Nabor, M. Mortada

The effects of anisotropic or directional permeability on the areal sweep efficiency and the flow capacity are examined. The paper points out the importance of taking directional permeability into consideration in planning a flood. It analyzes the two-dimensional flow pattern associated with the skewed line drive for a unit mobility ratio. The direct and staggered line drives are treated as special cases of the skewed line drive. Analytical expressions are developed for the areal sweep efficiency at breakthrough and the flow capacity. They are related to the spacing between like wells, the distance between a row of injectors and the nearest row of producers, and the degree of skewness of the line drive. The latter quantity is defined such that it is equal to zero for the direct line drive and equals one-half for the staggered line drive. The areal sweep efficiency and the flow capacity depend also on the orientation of the flood pattern with respect to the principal axis of anisotropy. The paper provides a simple method for determining the areal sweep efficiency and the flow capacity for a formation in which the permeability in the bedding plane is anisotropic. INTRODUCTION Directional or anisotropic permeability is manifested by the ability of the formation to conduct fluids more readily along certain preferred directions. This situation occurs in many producing formationsl-4 and is usually attributed to depositional features in which the sand grains are oriented in a preferred direction. In some cases it results from the formation of a major and a minor fracture system. 4 Directional permeability should be taken into account in many phases of the production and exploration activities. Recognizing its existence in the formation of interest and planning accordingly can lead to increased recovery and substantial savings. For instance, the areal sweep efficiency in a water flood depends to a great extent on the orientation of the flood pattern with respect to the principal axis of permeability. Anisotropic5 permeability is specified by the directions of its three principal axes and the permeability along each axis. The principal axes of permeability are mutually perpendicular. This paper deals with the areal sweep efficiency at breakthrough and the flow capacity for formations with aniscltropic permeability. The flood pattern considered consists of alternate rows of injecting and producing wells. The rows of wells are parallel and form a developed, skewed line drive which is illustrated in Fig. 1. The staggered and direct line drives are treated as special cases of the skewed line drive. PRELLMINARY CONSIDERATIONS The formations considered in this paper are assumed to be of uniform thickness and porosity. The permeability is uniformly anisotropic, with each principal axis of permeability maintaining

By R. L. Skaggs, D. T. Peterson

The effect of carbon in solid solution on the plastic behavior of thorium was studied by measuring the flow stress of Th-C alloys from 4.2" to 573°K and at several strain rates. Carbon was found to strengthen thorium primarily by increasing the thermally activated component of the flow stress. The strengthening due to carbon was directly proportional to the carbon content and decreased rapidly with increasing temperature up to 423" K. The flow stress also increased with increasing strain rate. The strengthening appears to be due to a strong short-range interaction between carbon atoms and dislocations. A yield point was observed in the Th-C alloys which increased with increasing carbon content. JTREVIOUS study of the mechanical properties of thorium has been confined largely to the measurement of the engineering properties. Work prior to 1956 has been summarized by Milko et al.1 who reported that additions of carbon to thorium sharply increased the room-temperature strength. In addition, the yield strength was observed to decrease rapidly over the temperature range from 25" to 500°C. In 1960, Klieven-eit2 measured the flow stress of thorium containing 400 ppm C. He found that over the temperature range from 78" to 470°K the flow stress was strongly dependent on temperature and rate of deformation. A drop in the load-elongation curve, or a yield point, was observed over most of the above temperature range. Above 470°K, the flow stress was nearly independent of temperature and strain rate. This strong temperature and strain rate dependence of flow stress is not generally observed in fcc metals. It is, in fact, more typical of the behavior reported for bcc metals. Bechtold,3 Wessel,4 and conrad5 have pointed out the striking difference between the commonly studied bcc metals and fcc metals in regard to the effect of temperature and strain rate on the flow stress. Zerwekh and scott6 studied the plastic deformation of thorium reported to contain 12 ppm C. They found that this material did not obey the Cottrell-Stokes law as expected for fcc metals. In addition, they found values of the activation volume smaller by an order of magnitude than expected for an fcc metal. They concluded that thorium was strengthened by a randomly dispersed solute. Thorium differs from many other fcc metals that have been studied extensively in that it shows a relatively high carbon solubility at room temperature. Mickleson and peterson7 report the solubility limit at room temperature to be 3500 ppm C. The lowest value reported is that of Smith and Honeycombe8 who report the limit to be 2000 ppm C at 350°C. The pres- ent investigation was a systematic study of the flow stress and yield point phenomenon of thorium over a broad range of carbon content, temperature, and strain rate. EXPERIMENTAL PROCEDURE The thorium used in this investigation was produced by the reduction of thorium tetrachloride with magnesium as described by Peterson et a1.' Chemical analysis of the original ingot after arc melting and electron beam melting is shown in Table I. Alloys were prepared by arc melting this thorium with high-purity spectrographic graphite. Threaded specimens with a gage length 0.252 in. diam by 1.6 in. long were used for the constant stress or creep measurements. These specimens were machined from rod which had been cold-rolled and swaged to % in. diam. Tensile specimens were prepared by swaging annealed 3/8 -in.-diam rod to 0.102 *0.001 in. The as-swaged wire was cut to lengths of 2 in., annealed, and the center 1-in. gage length elec-tropolished to 0.100 ±0.001 in. The specimens were gripped for a length of 3 in. at each end by a serrated four-jaw collet which was tightened by a tapered compression nut. No slipping occurred in the grips and negligible deformation was observed outside the 1-in. gage length. Both the creep and tensile specimens were annealed at 730°C under a vacuum of 1 x X Torr. The resulting structures consisted of equiaxed recrystallized grains with a grain size of 3200 grains per sq mm for the tensile specimens and 2200 grains per sq mm for the creep specimens. After the specimens were prepared, samples were analyzed for nitrogen, oxygen, and hydrogen. The results of these analyses are given in Table 11.

Jan 1, 1969

By G. B. Gould

Every mathematical statement of a measure of anything (as distinguished from a count') is followed by a qualification—always implied if not explicity stated—-that the statement is only an estimate, which may vary from absolute accuracy within certain limits, according to the accuracy of the measuring instrument. This basic fact, so obvious to engineers, is not so well understood by many that make use of laboratory coal tests. The tolerances that must be allowed in measurements of coal quality are very important in both their commercial and engineering applications—and the significance of these tolerances depends upon the use made of the laboratory coal test. Only in recent years has there been any intensive study of this question of tolerances as applied to measurements of coal quality. There is very little in our published literature on the subject, and what has been published is not readily comprehended by those who sell coal, or buy coal and use it, unless they happen to have made a study of the subject themselves. There is great need for more widespread understanding of what lias been discovered and how it applies to the use of laboratory coal tests. What follows is an attempt to condense into a few pages the results of the major investigations so far reported, and to explain their meaning in simple nontechnical language, with the addition of some recent illustrative experiments. Determining Deviations in Quality It is necessary at the outset to distinguish between: (I) the accuracy with which one test indicates the quality on one lot of coal; and (2) the accuracy with which a test of one lot of coal represents the average quality of a commercial product, of which the buyer will presumably receive a series of shipments—perhaps many thousands of tons. It is necessary to make this distinction because coal is a variable material—a heterogeneous mixture of a variety of components—the quality of which will vary from one lot to another. Therefore an absolutely accurate determination of the quality (if that were possible) of one lot of a coal is an

Jan 1, 1938

By G. B. Gould

Every mathematical statement of a measure of anything (as distinguished from a count') is followed by a qualification—always implied if not explicity stated—-that the statement is only an estimate, which may vary from absolute accuracy within certain limits, according to the accuracy of the measuring instrument. This basic fact, so obvious to engineers, is not so well understood by many that make use of laboratory coal tests. The tolerances that must be allowed in measurements of coal quality are very important in both their commercial and engineering applications—and the significance of these tolerances depends upon the use made of the laboratory coal test. Only in recent years has there been any intensive study of this question of tolerances as applied to measurements of coal quality. There is very little in our published literature on the subject, and what has been published is not readily comprehended by those who sell coal, or buy coal and use it, unless they happen to have made a study of the subject themselves. There is great need for more widespread understanding of what lias been discovered and how it applies to the use of laboratory coal tests. What follows is an attempt to condense into a few pages the results of the major investigations so far reported, and to explain their meaning in simple nontechnical language, with the addition of some recent illustrative experiments. Determining Deviations in Quality It is necessary at the outset to distinguish between: (I) the accuracy with which one test indicates the quality on one lot of coal; and (2) the accuracy with which a test of one lot of coal represents the average quality of a commercial product, of which the buyer will presumably receive a series of shipments—perhaps many thousands of tons. It is necessary to make this distinction because coal is a variable material—a heterogeneous mixture of a variety of components—the quality of which will vary from one lot to another. Therefore an absolutely accurate determination of the quality (if that were possible) of one lot of a coal is an

Jan 1, 1938

By V. Carlson, J. A. Rowe, E. O. Bennett

This report concerns the microbial flora found throughout the surface facilities of six water-injection systems in Texas and Oklahoma. Each system is described in detail and water quality data are presented when available. Pseudomonads and sulfate-reducers (Desulfovibrio) were the predominate organisms found in the systems. The genera Sphaerotilus, Bacillus, Achromobacter, Minococcus, Clostridium, Flavobacteriurn, and Sarcina were also isolated in significant numbers. Molds, iron bacteria, sulfur bacteria and soil bacteria were found less frequently in the systems. It is concluded that a complex microbiological flora exists in oilfield water-injection systems. A tendency for waters to deteriorate in passage through solid bed filters is also noted. INTRODUCTION The injection of water into a subsurface formation is often accompanied by numerous problems, the origin and nature of which are only partially understood. During the past few years, investigators have indicated that bacteria are involved in some of these problems even though very little specific information is available concerning this matter. It is surprising that investigators working on microbiological problems in water injection have never made a complete study of the flora of these systems. The logical course to pursue is that of identifying the microbial flora and then examining each organism individually and in combination with others to determine if they produce a particular problem. The objective of the present investigation is to identify the microbial flora in water-injection systems. In addition. the quality of the water in most of the systems was determined and some observations made on the extent to which passage through surface-handling facilities affects water quality. The information derived from this study can serve as a foundation for additional basic investigations. EXPERIMENTAL PROCEDURE The microbial flora of six water-injection systems were studied at 50 sampling points. These systems are experiencing problems which may be attributable to micro-organisms, and two were being treated with antimicrobial agents at the time of this study. The samples were obtained by inoculating 1.0 ml of the water onto the different media at the sampling point. The media used for primary isolation included nutrient agar, Strokes agar, M-10-E medium and semi-solid thioglycollate medium. The cultures were returned to the laboratory, incubated at 30°C until growth occured; then, pure cultures were obtained from individual colonies. The pure cultures were identified according to standardbacteriologica1 procedures. Those organisms not fitting identification criteria exactly, even though exhibiting close similarities to known species, are listed as unidentified organisms. Membrane filters were used to make water-quality tests at intervals before and after collection of the microbiological samples. The filter discs were 47 mm in diameter with a 0.45-micron particle-size retention, and the tests were done under closed conditions at 20) psi. Membrane flow rates are plotted as increasing filtrate volumes, and the total solids contents arcs given when available. An effort was made to sample at each place in a system where a change in environment occurred. In this manner, the effect of passage through filters, tanks, hay sections, etc., on the microbial flora and water quality could be observed. RESULTS SYSTEM A This waterflood system (Fig. 1) is about four years old and daily pumps approx&ately 110,000 bbl of oil-free sour water from a 3,700-ft deep reef. The

By AIME AIME

THE organization of the Board of Directors was effected Tuesday afternoon in ' executive session. Messrs. Bassett,. Buehler, Cates, Daveler, Eavenson, Eilers, Elliott, Evans, Guess, Lovejoy, Moulton, Rickard, Tally, Turner,, Wraith and Wright being pres¬ent. Messrs. Tally, Bassett, Rickard, Mathews and Cates were chosen to form the Executive Committe with Mr. Tally as Chairman, and Messrs. Krumb, Guess and Moulton for the Finance Committee with Mr. Krumb as chairman. Karl Eilers was reelected as Treasurer and H. Foster Bain as Secretary, after which the Directors adjourned for the usual dinner reported elsewhere. At the open session of the Board, at which in addition Messrs. Allen and McAuliffe were present, the work of organization was continued with appointments to committees as listed herewith. Various actions of the Executive and Finance Committees were confirmed, the budget for 1931 formally approved, the meeting dates for the year adopted, 34 new Members and Associates, 13 Junior Members and 31 Student Associates admitted by election and transfer and 29 resignations accepted.

Jan 1, 1931

By W. F. Kieschnick, and W. J. McGuire&apos, Eugene Harrison

This paper concerns the induction and extension of fractures into rock formations as involved in drilling, completing, and production stimulating operations on wells. Conclusions concerning formation breakdown are derived from (I) a review and extension of published analyses relating to mechanical theories of rock stress and the state of stress in the earth's crust and (2) a correlation of field data from fracturing operations. Conclusions concerning the mechanics of fracture extension, which indicate the relationship between fracture dimensions and rock properties, depth, and volume of injected fluid, are tentative and largely establish limits of relationships. These conclusions are derived from stress calculations, limited field data, and laboratory experimental studies. The experimental work involves the study of the stresses at the fracture boundaries and the geometry of pressurized fractures by means of photo-elastic modeling methods. Results of this investigation indicate that a large majority of pressure induced wellbore fractures are vertical, particularly in deeper wells; and variations in the pressures necessary to create and extend fractures can be explained largely on a basis of established rock properties. It is also shown that variations due to tectonic forces should usually be expected to be slight. Other results indicate that during that extension of fractures rather large fracture volumes are temporarily created by the parting of the formation. The purpose of this paper is to present the results of calculations and laboratory experiments concerning the mechanics of fracture induction and extension with a view to broadening existing knowledge relating to these phenomena. It is believed that continued progress in developing knowledge of this type is important to the further development of techniques for drilling and completing wells. When a wellbore is subjected to sufficiently high pressure, fracturing of the formation surrounding the well occurs. This phenomenon is a common cause of lost circulation during drilling operations, is a deliberately planned factor in most squeeze cementing operations, and is the controlling factor in many productivity stimulation operations. The exploitation of formation fracturing in recent years has been particularly impressive in sand-viscous fluid well stimulation treatments which were first introduced to the industry by the Stanolind research group1,2,3 several years ago. The borehole fracturing phenomenon has been analyzed mathematically by application of the principles defining the elastic and inelastic behavior of thick-walled cylinders of homogeneous, isotropic, and impermeable material.1,5,6 This basic argument is re-presented in this paper along with new considerations which strengthen it. In addition to these arguments relating to fracture induction some calculations and experimental results relating to the mechanics of fracture extension are presented. Since the nature of the pressure induced fracture at distances from the well is of great importance in considerations of well stimulation procedures there appears to be a need for an integrated approach which deals with the mechanics of fracture extension as well as with wellbore rupture.

Jan 1, 1955

By R. E. Boehler, N. C. Ludwig

At Buffington Station, Gary, Ind., Universal Atlas Cement operates fourteen 8 x 101/2 x 155-ft cement kilns in mill 6 and two 11 x 360-ft kilns in the Harbor plant. The No. 11 and 12 kilns in mill 6 are equipped with Manitowac recuperator sections. This report describes studies in measuring exterior shell temperatures on several of these kilns and the development of a traveling radiation pyrometer with certain novel features. Preliminary Work: At first various temperature-sensing devices were placed on the steel shell: 1) crayons with calibrated melting points, 2) colored paints with temperature-calibrated pigments, 3) aluminum paints with temperature-calibrated binders, and 4) metal-stem dial thermometers. The colored paints and aluminum paints failed to indicate the temperatures correctly. The crayons and thermometers did indicate fairly correct temperatures, but it proved impossible to apply enough of these on the shell to detect all the potential areas where hot spots might develop. Furthermore, considerable labor was required to apply these sensors and read the temperatures. Consequently no further work was done with these devices. Formation of Hot Spots: In the burning or clinker-ing zone of a cement kiln, the thickness of the protective coating and thickness of the brick govern the amount of heat transmitted to the steel kiln shell. Usually the protective coating consists of 4 to 8 in. of fused cement clinker. The formation of a hot spot is usually caused by loss of coating? that is, localized areas of the coating become thin or fall away from the refractory. This is generally caused by excessive temperature in the burning zone over a fairly long period of time. It may also be caused by a sudden thermal change in the burning zone. Variations in raw feed composition and in feed rate require changes in the fuel and air rates, and when these are not appropriately altered, conditions may develop in the kiln that will result in loss of coating. Luminescence on the kiln shell indicates that a hot spot has developed to a point that usually alters the refractory's thermal conductivity properties. When this thermal weakness zone occurs in the burning zone of the kiln, constant vigilance is required to protect it by maintaining proper coating. Even so, it has been the writers' experience that within a period of several days to about four weeks the hot spot usually recurs with greater severity. This necessitates shutting down the kiln and re-bricking the affected area. One of the prerequisites of a good burnerman is the ability to maintain a protective coating despite the many variables in operation. When he knows that it is getting thin or that an area has dropped off, he reduces the firing rate and kiln speed and brings feed into the affected area in an effort to rebuild the coating. But when powdered fuel is burned, the atmosphere of the kiln may prevent the burnerman's observing the condition of the coating closely at all times without taking off the fire. It is not considered good practice to do this frequently, as it imposes a thermal shock on the coating and upsets operation of the kiln. To help the burnerman scan the shell of the kiln along the burning zone, therefore, a radiation pyrometer, connected to a potentiometric recorder, was mounted on a slowly moving steel cable. The theory of operation, construction details, and adaptability of the radiation pyrometer are included in an excellent monograph' and also in a textbook.' Shell temperatures of the Atlas Cement kilns were measured with a Brown Instruments Div. low intermediate range Radiamatic unit, of range 200" to 1200°F, and a circular chart Electronik potentio-metric recorder, of range 500" to 1000°F. In Bulletin 59095M the supplier publishes standard calibration data (millivolts vs degrees Fahrenheit) for this radiation pyrometer, These data, however, apply only to flat surfaces having emissivities of unity. Calibration of Radiation Pyrometer for Use on Curved Surfaces: When applied to surface temperature measurements, the radiation pyrometer reading depends on the nature of the surface, the material of which it is composed, and also to some extent on the temperature of the surroundings. Although the present radiation pyrometer is designed to give a calibrated response under ideal (black body) conditions when used commercially, it must be calibrated empirically. The calibration procedure, given below, follows that described by Dike (Ref. 1, pp. 38-39). Calibration tests on plane and curved surfaces showed that the response of the radiation pyrometer was very

Jan 1, 1961

By N. De Nevers

A calculation method has been developed for carbonated water flooding. This method takes into account the effects of oil viscosity reduction and oil swelling, due to carbon dioxide transferred to the reservoir oil from the carbonated flood water. It shows the effect of changes in carbonation pressure and carbonated water slug size. The method is based on a Buckley-Leverett-type linear flow model. The mathematical approach is similar to that developed by Welge, et al, for enriched gas drives. However, it contains a very important advance over the previous papers in its solution to the problem of injection of a "slug" of carbonated water rather than continuous injection of fluid of constant composition. The method also shows the chromatographic nature of the transport of Co 2 in a carbonated water flood. Sample calculations indicate that viscosity reduction is the most important source of incremental oil recovery. Swelling contributes to a lesser extent. Raising the carbonation pressure increases recovery, but this effect diminishes as pressure is increased. increasing the slug size also increases recovery, but again this effect diminishes as slug size is increased. INTRODUCTION In recent years there has been considerable interest in carbonated water flooding as a secondary recovery method. Several carbonated waterflood field tests are now in progress.1'3 Yet there is no published method for predicting the results of a carbonated water flood. This paper presents a method for predicting the oil recovery for ; carbonated water flood. In a carbonated water flood. a carbonated water slug is injected into the reservoir. Then this slug is driven forward by injecting plain water behind it. The field operator must choose the carbonated water slug size and the CO2:water ratio in the slug. The calculation method presented here shows the effects of both slug size and CO2:water ratio on oil recovery. The calculation assumes that the addition of carbon dioxide to flood water increases the oil recovery by two mechanisms: 1. Carbon dioxide transfers from the flood water to the oil, thereby reducing the viscosity of the oil. At the lower viscosity, the oil flows more readily and is more easily swept out by flood water. 2. The carbon dioxide which transfers to the oil phase causes the oil phase to swell. Assume that the residual oil saturation (volume per cent of pore space) is the same whether the oil is carbonated or not. Then the same volume will contain less stock tank oil if it is carbonated, due to oil swelling. The plain water driving the carbonated water slug forward will strip the CO2 out of this residual oil., causing this oil to shrink. This swelling-followed-by-shrinking leaves behind a lower residual oil saturation than would have been left behind by a plain water flood. Since some oils swell to almost twice their original volume at reasonable CO2 pressures, this effect may be large. Two additional benefits sometimes claimed for carbonated water flooding are: 1. Adding CO2 to flood water greatly increases the water injectivity in some cases.1 (Injectivity = BPD of water injected/ft of sand x differential bottom-hole pressure, psi.) CO2 forms a weak acid in water; injecting carbonated water is, therefore, somewhat like acid treating. If a formation responds to acid treating, it will probably have a higher injectivity for carbonated water than for plain water. Although the increased injectivity increases the profitability of a carbonated water flood over a plain water flood, it does not change the production curve based on total fluid injected. Therefore, the method presented here does not show any effect of increased injectivity. 2. It is claimed in one patent4 that carbon dioxide will react with some constituents of the oil to form detergents. These detergents are expected to increase the oil recovery. There is little data available to support this idea. The calculation method presented here takes no credit for extra oil recovery due to the formation of any such detergents.

Jan 1, 1965

By M. A. Torcaso, P. Raimondi

To study mass transport in systems simulating oil recovery processes, different porous media were saturated with a mobile (carrier phase) and a stationary phase. Slugs of carrier phase containing a small amount of solute were displaced with pure carrier phase. By analogy to the chromato-graphic processes, the velocity of the solute can be predicted from a knowledge of the partition coefficient and the saturation provided that equilibrium between the two phases exists. Equilibrium was found to exist for different porous media, solutes and rates. The conditions were varied over the range normally encountered in the laboratory and in the field. The longitudinal dispersion of a solute undergoing interphase mass transfer was also investigated. INTRODUCTION The production of hydrocarbons by gas cycling, enriched gas drive and C02 or alcohol displacement involves, among other factors, relative motion between two phases and compounds, hereafter called solute, which are soluble in both phases. The solute is carried forward by the faster flowing phase at a lower velocity than the average velocity of that phase. Retardation of the solute is caused by chromatographic absorption and desorption in the slower flowing phase and by the degree of departure from equilibrium. At equilibrium the concentration of solute in the two phases can be related by the equation* CSw = K Cso ,............(1) where C, and C are the concentration of solute in the aqueous and oleic phases respectively and K is the equilibrium ratio, or partition coefficient. Displacement theories must contain an explicit assumption with regard to equilibrium, i.e., whether the compositions can be related by Eq. 1. The existance of equilibrium depends, in general, on the relative velocity between the phases. Unfortunately, other factors such as gravity segregation and viscous fingering, also depend on velocity. For this reason, whenever effects of rate on displacement were observed,2-4 it was practically impossible to discern what caused them — lack of equilibrium or the factors mentioned above. Equilibrium between phases has been the subject of extensive studies in fields such as extraction or chroma to graphy. It has received only small attention in flow through the type of porous media encountered in oil production.5 For this reason a method was developed which makes it possible to study the movement of a solute as it is affected by rate, type of porous media, partition coefficient and carrier phase, but in the absence of segregation or fingering. The information obtained enables one to determine when the assumption of equilibrium can be made. Briefly, the method consists of (1) saturating the core with a mobile and an immobile phase, (2) injecting a slug made up of the same fluid as the mobile phase and a small concentration of mutually soluble solute, (3) measuring the lag and the peak height of the slug at arrival and (4) correlating these variables with fluid properties such as partition coefficient and mixing constants of the medium. PROPOSED MECHANISM The principles of chromatography are combined with the equation of longitudinal mixing to predict the velocity of a solute slug compared to the bulk velocity and the peak height of a slug.6-9 The equation so obtained is valid under equilibrium conditions only. Therefore, a comparison between experimental and predicted results will give a measure of departure from equilibrium. This work was done with either the oleic or the aqueous phase being immobile. For simplicity, the following development is based on the case where the oleic phase is immobile. However, the treatment is the same in either case.

Jan 1, 1966