Search Documents

By Rossiter W. Raymond

The most important event in the history of mining in the United States was the discovery of gold in California, which led to the rapid development, not only of a new industry, but of a new empire. The settlement of the Pacific slope was making but slow progress under the inducements of a scanty commerce in furs and hides. The liberal grants offered by the government to agricultural pioneers in Oregon had attracted to the Willamette valley and other favored points in that territory a limited population; hot the government of the Mexicans in California was positively unfavorable to enterprise. It is true that gold was known to exist in the country. Rumors of it were current among the trappers, and placers had actually been worked in what is now Los Angeles County; but the padres of the neighboring mission of San Fernando had discouraged the business as demoralizing to their flocks; and the existence of the precious metal was not widely known; still less was its distribution northward suspected. The discovery of gold, January 19th) 1848, at Sutter's saw-mill, on the south fork of the American River, afterward famous as Coloma, in Eldorado County, was the beginning of the new epoch. This discovery was made by James Wilson Marshall, a partner in the mill, who observed in the debris washed down by the tail-race a glittering fragment of a mineral unknown to him. Suspecting it to be gold, he sought for further specimens, and obtained, in the course of a few days, several ounces of it. Tests subsequently

By J. C. Smock

The magnetic and polaric properties of magnetite, or magnetic iron ore: are fundamental facts in magnetism. The disturbing erect of this mineral upon the magnettic needle in land surveying must have been very early observed. The difficulties arising from this source are still almost insurmountable in making accurate surveys in districts where this ore occurs to any extent, either as a rock

The opening session* was held in the hall of the Franklin Institute, on Tuesday evening, June 20th, President Holley in the chair. The President introduced Mr. Franklin B. Gowen, who addressed the Institute as follows: Mr. President and Gentlemen : It has been made my pleasant duty to welcome you to the city of Philadelphia, but in so doing I shrink, with some natural reluctance, from appearing before an audience so much better fitted than I am to cope with the subjects which are usually presented to your association, and I can assure you that no trembling neophyte, first venturing within the groves of the Academy, or walking, with sacred reverence, in the gardens of Plato, felt more natural dread at the consequence of his own temerity than I do in appearing before an audience composed of geutlenlen belonging to a scientific society of which I am but too sensible that until now I have been a member only in name. In olden times, upon entering a temple sacred to the worship of any deity, it was customary to lay some offering upon the altar, and, in following this custom, it will be my duty, before entering upon the more pleasant part of the task which is assigned me, to call your attention to a subject with which I have been somewhat familiar, and which is entirely within the scope of your investigations, in the belief that I may throw out some hints or some suggestions of a process that may be fully developed by some of you, and eventually lead to the introduction of that which may be of service to the human race. The satirist who described the philosopher whom he met, as engaged in the process of extracting sunbeams from cucumbers,

By John B. Pearse

The works heretofore used in carrying out the Bessemer process have been constructed substantially as follows: The whole works or plant has been divided into three parts. 1. The division in which the metals for conversion are melted. 2. The division in which the metals are converted into steel and cast into ingots, or the converting department. 3. The engine and boiler department which supplies the hydraulic power required to move the machinery, and which furnishes the blast used for converting the crude metal into steel.

By James F. Lewis

The hematite iron ore mines east of the Hudson River are confined to a strip of country ten to fifteen miles wide, commencing on the south, near Fishkill, running northeast through Dutchess County, and striking the southeast corner of Columbia County, extending into Litchfield County, Conn., known there as the Salisbury mines, cropping out again at and near West Stockbridge, Lanesboro, and Cheshire, Mass., also in Bennington County, Vermont. The ores are found in beds or deposits and in veins, with a limestone ledge on one side and slate on the other, from six to thirty feet from the surface. Portions of it are quite solid, called rock ore, but

By L. G. Jones

The objective of this paper is to present a better method for determining reservoir reserves of gas from extended drawdown test data. Means are proposed for conducting the test with or without the requirement of prior build-up to a static reservoir pressure. The test is based upon solutions to unsteady-state gas flow developed in an earlier paper and is believed to be superior to previous methods for determining reserves from flow-test data. The limitatiom of such tests are discussed, with emphasis on the subject of test duration. The test can supply useful information if the mechanics and limitations of it are understood. This type of information can be particularly valuable in deciding whether or not to drill offsets to producing wells. INTRODUCTION In the life of any gas or oil field, the engineer in charge must decide if certain wells need stimulation and whether or not to drill more wells in the field. Well production tests can aid in making these decisions. For instance, the value of stimulation can be estimated by determining formation permeability and the amount of well damage from drawdown or build-up tests. These tests are freely used in the oil industry. A less well-known type of analysis uses data from extended drawdown tests for determining reservoir reserves. Into this category fits the "reservoir limit test" of Jones,'-3 and the "reservoir reserve test" presented in this paper. The chief value of these tests is their ability to "see" reservoir boundaries and thus prevent costly dry holes. Methods for determining reserves of both gas and oil reservoirs were proposed in the work of Jones. His proposals appear sound for single-phase liquid flow. The generality of his work on gas flow is open to question, however, since the analysis is based upon solutions to liquid-flow equations. The test for gas reserves presented here is derived from unsteady-state gas-flow solutions presented in a previous paper.' In deriving this reservoir reserve test, it has been assumed that all of the reservoir energy is provided by expansion of the reservoir fluid and that the reservoir is radially symmetrical. No mention has been made of the effects of deviations from Darcy's law such as quadratic (turbulent) flow and well damage. Since the terms representing these effects are function of rate only, since the rate is constant and since the data are plotted in a differential form, it seems reasonable that these terms will not affect the test results in any way. No attempt has been made to estimate the effects of water drive, faults, or other departures from an idealized and isolated reservoir. For a qualitative discussion of some of these effects, the reader is referred to a paper by Jones3

By Harold L. Overton, Leonard B. Lipson

The first paper in this series1 outlined practical methods for applying the theory of steady-state flow of an ideal Bingham plastic liquid through a circular pipe and axially through a stationary concentric annulus to the enginekring analysis of friction loss problems. In the present paper the properties of another useful rheological model, the pseudoplastic generalized Newtonian liquid, are considered. The terminology applied to the rheological models is derived from the following classification of rheological models. The term "generalized Newtonian" is a class designation applied to models which do not have yield points but which exhibit a dependence on shear rate. The term "pseudoplastic" applies to a sub-group of models within this class which exhibit a decrease in "viscosity" with increasing shear rate. The term "power model" refers to a mathematical model proposed for describing the behavior of the pseudoplastic liquid. This model is of interest since certain of the salt saturated, oil emulsion, and "low solids" drilling fluids, inverted emulsion- and oil-base drilling fluids, aqueous gels, gelled crudes and blocking agents employed in hydraulic fracturing operations are of this type. The present treatment considers the equations describing steady-state Poiseuille flow through a circular pipe and a stationary concentric annul us, and Couette flow between concentric rotating cylinders for a particular pseudoplastic model. As in the previous paper' it is demonstrated how these equations can be applied to engineering calculations of friction losses. Graphical procedures for recovering the various flow curves and cxample calculations are included. A method for simplifying turbulent-flow correlations is also proposed. PROPERTIES OF PSEUDOPLASTIC GENERALIZED NEWTONIAN LIQUID In contrast to the Bingham plastic the generalized Newtonian liquid is yield point-free. However, in the flow region the relation between stress and rate of shear is not linear as it is with ordinary Newtonian fluids. If the stress increases with shear rate at a less than linear rate the system is called pseudoplastic. While there is no yield point in this system, the Bingham advocates' coined the term pseudoplastic to distinguish this kind of behavior from Bingham plastic behavior. For the pseudoplastic a plot of flow rate vs friction loss or rpm vs torque, from pipe flow or rotational viscometer measurements, respectively, generally has three qualities: (1) it can be clearly extrapolated toward the origin; (2) there is very little slope at the origin. This accounts for the high apparent viscosities at low mixing and pumping rates; and (3) the flow curve possesses a continuous curvature which gradually diminishes as the shear increases. Rheological measurements on a pseudoplastic liquid can be easily misinterpreted as indicating Bingham plastic behavior. A typical case in point is illustrated in Fig. 1. These data were obtained from measurements on an oil-emulsion mud using a commercially available rotational-type viscometer (Model 35 Fann V-G meter). Applying the "two-point"* method3 gives a plastic viscosity of 17 cp and a yield point of 24 lb/100 ft2; at the same time weighting material settles very easily in this system. This illustrates that it is not possible to distinguish between the pseudoplastic and the Bingham plastic from only two data points. It is possible to determine qualitatively which model is more appropriate by a combination of dynamic and static measurements. Experience indicates that the simplest criterion for the Bingham plastic liquid is that the calculated dynamic yield point should be of the same magnitude as the measured gel strengths. For the oil-emulsion drilling fluid illustrated in Fig. 1, a value of 4 lb/100 ft2 was obtained for gel strengths measured

By L. L. Handy

The most common method of identifying hydrocarbon-bearing strata in a well that penetrates many different formations involves measurement and interpretation of the electrical properties of the formations as determined by electrical logs. Even though this method is used extensively, and even though in a great many instances it is capable of indicating presence of oil or gas, situations arise for which it is extremely difficult, if not impossible, to deduce the presence of hydrocarbons. These situations may involve the following. 1. A thin formation, bounded by highly resistive formations, in which it is impossible to obtain the actual resistivity of the uninvaded zone with existing logging devices. 2. A formation in which invasion has been so extensive that a value for the uninvaded zone resistivity cannot be obtained. 3. A very shaly formation in which the resistivity index, I, is lower than that usually associated with productive formations. 4. Laminated formations comprised of thin productive sands separated by thin shale streaks in which the individual sand and shale streaks are too thin to permit measurement of uninvaded-zone resistivity with existing logging devices. 5. Productive formations in which the water saturation is high. To extend the utility of electric log interpretation to identification of hydrocarbons in all types of formations, there is strong incentive to find a method not subject to these limitations. Some time ago, in connection with research on the wettability of reservoir rock, an investigation was conducted in which the resistivities of cores were measured shortly after they were removed from a core barrel, and again after they had been extracted and restored to their original oil and brine saturation.' The resistivities after extraction were generally lower. Other tests made on the cores indicated that they were more nearly water wet after they were extracted; thus, it was assumed that the observed changes in resistivities were due to a change in wettability of the cores. Other experiments'," have shown that resistivities of rock samples are sharply dependent on wettability. These experiments have shown that oil-wet samples are more resistive than water-wet samples. To obtain an understanding of how the wetting properties of the surfaces of core material affect electrical resistivity, a series of experiments was conducted. Two groups of core samples were prepared for testing. One group contained brine, but no residual oil. The other group was saturated with brine, flooded with oil to a low water saturation, then flooded with brine to a final residual oil saturation. Resistivity measurements were made on each group. Both groups were then flooded with the original brine to which a chemical had been added that renders sand and clay surfaces preferentially oil wet, a so-called reverse-wetting agent. very little change in resistivity was observed in cores containing only water. The group containing residual oil, however, showed resistivity increases of 100 to 200 per cent. These experiments showed that the resistivity of a core containing oil could be altered by changing wet-tability of the core. Moreover, the possibility was introduced that reverse-wetting agents might be employed as the basis for a logging method for identification of oil-bearing strata. Since behavior of a porous rock containing gas and water might be expected to be similar to that of a rock containing oil and water, such a method should also be applicable to identification of gas-bearing zones. In principle the wettability of the invaded zone could be reversed without altering conductivity of the interstitial water or the hydrocarbon saturation therein. Those strata showing significantly increased invaded-zone resistivities would, therefore, contain hydrocarbons; those with no significant change would be filled only with water. Addition of a reverse-wetting agent to a hydro carbon -bearing zone which is, by nature, already preferentially oil wet would not result in an enhancement of its resistivity. It is generally believed, however, that most hydrocarbon-bearing strata are preferentially water wet.

By F. W. Jessen, N. Lamont

The effect of an existing fracture or joint plane, which may exist in a rock, on the extension of a hydraulically induced fracture through the rock has been investigated in the laboratory. By use of a series of models made from various outcrop rocks, the width and orientation of existing fractures do not alter the extension or direction of the hydraulic fracture. The results, which conform quite well with the Grifith theory of failure, are illustrated by a series of photographs of the various models. INTRODUCTION The process of hydraulic fracturing has been widely used in the oil industry since its introduction in 1948 and has made possible the production of many reservoirs which would have been uneconomical prior to this process. Numerous studies of the mechanics of the process2-h nd of the effects of oriented fractures on recovery'.' have appeared in the literature. The increased recoveries led to attempts to predict the orientation of hydraulic fractures at the wellbore and to the development of methods to control this orientation. Field tests have been devised"." which indicate the validity of theoretical predictions of fracture orientation at the borehole. The direction of extension of hydraulic fractures from the borehole has not received much attention since most studies have predicted, at least implicitly, that a fracture will continue in the same plane unless a change of state of stress in the rock occurs. Since such predictions are based on the assumption of rock homogeneity, the effects of rock heterogeneities are left unknown. Most sedimentary rocks are composed of layers which reflect the changing depositional conditions of geological time. In addition, the more competent rocks frequently are fractured and jointed as the result of structural deformations or tectonic movements. Evidence of such joint systems exists in many surface outcrops, and it can be assumed that similar systems occur in many subsurface rocks, although the individual joints may be rather tightly closed due to overburden forces. This paper investigates the effect that an existing fracture in a rock would have on the direction or orientation of an extending or advancing hydraulic fracture when it intersects the existing fracture. The study was conducted on small rock models under triaxial stress conditions. The results indicate that an existing fracture will have little effect on the hydraulic fracture. THE MODEL STUDY The rock model represents an elemental portion of the earth containing an existing fracture which is located at some distance from a borehole from which hydraulic fracture is being extended. A generalized view of a model is presented in Fig. 1. The models were constructed from cement blocks and natural rocks. The various materials tested and some of the physical properties are listed in Table 1. The rocks were cut into rectangular blocks with dimensions of 11/2 X3Y2 X4 to 8 in. An initial slot was cut into the end face of each block along the longitudinal axis normal to the top face. The slot was filled with plastic aluminum along its outer edge to a depth of 1/16 in. A hole of 3/8-in. diameter was drilled into the sealed slot at the center of the end face to serve as an entry port for the fracturing fluid. Existing fractures of various types were simulated in the models. A "hairline" fracture with essentially zero width was created by cutting the model into two parts with a diamond saw and replacing the two parts together along the cut. A finite-width fracture was created by placing a layer of sand grains between the two faces of the cut just described (see Fig. 12). A large open fracture (similar, perhaps, to a vug) was created by cutting the model in two parts and removing a portion of the downstream block to a depth of % in. See Fig. 13. For a fourth

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system's behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell' have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine' have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.' A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By Persifor Frazer

There is always a difficulty in assigning limits to the scope of a science or art, and among the dificulties of that kind besetting practical or stratigraphical geology is that of determining what kinds of natural objects shall be considered as the goals of its search. The term mineral, in its broader sense, includes every solid, liquid, or

By C. H. Mathewson

MICROSCOPIC rnetallography has been exploited quite well enough to bring about a very general understanding that the typical metal or alloy is composed of minute crystalline particles blended into a coherent microstructural mosaic. One does not have to be a specialist in metallography to realize that the properties of such an aggregate are essentially a summation in appropriate form of individual effects derived from the shape, size, placement or orientation and cohesive characteristics of the component particles. It is clearly of great importance to consider the various forms of discontinuity which may occur at the boundaries between these crystalline particles. When a number of particles each possessing the same orderly arrangement of atoms are brought together into a close-fitting system of purely haphazard contacts, not unlike a handful of snowflakes compacted into a snowball, the particles are said to possess random orientation. There are no generally accepted views concerning the arrangement of the atomis or the constitutional reaction between atoms where one crystal meets another, but these contact regions are characterized by strength rather than weakness and it is customary to require the presence of many rather than few boundaries in preparing metal for useful service. Under certain circumstances all of the particles in a metal may be nearly alike in orientation. Other general tendencies of orientation may be associated with particular forms of mechanical and thermal treatment. In contrast with these cases of fortuitous or statistical diversity of orientation, we often find adjacent particles united along a plane which possesses a grouping of atoms belonging equally well to both structures. This, of course, determines a fixed relationship between the two orientations and it is always possible to derive one from the other by some form of rotation or reflection prescribed by the symmetry of the crystal structure under consideration. Particles united in this manner are known as twin crystals although the term refers to the form of association rather than the number of individuals concerned.

By R. W. Snyder, H. J. Ramey

This paper presents the results of applying the Buckley-Leverett' displacement theory to petroleum reservoirs consisting of a finite number of layers. The layers are assumed to communicate only in the wellbores, and the reservoir may be represented as a linear system. Most previous investigations of this nature were limited by assumptions and by inconsistent calculation techniques. This study improves on previous work by applying the Buckley-Lev-erett displacement theory to a noncommunicating layered system where permeability, porosity, initial saturation, residual saturation and relative permeability vary from layer to layer in a logical and consistent manner. Gravity and capillary-pressure effects are neglected. A modification of the Higgins-Leighton calculation method was used in this study. Waterflood predictions were made with all properties varying, and then with only permeability varying using several inability ratios. These results were compared with the Stiles and Dykstra-Parsons predictions. It is shown that the latter methods generally give poor values for the breakthrough recovery and pessimistic predictions for the performance after breakthrough. Similar results were obtained for a gas-displacement case. lNTRODUCTION Field experience with immiscible displacement usually shows constant producing conditions until breakthrough of the displacing fluid. Then oil production continues at increasing displacing-to-displaced fluid ratios until the economic limit is reached. Three different ideal mechanisms are known that will produce this behavior: (1) relative permeability effects as described by Buckley-Leverett frontal advance theory,' (2) vertical stratification as considered by Stiles,2 Dykstra and Parsons5 and others and (3) different path lengths involved in areal (two-dimensional) flow between wells as described by Dyes et al.4 Without question, a combination of these factors modified by formation heterogeneity and other known and unknown factors actually does control the behavior of real systems. This paper presents results of an investigation of certain factors that should affect performance but which have received little attention to date. In 1944, Law5 demonstrated that porosity and perme- ability are often found to have normal and logarithmic-normal distributions, respectively. throughout cored intervals in natural formations. This led to the concept of the noncommunicating, multilayered reservoir model for immiscible displacement. This model assumes that the reservoir is composed of a number of layers that communicate only at the wellbores. Each layer is individually homogeneous, but may be different from every other layer. Stiles' presented one of the earliest applications of this model to waterflood performance. In addition, Stiles assumed that the initial saturations and relative permeabilities were the same for each layer, porosity was the same. displacement was piston-like, fluids were incompressible and injection into each layer was proportional to that layer's permeability capacity (permeability-thickness product). The last assumption would be true if the mobility ratio for the displacement were unity.21 Dykstra and Parsons" used the same model as Stiles, but rigorously included mobility ratios other than unity for piston-like displacement. Dykstra and Parsons used their general result to produce charts for log-normal permeability distributions between layers. Similarly, Muskat6 Pub1ished analytical solutions for linear and exponential permeability distributions. In 1959, Roberts' described a scheme for calculating water-drive performance for the noncommunicating, layered reservoir model which considered two-phase flow in the displaced region. Roberts used the same model and assumed that the injection rate into a layer was proportional to that layer's permeability capacity, but that flood front locations could be evaluated from the Dykstra-Parsons results. These assumptions are inconsistent, and a material balance cannot be maintained except for a mobility ratio of unity. At the same time, Kufus and Lynch8 coupled Buckley-Leverett displacement theory with the layered model to provide an improvement of the Dykstra-Parsons method that was consistent. In 1960, Higgins and Leighton9 resented a numerical method for calculating waterflood performance also considering two-phase flow in the displaced region. The result was used to investigate variation in absolute permeability and oil viscosity. An excellent, detailed history of using the noncomrnunicating, layered reservoir model was presented by Nielsen.'" The preceding techniques (and many related ones) were similar in that differences in initial saturations, residual saturations and relative permeabilities from layer to layer were neglected. It is well known that the irreducible water saturation is an important function of absolute permeability. Calhoun11 showed that the irreducible water saturation

By A. M. Skov, L. F. Elkins

Inclusion of anisotropic permeability in mathematical analysis of pressure transients observed during development of the huge Spraberry field indicates a major fracture trend which is in good agreement with that observed by fluid-injection tests spread over a 12- by 17-mile area. Delineation of this trend is important in selecting a pattern of injection for the pending large-scale water flooding in this field. Determination of reservoir parameters yielding best agreement between calculated pressures and observed reservoir pressures in newly completed wells was made using an IBM 650 computer. INTRODUCTION The Spraberry field covering 400,000 acres is a tight sand of less than 1-md permeability cut by an extensive system of vertical fractures. Primary recovery dominated by capillary retention of oil in the fractured sand matrix blocks is less than 10 per cent of oil in place. Strong forces of capillary imbibition of water into the sand, coupled with water flow under dynamic pressure gradient, indicate considerable increase in oil recovery can be achieved through water flooding. Best results will occur if the pattern of water injection is selected to force the water flow across the grain of the major fracture system. Existence of an oriented vertical fracture system in the Spraberry, observed first in cores, was highlighted more recently by the 144-fold contrast in permeability along and at right angles to the major fracture trend required to match relative water breakthrough times in Humble Oil & Refining Co.'s waterflood test there. Spraberry operators since have conducted two gas-injection tracer tests for further areal confirmation of the fracture trend. Re-analysis of early reservoir pressure transients for evidence of anisotropic permeability has permitted many more local determinations of major fracture trend without resort to further field tests. This paper is limited to updating analysis of reservoir pressure transients to include anisotropic permeability as a test for orientation of the major fracture trend in the Spraberry. The reader is referred to Refs. 1 and 2 for information about general Spraberry reservoir performance and to Refs. 3 and 4 for information about significance of fracture orientation in selection of the injection-well pattern for water flooding the Spraberry. RESERVOIR PRESSURE DATA—DRIVER AREA During early development of the Spraberry Driver area, Sohio Petroleum Co. made the extra effort to measure the initial pressure in each of the 71 wells in a 5-mile-long area immediately after completion. Progressively greater reductions in pressure ranging up to 400 psi were observed throughout the six-month development period. Detailed data are presented in Ref. 1. Since the reservoir oil was undersaturated some 300 psi initially, early reservoir performance involving 55 new well pressures is subject to analysis as flow of a single compressible fluid in a porous media. Assumption of uniform permeability in all directions yielded good agreement between calculated pressures and observed pressures of these wells in the earlier study,' but subsequent, additional, mathematical development to include anisotropic permeability in the transient pressure considerations and present availability of electronic computers to perform the much more extensive arithmetical calculations now yield even better agreement. The previous analysis, assuming uniform permeability, consisted essentially of calculating pressure reduction expanding circularly around each producing well and summing these effects at the time and location of each newly completed well for comparison with the measured pressure reduction. Permeability, effective fluid and rock compressibility, and permeability X thickness were varied until the best match with measured pressures was obtained. The present analysis, assuming anisotropic

By W. S. Ayres

The boiler from which these stay-bolts have just been obtained was that of the locomotive Catasauqtia, Lehigh Valley Railroad, built at the company's shops, South Easton, Pa., in 1864. The iron is Lowmoor, and has been in use nearly ten years. The stay-bolts have become highly crystallized from some cause, and are very brittle, only requiring three blows of an ordinary machinist's hammer to break them. On the other hand, when annealed they are remarkably tough, as shown by the doubled sample. The majority of the broken bolts were found to be along the curve of the shell from a to c, Fig. 6. All of them were broken just inside the shell, and in no case next the fire-box. The fractures seem to have been of slow procedure, beginning at the top of the bolt, gradually working deeper and deeper until the pressure of steam overcame the tenacity of the remaining parts. Unequal contraction and expansion of the parts of the boiler seems in this case to be the only existing cause for a fracture. We cannot, from the data obtained, prove to an absolute certainty the manner in which the combined forces generated by the unequal contraction and expansion of the parts have arrived at their resultant work ; but we can, from these data and from existing circumstances, arrive at what seems to be a fair supposition, and thus at least anticipate the truth. The contraction and expansion in locomotive boilers has a large range between its extremes, and is very variable, both with reference to time and to their component parts. The shell and fire-box do not, perhaps, at any time, contract and expand at the same rate, nor proportionally; but, on the contrary, owing to the great changes in temperature that take place both in the atmosphere around the shell, and in the fire in the fire-box, their condition is often changed suddenly, and widely from that which they just previously had. This racking tendency produces strains in the different parts, and it is fair to suppose that the continual motion of the two

By J. B. Cheatham

Applications of the mathematical theory of plasticity promises to lead to the solution of many drilling and rock mechanics problems. ,Because of mathematical considerations, the inelastic behavior of rock has frequently been represented by a perfectly plastic model in conjunction with a yield criteria of the Coulomb or Mohr type. The totality of all stress states for which a solid ceases to behave elastically can be represented as a limit surface in stress space. ,Probing of such limit surfaces indicates details of strain hardening which are not provided by the standard triaxial testing procedure. Probing tests of the limit surfaces have been performed on Cardova Cream limestone to provide data for extending plasticity theory to cover situations in which consolidation and strain hardening are present. Test results indicate that this highly porous limestone undergoes a permanent volume decrease when it is subjected to hydrostatic pressures in excess of 3,500 psi. A virgin sample tested under a confining pressure of 1,500 psi has a yield strength of 1,700 psi; however, if the sample is subjected to a consolidation pressure of 5,000 psi, before testing at 1,500 psi, the yield strength is raised to 2,300 psi. Thus, both consolidation and strain hardening are important considerations in describing the mechanical behavior of this limestone. Tests conducted with the axis of the core having different orientations indicate that this rock is also anisotropic. Portions of the initial and subsequent limit surfaces are determined for samples loaded either perpendicular or parallel to the bedding planes. INTRODUCTION Previous experimental work in rock mechanics indicates that no mathematically tractable constitutive theory is inclusive enough to describe the mechanical behavior exhibited by all types of rocks under all conditions of stress and temperature. Indeed, the type of deformation encountered in a single type of rock is known to depend upon the stress and temperature conditions in the rock during deformation.l-5 Certain rocks and minerals, notably those minerals composed of ionic salts, have been shown to exhibit plastic deformation when tested under conditions of high confining pressure.2 Since the mathematical theory of plasticity provides simplifications over the theory of elasticity in certain types of problems, such as those in which limit analysis can be applied, it is of interest to know under what conditions plasticity theory may be applied to rock mechanics problems. The following factors determine the nature of the deformation a particular specimen will undergo: (1) the microscopic structure of the rock, i.e., the structure visible under an optical microscope, including number of phases, porosity, distribution of phases, (2) the mineralogical structure of the solid phases, (3) the conditions of stress and the rate of change of stress and (4) the temperature. Extensive experiments on Yule marble, Carthage marble, Solenhofen limestone and other calcerous rocks indicate that these relatively nonporous rocks deform plastically under certain conditions of loading, and creep under other conditions of loading.l-3 This study is concerned with the behavior of Cordova Cream limestone (Austin chalk) which is also composed almost entirely of calcite and thus has the same mineralogical composition but, because of a rather large porosity, it possesses a different microscopic structure. This investigation was undertaken to learn if Cordova Cream limestone deforms plastically despite the embrittling effect of the pore spaces, and to provide data which can be used to determine whether the mathematical theory of plasticity can describe the mechanical behavior of Cordova Cream limestone. A brief discussion is given of some of the underlying principles of plasticity theory concerning stress space and limit surfaces. Strain hardening is defined and a discussion is given of how the mechanical behavior of a material can depend on the previous history of the loading. This is followed by a description of the experimental procedures

By K. H. Ribe

Water often enters an oil reservoir during completion or workover operations on a well and forms a partial "water block" to oil production. A mathematical study of radial two-phase flow, neglecting capillary effects, has been employed to study the formation of such a water block and subsequent re-moval by production from the well. The effects of reduced oil permeability about the well on the well productivity were studied. The fluid saturation distributions about the well during formation and removal of the water block have also been computed. Several relative permeability relations and viscosity ratios were employed. If water has invaded the formation, its influence through relative permeability effects alone can cause the following. 1. Oil productivity will be depressed for extended periods after production is resumed and will build up only gradually as the water is removed. 2. Oil injected for treatment of water blocking will delay rather than promote restoration of full well productivity by enlarging the region invaded by water. Thus, unless the specific action of chemicals contained in the oil is needed, oil injection appears undesirable. INTRODUCTION During oilwell workover operations, water may enter the oil-bearing formation from the wellbore. When production is resumed, oil must flow through the region invaded by this water. The presence of this region can cause both well productivity and oil production rate to be low and oil to be produced with high water-oil ratio for some time after production is initiated. This situation is sometimes described as a water block. The introduction of water into the formation may result in other actions which also lead to reduction in well productivity and which are also usually included in the connotation of the broad term, water block. Often considered, for example, are the possibilities of clay swelling by contact with fresh water and the formation of emulsions with the formation oil. If it is suspected that such specific actions have taken place, remedial treatments are undertaken which usually involve the injection of chemicals in oil. Since the introduction of water, even in the absence of specific interactions with the formation or oil, will cause a temporary water block (which might be misinterpreted as evidence of a more severe situation), it is of importance to evaluate the magnitude and duration of this blocking which results purely from the reduction in relative permeability to oil in the vicinity of the wellbore. It is also of interest to evaluate the effect of oil injection on the productivity of a well blocked by water in this manner. Inasmuch as this unfavorable condition may persist for some time, it may lead to premature condemnation of a workover or premature abandonment of a potentially productive pay zone. A quantitative evaluation of the influence of water entry on the oil productivity through changes in relative permeability was made by solving a radial form of the Buckley-Leverett equation. The distribution of water saturation around the wellbore during the entry of water was calculated and was followed by a similar calculation of the saturation distribution during the period of resuming production. At any stage in the removal of the invading water, knowledge of the distribution of its saturation permitted calculating the attendant loss in oil productivity. The influence of the shape of the relative permeability relationships was also evaluated by carrying out the calculations for two hypothetical cases. Further, the effect of the oil-water viscosity ratio was examined by repeating the calculations, for several ratios of unity and greater, with the same relative permeabilities. Fi-nally, results are presented to show how the length of time a well must be swabbed to resume production depends on the length of time it has been subjected to invasion by water. STATEMENT OF THEORY Differential Equations Water is assumed to enter a producing formation which is initially at the connate-water saturation and contains no gas. The water and oil are treated as in-

By M. R. Tek, R. L. Nielsen

The scaling laws as formulated by Rapoport relate dynamically similar flow systems in porous media each involving two immiscible, incompressible fluids. A two-dimensional numerical technique for solving the differential equations describing systems of this type has been employed to assess the practical value of the scaling laws in light of the virtually unscalable nature of relative permeability and capillary pressure curves and boundary conditions. Two hypothetical systems — a gas reservoir subject to water drive and the laboratory scaled model of that reservoir — were investigated with emphasis placed on water coning near a production well. Comparison of the computed behavior of these par ticular systems shows that water coning in the reservoir would be more severe than one would expect from an experimental study of a laboratory model scaled within practical limits to the reservoir system. This paper also presents modifications of the scaling laws which are available for systems that can be described adequately in two-dimensional Cartesian coordinates. INTRODUCTION Present day digital computing equipment and methods of numerical analysis allow realistic and quantitative studies to be carried out for many two-phase flow systems in porous media. Before these tools became available the anticipated behavior of systems of this type could be inferred only from analytical solutions of simplified mathematical models or from experimental studies performed on laboratory models. To reproduce the behavior of a reservoir system on the laboratory scale, certain relationships must be satisfied between physical and geometric properties of the reservoir and laboratory systems. Where the reservoir fluids may be considered as two immiscible and incompressible phases, the necessary relationships have been formulated by Rapoportl and others.2-5 Rapoport's scaling laws follow from inspectional analysis of the differential equation describing phase saturation distribution in such systems. It will be recalled that these scaling laws presuppose three conditions: (1) the relative permeability curves must be identical for the model and prototype; (2) the capillary pressure curve (function of phase saturation) for the model must be linearly related to that of the prototype; and (3) boundary conditions imposed on the model must duplicate those existing at the boundaries of the prototype. These three requirements seldom if ever can be satisfied in scaling an actual reservoir to the laboratory system because: (1) The laboratory medium normally will be unconsolidated (glass beads or sand) while the reservoir usually is consolidated. Relative permeability and capillary pressure curves are usually quite different for consolidated and unconsolidated porous media. (2) The reservoir usually will be surrounded by a large aquifer which could be simulated in the laboratory only to a limited extent. (3) Wells present in the reservoir would scale to microscopic dimensions in the laboratory if geometric similarity is to be maintained. In view of these considerations, rigorous scaling of even a totally defined reservoir probably would never be possible. The purpose of this paper is to assess the practical value of the scaling laws in the light of the unscalable variables. This has been done by carrying out numerical solutions in two dimensions to the differential equations describing the flow of two immiscible, incompressible fluids in porous media for a field scale reservoir and a laboratory model of that reservoir. While both the reservoir and the laboratory model were purely fictional, each has been made as realistic and representative as possible. The field problem selected as the basis for the investigation was an inhomogeneous, layered gas reservoir initially at capillary gravitational equilibrium and subsequently produced in the presence of water drive. The laboratory model of this reservoir was designed to utilize oil and water in a glass bead pack. The numerical treatment employed was similar to that of Douglas, Peaceman and Rachford6 and it included both capillary and gravitational forces as

By J. P. Heller

The purpose of this discussion is to comment briefly on one of the points raised in the paper. The topic is that referred to in the Conclusion 2, in which the authors note a "disproportionate increase in mixing rate" as they proceed from a displacement performed at viscosity ratio 0.998 to one at 1.002. The recovery curve of this displacement, Fig.2A, shows a quite noticeable difference between the two runs in which displaced and displacing liquids were interchanged. This difference is ascribed by the authors to the change in viscosity ratio between the two experiments. This interpretation of the experimental results must certainly be based on rather precise measurements of the fluid viscosities or, at least, of the ratio of their viscosities. The ordinary Ostwald viscometer, employed in the usual way, is only just able to distinguish between liquids whose viscosities differ by two or three parts per thousand. The experimental refinements which then must have been necessary to measure with any precision a viscosity ratio of 1.002 would be worthy of mention. Indeed, one would have expected the authors to devote at least a paragraph to the display of this technique. For the same reason - to emphasize that the details of such unusually critical experimental work should be reported and discussed in the literature — I would have liked to see some comment by the authors on the requirements for temperature control during both the viscosity measurements and the displacement experiments. The temperature coefficient of viscosity of Soltrol C at room temperature has been measured in this laboratory at about - 1.9%/1°C. How much did the coefficients for the fluid mixtures used in the experiment differ from this and from each other Even if the two fluids were identical in this respect, it would appear that a temperature difference of .l°C between them in the flow system would produce a viscosity con- trast equal to that imposed by their composition. Reference to the experimental solutions of these problems would render Conclusion 2 more credible. Published theoretical work on the stability problem contains no indication of such extreme sensitivity to the viscosity ratio in the immediate neighborhood of unity. On the contrary, perturbation calculations such as those made by Chuoke1 and Perrine2 lead to initial finger growth rates proportional to the first power of the viscosity difference. The later instability growth calculated by Scheidegger,3 using the assumption of zero pressure gradient perpendicular to the average flow, is also simply proportional to the viscosity difference. The conflict regarding the sensitivity of the dependence of instability growth or mixing rate on the viscosity contrast needs explicit recognition in the literature. In particular, it would seem that the experimental work reported here should be supported by a more detailed description of the viscosity controls. Further discussion of possible alternative interpretations of the data would also be in order. For instance, is it possible that the unconsolidstted bead pack used changed its characteristics between the runs (due, perhaps, to the handling and inverting involved in the injection of the front) Or, are the data points of Fig. 2A averages from several runs performed in a series in which the two liquids alternated as displaced and displacing fluids? To summarize, it seems to me that the relationship asserted by Conclusion 2 in the paper needs the support of more detailed experimental evidence. REFERENCES 1. Chouke, R. L., van Mews, P. and van der Poel, C.: "The Instability of Slow Immiscible Viscous Liquid-Liquid Displacements in Porous Media", Trans., AIME (1959) Vol. 216, 188.