Search Documents

By J. C. Bokros, A. S. Schwartz

The irradiation-induced densification of pyrolytic carbon has been shown to take place by a simple first-order rate process. The values of the densification constant, kd, for isotropic carbons are 1.0 x 10-22, 1.4 X 10-22 and 2.8 X 10-22 (nut)-1 at irradiation temperatures of 520o, 670°, and 1000oC, respectively. The apparent activation energy obtained from the temperature dependence of kd for isotropic carbons is about 0.2 ev. Estimates of the maximum stresses that develop in isotropic carbons owing to shrinkage onto an unyielding support show that stresses in excess of 15,000 psi may be avoided during irradiatian at 600" and 1000°C if the initial coating density is kept above 1.5 and 1.9 g per cu cw, respectively. If a carbon coating on an unyielding support is not completely isotropic, stresses may continue to build up at high neutron exposures after densification is complete because of continued anisotropic distortion. Stresses in excess of 15,000 psi from this source may be avoided during irradiation at 1000°C if' the Bacon anisotropy factor is kept below 1.3. A variety of applications of carbon and graphite in nuclear technology has stimulated numerous studies of the effect of neutron irradiation on these materials One of the most interesting problems associated with the use of carbons and graphites in the high fast neutron fluxes of nuclear reactors arises as a consequence of the neutron-induced dimensional changes that occur during irradiation. Energetic neutrons displace car- bon atoms from their lattice positions; and, depending on the temperature, the interstitials and vacancies produced migrate to form clusters, recombine, or are annihilated at dislocations or crystallite boundaries. The competition of various sinks for the point defects results in structure-sensitive bulk-dimensional changes that vary with temperature in a complicated way. This paper deals with the neutron-induced dimensional changes in turbostratic carbons that are being used as coatings on nuclear fuel particles. Details of the structure of these carbons and the variations of their structure with deposition conditions have been reported.'-3 It has been shown in previous work that the neutron-

Jan 1, 1968

By J. Douglas, W. T. Ford, G. E. Henderson, I. F. Roebuck

An implicit numerical method is presented for simulating the differential and algebraic relations governing one-dimensional three-phase flow in porous media. The method is based upon compositional representation of the hydrocarbon system. Variable physical properties and water-oil capillary forces are included in the formulation; however, the effect of gravity is ignored. The effects of changing composition and mass transfer are considered through the use of known phase behavior concepts and correlations. The differential and algebraic equations and the numeric approximations to these equations are presented, and the computing algorithm is discussed in detail. This method is used to simulate a single-phase, two hydrocarbon-component gas displacement laboratory experiment. The method is also applied in simulating a solution gas drive and a gas injection problem using a reservoir oil represented as a nine-component mixture. For comparative purposes, the latter two problems also are examined with a one-dimensional volumetric model. The results of the simulation of the laboratory experiment demonstrate the relevance of this method to problems wherein the effects of mass transfer significantly enter the displacement process. The comparison of the results of the volumetric and compositional methods applied to a reservoir system illustrates not only the advantages and benefits derived from examining the hydrocarbon system composition, but also the necessity of this approach in application to a general reservoir problem. INTRODUCTION The advent of large memory, high-speed digital computers enabled reservoir engineers and mathematicians to pool their knowledge in the development of sophisticated techniques for the prediction of oil and gas reservoir performance. These techniques have progressed from volumetric and trend analysis calculations to the development of "mathematical models", or reservoir simulators. These latter models, in general, utilize finite difference approximations to the rather complex partial differential equations that mathematically describe the physics and thermodynamics of fluid flow in porous media. Many forms of this "model" type are available today, and much of the recent literature has been devoted to either improving known methods or developing new techniques for use in these prediction methods. The more advanced techniques in use at this time consist of the solution of differential and algebraic relations that describe a reservoir geometry in one, two or three dimensions and that incorporate reservoir hydrocarbon and water fluid systems that obey specified physical laws, having properties that might easily be measured in the laboratory. This type of model utilizes formation volume factors and other volumetric data, expressed directly as functions of pressure, to represent the reservoir fluid system. Experience gained with this type of simulator, commonly known as a "beta-factor model", has indicated that the approach has many outstanding defects, most of them due to the assumption that those properties measured in the laboratory completely describe the behavior of the fluid system throughout the life of the project under study. The effects of varying hydrocarbon composition and mass transfer between the phases of this hydrocarbon system must, by definition, be ignored. These defects, which may often be ignored in

Jan 1, 1970

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 S. H. Ward, R. A. Barker

THE Juniper Prospect is in Carleton County, N. B., at approximately 46" 31' north latitude, 67" 20' west longitude. During the summer of 1955 an area in west-central New Brunswick was selected for aerial electromagnetic prospecting with American Metal Co.'s helicopter-mounted equipment. Generalized, small-scale government geologic maps had suggested that this area was located along the western border of the northeast-trending granitic backbone of northwestern New Brunswick and might be underlain by the Ordovician-Silurian sediments and volcanics that are the host rocks for the recently discovered sul-fide deposits of the Bathurst-Newcastle district. The helicopter electromagnetic equipment located an anomalous conductive area, coincident with a small aeromagnetic closure, about five miles west of the town of Juniper. Subsequently 30 mining claims (total area about 6600 x 8000 ft) were staked and company personnel conducted geological, geophysical, and geochemical investigations on the property. The claims and immediately contiguous areas were mapped geologically, and an area of about 2000 x 4000 ft was soil-sampled and also surveyed by ground electromagnetic and gravity methods. Finally six diamond drillholes, totaling 3500 ft in length, were drilled to test the anomalous conditions discovered by the several exploration techniques used. The discovery of base metal sulfides during drilling was a technical success for the exploration methods, but the quantity of sulfides was too small to be profitably exploited.

Jan 1, 1959

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 R. P. Harvey, J. E. Walstrom, A. A. Brown

In view of the large number of wells directionally drilled from offshore platforms and urban drillsites, greater interest is being focused on the directional survey and on the uncertainty inherent in calculating the position of the bottom of the hole, or other benchmarks in the wellbore. Of the two methods presented here for evaluating this uncertainty, the first is analytical and uses the principles of probability and statistics. The second method uses Monte Carlo simulation to generate a pattern of the scatter of possible bottom-hole locations. Both methods are based on a model for which the usual assumption is made that the curvilinear axis of the wellbore is adequately represented by a large number, M, of linear segments joined end to end from the top to the bottom of the hole and closely following the course of the hole. The errors due to this approximation are assumed to be negligible and are not studied here. One form of the model is to consider that for a given hole there exists a perfect set of data describing it, of the above form. Each linear segment has an angle of inclination with the vertical 7, in radians, an azimuthal angle from true north xi in radians, and a length Si, in feet. The subscript variable, i, extends from 1 to M. Each of the M linear segments representing the axis of the wellbore has three components in the X, Y, and Z directions. These components are calculated from the equations: The bottom-hole coordinates as determined from the reference set of data for the wellbore are given by the following for i = 1 to M: X=SX, Y = Sy, Z=SZ; . (2) A directional survey instrument, lowered into the wellbore to the M segments, or stations, measures the values of the 7; and the 3,. In general, because of errors from various sources, the resulting angular measurements will not be precisely equal to the 7" and xi. Random variables l, and Ai are therefore considered with associated probability density functions fi(li) and gi(Ai), respectively, which represent the distributions of potential survey measurements for this hole. Since it is assumed that the values Si are accurately measured they are treated deterministically. The random variables corresponding to the surveyed bottom-hole coordinates are given, for i = 1 to M, by:

Jan 1, 1970

By H. R. Cooper

Beginning from the point of view that process control of ore flotation can be approached by on-line analysis of processing behavior, an equation was developed which has characteristics useful for a multi-variable feedback loop. The application of this equation to optimizing control with a data acquisition system will be discussed in Part 11. OBJECTIVE OF CONTROL The purpose of developing a feedback method for concentrator control is to provide a means for determining ore processing behavior in order that data describing the ore can be used to make appropriate changes in control conditions while the ore continues to flow through the flotation circuit. In absence of simple linear relationships between pairs of variables in the flotation process, attention is directed to the need for a generalized relationship which possibly would incorporate all pertinent variables. It would be possible to apply such an equation to the control of variables in a concentrator in real time according to operational requirements as ore gradually changes in quality and processing behavior. To meet this objective, it is necessary to define relationships describing flotation performance in terms of performance variables capable of being measured on-line, either with on-line instruments and analyzers or with manual methods lending themselves to quick determination of meaningful information. Relationships derived to meet this objective form a model for the flotation process. The evaluation of flotation parameters of the ore and the execution of control on the basis of the derived information while the ore is continuing through the process would constitute a feedback control system. Data for the following control and performance variables can be made available: 1) ore feed rate, 2) collector addition rate, 3) lime addition rate, 4) pulp density, 5) assay or multiple assays in feed, tails, or concentrate sufficient for material balancing, 6) particle size distribution in flotation feed, 7) froth height, and 8) frother addition rate. Other process variables, such as flotation air input and mechanical agitation intensity, are usually not considered as control variables and were not taken into account in this analysis. Other data, such as that describing the ore, would be useful for concentrator process control. However, information of this type (such as ore hardness, grain size of the mineral, chemical species of the mineral, and ore surface properties) being not presently measurable in the on-line sense, has the effect of causing disturbances in the system which can be detected only indirectly. The process control system being sought must recognize disturbances by their effect on measured performance followed by determining control changes to compensate for them in meeting the objectives of improved performance (perhaps higher throughput, maximum return, etc.). A suitable model for the process having capability of recognizing the effects of ore disturbances should provide the means for introducing proper compensation through control changes. This would be done by processing on-line data with the model to calculate coefficients reflecting ore characteristics. This procedure can only be applied in a typical flotation system when on-line data can be processed at faster rates than the usual rate at which the ore in flotation feed takes on significantly different characteristics. For example, material balance data can be accumulated in approximately 5-min intervals using on-line X-ray fluorescence analysis. The flotation system time constant might be 15 min including grinding. Ore changes having time constants longer than 1 hr could be accommodated using on-line feedback control by this procedure. FORM OF A FLOTATION MODEL A suitable form for mathematical expression of the flotation process is an equation involving rates, when each of the components entering into the separation can be treated independently by a rate equation. A first-order, kinetic rate equation can be written on

Jan 1, 1967

By R. G. Turner, M. G. Hubbard, A. E. Dukler

Gas phase hydrocarbons produced from underground reservoirs will, in many instances, have liquid phase material associated with them, the presence of which can affect the flowing characteristics of the well. Liquids can come from condensation of hydrocarbon gas (condensate) or from interstitial water in the reservoir matrix. In either case, the higher density liquid phase, being essentially discontinuous, must be transported to the surface by the gas. In the event the gas phase does not provide sufficient transport energy to lift the liquids out of the well, the liquids will accumulate in the wellbore. The accumulation of the liquid will impose an additional back pressure on the formation that can significantly affect the production capacity of the well. In low pressure wells the liquid may completely kill the well; and in the higher pressure wells there can occur a variable degree of slugging or churning of the liquids, which can affect calculations used in routine well tests. Specifically, the calculated bottom-hole pressures used in multirate backpressure tests will be erroneous if the well is not removing liquids on a continuous basis, and gas:liquid ratios observed during such a test may not be correct. Several authors1,3,8,14 have suggested methods to determine if the flow rate of a well is sufficient to remove liquid phase material. Vitter14 and Duggan' proposed that wellhead velocities observed in the field would be adequate for keeping wells unloaded. Jones8 and Dukler3 presented analytical treatments resulting in equations for calculating, from physical properties, the minimum necessary flow rate. An analysis of these studies indicates the existence of two proposed physical models for the removal of gas well liquids: (1) liquid film movement along the walls of the pipe and (2) liquid droplets entrained in the high velocity gas core. Although there probably is a continuous exchange of liquid between the gas core and the film, they will be treated separately for the purposes of this study. The development and comparison of these separate models with experimental data will permit the determination of which, if either, is the controlling mechanism for the removal of liquids from gas wells. The Continuous Film Model Liquid phase accumulation on the walls of a conduit during two-phase gas/liquid flow is inevitable due to the impingement of entrained liquid drops and the condensation of vapors. The movement of the liquid on the wall is therefore of interest in the analysis of liquid removal from gas wells. If the annular liquid film must be moved upward along the walls in order to keep a gas well from loading up, then the minimum gas flow rate necessary to accomplish this is of primary interest. The analysis technique used follows Dukler2 and Hewitt5 and involves describing the profile of the velocity of a liquid film moving upward on the inside of a tube. The minimum rate of gas flow required to move the film upward is then calculated.

Jan 1, 1970

By H. J. Modi, D. W. Fuerstenau

pH effectively regulates the flotation of corundum through its control of the surface charge. Since collector ions function as counter ions in the double layer, the collector must be anionic when corundum is positively charged and cationic when it is negatively charged. Even a long-chained xanthate will float positively charged corundum. Association of the hydrocarbon chains of the adsorbed collector ions gives rise to increased adsorption and increased rate of flotation. Because of this phenomenon, neutral organic molecules aid flotation although the unionized collector cannot promote flotation independently. Na2SO4 was found to depress the flotation of corundum with animic collectors through ionic competition whereas the same salt activates corundum for flotation with cationic collectors. A close connection between ion association, electrokinetic behavior, and floatability is suggested. As part of the increasing interest in the physical chemistry of flotation, strong experimental evidence has been obtained which indicates that organic ions are adsorbed as counter ions by non-metallic oxide minerals mainly because of the electrical forces which exist at interfaces(1-3). Thus flotation of such minerals without activation should then be appreciable only when the mineral and the collector are oppositely charged. Since only limited work has been reported on the effect of surface charge upon flotation behavior(3) this general hypothesis on the flotation of insoluble oxide minerals has not been tested thoroughly. Investigating the flotation of a mineral for which extensive electrokinetic data exist yields a way to test this hypothesis more extensively. The electrokinetic behavior of corundum in solutions containing both organic and inorganic electrolytes has been studied quite thoroughly' Electro-kinetic measurements have shown that in the case of corundum, simply by controlling pH the solid surface can be charged either positively or negatively. For this reason corundum was chosen for study. The object of this paper is to present a study of the flotation of corundum and to explain the results in the light of our knowledge of adsorption mechanisms gained from streaming potential investigations. This paper is concerned mainly with the effect of surface charge on flotation with anionic and cationic collectors, but in addition such flotation phenomena as collector hydrolysis, Van der Waals attraction between hydrocarbon chains, ionic competition, and activation are discussed. While the present paper is confined to a reporting of the data obtained and some of the more immediate conclusions, it will be seen that the results, in combination with those on other systems, do promise a

Jan 1, 1961

By S. H. Ward

Since the advent of the first airborne electromagnetic system, it has been evident that such systems were inherently limited to shallow depths of exploration of the orderof 100 to 200 feet. Hence in 1948 the author and his associates commenced a study of the feasibility of applying natural audio frequency magnetic fields to investigation of the electrical properties of the subsurface in order to provide a system capable of overcoming this limitation. This research culminated, in the summer of 1958, in the development of a new airborne prospecting device — AFMAG. This device has detected surface massive sulphide bodies while flown at heights as great as 3000 feet. Hence, by inference, it is deduced that it could find similar sulphide bodies if they were buried 2500 feet below surface and the aircraft flown at 500 feet. Thus a considerable improvement over conventional airborne electro-magnetic systems is obtainable with AFMAG. Airborne exploration in mountainous terrain involving deep valley fill now becomes possible. Examples are presented of the results of actual airborne AFMAG test surveys over known deposits of massive sulphide mineralization. Interpretation of the data is discussed with particular reference to the knowledge of size, dip, depth, depth extent, and conductivity of causative bodies as deduced from the flight records. The possibility of mapping major geologic structures, for both mining and oil exploration purposes, is discussed with the aid of examples of airborne AFMAG data. INTRODUCTION Electrical prospecting with natural audio frequency electromagnetic waves as a source has, within the last three years, become a practical reality. Both airborne and ground survey instruments under the trademark "AFMAG" (for audio frequency magnetics) have been used successfully in the search for conductive mineral deposits. It is the object of this paper to discuss the principles of airborne AFMAG and to illustrate its application by means of examples obtained over known geological situations. WHY AFMAG? In view of the large number of good airborne electromagnetic systems available for purchase or contract surveying, one may well ask, "Why another new system?". Primarily the answer lies in the fact that AFMAG does not suffer from the very limited depth of exploration available with, and inherent in, any conventional airborne electromagnetic system. AFMAG is not conventional insofar as an artificial transmitter is not required. Only a receiving and detecting device is employed. The electromagnetic fields measured are those from millions of electrical discharges in the atmosphere (e.g;..- lightning strokes). Their random distribution in time and space make these fields difficult to employ and it is only with recent advances in the electronics of geophysical instrumentation that this has become possible. A second reason for AFMAG and one of its

Jan 1, 1961

Jan 1, 1963

By V. N. Cox, C. H. Grant

Both the underwater method and the rock cratering method contribute useful information in evaluating and comparing new explosive compositions. Results indicate that metallized explosive systems which contain no high explosives can be formulated to offer more energy and to do potentially more work than other slurry explosive systems. Commercial field trials are presently being conducted with a metallized explosive that does not contain any high explosive. Although there have been only a few commercial blasts to date, the results are encouraging. The Dow Chemical Co. has been engaged in explosives research since 1953. The incentive for this research was based on the desire to find additional markets and uses for magnesium metal of which DOW is a heavy producer. The thermodynamic calculations for explosives incorporating light metals indicate energies approximately twice that of nonmetallized explosives currently being used. Early in 1960, additional patent rights were acquired from independent inventors and a development station at Virginia, Minn., was set up to conduct a commercial testing and market evaluation program with this new product. The actual manufacturing of the explosive was to take place in the boreholes. Due to the complexities involved in controlling the manufacture of the explosive in the borehole, performance was sometimes in doubt. Furthermore, experience indicated that the metallized explosive systems could not be evaluated by the standard laboratory methods then being used, and it was necessary for us to look for more suitable research methods. Tests were conducted on several known methods for explosive evaluation. The most promising were the Navy underwater method for energy evaluation and a modified Livingston crater method intended to compare explosive performance in rock. After shooting over 2000 charges of 8 to 60 lb under water and in taconites and granites with many types of explosives, a metallized explosive was developed which is now being manufactured at a pilot plant in Virginia, Minn. and is being evaluated commercially. UNDERWATER METHOD An underwater method of evaluating explosives has a number of advantages. A permanent site can be established and operated indefinitely. A large number of charges may be fired in a short length of time. The medium in which the explosive is shot is exactly the same for successive shots. Data are recorded electronically minimizing human judgment. A disadvantage of the underwater method is the lack of confinement provided by the water medium. Many explosive compositions require a greater degree of confinement to develop good energy release. Although the relative significance of the two components of energy release that are usually measured, the shock energy and the bubble energy, and their relationship to breaking rock are not clearly understood, the underwater method is a valuable research tool for screening explosive compositions and gaining an insight into the manner in which explosives release their energy. Fig. 1 shows schematically the underwater apparatus. The piezoelectric gages are calibrated tourmaline crystals which convert the shock wave to an electrical impulse which is transmitted to an oscilloscope and recorded on Polaroid film. The shock traces are picked up by scopes 1 and 2. Scope 1 is timed to pick up the full trace whereas scope 2 is timed to cover the first increment of the shock trace expanded, and therefore allows full inspection of the propagation and peak pressures. Scope 3 records the bubble trace. When the explosion occurs, a bubble of gaseous detonation products is formed. The outward momentum of the water keeps expanding the bubble even after

Jan 1, 1963

By C. E. Melbye

A serious problem confronted the city of Gunnison early in 1958 in that, for a few months during each spring runoff, the water supply derived from the Gunnison River became polluted to an unsafe point. The city management had the choice of building a new filter plant to replace the old unworkable one or drilling wells to develop a ground water supply. The cost of a new filter plant was estimated at $180,000 compared to $50,000 for 6 wells, pumps, and additional water lines, hence their decision to attempt the development of sufficient ground water. Six wells at 300 gpm per well were estimated as necessary to handle the peak usage of 1800 gpm, compared to only 600 gpm winter usage. Although the water table is only about 7 feet below the surface, existing private wells drilled at random locations were producing much less than 300 gpm. The problem was to find sufficient quantities from about 6 wells. Since the presence of an ancient bedrock channel was suspected, which would be expected to localize greater quantities of water, the city of Gunnison authorized a resistivity survey to investigate bedrock topography below unconsolidated gravels. Survey results were highly satisfactory and led to the preparation of a fairly detailed bedrock contour map, Fig. 1. This map was used in selection of all well sites. The city of Gunnison is situated in a flat, broad, Quaternary gravel-filled valley (terrain dips south at 20 feet per mile) at the junction of the Gunnison River and Tomichi Creek. The valley, about 2-4 miles wide, is bordered by hills of Tertiary volcanics to the northwest, Cretaceous and Jurassic sediments to the north and northeast, and pre-Cambrian granite and gneiss to the south. One previous deep test well at the west end of the city hit pre-Cambrian granite at 496 feet, which, allowing for 40 feet of surface gravels, makes the older sedimentary section 456 feet thick. This section is probably all Jurassic Morrison formation, consisting of the Saltwash sandstone overlain by the much thicker section of Brushy Basin shales and mudstones. According to the driller's logs, each well was bottomed in what is believed to be Brushy Basin shale. FIELD PROCEDURE A standard Gish-Rooney type of resistivity apparatus was employed, which used 45-volt dry batteries as a direct-current power source. Applied voltages could be selected from 45 volts to a 270-volt maximum. Copper-coated steel stakes were employed as current electrodes and porous porcelain pots as potential electrodes to avoid possible polarization effects with direct current. The expanding Wenner electrode configuration was utilized for this investigation and measurements were made at 20 40, 60, 80, and 100-foot spacings,

Jan 1, 1961

By Thomas C. Wilder

Jan 1, 1965

By J. D. Rice, S. C. Pitzer, C. E. Thomas

Interpretation of pressure build-up data obtained in the conventional manner has often been difficult because of the deviation from theoretical behavior. Major causes of this deviation have been attributed to damage and afterflow, and to fluid redistribution in the wellbore which, in extreme cases, can result in a pressure hump in the early portion of the build-up curve. Theoretical investigations show that bottom-hole pressure is definitely influenced by phase redistribution in the tubing column during surface shut-in tests, and that the magnitude of the effect is sensitive to the producing gas-oil ratio and stabilized rate of flow in the well. For field experiments a wire-line tubing packer, which can be run in the tubing against a stabilized flow rate, was developed for bottom-hole shut-in tests. By use of this bottom-hole .shut-in method, pressure humps previously observed in surface shut-in tests were completely eliminated and the effects of afterflow minimized. Analog and digital computer studies have been made to obtain theoretical curves for comparison with field curves, and remarkable agreement between the results of bottom-hole shut-in tests and theoretical curves has been obtained. INTRODUCTION Pressure build-up data from shut-in wells have been used by the petroleum industry to determine the permeability of the formation, to estimate wellbore damage, to evaluate the static reservoir pressure and to speculate reservoir volume. Calculation of these factors is based on methods of analysis developed for theoretical systems whose build-up curves have a charac- teristic shape when the wells are shut-in at the sand faCe.1,2,3,5,8,10 However, field curves obtained by conventional surface shut-in methods do not always exhibit this characteristic. Such factors as stratification, rock heterogeneities and irregular reservoir geometry can cause the character of a build-up curve to deviate from that predicted by theory for a simple system. In addition, all field build-up data are affected by the methods utilized in obtaining the measurements. Conventional surface shut-in techniques allow two effects to occur which contribute directly to the manner in which pressure builds up at the sand face. The first, afterflow, has been recognized for some time and methods are available for estimating the magnitude of its effect and for correcting its influence.' The second, phase segregation in the tubing column after shut-in, has been reported only recently and appears to exert a considerable influence on the character of field build-up curves.!' In extreme cases, phase segregation can produce "pressure humps" in the early portion of the build-up curve. Such curves have been considered anomalous and have defied analysis. It is the purpose of this paper to discuss the effect of phase segregation in the tubing string during build-up and to present a bottom-hole shut-in technique for obtaining field build-up data which minimizes the influence of afterflow and eliminates phase segregation effects. It will be shown that reliable reservoir information can be calculated from build-up measurements obtained using this method, even though conventional surface shut-in tests on the same wells yield anomalous data. DISCUSSION OF PHASE REDISTRIBUTION EFFECTS Some of the effects of phase redistribution on pressure build-up curves have been described previously by Matthews and Stegemeier" who have presented evidence that phase segregation is responsible for the pres-

By R. F. Shurtz

The problem of estimating the precision of systematic samples from a mineral deposit is attacked by interpolating the quality, or other attribute measured, by using Fourier approximation. Such approximation using trigonometric functions yields a particularly simple expression for the variance of the mean value of the attribute sampled. The increasing volume of publication in which the science of statistics has been applied to the problem of predicting the quality of an ore deposit from sets of samples taken from it might lead one to believe that engineering judgment can be replaced by statistics. There is, in fact, an implicit application of the principle of indifference or of "no advance knowledge" which has been discredited in modern statistical theory and practice. These pages will show that it is not always advantageous, efficient, or safe to apply this principle. This problem will be dealt with by using series of orthonormal functions to describe the distribution of the quality of ore in space as contrasted with the more commonly used distribution of quality as a frequency function. Mining engineers are rarely well acquainted with the use of such series or of the functions comprising them. These uses are commonplace in the fields of communication and information theory and it is not surprising that applications taken from these fields should be found appropriate for the analysis of the information contained in a set of samples from an ore deposit. It is hoped that presentation of the following basic, but incomplete, results will generate wider interest in this method and a fuller development of it as well as of the necessary computational techniques. FUNDAMENTAL PREMISES Random variables in the strictest sense are not very common in mine sampling. The sets of samples are almost always taken systematically and are, hence, not subject to all the laws governing random events. Therefore, unless there is indeed random experimental error in the determinations of quality, the treatment of systematic samples of n members as though they were random samples of n members often leads to calculation of confidence limits for the mean quality which are larger than they should be with the result that unnecessary additional sampling is done. The reason for this will emerge in the following discussion leading to Eq. 24 and its generalization. In order to inquire into such questions on a rigorous basis and at the same time avoid arithmetic complications that would obscure the concepts involved, a one dimensional distribution of quality will be assumed. The quality, Q, in suitable units is expressed as a function of the distance, x, along some straight traverse of the deposit. Such a function

Jan 1, 1967

By G. T. Engel, W. G. Gruzensky

Jan 1, 1960

By J. W. Givens, P. B. Crawford

A potentiometric model study has been made of the effect of vertical fractures existing in the matrix of the reservoir on the flooding or cycling performance. Fractures can have unusual flow characteristics. Fluid entering one side of a fracture can emerge on either the same or opposite side of the fracture, depending on the particular streamline. Fracture orientation has a great influence on the sweep efficiency. However, sustained fluid injection may still permit large areas to be swept. INTRODUCTION One can find considerable information available in the literature on the effect of fractures originating at the wellbore on sweep efficiencies of water-flooding or gas cycling programs.1-5 However, there arc few, if any, quantitative data reported in the literature on the effect of vertical or horizontal fractures existing out in the reservoir matrix on secondary recovery performance. Several papers have been presented on the effect of vertical and horizontal fractures on the productivity or conductivity of waterflood and gas injection patterns, but in all cases the fractures initiated at the well and simulated commercial fractures.6,7 It is believed that natural fractures may exist out in the reservoir matrix and would effect the displacement performance. The purpose of this report was to study the effect of a few isolated vertical fractures existing out in the reservoir matrix on the performance of waterflood or gas injection patterns. DESCRIPTION AND EQUIPMENT Since the potentiometric model was described by Lee8 it has been widely used to study numerous types of fluid displacement problems. The potentio- metric model can be used for fluid displacement studies when the following assumptions are valid: (a) steady-state conditions exist; (b) the mobilicy ratio is one; and (c) the capillary and gravitational effects can be neglected. Five-spot and direct line drive square patterns were studied using a copper strip of the desired length and orientation to simulate a vertical fracture in the reservoir matrix. The 20 x 20-in. model was considered to represent only one element of an infinite array of similar patterns. DISCUSSION AND RESULTS FIVE-SPOT PATTERN — ONE FRACTURE BETWEEN INJECTION AND PRODUCING WELLS Fig. 1 shows a quadrant of a five-spot pattem with the vertical fracture existing out in the matrix along a line connecting the injection and producing wells. Length of the fracture was equal to 35.4 per cent of the distance between injection and producing wells. In this particular pattern it is noted that fluid breaks through from the injection well to the producing well when the dimensionless time is equivalent to 26.9. Dimensionless time is defined as volume of fluid injected divided by the volume of displaceable fluid in the pattem expressed as a per cent. The cross-hatched areas of this

Jan 1, 1967

By Don K. Smith

The Reclamation Program for the Northern Zone of the Indus Plains in West Pakistan involves the construction of 30,000 irrigation wells to serve 20 million acres with an annual pumpage of nearly 40 million acre-ft. Detailed hydrologic studies indicated that mild steel would be satisfactory for casing; however, within two years after construction, about 10% of the wells began to fail due to encrustation and corrosion of the casing. A search to find a noncor-rosive substitute for the steel casing disclosed that only fiberglass-reinforced plastic pipe satisfied all requirements. The results of initial use of fiberglass have been highly satisfactory and have resulted in a lower over-all well cost. Advances in technology pomise further reduction in prices and improvement in performance. The northern zone of the Indus Basin comprises a vast area of essentially level land in West Pakistan which is irrigated by means of a complex system of canals that distribute more water to more land than any other irrigation system in the world. In spite of the development of this extensive irrigation system, production of food and fiber has not met the requirements of the country. The major reasons that crop production has not reached anticipated levels are the lack of adequate irrigation supplies and effective subsurface drainage. As a result, yields per acre are among the lowest in the world for areas where irrigation is practiced, and many thousands of acres of valuable crop land have been removed from production or seriously affected by salinity and waterlogging. To rectify the problems of irrigation supply and drainage the Government of Pakistan has undertaken a program of Salinity Control and Reclamation which is being administered by the Water and Power Development Authority. The Authority retained Tipton and Kalmbach, Inc., of Denver, Colo., as engineers for the program in the northern zone of the Indus Basin. This program consists of the construction of about 30,000 irrigation wells to provide drainage for waterlogged lands and full irrigation supplies, including leaching requirements, for all cultivated lands. The program will ultimately ensure the productivity of 20 million acres and involve an annual pumpage from groundwater storage of nearly 40 million acre-ft of

Jan 1, 1970

By H. J. Gerritsen

From presently available components a portable, rugged, reliable apparatus can be built which will be able to detect methane concentrations of 0.1% and lower in air. Sensitivity and design considerations of the apparatus are given. Ways to build even simpler and cheaper apparatus in the future are indicated. Traditionally the problem of detecting and monitoring the presence of dangerous gases present in small concentrations in air has been a difficult one. Whether canaries,' or sophisticated chemical or thermal conductivity methods are used, they seem to be either slow, cumbersome, or otherwise not completely satis factory. A seemingly natural way of gas detection is to measure the strength of absorption undergone by light of a wavelength that the gas to be detected absorbs. This method has been difficult to apply in view of the low power output over a narrow wavelength range in the infrared from traditional glow bar light sources. Besides, the apparatus to select such a narrow band of wavelengths is quite costly and complex. The use of a laser solves these problems simultaneously. Although a variety of lasers are available to detect a variety of gases (such as carbon monoxide or methane) whose early detection would be desirable, this paper deals exclusively with the use of lasers for methane detection. ABSORPTION SPECTRUM OF METHANE IN AIR The absorption spectrum of pure methane is very well known.2-4 There are four fundamental vibrations, all occurring in the infrared of which only two absorb strongly. They occur at 3.3 µ and 7.7 µ. The absorption spectrum is quite complicated due to the fact that the methane molecule can rotate in a variety of rotational quantum states. These rotational states are furthermore split by anharmonic forces such as the Coriolis force.4 When methane is present in air at atmospheric pressure, much of the very fine structure due to the anharmonic force is washed out since the pressure broadening is of the same magnitude as the splitting due to anharmonicity. For most of the spectrum at 3.3 µ, the rotational lines are typically separated by 0.01µ while the individual linewidth is 0.0001 µ. It is thus very important that the laser emit at a wavelength that coincides with one of these rotational lines. Before elaborating on the point that such a laser exists and how it performs, it is necessary for later discussion to point out that in the

Jan 1, 1967