Search Documents

By B. A. Eaton

The subject of many discussions and technical papers in the last 20 years has been the prediction of the well-bore pressure gradients that are required to induce or extend fractures in subsurface formations. The subject merits this attention because of the frequently recurring problems that arise from an inability to predict fracture pressure gradients. Encountered in several common types of operations in the oil industry are problems associated with the prediction of formation fracture pressure gradients. When wells are being drilled in both new and old fields, lost circulation is often a very troublesome and expensive problem. Complete loss of circulation has been disastrous in some cases. Many times, such disasters could have been avoided if techniques for calculating fracture pressure gradient had been employed in the well plans, and if casing strings had been set, and mud weight plans had been followed accordingly. In areas of abnormally pressured formations, the prediction of fracture gradients during the well-planning stage is extremely important. In fact, it is as important as the prediction of formation pressure gradients, which has received a great deal of attention in recent years. There are several published methods used to determine fracture pressure gradients. However, none of these methods appears to be general enough to be used with much reliability in all areas. In 1957, Hub-bert and Willis published a classical paper that included the development of an equation used to predict the fracture extension pressure gradient in areas of incipient normal fau1ting.l Overburden stress gradient, formation pore pressure gradient and Poisson's ratio of rocks were the independent variables that were shown to control fracture pressure gradient, the dependent variable. In 1967, Matthews and Kelly published another fracture pressure gradient equation that is different from that of Hubbert and Willis in that a variable "matrix stress coefficient" concept was utilized.3 Later the same year, Costley wrote about a similar idea.5 Goldsmith and Wilson used a least-squares curve-fitting technique and field data from the Gulf Coast area to correlate fracture pressure gradient with formation pore pressure gradient and formation depth.' They noted that the fracture pressure gradient increased with increasing depth while the pore pressure gradient remained constant. In each of these cases, the problem for which a solution was sought was to determine the bottom-hole pressure gradient required to initiate or extend a fracture. Results of the previous work show that fracture pressure gradient is a function primarily of overburden stress gradient, pore pressure gradient, and the ratio of horizontal to vertical stress. There is argument for a fourth variable in that in many cases breakdown fracture pressure gradient is greater than the fracture extension pressure gradient. However, if the fracturing fluid is able to penetrate the formation through the pores or existing cracks, there is very little

Jan 1, 1970

By K. E. Brown, J. P. Brill, A. R. Hagedorn, T. C. Doerr

A recent correlation for multiphase vertical flow by Hagedorn and Brown1,7 has been used to examine the effect of numerous variables such as pipe size, gas-liquid ratio, liquid flow rate, liquid viscosity, liquid surface tension, oil API gravity, water cut, etc., on flowing pressure gradients. Most previous correlations were derived by using a form of the general energy equation which assumed negligible changes in kinetic energy. Under certain flow conditions this can lead to serious errors in predicting pressure traverses. Examples of traverses calculated both with and without the kinetic energy term are included here. The Chew and Connully correlation for calculating changes in viscosity with changing pressure and temperature has been extended so that the effect of a changing viscosity with depth on pressure gradients can also be predicted. Examples illustrating the effect of tubing size show the importance of this correlation in determining pipe sizes for initial well completions. Limitations on production rates in small pipe sizes from excessive friction losses can be predicted with greater accuracy. The combined use of a vertical-flow correlation with a horizontal-flow correlation provides a method for determining the maximum flow rate possible from both flowing and gas-lift wells. Examples are included to illustrate application of the combined use of vertical and horizontal-flow correlations in optimizing production rates and minimizing excessive costs from over injection of gas in gas-lift wells. INTRODUCTION The problems of multiphase vertical and horizontal flow are immediately encountered in any producing well. Fig. 1' shows three stages of production: flow in porous medium, vertical flow and horizontal flow. All three stages affect the production rate from a well. If we are confronted with a particular well and consider only the variables involved in the vertical and horizontal stages of Fig. 1, we should be capable of calculating the flowing bottom-hole pressure necessary to produce a particular flow rate. A change in any of the variables results in a new solution to the problem. Recent approaches in developing horizontal and vertical-flow correlations have improved to the extent that the solution can now be obtained with more confidence than previously. VERTICAL AND HORIZONTAL FLOW CORRELATIONS Since Poettmann and Carpenter's' original work in multiphase vertical flow, there have been many additional vertical-flow correlations developed which have greatly improved the accuracy of pressure-loss calculations. These

Jan 1, 1967

By L. H. DeWald

During the examination of various commercial brands of copper wire bars, which had exhibited different degrees of adaptability for being drawn into fine gages of wire by the present-day high speed machine, it was found that each brand had a more or less characteristic crystal macrostructure. The extremes of these macrostructures are similar FIG. 1.-EFFECT OF METAL TEMPERATURE IN CASTING ON CRYSTAL MACROSTRUCTURE OF BAR. Reading from left to right, metal temperatures were 1166" C., 1149" C., 1132" C.. 1121" C., 1093" C. to the maximum variation shown in the photomacrographs reproduced in Fig. 1. As the various brands did not differ materially in chemical composition, it was concluded that the differences in structure were the result of variations in casting technique.

Jan 1, 1931

By H. A. Domian, H. Biloni, G. F. Bolling

Jan 1, 1965

By B. E. Warren

Jan 1, 1965

The income of the Institute is derived mainly from dues, advertising in MINING AND METALLURGY and sale of publications. These sources are fortunately supplemented by the interest from invested funds now amounting to over $500,000, a considerable portion constituting endowments for especial purposes. Aside from the medal funds the principal endowments are the James Douglas Fund for support of the Library, the Rocky Mountain Fund and the Seeley W. Mudd Memorial Fund. The income from the two last named is available for the support of research and a variety of other purposes. It is allocated by the Board of Directors upon recommendation of standing committees in each case. An active effort is being made to increase the endowments of the Institute, it being felt that unusual opportunities exist for research and public and professional service through the organization.

Jan 1, 1939

By James C. M. Li, J. Weertman

James C. M. Li (United States Steel Covp.)—The author has discovered a single analytical relation between the minimum, or the steady-state, creep rate and the applied stress confirming a statement made earlier by orn" that such relation should exist and should reduce to the power law at low-stress levels and to the exponential law at high-stress levels. However, as pointed out by the author, all the available theories of creep are not adequate to explain such a single relation. Inspired by such a statement, I have made a suggestion concerning a possible explanation13 of the new relation. This is based on the concept that the rate-controlling step is the cooperative motion of a short dislocation segment. The nature of and the reason for the cooperation have been discussed13 as well as the reason why the motion of the dislocation could be controlled by the mobility of the short segments. A natural short segment has been proposed13 to be the node formed at the intersection of two dislocations. The reason why the node may contain only a few atoms has not been discussed and the application of the proposed model to all the metals and alloys reported in this paper has not been attempted. It seems appropriate to include these in this discussion. The model was partly inspired by the substructure recently observed in the electron microscope on creep specimens. One such example is shown in Fig. 48 in Darken's Campbell Memorial Lecture. The substructure is a collection of dislocation intersections and irregular networks. No pileups have ever been observed in such specimens. A closer look at the network indicates that one set of dislocations is more straight than the other set sug- gesting that the other set may be mobile in the plane of the net. Since, for the network to survive the creep deformation, some dislocation reaction has to take place at the intersection of the two sets of dislocations, the mobile dislocation, even if it can slip in the plane of the net, will have difficulty in moving at these intersections. To illustrate, an idealized network is shown in Fig. 9. The thin lines represent the set of dislocations which are mobile and can slip in the plane of the net except at the nodes which are formed at the intersection of the two sets of dislocations. The other set of dislocations, represented by the thick lines, may have a Burgers vector which is not parallel to the plane of the net, and therefore may be immobile. Because of this nonparallel Burgers vector, the Burgers vector of the node, which is avectorial sum of the Burgers vectors of the two sets of dislocations, will also be nonparallel to the plane of the net. Therefore, the node cannot slip in the plane of the net and the set of mobile dislocations will have to drag the node along as indicated in Fig. 9. Now the problem is how many atoms could be involved in the nodes. Exact calculation of the length of the node is a difficult task especially since the length should be a steady-state length and not an equilibrium length. An estimate of the equilibrium length can be made by balancing the line tension of all the dislocations. The results would be a small fraction of the length 1 between nodes as shown in Fig. 10. In some instances such a calculation shows that the node will not be formed. But, since the calculation is only approximate, it does not rule out the possibility of the formation of a short node. In other instances the equilibrium length of the node is an appreciable fraction of I. Such a node would

Jan 1, 1963

By E. C. Olsen

This paper describes the unusual hazardous conditions encountered in mining thick, inclined coal beds under heavy cover in Carbon County, Utah. This description includes heavy roof conditions, sloughing ribs, steep pitches, and severe bouncing experience. Emphasis is placed on the top bench system of mining and development in the top half of a 16-ft seam, under the roof rock, which is bolted, and the recovery of the bottom coal on retreat; success of bolted roof and ribs; and shuttle car haulage on pitches up to 25%. The use of protective cages for operators of face equipment, the guarding of belt conveyor systems to prevent injury associated with deliberate entry, and the positive lockout of electrical equipment is described and illustrated. Most coal mining areas of the Western United States are characterized either by thick beds, steep pitches or heavy cover. Individually, each of these may present inherent safety hazards that in- fluence mining. When all are met at one location, the combined effect is to make the safe operation of the mining property a challenge to management. These combined conditions are met in the coal mines of the Columbia-Geneva Steel Div. of the United States Steel Corp., in southeastern Utah. This coal bed measures up to 16 ft in thickness, the pitch is as much as 25% and mining is under a cover approaching 3000 ft. The Sunnyside coal bed is found in the Book Cliff Mountains, a great escarpment extending from west of Castle Gate, Utah, some 180 miles to Grand Junction, Colo. West of Green River, Utah, the Book Cliffs swing sharply to the north because of the domelike uplift of the San Rafael Swell to the southwest. It is in this northerly extension of the Book Cliffs that the Sunnyside coal in the Blackhawk Formation of the Mesa Verde Group, Upper Cretaceous Age is found. The strata are transected by numerous faults with from a few inches to as much as a hundred feet of displacement. All of these faults have had some effect on the mining operation, either in roof control or, in the case of major faulting, the pattern of mine development. MINING PRACTICES Mining has generally been confined to the Lower Sunnyside Seam which measures from 7 to 16 ft in

Jan 1, 1963

By Samuel M. Shallcross

The principal function of the Bell limestone mine, of the American Lime & Stone Co., at Bellefonte, Pa., is to supply raw material to the company's modern rotary kiln plant at Bellefonte. Because mining costs are naturally higher than many open quarries with more favorably lying strata, it cannot compete, except to a limited extent, in the markets for road material and fluxing stone. In the broad sense, mining does have an advantage over the average open quarry when the end products are designed for and controlled by process and final inspection for the exacting demands of the industrial and chemical markets. Approximately 80 per cent of the output of the Bell mine is sold in the industrial and chemical markets. And as higher standards are set each year for construction and agricultural work, the remainder of the mine output goes to these two classes. Mining produces more fines than open quarries, and these fines are difficult to market; on the other hand, owing to absence of contaminating clay and impure rock, a greater percentage of the output of a mine is marketable. Because of the great production of fines, rotary kilns for lime burning are necessary when the stone must be mined. In 1921 the American Lime & Stone Co. operated 19 quarries and 72 lime shaft kilns spread over an area of 40 miles between Bellefonte and Altoona, Pa. At that time the lime-production capacity of the 72 kilns amounted to approximately 700 tons per day. After 15 years the company operates three rotary kilns and six shaft kilns, with a combined output of approximately 500 tons of lime per day, concentrated in its single Bellefonte plant (Fig. 1). Its only stone-producing plants are now the limestone mine at Bellefonte and a construction-stone open quarry at Union Furnace, near Tyrone, Pennsylvania. Since the middle of 1922, the management of the American Lime R. Stone Co. has been under the direction and ownership of the Warner Company, of Philadelphia. To the Warner Company's experience, gained through successful operation of lime plants in eastern Pennsylvania for many years, and the successful solution of material-handling problems

Jan 1, 1937

By D. N. Dietz

Jan 1, 1969

By Charles Hardy

?POWDER METALLURGY? is the production of semiformed or fully formed metal products by compressing metal powders. It had its beginnings in the fabrication of tungsten and molybdenum bars and wire by compression, sintering, and swaging. Powder metallurgy, although still in its infancy with other metals, has shown remarkable progress in recent years. Production of Metal Powders-The first step in any powder-metallurgy process is the production of metal particles of suitable shape and size. Mechanical processes for producing metal powders, as by crushing, stamping, or grinding, are usually not suitable, as the resulting particles may cause laminations and decreased strength in the finished articles. The principal metal powders now in use on a large scale are produced by electrolysis or by reduction of metal oxides by gases.

Jan 1, 1936

By L. A. Panek, F. D. Wang, M. C. Sun

A method of stability analysis of underground opening by the limiting equilibrium approach is presented. The Coulomb shear failure criterion is employed to find possible fracture surfaces and to evaluate the shear strengths along them. The normal stresses required for the calculation of shear strength, and the shear stresses for the calculation of total shear force along the fracture surfaces, are obtained by the finite-element method of plane-strain analysis. The factor of safety against Coulomb shear failure is calculated as the ratio of the shear strength to the shear force along a fracture surface, similar to the practice in slope stability analysis. Four shapes of openings are analyzed, with three ratios of horizontal lithostatic stress and four coefficients of internal friction. The critical fracture surfaces found by this method of analysis are consistent with the arching of roofs and slabbing of ribs commonly observed.

Jan 1, 1972

By David R. Mitchell, Ernest M. Spokes

There is evidence that study of minerals now considered to have susceptibilities too low for magnetic separation should be continued. Present concepts may be false. INFORMATION on magnetic properties of minerals other than magnetite is scanty and sometimes erroneous. Often there is no information at all. W. R. Crane's table of tractive forces published in 1901 is still the best guide in the field, but Crane's data can be used as a guide only, since he ignored such important factors as impurities in the minerals, size, and packing effect. During the 1930's work by the U. S. Bureau of Mines added greatly to knowledge of the behavior of certain minerals in a magnetic field. The USBM bulletin gave for the first time complete magnetic characteristics for magnetite and fragmentary data for a few other minerals. The work reported here covers a study of the magnetic-susceptibility relationships of the wolframite series of minerals and sphalerites as related to their chemical compositions. Theories are proposed to explain divergences found from commonly accepted beliefs. Gouy Method for Determining Magnetic Susceptibility: The Gouy principle has been adequately described in the literature1-3 and applied in various ways.'- It may be summarized by stating that a body of material (gas, liquid, solid, or powder mixed with gas or liquid) of cylindrical shape changes in apparent weight if suspended with one end in a zone of higher magnetic field strength than is found in that portion of the field surrounding the other end of the cylinder. This may be expressed as K2 - K1 A (H2 - H2) - mg. [1] In determining the susceptibility of mineral specimens by the Gouy method, the mineral is reduced to a fine powder and packed in a glass weighing tube. Precise measurement of the bore and depth of the tube is made unnecessary by a calibration procedure. Using the basic expression of Eq. 1, Hoare"

Jan 1, 1959

By H. J. de Haan, J. van Lookeren

In view of their size and favorable characteristics, Shell's heavy-oil fields on the eastern coast of Lake Maracaibo, Venezuela — known as the "Bolivar Coast" — are obvious candidates for the application of thermal recovery methods. During the last decade various methods have been tested, such as steam drive,', steam soak and in-situ combustion.3 Steam soak was discovered as a promising production method rather accidentally in 1969, during early steam drive testing in the Mene Grande Tar Sands1 When steam erupted at the surface due to breakdown of the overburden, the injection wells were back-flowed to relieve the reservoir pressure. This resulted in high oil production rates, all the more impressive because the reservoir is unproducible by primary means. It was concluded that injection of limited amounts of steam might be a very effective method for stimulation of heavy-oil wells. To obtain information that would be more widely applicable, testing of the process was continued in the eastern part of the Tia Juana field (Fig. 1). Here even more favorable results were obtained than in Mene Grande. This created sufficient confidence to prompt a large-scale project, which was considered essential for an adequate evaluation of ultimate recovery, steam requirements, mechanical problems and economic prospects. To obtain an early answer on the ultimate recovery a fairly depleted area was selected in the western part of the same field ("A" in Fig. 1). The first steam was injected in March, 1964. This paper describes the performance until the end of 1966. The production mechanism and well inflow performance are discussed in some detail. In view of the favorable results obtained so far, steam soak operations have expanded considerably since the start of this project. Presently, five other projects are under way, of which two are about the same size as the one under discussion, with further extensions being programmed. Description of the Project Reservoir Data Shell's part of the Main Tia Juana reservoir has a proven area of 8,800 acres and an average net oil sand thickness of 130 ft. It contained an initial quantity of 2,750 million bbl of stock tank oil (STOIIP), with an API gravity ranging from 10 to 15 API and a viscosity of 100 to 10,000 cp (at 122F). As a geological unit it is known locally as the Lower Lagun-illas member of the Lagunillas formation, which is early Miocene in age. The sand bodies are believed to be distributary channel sands deposited in a lower coastal plain. Most of the channel sands have lateral and some vertical communication. The field was developed mainly in the period from 1936 to 1950. The ultimate primary recovery is estimated at 18.1 percent of STOIIP, of which 12.5 percent STOIIP had been produced by the end of 1963 from some 540 wells. By then the average pressure had decreased from an initial value of 960 to 550

Jan 1, 1970

By James Barr

Geology and Mineralogy1 TENNESSEE phosphates are commercially divided into three varieties: Brown, blue and white. The first two only are now of commercial importance. The white phosphates of Perry county are not being mined at present owing to their pockety occurrence. The white phosphate is Fig. I.-DEPOSITS OF 75 PER CENT. BROWN PLATE ROCK UNDERLAIN BY SAND. the highest grade of the three varieties, running often more than 90 per cent. tricalcium phosphate. The brown rock is found in three varieties, grading from the dense almost vitreous gray rock, 76 to 82 per cent. tricalcium phosphate, down

Jan 9, 1914

By Waldemar Lindgren

THE following pages present data accepted by many geochemists and geologists regarding the succession of minerals and the temperatures of formation in ore deposits affiliated with igneous rocks. They also present the individual views of the author on the composition of the solutions that were active and on the mode of deposition. He does not regard these views as speculations or hypotheses but as forming a well motivated and consistent theory. It is not certain, however, that all investigators will subscribe to them.

Jan 1, 1936

By E. H. Herron, J. H. Perry

Mathematical simulation of reservoir behavior may be used to help understand reservoir processes and to predict reservoir behavior, thereby leading to the most economically desirable form of exploitation. In addition, simulation can be used as a tool for reservoir description to learn more about the physical nature of the reservoir and the mode of primary recovery. This use is essential in most reservoir studies and represents one of the more significant applications of simulation. Prediction of reservoir behavior provides information concerning displacement efficiency, optimum well locations, and the comparison of alternative producing processes. Since oil, water and gas all commonly occur in many reservoirs, a simulation that takes into account the simultaneous flow of three immiscible phases is a prerequisite for obtaining such information. Furthermore, a simulation that describes this flow in two dimensions provides the capability of evaluating combination gas- and water-drive reservoirs with either an areal or a cross-sectional description. Many reservoir problems require such a simulation. The objective of the research described here was to develop a mathematical model and associated computer program for accurate, efficient, and economical prediction of the reservoir flow of three phases in two-dimensional geometry. Several authors have discussed multidimensional reservoir simulation of two phases.l-3 Gottfried et al.' presented an analy- sis of three-phase simulation in one dimension. Fagin and Stewart" and Garrett6 have discussed the simulation of three-phase flow in two dimensions, but their solution techniques were different from the method discussed in this paper. In particular, previous models and computer programs have usually involved either fewer capabilities or inefficient and less rigorous descriptions resulting in limited applicability. In the following sections, we present a description of the model, an evaluation of the method, and the application to reservoir problems. Additional and supporting material is presented in the Appendix. The Three-Phase Model General Description The three-phase model consists of the mathematical description of the flow of oil, water and gas in a reservoir. Conceptually this description considers flow in all three space dimensions, but in the subject simulator it is assumed that no flow occurs in the third dimension. The model includes the effects on reservoir behavior of fluid and rock compressibility, viscosity, gravity, capillary pressure, relative permeability and gas solubility. Within the two-dimensional restriction, the reservoir considered may be completely heterogeneous and anisotropic. The differential equations governing multiphase, multidimensional flow in a reservoir are too complex to be solved analytically, so they are approximated by a system of finite difference equations and solved

Jan 1, 1970

Executive JOHN M. LOVEJOY, Chairman ERLE V. DAVELER G. B. WATERHOUSE WISER JUDSON WILLIAM WRAITH J. V. W. REYNDERS, Consultant Finance HENRY KRUM, Chairman PAUL D. MERICA R. M. ROOSEVELT H. G. MOULTON, Consultant Admissions. RALPH M. ROOSEVELT, Chairman GEORGE D. BARRON, Vice-chairman JOHN L. CHRISTI PHILIP KRAFT W. A. SINSHEIMER

Jan 1, 1936

By W. R. Hibbard

Evidence is presented which confirms previous findings that models of solution strengthening depending solely on lattice parameter changes are incomplete. Direct evidence for the Suzuki interaction of chemical effects was obtained from alloys with constant lattice parameter, for which solution strengthening was nearly independent of temperature from 300° to 20°K. SOLUTION hardening principles have been reviewed recently by Parker and Hazlett.' In addition, the interactions of dislocations and solute atoms were reviewed by Cottrell.' It is not proposed to review this information again here, but to refer only to those details which are pertinent to the developments which have occurred since these papers were published. Literature Survey To set the stage, it has been known for some time that solid solutions are stronger than the pure parent metals, and that the strength increases with increasing amounts of alloying element. One of the most extensive earlier investigations was that of Lacy and Gensamer,3 who reported that the solid solution strengthening of iron is a power function of the composition, that the relative strengthening is a function of the solubility limit, and that the rate of strain hardening is a function of the amount of solution hardening. Thus, the slope of the plastic portion of the stress-strain curve increases with increasing yield strength. An approach to the development of principles resulted from work by Gulyeav,4 who reported that the strengthening effect of solid solution alloys is periodic in character when plotted as a function of the atomic number of the solute. A more quantitative approach was to relate the amount of solution hardening to the lattice parameter of a solid solution. This was done quite early by Norbury,5 and later by Brick, Martin, and Angier,6 Frye and Hume-Rothery,7 and French and Hibbard. This relationship is only approximate in character. The next development was to relate the amount of solution hardening to lattice stress, as proposed by Frye and co-workers.9 They calculated what they called ionic overlap from the lattice parameter changes represented as a change in volume, using Bridgman's hydrostatic compression equation. The resulting correlation was somewhat improved as compared to the simple correlation with lattice parameter. The next step was carried out by Dorn, Pietrow-sky, and Tietz,10 who suggested an effect of valence differences in addition to the effect of lattice parameter. Their proposed parameter shows flow stress as a single valued function of the concentration and a relationship involving these two atomic variables. In addition, they contributed new data on the temperature dependence of solution hardening which indicated that it was a strong function of temperature for aluminum alloys. Allen, Schofield, and Tate11 carried the valency effect even further by reporting data indicating that with a given solvent element the stress-strain curve is a function of electron density, entirely neglecting the size factor, which varied by a factor of two. This evidence led to a quantum mechanical calculation by Cottrell, Hunter, and Nabarro,12 based on the formation of a dipole around an edge disloca-

Jan 1, 1959

By Robert G. Nisle

The classic theory of pressure build-up in shut-in oil wells as developed by Horner and van Everdingen is based on two-dimensional radial symmetry in the well-reservoir system. Such symmetry does not exist in the case of a well which parfiially penetrates the producing formation. As a result of this lack of symmetry the use of the classic theory in such cases become questionable. In this paper the mathematical theory has been extended to include the case of partially penetrating wells. Numerical solutions illustrating the up-dication of the equations are presented. The effect of partial penetration on pressure build-up is shown by a comparison of synthetic pressure build-up curves derived from the numerical solution of the equation for partially penetrating wells for various degrees of penetration. It is shown that partial penetration is detectable from the characteristic shape of the pressure build-up curve and that formation productivity may be calculated from the pressure buildup data in a manner identical to that described by the classic theory. INTRO DUCTION The use of pressure build-up data on shut-in oil wells is a well-established technique for the measurement of reservoir productivity. The work of Homer' and van Everdingen' is well known. Since the publication of their pioneering work, others"' have elaborated and extended the technique. These authors have dealt exclusively with the case of wells that completely penetrate he producing formation. Two-dimensional radial symmetry is achieved, thereby, and the problem is greatly simplified. This ideal situation is seldom encountered in practice. But, in spite of this, the technique is often employed anyway. The question, therefore, arises: how much confidence can be placed in the results of such calculations when the data are taken on wells that only partially penetrate the producing formation? It is the purpose of this work to answer this question. Only the single-phase fluid case will be considered. van Everdigen' has presented a rigorous solution of the two-dimensional case, but because of mathematical complexity, this solution has not proved to be particularly useful. Horner' utilized Kelvin's point source solution in two dimensions and showed that it is a sufficiently good approximation for oil reservoir studies. The geometry of the well and the producing formation is illustrated in Fig. l. Geometrically the formation may be described as an infinite slab. The infinite slab is defined as the solid (the porous medium in this case) bounded by two parallel planes. Thus, the slab extends to very great distances in both the positive and negative directions of two of the coordinate axes, but is finite in