Search Documents

By H. C. H. Darley

The principal causes of unstable boreholes* have been known for many years. For example, in a paper published in 1938, Halbouty and Kaldenbachl listed nearly all the classes of troublesome shales that we know today. Much progress has been made since then, but we still have difficulty in identifying the cause of hole problems and consequently in selecting the right technique for combatting them. This subject was discussed at length by Kelly2 in a recent paper and much valuable information was provided. In the Shell Laboratory it was felt that there was a need for a better understanding of the basic mechanisms underlying the various forms of hole instability, and accordingly, a study was undertaken. This entailed reviewing the physics and chemistry of shales, developing tests for characterizing shales, and constructing a model to study borehole stability. Full simulation of underground conditions was not attempted; for example, we used reconstituted shale specimens, and often simple solutions or suspensions rather than practical drilling muds. Such procedures are necessary in order to establish mechanisms, but caution must be exercised in applying the results to the more complex conditions in the field. Part I — Compaction Swelling and Dispersion of Shales Compaction of Sediments Shales are formed by the compaction of sediments. Water is expressed as sediments are buried by subsequent layers, and the degree of compaction is proportional to the depth of burial, provided the water is able to escape easily to permeable formations. The younger sediments, when recovered at the surface, soften and disperse when mixed with water; the older shales have undergone diagenesis (alteration of the clay minerals, secondary cementation, etc.), and remain hard and unaffected by water. In this paper we shall apply the term "shale", as the drilling engineer does, to everything from clays, which are highly reactive to water, to completely lithified materials such as claystones and slates, which are completely inert. However, because these materials behave quite differently when encountered in the drilling well, it is desirable to develop a simple means of characterizing them. The feature that distinguishes the various shales is their dispersibility in water —soft clays disperse readily, harder shales disperse slowly when agitated, the lithified materials will not disperse at all unless milled. In the dispersibility test used in our study, a weighed amount of dry shale particles in the 12-to-100-mesh-size range was placed in a tube filled with distilled water and allowed to soak for 16 hours. The

Jan 1, 1970

By J. E. Forciea, O. E. Palasvirta, L. G. Hendrickson

AH pilot and commercial plants working with Mesabi Range taconite employ wet magnetic separation. Progress is being made with a dry magnetic process, but this has not yet been applied to taconite except in the laboratory and pilot plants. Flowsheets in all five plants are similar, including multiple stages of grinding, each followed by magnetic separation, which rejects a finished tailing. The flowsheet in Fig. 1 is used in the preliminary and commercial plants of the Erie and Reserve mining companies and in Unit C of Oliver's Pilotac plant. At Erie the desliming operation is omitted. In the Reserve plants two ball mills are used with each rod mill. Pilotac uses flowsheet 2 in Unit A-B and has also tested flowsheet 3. Principal differences in the three flowsheets are the number of grinding stages and the type of roughing. The terms open-circuit roughing and closed-circuit roughing are best defined by the diagrams. Special Test Data Given in Table I: For the Erie and Reserve commercial plants, Table I reflects routine conditions. Pilotac, the experimental plant, has been run with several flowsheets, including those followed at Erie and Reserve's commercial operations. Pilotac data for these last two flowsheets are similar to commercial plant results and are not repeated in Table I; statistics are given, however, for special test conditions. using other flowsheets. when the objective was to 'to extra-grade concentrate, Variations Among Taconites: Although crude ore tributary to these plants is generally classified as taconite, there are important differences in mineralogy and fineness of dissemination.

Jan 1, 1959

By W. A. Wood, H. M. Bendler

Cylindrical fatigue specimens of OFHC* copper were subjected to alternating torsion while under axial tension. The superimposed tension strongly influenced the fatigue life of the specimens as well as causing characteristic changes in the structure of the fissures that form and develop throughout the life of the specimen. A further effect of the superimposed tension is a continuous extension of the specimen throughout its life in a manner formally similar to creep at elevated temperatures. RECENTLY, studies of the structural changes leading to fracture of metals by fatigue have been made by a technique which reveals, in section, the lattice distortion arising from slip movements inside the metal and the disturbances built up by these slip movements at the metal surface.1'2 According to this approach, fatigue develops as follows: a) fine slip movements concentrate in bands, producing little strain-hardening and thus permitting abnormal totals of plastic strain before fracture and, therefore, abnormal amounts of lattice distortion along each slip band; b) both forward and backward slip movements occur in each band which, according to their distribution across the band, may build up surface disturbances that include sharp notches; c) these notches early in the life of a specimen penetrate the distorted slip zones and transform them into deep, fine fissures; d) the fissures on continued cycling locally saturate the metal and cause local collapse in the form of a microcrack. Thus, the significant consequence of fatigue or cyclic

Jan 1, 1962

By J. W. Mills

In 1956 the discovery of nickel on Jumbo Mountain, Snohomish County, Washington, focused attention on this part of the Cascade Range, far more renowned for its timber than for its mineral resources. Hand specimens assaying as high as 13 pct Ni encouraged Discovery Mines Inc. of Mount Vernon, Wash., to stake 12 claims in Township 31 North, Range 9 East of the Willamette Meridian. The property was leased to Climax Molybdenum Co. in January 1957. Trenching, sampling, and geologic studies were carried out during the late summer of 1957 when most of the snow had disappeared. Although the deposits are only about four miles south of the logging town of Darrington and about two miles. as the crow flies, southeast of the end of a good gravel road from that town, all camp supplies and equipment were air-dropped, owing to the extremely rugged nature of the terrain. GEOLOGY The rocks of Jumbo Mountain comprise a series of folded Tertiary sediments forming a belt one half to three quarters of a mile wide, trending northwest about parallel to the strike of the beds. This belt is sandwiched between two extensive in-trusives—gabbro on the northeast and quartz dio-rite on the southwest—and has been intruded by dikes of fresh and serpentinized dunite. Shear zones follow many of the dunite-sedimentary rock contacts, and in these the nickel mineralization is concentrated. A few transverse faults offset the sediments and ultrabasics. Sedimentary Rocks: Two kinds of quartzite make up the greater part of the sedimentary section. One is dark purplish gray, very fine-grained, and thin-bedded; the other is light gray to white, medium-grained, and well-bedded. Locally the light-colored quart lites contain grit and conglomerate horizons, rich in pebbles of white chert and very fine-grained quartzite. Intercalations of dark argillites near the summit of the mountain give many of the cliffs a striking banded appearance. Dark gray to black, aphanitic, thin-bedded argillites occur throughout the entire sedimentary belt. Beds vary from a few inches to a few tens of feet thick and some, near the summit, contain numerous fossile (willow?) leaves. Only the thicker beds were mapped. The sediments are quite fresh on the northeastern side of the belt and become more metamorphosed (hornfelsed) toward the southwestern side. Brown biotite is present in the metasediments in increasing amounts as the quartz diorite is approached. Near the intrusion, porphyroblasts of plagioclase, up to 1 in. across, testify to the soaking of sedimentary rocks by solutions emanating from the quartz diorite. Most of the sediments mapped lie along the southwestern limb of an anticline, the axis of which follows the crest of the mountain and plunges at a low angle to the northwest.

Jan 1, 1961

By V. J. Parks, Leonard Obert, A. J. Durelli

The state of stress in the region where core discing initiates has been investigated through the use of three-dimensional photoelastic models and the results of this study have been compared with those of a previous investigation of rock models. This comparison showed that discing initiates at a point of maximum shear stress, but that the magnitude of the shear stress, as determined photoelastically, is much larger than the shear strength of the rock as deter-mined from conventional triaxial testing. Moreover the shear stress required to produce discing is not constant, but depends on the applied stress field. Possible explanations for these effects are included in this article. In both exploration core drilling and overcore drilling to determine the stress in rock, occasional sections of the core obtained are broken into a series of discs.I* This discing phenomenon has been investigated in the laboratory to ascertain the state of stress required to produce it. The objective of the study was to obtain information that would make it possible to approximate the state of stress in subsurface rock from which disced cores are obtained. Referred to as the discing test, this laboratory study was conducted to determine the applied axial and radial loads required to produce discing in cylinders of rock from several sources. The dimensions of these cylinders are given in Fig. 1. Discing was produced by first applying an axial load to the cylinder and then diamond drilling the core along its axis. The radial load was then incrementally increased until discing initiated. The discs fractured or ruptured* on surfaces approximately normal to the axis of the core (planes A, B, C, D, E, Fig. 2). In many instances, the bottommost disc fractured before the drill had advanced past the points F and F', cre- ating a mush room-shaped fracture surface FEF'. Alternatively, it was also found that the same shaped fracture surface could be created by first drilling in an unloaded core to some point, say F and F', and then triaxially loading the core until fracture occurred on the plane FEF'. Jaeger and cook4 produced similar fracture surfaces in specimens subjected only to a radial biaxial load. From these laboratory studies and observations made in the field, it can be inferred that fracture or rupture initiates in the proximity of the lowest point in the kerf area, (at F and F', Fig. 2) and that the fracture or rupture surface progresses inward from this point. The experiments described above, however, do not provide any information as to the state of stress at the point at which fracture or rupture initiates. In the investigation discussed in the following sections, a three-dimensional photoelastic analysis was conducted to determine the stress distribution in the proximity of the kerf area of a triaxially loaded plastic (epoxy) cylinder containing a short section of core.

Jan 1, 1969

Jan 1, 1969

By G. E. Perry, R. D. Seba

The steam soak process is the most widely applied and most successful thermal supplemental recovery process in use today. This process, which consists of injection of steam in various quantities into a well that is subsequently used as a producer, can be repeatedly applied to the same well. The treatment, if successful, will increase the production rate of the well for a given period. For this reason, the steam soak process should be considered primarily as an acceleration method rather than as a method for producing oil that could not be recovered by primary means. However, application of this process to reservoirs containing viscous crude may produce oil that otherwise would have zero present value or would be produced at such low rates by primary methods that it would not be profitable. A major necessity for the successful application of the steam soak process is a source of natural energy. This source of energy may exist in one or more of the following forms: expansion of fluids, compaction, gravity drainage. This study was undertaken to determine the effect of gravity drainage on steam soak production and recovery and to develop a method for predicting production rates, ultimate recovery, and heat requirements when applying this process to a reservoir in which the primary driving mechanism is gravity drainage. (For additional discussion of some of the factors that affect recovery of viscous oils by the steam soak process, see Ref. 5.) Model The model developed for analyzing the performance of the steam soak process in gravity drainage reservoirs is a modified form of that developed by Matthews and Lefkovits.1,2 The model used in this study consists of two concentric cylindrical volumes of productive formation with cap and base rock of impermeable material. The inner cylinder, extending from the radius of the wellbore to some predetermined radius, is referred to as the hot zone; and the outer cylinder extending from the exterior boundary of the hot zone to the outer drainage boundary will be referred to as the partially heated zone. Fig. 1 is a schematic cross-section of the model, indicating its thermal properties. Horizontally the reservoir is divided into two regions, the upper region having a high gas saturation and residual liquid saturation, and the lower region containing a high liquid saturation. The boundary between these two regions is referred to as the free surface, which drops and changes shape as fluid is produced from the reservoir. The model was used to analyze the steam soak process by applying the following assumptions: 1. The horizontal formation is composed of one or more homogeneous noninterconnected layers containing incompressible fluids. Each layer is treated separately. 2. The hot zone is heated by the initial soak and is maintained at the steam temperature corresponding to atmospheric pressure by resoaking. 3. The height of the free surface at the producing well is at all times zero (producing well is pumped off). 4. The position of the free surface at the start of

Jan 1, 1970

By J. C. Martin

Equations for three-phase, three-dimensional, compressible flow (including capillarity) are reduced to two-dimensional relations by a partial integration. This reduction allows three-dimensional flow problems to be treated with mathematics for only two spatial dimensions. The results can be used to formulate flow equations for two-dimensional reservoir simulators in which the effects of capillarity and fluid segregation in the third dimension are represented. Such reservoir simulators would retain many of the aduantages of two-dimensional simulators while simulating three-dimensional effects. The principal restriction of the method is that the thickness of the reservoir should be small, compared to the distance across the reservoir. INTRODUCTION In recent years, computers have been used to calculate the performances of many reservoirs. Most of the detailed calculations, however, are based on finite difference solutions of the flow equations, and present day computers are seldom able to handle a sufficient number of cells to produce entirely satisfactory solutions, even for reservoirs represented by two-dimensional arrays of cells. The situation becomes much worse when one wishes to approximate the reservoir by a three-dimensional array. A great saving in computation or a more detailed solution can be obtained for many reservoirs by using the partial integration of the equations of flow presented in this paper. The integration reduces the three-dimensional equations to two-dimensional relations; and, for studies of two-dimensional flows in vertical cross-sections, the equations are reduced to one-dimensional relations. Most reservoir performance calculations currently are based on one- or two-dimensional flow re1ations.l-3 In some cases flow in the third dimension is approximated by assuming a particular type of vertical saturation distribution, such as gravity segregation.l.2. The relations developed in this paper approach those for segregated flow as the capillary pressures approach zero, and they approach those for uniformly distributed saturations as the capillary pressures are increased. For this analysis, the ratio of the reservoir's thickness to the maximum distance across it must be small. The capillary pressures between the oil and water should also be small compared to the maximum pressure difference in the reservoir. This requirement is met by most reservoirs. It is assumed that the capillary-pressure curves are well defined, whether or not hysteresis effects are included. Also, the reservoir must have sufficient vertical permeability to allow the fluids to segregate. The results presented here provide a firm theoretical foundation to Coats' et aL5 assumption of vertical equilibrium and extend the relations to three-phase flow. Coats' assumption of vertical equilibrium, which he verified by calculations and experiment, is developed here mathematically from basic flow equations. DISCUSSION SATURATION AND PRESSURE DISTRIBUTIONS Appendix A presents a mathematical analysis of fluid flow in reservoirs where the ratio of thickness to maximum distance across the reservoir is small. The results indicate: (1) that the fluids along any line perpendicular to such a reservoir's upper surface are in static capillary equilibrium (vertical equilibrium); (2) that, to a first approximation, the fluid pressures and properties are functions of only areal position in the reservoir and time; and (3) that hydrostatic pressure gradients exist along any line perpendicular to the reservoir's upper surface. The results might be expected after studying several physical considerations. First, no flow is allowed normal to the upper and lower reservoir boundaries, which are relatively close together.

Jan 1, 1969

By R. N. Edmondson, L. Mair, M. Tenenbaum

It is the policy of The Metallurgical Society to provide, in the TRANSACTIONS OF THE METALLURGICAL SOCIETY OF AIME, a prompt and accurate medium for publication of reports of significant new research in both the scientific and engineering aspects of metallurgy. Papers submitted to the Society are judged by competent reviewers according to established criteria for technical merit, accuracy, originality, and presentation. STAFF: Editor, Dr. Gerhord Derge; Editorial Assistant, M. A. Redmyrski; Staff Coordinator, James J. Burke. THE METALLURGICAL SOCIETY: W. R. Hibbard, Jr., President; J. C. Kinnear, Jr., Past-President; John Chipmon, Vice President; T. D. Jones, Treosurer; R. W. Shearman, Secretary. PUBLICATIONS COMMITTEE OF THE METALLURGICAL SOCIETY: A. w. Thornton, Chairman; Dennis Carney, John Chipman, J. C. Fulton I. H. Hollomon H. H. Kellogg, A. E. Lee, R. Maddin, J. D. Sullivan, and D. Swan. Review and selection of manuscripts is under the supenision of the following Publications Committees: Extractive Metallurgy Division: T. 0. King Chairman; S. J. Dickenson, R. P. Hermsdorf R. Schuhmann. Jr.. and P. T. Stroup. Institute of Metals Division: Robert Maddin, Chairman F. Lincoln Vogel Jr., Vice-chairman, W. A. Baikofen, R. W. Boliuffi J. H. Bechtold M. B. Beyer. R S. Busk, C. G. Gaetzel, Paul Gordon C. D. Groham, Jr., Robert F. Hehe-mann,' R. I. Joffee, F. c. Joumot Jr lrving R. Kramer A. J. Lena Harold Morglin, C. J. McHargue, Russell A: Meussner, J. B. Newkirk H. W. Paxton, G. E. Pellissier, I. S, Servi, Alon U. Seybolt, A. J. Shaler, and Jack Washburn. Iron and Steel Division: J. C. Fulton, Chairman; T. B King, F. C. Langenberg, W. 0.. Philbrook, and T. 0. Winkler. Published bimonthly by the American Institute of Mining Metallurgical and Petroleum Engineers Inc., 29 West 39th Street, New York 18, N. Y. Telephon: pEnnsvlvanio 6-9220. Subscription 522 Per year for non-AIME members in United States and North South and Central America; $25 foreign;, $5 for AIME members. Single copies, 54.00, single copies foreign, $5.00. The AIME is not responsible for any statement made or opinion expressed In its publications . . . Copyright 1958 by the American lnstltute of Mining Metallurgical and Petroleum Engineers Inc. . . . Registered cable address AIME New York , . . lndexed in Engineering Index Industrial Arts Index and Chemical Abstrats . . . Application for second class mail privileges ii pending at New York, N. Y. Volume 212, Number 4.

Jan 1, 1959

By D. W. Smith, Robert F. Mehl

It has been shown by numerous studies that the orientations of new metal crystals are determined by the orientations of the crystals in the original matrix, whether these new crystals are formed by recrystalliza-tion, by precipitation from solid solutions, or by transitions in the solid state. The orientations of crystals formed by recrystallization of cold-worked metals that bear preferred orientations are occasionally simply related to the initial preferred orientation; often, however, the relationship is not simple, but complicated and obscure in origin.1 The orientations of crystals formed by precipitation from solid solutions have been shown by many studies in this and other countries to be determined by the orientation of the crystals of the solid solution stable at the higher temperatures, though here again complicated cases have been found and no general theory is available.2 Related to this latter class are the relationships found in crystals formed by suppressed transforma-

Jan 1, 1935

By D. W. Smith, Robert F. Mehl

It has been shown by numerous studies that the orientations of new metal crystals are determined by the orientations of the crystals in the original matrix, whether these new crystals are formed by recrystalliza-tion, by precipitation from solid solutions, or by transitions in the solid state. The orientations of crystals formed by recrystallization of cold-worked metals that bear preferred orientations are occasionally simply related to the initial preferred orientation; often, however, the relationship is not simple, but complicated and obscure in origin.1 The orientations of crystals formed by precipitation from solid solutions have been shown by many studies in this and other countries to be determined by the orientation of the crystals of the solid solution stable at the higher temperatures, though here again complicated cases have been found and no general theory is available.2 Related to this latter class are the relationships found in crystals formed by suppressed transforma-

Jan 1, 1935

By E. F. Johnson, V. O. Naumann, D. P. Bossler

A method is presented for calculating individual gas and oil or water and oil relative permeabilities from data obtained during a gas drive or a waterflood experiment performed on a linear porous body. The method has been tested and found both rapid and reliable for normal-sized core samples. INTRODUCTION Individual oil and gas or oil and water relative permeabilities are required for a number of reservoir engineering applications. Chief among these is the evaluation of oil displacement under conditions where gravitational effects are significant, such as a water drive or crestal gas injection in a steeply dipping oil reservoir. Numerous proposed methods of obtaining relative permeability data on reservoir core samples have been too tedious and time consuming for practical use, or have yielded questionable and sometimes inconsistent results. A method has been developed by which the individual relative permeability curves can be calculated from data collected during a displacement test. The method is based on sound theoretical considerations. Using this method, with a properly designed experimental procedure, relative permeability curves can be obtained using core samples of normal size (i.e., 2 to 3 in. in length and 1 to 2 in. in diameter) within a few days after receipt of the core. In a recent publication D. A. Efros1 describes an approach to the calculation of individual relative permeabilities that is based on the same theoretical considerations. We believe the approach described in the present paper is more adaptable to practical application than the method implied by Efros. In addition, comparisons with independently determined relative permeabilities are furnished to substantiate the reliability of the new method. DERIVATION Previously the theory of Buckley and Leverett as extended by Welge2 as been used to calculate the ratio of relative permeabilities. In the derivation which follows, this theory is further extended to permit calculation of the individual relative permeabilities. The theory assumes two conditions which must be achieved before the method is applicable. They are that the flow velocity be high enough to achieve what has been termed stabilized displacement,3 and that the flow velocity is constant at all cross sections of the linear porous body. In stabilized displacement the flowing pressure gradient is high compared with the capillary pressure difference between the flowing phases. The high pressure drop insures that the portion of the core in which capillary effects predominate will be compressed to a negligibly small fraction of the total pore space. The assumption of constant flow velocity at all

By W. C. Peters, M. W. Pullen, C. A. Bays

During the past three and a half years considerable improvement in the techniques of solution extraction of salt has been made by the use of wells which are cross-connected by hydraulic fracture at the base of the salt. Fracture is initiated by the application of hydraulic pressure. Underreaming and perforating techniques are used to control initiation. The fracture is then extended by continuous pumping until a pressure communication is obtained. Pressure communication is followed by washing-through to obtain a low-pressure connection. Well construction using non-shrinking cements that will bond with the salt and with associated strata is essential to the process. The connecting procedure described permits solution to attack the salt deposit from its base upward, thereby permitting greater recovery, fewer well repairs, better saturation of brine, and overall economy of development and operation. The fracture method of cross-connection requires special pumping of the fluid and is limited by present experience to initial well spacings up to about 1200 feet. The method has been used in bedded salts and is applicable for wells in salt domes. Since spring of 1956, new developnents in salt extraction wells have followed procedures substantially different from those in prior use. Although multiple well systems, connected either by the natural consequence of the solution process or by undercutting under protective pads, have been used previously, the past three and a half years have brought about new designs and new operation procedures based on cross-connection of the wells using hydraulic fracture. To date these techniques have been applied largely in the bedded salt deposits of Michigan, Ontario, Ohio, West Virginia, and New York. To the writers' knowledge they have not been used in salt dome wells. The methods have not yet been in use for a sufficient period to permit any statistical comparison with the older methods. It is evident from experience to date that the methods result in higher recoveries, lower capital and operating costs, better saturation of produced brine, reduced well maintenance and repair, and in considerably fewer production units to insure

Jan 1, 1961

By P. K. Leung, S. R. Allada, P. M. Dranchuk, D. Quon

The alternating direction explicit procedure (ADEP) makes use of the boundary conditions to reduce multi-dimensional problems to a series of one-dimensional problems. The method, previously applied to reservoirs containing only an under-saturated oil, has now been extended to cover the case of natural gas reservoirs. l his involves solving a non-linear partial differential equation, application of the procedure is straightforward and no calcuIational problems were encountered. INTRODUCTION There has been a growing interest in formulating mathematical models of petroleum and natural gas reservoirs — models which permit the engineer to examine and evaluate the physical and economic consequences of various alternative production policies. The tremendous reduction in the cost of solving such models in recent years has made possible their use as an almost routine management tool. This reduction is the result not only of improved computer hardware but also of the development of more efficient mathematical techniques. The present paper is concerned with the application of a recently proposed numerical method (ADEP) to two-dimensional gas reservoirs. STATEMENT OF THE PROBLEM Given a two-dimensional Region R, bounded by a closed Curve C (Fig. 1) such that the behavior on Curve C is known, the differential mass balance for each fluid phase in R, neglecting gravitational effects and assuming Darcy's law for fluid transport, can be written as:

Jan 1, 1967

By N. A. D. Parlee, M. M. A. El-Naggar, G. 8. Horsley

The apparatus was of a modified Sieverts type1&apos;2 with a reaction tube designed to function also as the cell component of an oxygen pressure gage. The reaction tube assembly is shown in Fig. 1 and described in Table I. The main body is a 1-in.-OD 15 mol pct calcia-stabilized zirconia tube. The ground-glass joint, g, the male part of which carried a prism for optical pyrometer use, was connected to the remainder of the Sieverts apparatus. The calibration of the cell (at its 1000°K operating temperature) for internal oxygen pressures (PO2in) in the range 0.2 to 10 9 atm was made against internal air pressures measured by mercury manometer and McLeod gage. At very low po2in values down to l<Tn atm the performance was checked by comparing pQ calculated from electromotive force measurement when the apparatus was filled with a) pure CO gas and b) a CO-CO2 mixture of the ratio CO-CO2 = 1:0.398, with Po, values calculated from the equilibrium constants of the following reactions at 1000°K:

Jan 1, 1968

[Executive J. V. W. REYNDERS, Chairman WILLIAM H. BASSETT E. DEGOLYER I-I. C. BELLINGER RICHARD PETERS, JR. Finance GEORGE OTIS SMITH, Chairman H. A. GUESS F. JULIUS FOHS Admissions GEORGE D. BARRON, Chairman DAVID LEVINGER P. A. MOSMAN H. K. MASTERS J. A, SINGMASTER Alternates F. G. BREYER C. W. HAMILTON F. E. PIERCE]

Jan 1, 1928

By W. A. Burns

Vertical flow is an important mechanism in many petroleum reservoirs. Yet no adequate method has heretofore been proposed for determining the in-situ vertical permeability in a reasonable length of time. Core measurements of vertical permeability do not necessarily represent effective in-situ values. Cores may be representative of only localized pockets of high or low permeability, and total reliance on cores might cause one to overlook the presence of vertical fractures altogether. Conventional single-well tests, interference tests, and pulse tests commonly are used to determine horizontal transmissibility, khh/u, but generally are not designed to measure vertical trans-missibility or other groups containing vertical permeability, ku. A new well test described here provides a measure of in-situ vertical permeability as well as the previously measurable horizontal permeability. This "vertical test" determines the vertical diffusivity, kv/gcu, and horizontal transmissibility, khd1/u (where dl is the length of injection or production perforations). The test is basically a short-term interference or pulse test normally lasting less than a day. Conventional down-hole equipment is used between two isolated sets of perforations in the same well, as shown in Fig. 1. A "layer-cake" description of a reservoir can be obtained by running two or more vertical tests in the same well, each over different vertical intervals. Each test can be analyzed for the effective vertical and horizontal properties of the interval tested. In addition, a vertical well test can determine the effective- ness of a shale barrier or tight zone by testing across it. Field use to date of vertical testing has confirmed its utility in reservoir description. Mathematical Basis Interpretation of vertical tests for ku and kh is based on comparison of measured pressure response with computer-generated response. Pressure behavior during a vertical test in a homogeneous, anisotropic, infinite-acting reservoir (infinite in radial extent and bounded in both directions vertically) is given by the diffusion equation for slightly compressible fluids with negligible gravity effects.l A schematic diagram of the problem is presented in Fig. 2, and an analytic solution, based on the image

Jan 1, 1970

By A. E. Anderson, H. J. Hilll

Laboratory data for 155 field muds and 77 .shales have been used to develop correlations to estimate the net streaming potential component of the SP. Analysis of these data and comparison with field tests show that streaming potential corrections must be made to obtain an accurate estimate of the electrochemical potential for use in the determinntion of formation water resi-iiviy. The self-potential curve is related to formation water resistivity, and considerable effort within the oil industry has been directed toward clarifying this relationship. The Schlumbergers and Leonardon1,2 recognized that two types of potentials, electrochemical and electro-kinetic, resulted in the net potential variation in the borehole. Most of the later papers have dealt with the electrochemical component and its controlling parameters. The investigations of Wyllie,2" Salisch&apos; and Moore5 provide the most extensive sources of data on streaming potentials across mud cakes. These authors observed streaming potentials of sufficient magnitude to be of critical importance in SP interpretation, if it is assumed that streaming potentials develop only opposite permeable beds. Streaming potentials across shales have been reported by Schenck and by Gondouin, et al.7 Since SP magnitude is the difference in potential opposite shale and sand beds, these authors suggested that the net streaming potential component may be negligible. The present investigation was initiated to evaluate empirically some of the parameters controlling the magnitude of the net borehole streaming potential and to determine whether useful relations between mud properties and net streaming potential could be established. Field muds were used in this investigation. These muds were characterized by the amount and type of chemical treatment utilized for five days before Sampling. Filter loss, mud resistivity, mud density and filtrate resistivity were determined for each mud with conventional methods and apparatus. Mud filtrate pH was estimated with Hydrion pH paper for the appropriate range. Streaming PotentIals Across MuD Cakes Commercial calomel-saturated KC1 electrodes are subject to potential changes attributed to contamination of the very small fiber-tip bridge by particles of clay from the drilling mud.3 To minimize this difficulty, a large diameter bridge (Fig. 1) was constructed with saturated KC1-agar gel and a filter disc of microporous porcelain. Experiments showed that streaming potentials across this bridge are less than 0.2 mv at pressure dif-

By Frank Rieber

THE study of earthquakes long ago developed the fact that by studying the travel times of the various groups of waves from the same earthquake, as received on seismographs at varying distances, major discontinuities in the earth&apos;s crust could be reasonably well defined. Using this original discovery, and adapting the mathematical methods originally devised in seismology for determining the depth below the surface at which such discontinuities occur, other investigators were able to develop a new branch of seismology which promises to be extremely useful in extending the boundaries of our information on local conditions below the earth&apos;s surface. As an outstanding example of the applications of this work, we have the so-called "seismograph" exploration of the Gulf Coastal region, which has resulted in the discovery of a large number of salt domes, from many of which oil fields may be developed. The original seismological investigations, from which different major discontinuities in the earth&apos;s crust were deduced, as well as these later successful experiments in locating salt-dome structures, involve relatively simple problems. In each case we are dealing with surface materials in which elastic waves will be propagated with a fairly definite velocity. Under this material we have something quite different, in which the elastic waves will be propagated with a much higher velocity. The transition from one medium to the other is likewise relatively abrupt. For example, in the Gulf Coastal region the sedimentary formations have a velocity of 6000 or 7000 ft. per sec., while the velocity of elastic waves in the cap rock of a salt dome may exceed 15,000 ft. per sec. Under such conditions, delineation of the formation is a relatively simple matter, and, furthermore, there is very little opportunity for ambiguous results.

Jan 1, 1928

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.&apos; 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