Search Documents

By R. D. Carter

Methods are presented for calculating the performance of multiwell gas reservoirs. These methods account for two-dimensional, unsteady-state flow of a non-ideal gas through a heterogeneous formation. The performance can be obtained for gas reservoirs with known physical properties and configuration subjected to specified development and production schedules. These methods may prove helpful in optimizing the exploitation of gas reservoirs. The numerical solrrtion procedures given for this non-linear problem are extensions of known techniques for the linear diffusion equation. Two of these procedures are based on a new difference equation given by Saul'ev. Limitations and advantages of these procedures are discussed. An example solution is presented and is used to compare the calculation methods. INTRODUCTION Li planning the operation of a multiwell gas reservoir, it is desirable to predict its performance for each alternate development and producing schedule under consideration. This paper presents three numerical procedures that can be used to calculate such gas reservoir performance using a digital computer. These methods account for two-dimensional transient flow of a non-ideal gas within the reservoir. Extension of these methods to three-dimensional flow is possible. The paper is divided into three parts. The first describes the problem for which a solution is wanted and three numerical procedures that can be used to obtain solutions. The second discusses the information necessary to set up a problem. This part will be of greatest interest to the field engineer. The third gives data and results for an example problem. PROBLEM AND SOLUTION METHODS In the following development, a gas reservoir is regarded as a sealed, porous, permeable, heterogeneous body in which a single, non-ideal gas phase is flowing. Flow is permitted in the x- and y-horizontal directions, but is considered to be negligible in the vertical or z-direction. The reservoir can be tapped by one or more wells. Methods to be described provide a basis for calculating the complete performance of the reservoir, including rates and cumulative production from the wells, and pressure distribution and decline in the reservoir. If the isothermal flow of a non-ideal gas obeying Darcy's Law occurs in two dimensions in a region R, the governing equation in this region is:

Jan 1, 1967

By G. Y. Chin, W. F. Hosford

It has long been recognized that, for a crystal to undergo an arbitrary shape change by crystallographic shear, at least five independent shear systems must be activated. Taylor1 analyzed the deformation of fcc crystals by {111}<110> slip and showed that the principle of minimum work could be used to determine which combination of slip systems might act to produce a given shape change. Later Bishop and Hi112 attacked the same problem by first identifying the state of stress which could physically activate five or more {111)(110) systems simultaneously and then invoking the principle of maximum virtual work to determine which combination would act for a given shape change. Taylor&apos;s general approach has been used to calculate the active slip systems and the relative strengths during axisymmetric flow of cubic crystals with other deformation modes, including {112}<111> slip3 and {111}<112> twinning.&apos; It is the purpose of this note to identify the stress states capable of activating simultaneously five or more of these systems. {111}<112> Twinning. The {111}<112> systems are defined in Table I and illustrated in Fig. 1. The sense of the shears are those for twinning in fcc crystals. The shear stress on any system may be written in terms of the stress components expressed relative to the 1,2,3 cube axes. For example, the stress For five or more of the systems ij to be activated simultaneously, it is necessary that Tij/T be equal to unity in at least five of the expressions [3], while Tij < T for the remainder. (These restrictions assume that the same shear stress T is required to activate each system and that the shear stress on any system cannot exceed that for its activation.) The discrete values of the stress components A through H which can satisfy these restrictions are listed in Table II together with the values of Tij/T for all systems. It may be noted that Tij/ T is not necessarily zero on all nonactive systems as found by Bishop and Hill for <110> {111} slip. In fact, for the case of twinning Tij/T may reach negative values as large as -2 for shear stress in the reverse direction, Table 11. This means that, if those stress states are to be valid, the shear stress to activate reverse twinning must be greater than twice that for forward twinning. Such stress states are probably valid in view of the directionality of the twinning process. There has been no documentation of reverse twinning thus far. It may be noted that the stress states fall into four classes, which may simultaneously activate eight, six, six, and five twin systems respectively. These solutions may be applied to the {112}(111)

Jan 1, 1970

By A. J. Miller, G. R. Brophy

THE Grossmann principle1 for the calculation of hardenability of steel from composition is attractive because of its simplicity. It postulates that the hardenability of a steel for any particular grain size may be expressed in terms of an ideal critical diameter, which is the product of an ideal critical diameter for a "pure" iron-carbon base multiplied by a succession of factors for each element of composition contained, including the incidentals. It infers that the multiplying factor for a given quantity of an element is constant for all combinations of composition. Grossmann has stated that in the great majority of cases the experimental values of hardenability are found to be well within 10 per cent of the calculated values. The validity of the scheme has been confirmed in principle by Crafts and Lamont2,3 and by Kramer et al .4 Their results differ from the original in the magnitude of the effect of some elements but, since they accepted the product concept and followed its direction, it is to be expected that similar results would be obtained. Crafts and Lamont, in their second paper, recognized that the factor relations for manganese and nickel were not continuously linear, but beyond certain critical amounts increase sharply. Austin, Van Note and Prater5 demonstrated nonlinear relations between hardenability and alloy contents. Their curves for chromium, silicon, manganese and nickel were all convex upward to show decreasing incremental effects. They concluded that Grossmann&apos;s factors underevaluate the hardening effect of low percentage additions of these elements, and overevaluate the larger additions. [ ] More recently W. Steven6 has found that the factor method fails when applied to complex steels containing two or three principal alloying elements other than carbon. Multiplying factors derived from his results bear a complex relation to composition and are not constant for a given amount of an element as the base composition varies. C. A. Liedholm7 found that, using either the Grossmann or the Crafts-Lamont multiplying factor values, not only were

Jan 1, 1945

For the third year in a row total value of the Nation&apos;s mineral production broke all previous records. A new peak of $20.4 billion was reached in 1964. So says the U. S. Bureau of Mines. The Bureau&apos;s preliminary estimate for 1964 is $800 million higher than the value reported for 1963, and $1.6 billion greater than the 1962 figure. The unprecedented demand for minerals, has stimulated development of new mines and expansion of processing capacity in many parts of the country.

Jan 2, 1965

By W. David Leyshon

This paper compares the general characteristics of Florida phosphate rock for the manufacture of fertilizers and phosphoric acid. The subjects covered are relative reactivity for superphosphate manufacture and various properties of the rock for use in phosphoric acid production, including gypsum filterability, sulfuric acid consumption, foaming, relative cor-rosiveness, chemical composition, and others. Comprison will include rocks of North Africa, the western U.S., and other locations. This paper reviews some of the important properties of phosphate rock, giving particular attention to a comparison of Florida phosphate with some of the major phosphates produced in the world. It also discusses some of the techniques available for evaluatine Phosphate. rock for phosphoric. acid and fertilizer use. The use of Florida rock outside the U.S. has been increasing steadily, evidence in itself of certain desirable properties of this rock. Table I is a list of the major producers and exporters of phosphate rock as of 1965. Florida, meaning in this case Central Florida, was the largest producer of phosphate rock and the largest exporter of rock. Of a total world production estimated at 60 million tons, about 50 million tons were produced during 1965 by the areas listed in Table I. One major producing area omitted from Table I and from subsequent evaluations in this paper comprises Nauru, Ocean, and Christmas Islands, whose output totaled about three million tons in 1965. This phosphate is marketed entirely by the British Phosphate Commissioners, primarily in Australia and New Zealand, and does not get into the world market elsewhere to a significant degree at the present time. South Africa, Egypt, Israel, and North Vietnam also produce lesser quantities of rock for the international market, and these materials are also excluded from this review. The phosphate industry is now seeing the introduction of several new areas that will be significant factors in the rock market in the future. The North Carolina and North Florida areas each have capacity at present or in the near future for producing about three million tons per year. The Djebel Onk mine in Algeria has a potential of about 800,000 tons per year. In addition, announcements have been made of significant development of deposits in Peru and in the Spanish Sahara, and sizable expansions are planned in several of the existing areas that are now of minor importance. Competition for the rock market appears to be increasing so that more attention will probably be paid to evaluating phosphates available from different areas. It is not possible here to go into all the properties of even a few of the materials listed in Table I. However, it is practical to discuss some comparative properties important for making phosphoric acid and fertilizers, and to touch on the techniques used for the evaluation of various phosphate rocks. Table II gives typical analyses of the major phosphate rocks available in the world as of 1965. There is, of course, an inherent danger in labeling any analysis as typical, but it appears that variations between two samples of rock within one given area are generally less than between two samples from different areas. The analyses listed were compiled from several years of data so they may not be representative of the present output of any particular area. The range of grades represented by these analyses goes from 66 BPL for Gafsa rock to 85 BPL for Kola

Jan 1, 1968

By H. Dykstra

Calculated wellbore pressures were obtained for parameters of radilcs ratio and permeability. In all cases bur two, after-production was allowed to occur for one day. The calculated pressure build-up data were compared with actual pressure build-up data from a condensate well. The paper discusses the comparison, gives reasons for, inability to match calculated (2nd actual data, and presents conc1usions derived from the comparison. INTRODUCTION Considerable information has been published on the decline in wellbore pressure for oil and gas fields resulting from production of fluids. Considerable information also has been published on pressure build-up for oil fields, but relatively little has been published for gas and condensate fields. Perrine&apos; summarized the equations derived for pressure build-up in oil wells and showed how the different methods could be applied. On the other hand, very little theoretical information has been published on pressure build-up in gas or gas-condensate wells.2-4 Tracy,&apos; by pointing out the similarity between the equation describing oil wellbore pressure decline and the equation describing gas wellbore pressure decline, showed how methods of pressure build-up analysis could be applied to a gas well. In analyzing pressure build-up data for gas or gas-condensate wells, it would be very desirable to compare actual pressure build-up data with calculated pressure build-up data. Such a comparison would lead to a better feeling for the quantitative picture of gas-well build-up curves and would help to establish a degree of confidence in a given build-up curve. This paper discusses the comparison and presents conclusions made from it. Calculated wellbore pressures are given for parameters of radius ratio and permeability. In all cases but two, after-production was allowed to occur for one day. The work was done in an attempt to evaluate the reserves of a gas-condensate field. As will be discussed later, it was not possible to make an evaluation. METHOD OF CALCULATION The method used for solving the depletion and pressure build-up behavior for a gas-condensate well was that developed by West, Garvin and Sheldon.5 or the IBM 704 computer program used, it was assumed that liquid condensing out of solution would have only a minor effect on the flow behavior. Therefore, the field was treated as having single-phase flow, with the gas saturation S, remaining constant at unity minus the connate-water saturation Sw. As will be discussed later, this may not have been a good assumption. It also was assumed that oil production could be considered as equivalent gas production by using a conversion factor based on an analysis of the liquid produced. An outer boundary condition of no-flow was assumed. BASIC DATA The data required for the study were reservoir properties, fluid properties and production data. Reservoir properties included the following: thickness h, 179 ft; porosity, 0.194; connate-water saturation, 0.43; gas saturation, 0.57; permeabilities k, 1.0, 0.5, 0.25 and 0.15 md; radius ratios re/rw, 800, 1,000 and 2,000; wellbore radius rw, 5 in.; initial reservoir pressure P, 6,529 psia; and reservoir temperature, 272°F. The selection of permeabilities was based on core-analysis data showing an average air permeability of about 2 md. With a relatively high connate-water saturation, the average effective gas permeability would be about one-half, or less, of the air permeability. The selection of radius ratios was based partly on the observed rapid decline in wellbore pressure during production and partly on the fact that the well was believed to be located in a relatively small fault block. Reservoir fluid properties were obtained from a reservoir fluid study on a recombined sample of gas and condensate. Gas formation-volume factors were calculated from pressure-volume measurements made at reservoir temperature. Gas viscosities were calculated by the method of Carr, et al,G from an analysis of the re-combined sample. These data are shown as a function of pressure in Table 1. The production schedule from the gas-condensate

By R. P. Smith

IN 1937 the employees of the Hollinger Consolidated Gold Mines Ltd., at Timmins, Ont., Canada, approached organized medicine for a plan to provide themselves and their families with a complete medical service including hospitalization. This approach was made because of two conditions: (1) the necessity for some organization whereby contributions could be made in instalments to offset the sudden and large expense connected with illness for which the average employee is not prepared, and (2) the need for a more cohesive quality of service to replace a contract practice then in vogue, which did not take care of the employee&apos;s family or dependents. In different parts of the United States and Canada similar plans were in force, not only for employed groups but for numbers of people in circumscribed areas. In Europe different countries had adopted various forms of health insurance plans. It then was no surprise that these employees as a group should ask for some similar plan. With the cooperation of the management, a plan was instituted. The plan has one object only-the provision of medical services, including the accessory services such as hospitalization, X-rays, nursing, and necessary drugs and medicines. It is not concerned with cash benefits, sickness and accident payments, or any form of monetary return as the result of illness, but is simply an organization to provide medical services and pay for them from a central fund. It is a nonprofit organization and the only shareholders are the employees and their families who contribute and who in turn receive services. The Hollinger plan is limited to the employees of the Hollinger mine, who live in Timmins and the immediate vicinity. At first, the membership was voluntary. In the third year, it was made a condition of employment, with the proviso that at the end of any Association financial year an employee may withdraw if he wishes to do so. After withdrawal, he cannot again become a member until the end of the next Association year, and then only by the payment of the back dues. This step was taken to ensure continuity, and to prevent employees from withdrawing when well and reentering when faced with impending illness or operations. The employee himself is subject to medical examination on application for employment, and if accepted is subject to further examination by a specialist in chest diseases under the Workmen&apos;s Compensation Board of Ontario. Thus it is possible to exclude those who have heart disease, tuberculosis, or other disabling conditions when first employed. Wives and families of employees are not subject to medical examination but other dependents of an employee-mother, father, or other relative-who are declared by the employee to be wholly dependent upon him must be examined. The employee

Jan 1, 1944

By C. S. Matthews, P. Hazebroek, F. Brons

A method has been developed for calculating the average pressure in a bounded reservoir. The reservoir is first divided into the individual drainage volumes of each well, by using the criterion that at steady state each individual drainage volume is proportional to a well&apos;s production rate. The average pressure in each drainage volume is then calculated by a method developed in the report. By volumetrically averaging these individual drainage volume pressures, the average pressure in the entire reservoir is obtained. To calculate the average pressure in each drainage volume, a correction is applied to the ordinary extrapolated pressure, i.e., the pressure obtained by extrapolating to infinite time the linear portion of the graph of closed-in pressure versus log[?t/(t + ?t)], where At is the closed-in time and t the production time. The correction, which is a function of the production time, is presented in graphical form for different shapes of the drainage area (horizontal cross section of the drainage volume). INTRODUCTION It is important to be able to find the volumetric average pressure in a reservoir so that the size of the reservoir may be determined from material balance calculations. It is also desirable to be able to find the approximate distribution of pressure within a reservoir for detection of fluid movement. &apos;The purpose of this paper is to present a method for calculating both the average reservoir pressure and the approximate distribution of pressure within a bounded reservoir — that is, a reservoir with no water drive. In reservoirs where the pressure builds up rapidly after wells are shut in, the determination of average pressure generally poses little problem, for one often need only average the final buildup pressures. It is when pressure buildup is slow that difficulties arise. For practical and economical reasons, the time allowable for closing in wells is limited. If at the maximum allowable closed-in time the pressure has not reached a con- stant value (and this is more often the case than is generally realized), calculation of average pressure presents difficulties. One method which has been used in the past for obtaining an "average" pressure involves the extrapolation to infinite closed-in time of curves of pressure as a function of closed-in time. Contour maps of these extrapolated pressures are used in arriving at an average pressure. In some cases wells are shut in for arbitrary Lengths of time, and contours of these arbitrary shut-in pressures are used for calculation. All these methods are somewhat arbitrary and suffer from lack of fundamental theoretical basis. For example, no precise meaning can be given either to the extrapolated or to the arbitrary shut-in pressure. Furthermore, there seems to be no reason why either set of pressures should be contourable. It was because of the vagueness of these old methods and their lack of a sound theoretical basis that, the present method for determining average reservoir pressure was developed. SYMBOLS A = drainage area of a well, cm&apos; A, = drainage area of well j, cm2 A , = total drainage area of all wells in a reservoir, cm2 c = coefficient of compressibility of fluid at reservoir temperature, atm-1 f = hydrocarbin-filled porosity, fraction h = net thickness of formation, cm k = permeability of formation, darcy p = pressure at a well, atm p0 = initial reservoir pressure, atm p* = pressure obtained by straight-line extrapolation of linear portion of plot of p versus log [At/(t+ At)] p = volumetric average pressure inside drainage volume of a well q = volume flow rate at a well, cc/sec a: prevailing reservoir conditions q, = total volumetric flow rate from a reservoir, cc/sec at prevailing reservoir conditions rw = radius of well bore, cm T = kt/fµcA t = corrected time of production of a well, seconds, defined by

Jan 1, 1955

By Alfred Weiss

The Daniel C. Jackling Award, established in 1953, is presented for significant contributions to technical programs in mining, geology, and geophysics. This year&apos;s recipient, Alfred Weiss, was cited "For his pioneering achievements in applying systems analysis, operations research, and computer methods in exploration, mining, and environmental control; and in disseminating the new technology to the worldwide mining community; and for his lecture, “New Concepts and Tools in Decision-Making.”

Jan 1, 1979

By Alan Lawley

AS part of a programme of experimental work on single crystals of rhenium and ruthenium, it was found necessary to have the crystallographic angles for these elements. The angles given in Table I were determined to the nearest minute of arc, by Univac Computer. The angle $ between (HKIL) and (hkil) is given by the formula: The angles given for rhenium should be particularly useful in view of the increasing interest in this element. For completeness, angular values are also given for cadmium zinc and magnesium; these form a more complete set of values than those given by Salkovitz.&apos; The assistance of Professor J. Hobstetter and Mr. J. Crozier in the programming of this data is acknowledged. The work is sponsored by the Office of Naval Research under contract No. NONR 551 (19).

Jan 1, 1961

By Raymond Mancha

The use of blower fans and tubina- for auxiliary face ventilation is a somewhat controversial subject. Most authorities agree upon the advantages of this system over the maintenance of line brattice for dead-end ventilation. They also recognize the system&apos;s disadvantages. However, there is considerable disagreement as to the proper relation of the advantages to the disadvantages. It is the purpose of this report to present the opinions of four authorities of diversified occupation and permit the reader to form his own opinions on the subject. ORsErvers Carel Robinson is a mining engineer consultant of wide repute, who has a background of operating experience. Thomas G. Fear is a mine operator whose present position is that of chief engineer for the Hanna Coal Co. at St. Clairsville, Ohio. E. H. Denny is chief of the Coal Mine Inspection Division of the U. S. Bureau of Mines, and is a very able engineer whose past experience well qualifies him to appraise this problem. Cloyd M. Smith, editor of the periodical Mechanization, is a ventilation consultant. For a number of years he was a member of the Mining Department of the University of Illinois and has studied mine ventilation for many years. It is Mr. Robinson,s opinion that closely spaced breakthroughs, although conducive to low gas accumulation past the last open crosscut during the driving of entries, are difficult to maintain and are constant Sources of fugitive air. Jute or canvas brattice is suitable for entry driving with

Jan 1, 1944

By Raymond Mancha

The use of blower fans and tubina- for auxiliary face ventilation is a somewhat controversial subject. Most authorities agree upon the advantages of this system over the maintenance of line brattice for dead-end ventilation. They also recognize the system&apos;s disadvantages. However, there is considerable disagreement as to the proper relation of the advantages to the disadvantages. It is the purpose of this report to present the opinions of four authorities of diversified occupation and permit the reader to form his own opinions on the subject. ORsErvers Carel Robinson is a mining engineer consultant of wide repute, who has a background of operating experience. Thomas G. Fear is a mine operator whose present position is that of chief engineer for the Hanna Coal Co. at St. Clairsville, Ohio. E. H. Denny is chief of the Coal Mine Inspection Division of the U. S. Bureau of Mines, and is a very able engineer whose past experience well qualifies him to appraise this problem. Cloyd M. Smith, editor of the periodical Mechanization, is a ventilation consultant. For a number of years he was a member of the Mining Department of the University of Illinois and has studied mine ventilation for many years. It is Mr. Robinson,s opinion that closely spaced breakthroughs, although conducive to low gas accumulation past the last open crosscut during the driving of entries, are difficult to maintain and are constant Sources of fugitive air. Jute or canvas brattice is suitable for entry driving with

Jan 1, 1944

By A. Ahmadieh, J. E. Dorn, Jack Mitchell

A detailed investigation of the disloccrtion mechanisms controlling prismatic, slip in a solid solutions of magnesium containing up to 15.9 at. pel Li revealed that low-temperature prismatic slip is controlled by the Peicrl&apos;s mechanism for all alloy compositions. At higher temperatures slip in the 7.9 at. pet Li alloy was found to follow the diclates of the cross-slip mechanism from basal to prism planes, while the higlier-lilhium alloys exhibited an athermal dejormation process which might be due COARSE-GRAINED polycrystalline aggregates of distilled magnesium exhibit very limited ductility at low temperatures.&apos; The ductility of polycrystalline magnesium, however, is much improved by solid-solution additions in excess of about 8 at. pct to shot-range ordering. The major effect of lithium additions on the low-temperature Shear] stress. for prismatic slip appears to be due to its efftct in lowering the Peierl&apos;s stress whereas (1 increase in the shear stress with increased lithium conlent at higher temperatures might he due to short-range order strengthening. Additions rip to 7.9 (at. pel Li appear to have only minor effects on the slackiug-jaidt energy. Li.2 Although the relative amounts of basal slip and twinning are not materially altered by such alloying, a substantial change takes place in the relative amounts of prismatic slip. Whereas prismatic slip is negligible and limited to regions of high stress concentrations in the vicinity of the grain boundaries in pure magnesium, extensive prismatic slip over entire grains occurs in the higher lithium content alloys. Furthermore, as the lithium content is increased above 8 at. pct, the flow stress for polycrystalline aggregates decreases. Comparison of the single-crystal data on pure magnesium for basal slip by Sheely and Nash3 and those for prismatic slip by Flynn, Mote, and Dorn4 indicates that the critical resolved shear stress for prismatic slip is much above that for basal slip.

Jan 1, 1965

By E. P. Abrahamson, S. L. Lopata

A general relalionship between the amount of glide shear (due to slip) and the macroscopic shape change has been developed. Since the deformation can be large, finite strain analysis is employed. In this treatnzent, the shape change is expressed by a deforlnation gradient matrix, F - [axi/aXj/ where X and x refer lo the initial and final positions of a particle. This matrix is readily evaluated in termas of the amounl of glide shear, a , and the direction cosines of the slip-plane normal and the slip direction for a single active slip systena. In the case of deforma tion from slip on several systems, the product F(a) of the several deformation gradient matrices is first calculated. Then, by assuming that the final configuration is reached by a long series of small shears of magnitude a, occurring more or less alternately in the several slip sys- IN problems of crystal plasticity, it is often necessary to relate the amount of glide in the operating slip systems to the macroscopic strain components. This relationship is relatively simple when only small strains are considered. In this case the separate strain contributions from several slip systems are additive. Such a procedure is incorrect in the case of large plastic deformation. To the authors&apos; knowledge, no general treatment of the latter problem has appeared in the literature. In specific instances, Mark, Polanyi, and schmid3 have derived the relationship between glide and axial elongation during tensile pulling of a single crystal when only one slip system is active. A similar relationship for duplex slip was worked out by v. Göler and Sachs.4 Taylor and Elam1 have also studied in detail problems of large plastic deformation. Due to uncertainties as to the slip systems at the time, ! terns, the final deformation gradient matrix can be obtained. Mathematically, the problem is reduced to finding the limit of F(aN as N — m while a — 0 in such a way that the Product Na = a finite conslant designatzng the accumulated amount of shear. It turns out that this limit is simply eafF1, where F, is the matrix whose elements are the coefficients of a in the F(a) twatrzx. Application is made of the present treatment to a fcc crystal of Permalloy compressed on the (110) plane and contrained to elongate in the [112] directzon. In addition, it is shown that one rnay readily obtain from the general analyszs the well-known formulas relating elongation, amount of glide shear, and amount of lattice rotation for crystals deforming by single and double slip under tension. however, they were mainly concerned with proving that slip occurs on {111}(110) systems in fcc metals. The present general treatment arose out of recent investigations of magnetic anisotropy in cold-worked Fe-Ni alloys6 as well as some related problems in strength anisotropy of single crystals.&apos; Detailed application of the general analysis will be provided in the case of a Permalloy single crystal compressed on the (110) plane and constrained to elongate only in the [112] direction. It will be seen that, as expected, the small strain approximation leads to significant errors after moderate straining. It will also be shown that the present general treatment yields the formulas of Mark et al. and of v. Goler and Sachs in those specific cases. GENERAL CONSIDERATIONS Consider a homogeneous deformation in which a material point initially at moves to (xl,xz,x3), both positions being referred to the same set of Cartesian axes. Such a deformation can be

Jan 1, 1967

By Henry J. Welge

The approximate methods which are now in use for calculating oil displacement from reservoirs by gas-cycling or gravity-drainage at constant gas pressure, or by water flooding, make use of fundamental relationships derived by Leverett&apos; and Buckley and Leverett.2 The mathematical equations needed are derived by applying Darcy&apos;s law to the flowing phases, and by material balance considerations. ln general, any treatment of this type gives. for any particular exploitation time considered, a plot of oil saturation against distance in the reser-voir. The oil recovery must then be obtained by integrating in some manner the area under the plot. A useful analytical method has been derived for computing the average saturation, and hence the oil recovery. Use of this method simplifies the calculations because it makes unnecessary any numerical integrations. and even the saturation distribution plots are not needed. A further advantage of the method is that knowledge of the relative permeabilities is required only for a limited and intermediate saturation range. In both the Buckley and Leverett method and the method discussed here, a linear sand section is assumed, and in the case of gas drive the gas pressure is assumed sufficiently constant both with respect to reservoir position and time so that changes in gas density. solubility, or reservoir volume factor are negligiklle. Thus, the exploitation contemplates oil displacement as by an immiscible phase. Examples are given to illustrate how the new method can be used. INTRODUCTION AND THEORETICAL BACKGROUND This paper treats a simplified method for computing oil recovery when the oil is displaced, from the reservoir sand by a fluid which, within the limits of accuracy desired, can be assumed to be incompressible and immiscible with the oil. The method makes use of two basic relations originally developed for the case of water displacing oil. However, the case of gas displacing oil saturated with gas at a constant (or nearly constant) pressure may also be considered a displacement by an immiscible fluid. This is possible for the reason that the concentration of the gas in the oil never changes if tile pressure is fixed. consequently, any additional free gas must remain undissolved in the oil, and so must act essentially as an immiscible phase. For convenience, the cape of gas drive will be considered first. Both of the two basic relations needed refer to a linear reservoir, or one which may be somewhat idealized so that a constant cross-section is exposed to fluid flow. Further, the displacing phase is injected at one end, or face; of the reservoir, while the produced fluids are understood to emerge at one point only. which is the opposite face of the reservoir. As indicated in Fig. 1, the reservoir may he inclined by some angle 8. The first of the required equations may be derived by applying Darcy&apos;s law to both displacing and displaced phases, and subtracting the two equations. The result for water drive has been given previously, as Equation (13) of the reference by

Jan 1, 1952

By L. L. Handy

Two equations are compared for calculating imbibi-tion rates of water into porous media. A linear relation was found to apply when the square of the volume of water imbibed was plotted against time for water imbibing in natural rock materials. The effect of adsorbed organic materials on surface wettability was determined from measurements of rate of imbibition, water saturation and water permeability. INTRODUCTION An understanding of the mechanisms by which water displaces oil from porous media requires an understanding of the role of capillary forces in the displacement process. The capillary forces are determined by properties of the water-oil-solid surfaces. Consequently, some measure of the surface properties of reservoir rocks is necessary if displacement behavior in oil reservoirs is to be predicted from laboratory measurements. Several experimental procedures have been proposed to assign quantitative wettability indexes to reservoir rock surfaces. The most recent of these proposals are those of Bobek, Mattox and Denekas&apos; and Amott.&apos; These indexes are designed to show a continuous variation from the preferential oil-wet to the preferential water-wet systems. They require measuring some property of the rock which is a function of surface wettability. The quantities are measured on unaltered core material and compared with values obtained for known oil-wet and water-wet extremes on the same material. These methods are useful but are semi-empirical in nature. They have the disadvantage that the measured quantities may be functions of other variables in addition to surface wettability. Some understanding of the behavior of water-oil-rock systems may be obtained from studies of the simpler water-air-rock systems. Surfaces of a given wettability are more easily maintained for the latter. Also, certain approximations can be made for the flow of water and air which are not justified for water-oil flow. The theoretical and experimental work discussed in this paper is part of a study of the role of capillarity in the dynamic displacement of air by water. The expec- tation is that observations for the water-air system will lead to generalizations which are applicable to water-oil systems. THEORY Imbibition is defined as that process by which one fluid displaces another from a porous medium as a result of capillary forces only. Differential equations have been derived to describe the displacement of oil by water during imbibition, but these equations have not been solved.".&apos; With certain assumptions, however, solutions are possible for water-air imbibition.

By H. I-Lieh Huang, O. D. Sherby, J. E. Dorn

RECENT investigations1-4 have suggested that the total plastic strain, C, for high temperature creep under a given stress can be correlated by means of a temperature-compensated time, te- where t is the duration of test, e is the base of natural logarithms, aH is the activation energy for creep, R is the gas constant, and T is the absolute temperature. A typical set of data on which this deduction was based is given in Fig. la. Inasmuch as the e-in t curves for the various temperatures are identical except for their displacement along the In t axis where f and are appropriate functions. If the increase in creep rate with an increase in temperature is exclusively ascribable to a single process involv- ing thermal activation, (T) should equal RT and therefore it = constant [2] where ?= te Assuming the validity of this relationship, ? = ?2 at the same strain, or where the subscripts refer to two alternate test conditions.

Jan 1, 1957

By A. J. Teplitz, J. K. Rodgers, E. M. Pohoriles

This note is presented in order to call attention to a newly developed laboratory apparatus that appears to have interesting possibilities for the study of lost circulation of the pressure-parting type, for a more realistic evaluation of sealing materials for this type of loss, and for investigation of other aspects of hydraulic fracturing. A study of lost circulation, conducted jointly by operating divisions of the Gulf Oil Corp. and the Gulf Research & Development Co., has been in progress for the past few years. Earlier phases of the study dealt with mechanical factors contributing to loss of returns in laminar type sediments.&apos; One aspect of this type of loss which has been widely discussed in the literature and elsewhere concerns the geometry of the fractures into which such losses occur. This subject is not merely a philosophical one, but has an important bearing on procedures for combating lost circulation difficulties and on other well operations. The further course of the collaborative studies previously mentioned has produced some information on this subject which, though limited in amount, is believed to be significant. In several typical cases of lost circulation in laminated shale and sand formations, the wells have been surveyed with instruments2,3 capable of detecting the points of loss with considerable accuracy. These wells were located in the Gulf Coast of Texas and Louisiana, in the Texas Panhandle, and in the Rocky Mountain area. In each of these instances the loss was shown to be occurring at some distance removed from the bottom of the hole and sharply confined to a very narrow interval, less than 1 ft in thickness. It is rather difficult to imagine that fractures of a vertical nature&apos; would be of the limited extent indicated by the surveys in these wells. Furthermore, the high rates of loss at relatively low pressure differentials observed in some of these instances together with the narrow intervals involved would appear to require the radial flow pattern which would be associated with a more or less horizontal fracture plane. In view of the foregoing observations and records of successful squeeze cementing operations, in attempting to construct a laboratory apparatus For simulating pressure induced lost circulation, it was decided that the design should be based on a horizontal type of fracture. It was realized that5,6 Vhe production of horizontal fractures in a formation with non-penetrating fluid requires some mechanism for transmitting vertical components of the hydraulic forces to the plane of parting. This was accomplished, as will be evident in the description of the apparatus that follows, by providing recesses in the wall of the borehole which are believed to be representative of wellbore configuration found in Gulf Coast wells. The possible influence of such borehole irregularities on both pressure fracturing and the treatment of lost circulation appears to have been neglected in many studies of these subjects. It is believed that the present apparatus, incorporating features that take into account some of the factors mentioned, will provide means of gaining a further understanding of pressure induced fracturing and possibly for more realistically evaluating the numerous sealing agents, particularly the non-time-setting materials employed in treating the pressure-parting type of loss. DESCRIPTION OF APPARATUS The test cell used in the laboratory investigations is shown in Fig. 1. A schematic diagram of the entire system is shown in Fig. 2.

Jan 1, 1956

By R. S. Silo

A new process for leaching oxidized copper ore (CuO) is presented. Its principle advantage is that it requires low cost raw materials. For large operations the only materials used are sodium chloride (NaCl) and calcium carbonate (CaCO3). In those operations of medium size, sulfur, which is abundant in the North of Chile, is used. However the process may also be applied to small copper plants provided low cost sulfuric acid is available as a substitute for sulfur. That is, it may be applied to those small capacity plants where large additional installations for the regeneration of the leaching solution are not needed. The method is comprised of six stages, namely: Crushing, Leaching, Washing of the Solution, Precipitation, Regeneration of the Leaching Solution, and Fusion. A flow diagram for this process in a large or medium-sized plant is given. LEACHING This process requires hydrochloric acid as the leaching solution for treating copper ores. This is easily produced from sulfuric acid and sodium chloride as follows: A small excess of NaCl is advisable to avoid the presence of free sulfuric acid. The experiments reported here were performed with 36 grams of HC1 per liter of solution to obtain a free acidity equivalent to 50 gpl of H2SO4 (the usual leaching solution concentration). Hydrochloric acid leaching may be done by either of the usual methods; i.e. percolation or agitation. The method selected depends on the characteristics of the ore. No matter which method is chosen, low pressure air can be used either as an air lift for circulation of the solution, or for agitation of the pulp with "Pachuca" tanks. In either case, there are no problems with respect to oxidation. Leaching by agitation produces good recovery; but at times, the problem of high consumption of leaching solution or the overproduction of grinding slimes which interfere with contact between phases negates this advantage. The percolation method was used in this study. Percolation took place in four 15-liter containers, 11 1/2 in. diam. by 12 in. height, having hemispheric bottoms. Each container was fitted with a wooden bottom having 1/4 in- holes. The air lift was a 3/4 in. plastic tube 10 in. high. Air at 6 psi was supplied from a rotating laboratory compressor. 10 kg of ore were charged to the container and counter-current leaching was used. The leaching solutions were fed through the center pipe in order to use the ascending solutions to reject the air in the charge. The solution circulated through the air lift in order to facilitate the movement according to the increase in density produced by the reaction between the phases present. Specific observations made in this report cannot be generalized because the leaching characteristics are dependent on the nature of the ore used, the method of leaching, and the production of primary or secondary slimes while crushing or leaching.

Jan 1, 1967

By G. Richards Gwinn

Graphite is one .of the so-called minor nonmetals that have become of major importance during the present conflict. Requirements of the domestic industry for graphite are relatively small and uses have changed so greatly in the past 25 years that different types, with few exceptions, are more nearly interchangeable in use. However, the quality of the graphite required for crucibles is exacting and under present manufacturing practices only a small part of the world production can be utilized economically by this industry. Graphite, for the domestic crucible industry, was supplied almost entirely by imports from Madagascar and Ceylon before the present war. The quantity required has been increased by the war and the greatest demand is for Madagascar flake. The fall of France to Germany in 1941 and the Japanese expansion southward threatened our supply of crucible-grade graphite, therefore the domestic industry was revived to supplement imports and to act as insurance in the event of complete stoppage of imports. Two plants were constructed for the Government, one at Ashland, Alabama, called the Crucible Flake Graphite Co., operated by Haile Mines, Inc., and one near West Chester, Pennsylvania, operated by the Benjamin Franklin Flake Graphite Co. Loans were made by the Reconstruction Finance Corporation to the Alabama Flake Graphite Co. at Ashland, the Ceylon Graphite Co., at Goodwater, Ala., and the Southwestern Graphite Co., at Burnet, Texas. However, by the end of 1943 the import situation had improved and it was no longer necessary to continue Government assistance to mine and stock-pile the high-cost domestic material. Graphite is found in virtually every country of the world. Most occurrences are of no economic importance because of their poor quality, prohibitive mining and milling costs, inaccessibility, and poor location with reference to the industrial centers where graphite is used. The principal commercial deposits are found in Ceylon, Madagascar, Chosen (Korea), United States, U. S. S. R., Canada, Germany, and Mexico. Graphite is a typical international mineral in that it is utilized mainly in countries noted for industrial activities, whereas the two islands, Ceylon and Madagascar, which produce the major part of the world&apos;s high-grade graphite, are virtually devoid of any manufacturing establishments requiring its use. The industry thus is violently disturbed by wars. Types oF Graphite Used in Crucibles Graphite of various kinds and from many countries has been used for making crucibles. Passau (Bohemian) flake was first used in Europe, but by the middle of the seventeenth century crucible formulas were established largely on Ceylon graphite and German Klingenberg clay. Ceylon graphite is granular and the grains are subdivided into more or less

Jan 1, 1948