Search Documents

By W. B. Wise, H. D. Attra, W. M. Black

The purpose of this paper is to illustrate a comparatively new approach for solving a problem that has plagued oil producers for many years—how to make the most money with available field production capacity. In a very basic operation such as a top-allowable one-well field, a computer and staff of engineers are not required to determine that the best operation is to produce that well at its maximum for as long a time as it will produce. However, as the problem becomes more complex and the number of wells and number of reservoirs increase, it probably will be a wise choice to devote some engineering thought and calculation time to determine the best way to operate that field on both an ultimate and day-to-day basis. Continued developments in the field broadly referred to as "operations research" have led to a number of useful engineering analyses for optimizing field operational problems. One of these techniques, termed linear programing, is parlicularly interesting because of its potential applications to the study of complex oilfield operating problems. Linear programing is concerned with the problem of planning a complex set of interde~endenl activities in the best possible (optimal) fashion. A field application of linear programing is being used to maximize daily income from a multi-reservoir producing operation. The initial results of this study indicate that the linear programing technique (I) can be applied effectively to field problems, (2) can be used to detect unsuspected avenues for improving operations and (3) can effect a significantly improved means of controlling field operations through the use of the highspeed computer. INTRODUCTION The usual approach for studying complex field operational problems has been to look in detail at several basic plans suggested by judgement and experience and, then, to choose the best plan from those investigated. This approach is entirely satisfactory provided: (1) that judgment and experience alone will reliably pick the best or near the best program as one of the plans to be investigated; and (2), and most important, that the factors which control field production will not change. This paper deals with the large group of field operations which are at least one-step removed from a satisfactory solution using the conventional method of analysis. Problems such as pipeline- and gathering-system design, transportation studies, large field development projects and multi-reservoir producing operations fall in this latter category. To find the best program for handling these problems when the best program is not obvious, a complete optimization study is required, and high-speed computers are required to continually update the program once it is selected. Operations of this type have an extremely large number of possible solutions. However, only one will result in maximum economic return and, more important, this one program must be re-evaluated and adjusted constantly to accomodate such changing field factors as the number of producing days, the price of crude, the discovery of new reserves, declining reservoir pressure with subsequent decline in production and a host of other factors which characterize specific operations. One of the most widely used, direct methods of optimization is a numerical procedure termed "linear programing", which is designed to investigate systematically each of the possible solutions involved in a given problem, determining rapidly and uniquely the optimum answer. As implied by the name, linear programing requires that all field variables be either defined or expressed in a series of interdependent linear equalities or inequalities. These relationships are solved simultaneously, using a form of matrix solutions to determine the optimum values for each producing variable that will maximize net profits. The solution techniques used in a recent field case study is termed the "simplex method" and is discussed in full detail in Ref. 1.

By James Humphrey

There is a natural resistance to change in any form. When it comes to implementing a game-changing technology like autonomous haulage on an operating mine site, the resistance can at times appear insurmountable. Whether based on real information or misguided preconceptions, it?s essential that we don?t underestimate the effort it will take to overcome this resistance, while clearly recognizing when it?s worth the effort. These realizations did not happen overnight for sites where autonomous haulage has been implemented. Nor did they happen overnight for manufacturers of autonomous mining equipment. Caterpillar, for example, began building toward autonomy more than 30 years ago and had a truck operating in the mid-1990s. Lessons learned during those early years led to the realization that there was much more to be learned ? not just about the machines themselves, but about the systems and knowledge that must be in place for a successful implementation. Working closely with customers, initiating real-life demonstrations on mine sites, having discussions with regulatory agencies and other stakeholders ? these and many other activities were necessary before autonomy could be fully and successfully launched. After decades of effort, Caterpillar?s first official commercial autonomous mine site went to work in 2011 in Farmington, NM.

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.' ' Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.'" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-

By J. L. Fitch, B. G. Hurd

The presence of gypsum in samples subjected to standard core analysis introduces serious errors in the measurement of water saturation and porosity. The magnitude of these errors, depending upon the type of analysis procedure used, has been evaluated experimentally and theoretically. Water-saturated pore volumes determined by vacuum distillation or retort procedures may be too high by as much as 48 per cent of the volume of gypsum present. The maximum error in this value determined by Dean-Stark extraction with toluene is 36 per cent of the gypsum volume. If pore space is measured by a Boyle's law type porosimeter the maximum erorr due to gypsum is 38 per cent of the original gypsum volume present. The conscientious core analyst may eliminate the errors caused by gypsum by either of two methods. He may apply suitable corrections for the water of crystallization removed from gypsum, or take special precautions to prevent dehydration of gypsum during core cleaning and alzalysis procedures. These alternatives are discussed and methods for obtaining reliable core analysis data are proposed. Gypsum (CaSO4 . 2H2O) occurs as a common mineral in many of the world's sedimentary basins. In addition to its presence as primary bedded evaporitic deposits, it is often found as secondary replacement zones and pore-fillings in host rocks. That it may cause false indications in porosity, both on neutron logs and in standard core analysis, has been generally recognized. False porosity interpretations from neutron logs result because the chemically bound water in gypsum attenuates neutrons as effectively as free water. Errors in core analysis result when gypsum is dehydrated prior to or during the analysis. Vacuum retort procedures used by many laboratories to determine fluid saturation and porosity effect a complete removal of water of crystallization from gypsum. The crystal water, distilled and collected with free water from the sample, is counted as water-sat11rated pore space. Dehydration of gypsum results in a decrease in volume occupied by calcium sulfate and a corresponding increase in pore volume. Thus, if gypsum is dehydrated by sample cleaning and drying procedures, subsequent porosity measurements with gas expansion or compression instruments will be high. The qualitative effect of gypsum on porosity determinations by retort procedures was discussed in a recent article by Bynum and Koepl3. They concluded that since false porosity was interpreted as being water saturated, determinations of oil content were not affected. Others in the industry, however, have pointed out that porosity values are often used with values of water saturation determined independently, e.g., from electric resistivity logs in making reserve estimates. Values of water saturation calculated from the formation resistivity using the Archie equation, S4= are actually low if values of are high due to gypsum. In this case, false porosity values due to gypsum could cause serious errors in calculations of oil-in-place. This investigation was made to determine quantitatively the effect of gypsum on porosity values determined by different core analysis procedures. EXPERIMENTAL PROCEDURES Dehydration of Gypsum The extent and rate of removal of water of crystallization from gypsum under conditions encountered in core analysis procedures were studied. Taken for the experiments were three types of relatively pure, naturally occurring gypsum, viz. selenite, satin spar and massive gypsum. Although all samples probably contained trace chemical impurities, each type lost the theoretical amount of water upon complete dehydration. An automatic recording thermobalance, similar to one described by Mauer2, was used to obtain dehydration rate data for samples heated in dry air. Water

By Michael S. Wieland

This report investigates whether it is tractable and thereupon reasonable to remove the booster influences when ranking the toxic fumes of non-cap-sensitive mining explosives. The hazard potential is represented by relative fume toxicity RFT, the resultant influence computed from a formulated rule with a restricted set (tally) of fume components and multiplying constants. The RFT result is compared to the rule criterion that represents the worst case tolerable toxicity stipulated by regulations or otherwise. The typical nonstoichiometric booster composition can render notable unwanted influence on the RFT ranking results. Historically, the remedy was to shoot the booster separately and reduce the trial concentrations by the relevant subtractions. Unreal negative concentration results were noted on rare occasions, revealing the faulty nature of the reduction procedure. With thermodynamic reaction chemistry codes TDRC, on the other hand, the booster ingredients can be retained in or removed from the reactant tabulation, yielding reduction factors taken as the concentration ratios for the two circumstances. The underlying presumption regarding reduction factors is that the ratios of component concentrations with and without the booster ingredients when resolved theoretically ought to equal those rendered experimentally, when tractable. Upon rearranging the presumed relationship, the unadjusted fume concentrations divided by theoretical reduction factors yield the wanted readjusted fume concentrations, thereby removing the booster influence without generating unreal negative concentrations. The work principle from thermodynamics restricted by the constraint of zero-net interaction loss is utilized in conjunction with the TDRC to resolve the wanted (Z-state) fume spectrum. This circumvents the recognized difficulty with traditional hydrodynamic detonation theory that is related to the nonideal reaction characteristics of typical mining explosives. The tacit characteristic of the Z-state resolved from the non-trajectory work principle is that it can not be remedied or tampered with, unlike the traditional trajectory techniques where the reaction products can be readjusted by changing formula representations, nullifying specific reactions or otherwise manipulating the thermodynamic trajectory process. Though the work principle theory is an unfinished formulation, the Z-state results prove worthy enough for rendering trends and comparisons and reaching otherwise useful conclusions. Further refinements should yield more reliable fume spectra and indicate explosive composition readjustments that would reduce toxic effects within underground mines and the wider world environment.

By H. H. Rachford, D. W. Peaceman

A system of partial differential equations describing miscible displacement of fluids in porous media is derived The system takes into account the influence of gravity, the spatial distribution of permeability, diffusion, and fluid viscosities and densities. A numerical procedure for approximating solutions to the differential systems has been tested for a horizontal two-dimensional geometry. In end-to-end displacements of oil with less-viscous solvents, the numerical solutions exhibited fingering qualitatively similar to that observed in laboratory models. Small random spatial variations in permeability about the mean value are sufficient to initiate fingering. Quantitative comparisons of computed results with laboratory data show good agreement. INTRODUCTION Miscible flooding of oil by solvent is receiving increasing consideration for field use.1-4 Unfortunately, such floods potentially present severe problems in loss of recovery through the by-passing of oil by expensive solvent. Consequently, their economic evaluation requires sound techniques for predicting recovery. The purposes of this work are to present a finite-difference method for calculating the multidimensional displacement of oil by solvent and to investigate the validity of the method by comparing results of calculations with data from displacements in laboratory models.5 The formulation of the method to simulate the model experiments treats a case of limited scope in the description of solvent flooding in the reservoir. The model experiments were carried out with fluids that were assumed to form an incompressible, ideal, two-component system with constant diffusivity. Establishing the validity of the method even in its present form provides a major step toward the goal of quantitatively evaluating individual solvent flooding projects. First, it demonstrates the fea- sibility of calculating the course of displacements which are dominated by an inherent macroscopic instability, i.e., viscous fingering. In addition, the method provides almost the only practical means of examining the effect of the size and extent of reservoir inhomogeneities on the development and propagation of the fingers, taking into account the important influences of diffusion and gravitational segregation. It might be inferred that such studies cannot yield results of practical value because in displacements dominated by macroscopic instability the pattern of finger development and the resulting performance should depend critically on small variations in the geometric distribution of reservoir inhomogeneities. The present method offers the capability inherent in computational techniques of predicting performance reproducibly with any arbitrary distribution of reservoir properties and, as such, provides a means of evaluating the sensitivity of behavior to uncertainties in reservoir definition and, thus, assessing the reliability of prediction of performance. THE PHYSICAL PROBLEM DATA USED FOR QUANTITATIVE TESTS OF CALCULATIONS In this work, a calculation for treating the miscible-displacement process is tested by comparing calculated results with experimental data. The experiments chosen for comparison have been described in detail,5 and the quantitative features are summarized in a later section. Briefly, oil was flooded by solvent of equal density from a thin rectangular channel in Lucite packed with uniform Ottawa sand. Under these conditions, two-dimensional geometry is considered adequate to represent the process. Inasmuch as the present work is concerned with testing of the method by comparison with specific two-dimensional data, the system of equations presented in the succeeding section will be oriented to defining the physical system specifically. THE DIFFERENTIAL EQUATIONS Suppose x, y to be a Cartesian co-ordinate system and define h(x, y) to be the height of a point above a horizontal reference plane. It will be assumed

By H. N. Hall

Various factors must be considered in an engineering evaluation of gravity-drainage reservoirs. Among these are: (1) the effect of producing rate on total oil recovery; (2) the effect upon well productivity and ultimate recovery of the pressure level maintained during the producing life of the reservoir; (3) the economic advantage of full or partial pressure maintenance; and (4) estimate of the rate of gas production and injection and the possible purchase of gas under conditions of full pressure maintenance to ascertain compressor facilities needed. All of these factors can be evaluated only when a reliable method is employed for determining reservoir performance in gravity-drainage reservoirs. The purpose of this paper is to present a general method for calculating the performance of a gravity-drainage reservoir. This method is applicable for conditions of complete pressure maintenance, partial pressure maintenance and normal pressure depletion. Provisions are made to take into account variations throughout the reservoir of reservoir configuration, changes in permeability and fluid composition. Based on the method presented in this paper, an IBM 650 computer program has been developed. The past performance of an actual gravity-drainage reservoir producing under conditions of declining pressure and no gas injection was duplicated using this program. INTRODUCTION In tilted reservoirs the production of oil is influenced by drainage of oil from upstructure to downstructure locations. When this downstructure drainage of oil is sufficient to cause effective segregation of the gas and oil in a reservoir, the reservoir is usually classified as a segregation drive or gravity-drainage reservoir. (Discussion will be restricted to gravity-drainage reservoirs which have no encroachment of edge water.) The important feature in gravity-drainage reservoirs is the density difference between reservoir oil and gas. These phases tend to segregate in the reservoir with the result that in the gas cap the oil saturation is maintained at a higher level by drainage of oil from the gas-cap area. Oil can be produced from the oil zone at a low gas-oil ratio and reservoir energy is thereby conserved. The standard material balance in not adequate for predicting gravity-dramage reservoir performance because it does not take into account the difference in saturation above and below the gas-oil contact. Several authors'.' have presented methods for calculating the performance of gravity-drainage reservoirs in which reservoir pressure is maintained constant by gas injection into the gas cap. Using some simplifying assumptions, these methods can be employed with a desk calculator to give acceptable results. The problem of predicting the performance of gravity-drainage reservoirs under the conditions of declining reservoir pressure is many time more complex than that of constant pressure. fierefore, attempis to develop a method suitable for desk calculation have required excessively simplified assumptions. In the past several years, highspeed digital computers have become more widely available for reservoir engineering problems. These corn puters are well suited to problems such as the prediction of the performance of gravity-drainage reservoirs with pressure decline. Many of the simplifying assumptions necessary for hand computation can be eliminated so that a realistic approach to the gravity-drainage process can be made. CONCEPTUAL PICTURE OF OIL MOVEMENT IN GRAVITY-DRAINAGE RESERVOIRS Before attempting to develop an analytical treatment for conditions occuring in a gravity-drainage reservoir, a concept should be formed concerning the movement of fluids in the reservoir as oil is produced. A review of the literature'.' shows that it is customary to classify gravity-drainage operations into two categories—(1) with complete pressure maintenance, and (2) with declining pressure. The same line of reasoning will be followed in presenting the concept of the movement of fluids in the reservoir because it is easier to visualize the movement of fluids under conditions of complete pressure maintenance. After discussing complete pressure maintenance, an analogy will be made between that and the case of declining pressure. It should be kept in mind throughout that the final aim for the problem of solving gravity-drainage performance with digital computers will be to develop a general program for any kind of gravity-drainage reservoir. COMPLETE PRESSURE MAINTENANCE One feature which is generally common in gravity-drainage reservoirs is a gas cap located at the top of the structure. This is shown in Fig. 1(a). Fig. l(b) shows oil saturations that might occur through the reservoir. In the gas cap, oil saturation is lower than

By H. D. Outmans

The mechanics of filtration are described by a theoretical-empirical nonlinear diffusion equation which, under certain circumstances, may,be linearized and then solved explicitly. For filtration under static conditions linearization leads to a boundary value problem analogous to a heat flow problem with a known explicit solution. The corresponding solution for static filtration is compared with published experimental data. Under dynamic conditions filtration takes place in three successive stages. During each of these the rate of filtration and/or the thickness of the filter cake are different functions of the filtration time. The proposed mechanism explains the observed high resistance of dynamically deposited filter cakes against erosion and also the connection between filtration rate and viscosity of the drilling fluid. Several of the quantities governing dynamic filtration have no counterpart in the static filtration mechanism. The static filtration rate is, therefore, not a reliable measure for the dynamic rate and vice versa. INTRODUCTION Filtration under simulated borehole conditions, i.e.. either from a static suspension or from a suspension flowing parallel to the filter cake (dynamic filtration), has been the subject of several laboratory studies.1-7 This experimental and theoretical work showed that many aspects of the filtration mechanism cannot be explained by an elementary filtration theory based on the assumption that the filter cake is incompressible. For a more satisfactory theory compressibility should be taken into account. This has been done. at least to some extent, in the filtration theory developed for the chemical industry. Surveying the literature (see Ref. 8 for a bibliography) it becomes apparent that, although several properties of filter cakes deposited from many different suspensions have been measured, including compressibility, the filtration equations are essentially empirical in nature. No attempt has been made, for instance, to develop filtration theory along the lines of consolidation theory.12 This theory, upon which a highly successful development in soil mechanics rests, would appear to be an excellent starting point for a filtration theory since the compressibility concept is an essential part of it. As we will use this theory in the following development it is well to state the four assumptions on which it is based: 1. The fluid flow through a compressible porous medium is governed by Darcy's law. 2. The rate of change in fluid content of an element of the porous medium is proportional to the rate of change of solid pressure between the particles. 3. The solid particles are incompressible within the range of pressures considered. 4. The total pressure on a surface normal to the line of flow is equal to the sum of the fluid pressure and the pressure between the particles (solid pressure) at that surface. To calculate the rate of filtration and other quantities of interest it is necessary to know the filter cake compressibility and the permeability as a function of the solid pressure. It has not been possible to calculate these functions from theoretical considerations and they will therefore be introduced in this paper as empirical expressions. Both are determined in a compression cell where the cake is subjected to a mechanical load. Permeability and compressibility are measured after the solid pressure has stabilized, i.e., after the excess fluid pressure has been dissipated and the uniform pressure is transmitted entirely by the solids. The permeability and compressibility thus determined are not the same as during filtration in the borehole, under either static or dynamic conditions, because then, due to frictional drag, the solid pressure and hence also the permeability and compressibility vary along a line normal to the direction of the mean flow and can only be defined locally. DERIVATION OF THE FILTRATION EQUATION FOR COMPRESSIBLE FILTER CAKES As the filter cake in the borehole is thin compared to the radius of the hole the filtration may be considered as linear. Taking the x-axis normal to the wall, with its origin at the formation, we have, according to Darcy's law,

By F. M. Stewart

A paper by Hurst and van Ever-dingen in 1949 led to the practical solution of many nonsteady-state flow problems.' Subsequently, applications of this material have been discussed by several authors. The purpose of this note is to outline a simplified direct procedure for water-influx pressure predictions for pressures above the bubble point. The prcposed procedure is unique in that after the "time-conversion factor" and the "influx constant" are obtained, predictions are made without trial and error. The elimination of trial and error is made possible by expressing certain PVT characteristics as explicit functions of pressure. The procedure is applicable to fields producing at pressures above the bubble point and under radial water drive from an infinite aquifer. In the usual applications it is required to predict (1) withdrawals given future pressure history or (2) reservoir pressure, given future withdrawals or production. The theory, assumptions and conditions under which the nonsteady-state solutions hold are covered in the references and are not repeated. The simplified direct procedure outlined here has been developed from published solutions for the so-called "pressure case". A similar direct procedure for the same special conditions has been developed for the "rate case"; however, in this instance the procedure is lengthy and offers no advantage over the trial-and-error solution. Consequently, this paper covers only the procedure based on the pressure case. DISCUSSION Under the pressure case prediction of withdrawals from a water-drive field, given the future pressure history, involves the use of a nonsteady-state water-influx term together with the material balance equation. The solution of a particular problem requires the following steps: (1) division of past pressure and production data into equal time intervals; (2) computation of cumulative water influx vs pressure drop for each time interval of past history using the material balance equation; (3) determination of the best values for the "time-conversion factor" and the "influx constant" by equating the nonsteady-state influx term with the influx computed in Step 2 (this requires a trial-and-error procedure); (4) prediction of influx for a given pressure drop using the influx term and the computed factors from Step 3; and (5) calculation of fluids produced using the predicted influx in the material balance equation. A forecast of the water-oil ratio or water production rate is required to complete the calculation. It is assumed that a reliable estimate of original oil in place is available, as well as the production and pressure history, and data on reservoir fluid characteristics. If the original oil in place is unknown, it may be computed simultaneously with the tinie-conversion influx factors2,3. However, the oil in place so computed may be considerably in error3; consequently, it is usually preferable to use oil in place estimated from pore-volume data unless there is considerable doubt as to its accuracy. The Material Balance Equation Above the Bubble Point Above the bubble point the simpli-

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 H. R. Froning, R. O. Leach

In wertability alteration flooding, a chemical agent is rnoved through a reservoir by the flood water to increase oil recovery by decreasing the degree of wetting of the rock by the oil. Substantial amounts of the chemical may be lost during movement through the reservoir. The extent of the loss, and therefore the economics of the process, depends in some cases on factors which are difficult to reproduce in the laboratory. Therefore, a short-duration, low-cost field test method is needed to permit evaluation of chemical requirements under actual field conditions. This paper describes a small scale rest conducted at a single well for measuring chemical requirements, thereby giving a more reliable evaluation of this important factor in the applicability and economics of the process. In the rest a small wafer slug containing the chemical agent and a nonadsorbed tracer is displaced into the reservoir by a known volume of wafer. The well is then placed on producrion. Chemical loss per barrel of pore volume contacted is calculated from [he fractional recoveries of the agent tested and the nonadsorbed tracer. The method has been used to determine within the actual reservoirs the chemical requirements for both a sandstone and a dolomite reservoir. Several chemical agents are potentially available for wettability alteration flooding, although none is universally applicable. For some applications of the method, chemical costs per barrel of additional oil recovered can be substantially less than one dollar. INTRODUCTION Wettability alteration flooding provides a means of increasing oil recovery from reservoirs by decreasing the degree of wetting of the rock by the oil and increasing the displacement efficiency of the flood water. Earlier studies demonstrated a relationship between oil recovery during waterflooding and the degree of wetting of a rock surface by an oil. The application of wettability alteration flooding to the Harrisburg field of Nebraska provided a field test' of this recovery process. Subsequently, additional laboratory and field tests have developed additional procedures for evaluating wettability alteration flooding, and have indicated where the process may be applicable. Applicability of this process to specific reservoirs is determined by a progression of tests to determine sus- ceptibility of the reservoir to alteration of its wettability, to indicate the degree of recovery improvement and to estimate the amount of chemical required to process the reservoir. The economics of applying improved oil recovery processes depends not only upon the degree of improvement in oil recovery achievable by the process but also upon the process costs and the timing of the income and the investment. Emphasis in this paper is on the expenditure aspects of the process. The work reported in this paper indicates that the chemical investments required for wettability alteration flooding are substantial. For evaluating the economics of a potential flooding application it is imperative that a sound estimate of the chemical requirements be made for the reservoir. Generally, true reservoir conditions are not adequately simulated in laboratory chemical propagation tests. Because of wide well spacings, many years might be required to obtain chemical propagation data from conventional pilots or inter-well tests. Consequently, n short-duration, low-cost method is needed to determine chemical requirements in the field. The potential applicability of wettability alteration flooding is discussed, as well as the economics of wettability alteration with respect to the inherent and imposed restrictions on the timing of income and investments. DETERMINATION OF CHEMICAL REQUIREMENTS In the process of moving a chemical bank through reservoir rock, some of the chemical agent lags too far behind the flood front to be effective or is otherwise lost to the reservoir system. The extent to which these losses occur overshadows the reductions in chemical concentration due to diffusion and to mixing with the reservoir fluids. Experience indicates that almost without exception, chemicals which induce a wetting change undergo either sorption reactions or chemical reactions with mineral constituents of the pore surfaces. Other reactions may occur between the added chemical and the reservoir oil and water. Even in limiting consideration to reactions of the relatively inexpensive inorganic salts, bases and acids, the reactions may be exceedingly complex. Reservoir pore surfaces consist of more than silica in sandstone reservoirs, and more than calcite or dolomite in carbonate reservoirs. Many mineral species are present, each exhibiting specific tendencies to react with an injected chemical. The reactions which occur can consume enough of the agent to have an important effect on economics. These reactions can cause a change in pH of the chemical bank, or may remove some of the active chemical by precipitation, adsorption or reaction to form a new chemical which may or may not be effective in

By Edward A. Barrett, Launa Mallett, Charles Vaught, Kathleen M. Kowalski, Michael J. Brnich, Lynn L. Rethi, Robert Randolph, William J. Wiehagen

The goal of the Pittsburgh Research Laboratory's training research is to explore the relationship between valid instruction and improved worker performance. Studies summarized in this paper include: (a) the development and assessment of training techniques that strengthen miners' ability to act competently in emergencies; (b) investigations of classroom simulations that enhance the perceptual, judgment and decision making skills of workers confronted with mine hazards; and (c) field studies of the linkages between financial investments in occupation at skills training and measurable outcomes.' The authors maintain that better ties between training and performance will aid long-term human resource development. The best way to establish this connection is through the use of objective and reliable data.

By William Gleason

Despite a global slowdown in most mining sectors, Caterpillar, the world?s largest mining equipment maker continues to forge ahead with new, and improved offerings to the global mining industry. On Nov. 10, as part of its Loading and Hauling Mining Focus Session at its Tinaja Hills proving grounds south of Tucson, AZ, Caterpillar unveiled to customers and media its two latest mining products, the 994K wheel loader and 6020B hydraulic shovel. The Caterpillar 994K wheel loader builds on a nearly 25-year history of the 994 model but with one noticeable difference ? the 994K is bigger ? much bigger. Weighing in at 239,142 kg (527,218 lb) the 994 K is up to 20 percent larger and comes with up to 30 percent increase in payload capacity. This behemoth of a machine is currently being tested at Caterpillar?s Tucson Proving Grounds, Anglo Gold Ashanti?s Cripple Creek and Victor gold mine in Colorado and Taseko?s Gibraltor copper and molybdenum mine in British Columbia, Canada.

By M. K. Larson, J. P. Dunford, E. G. Zahl, T. M. Brady

Researchers at the Spokane Research Laboratory of the National Institute for Occupational Safety and Health, Spokane, WA, have collaborated with three Western longwall coal mines in an on going effort to develop technologies that will aid in providing safe and stable working areas. The goal of the research described here is to develop a stress monitoring system that will provide immediate information to mine managers for making daily safety decisions as areas of poor ground are mined through. Initial work has focused on answering preliminary questions regarding the reliability and use of stress change patterns. Research is concentrated on monitoring horizontal stress because horizontal stress is transmitted over long distances through stiff strata, thus allowing an extended length of entry to be monitored. This paper presents an explanation of the concept, key results from field tests at two mine sites, and a proposed process for implementing a monitoring system. System design layouts, instrument use, data collection and interpretation methods, and processes to present findings to mine staff are described. Additional validation and correlation with actual failure mechanisms are required before this approach can be recommended at ongoing operations. However, initial results indicate that this approach shows promise in mines prone to bumps and roof falls associated with large stress changes. In particular, analyses of measured changes in horizontal stresses appear useful.

By J. J. Bodmer

Although the similarity between puzzolana, or trass, and blastfurnace slag, as seen by comparison of the analyses, is a well-known fact, blast-furnace slag has not been used commercially as a substitute for those cementing materials. The reason, the writer apprehends, lies in the fact that unless such slag is disintegrated or subdivided by rolls, the process must either be too costly, or the material is not in a fit and proper condition for the purpose. In order to produce a reliable slag cement, the slag must be ground together with the lime into an impalpable powder. The subdivided slag must, therefore, be perfectly dry, and, at the same time, friable. The stronger the hydraulic properties of the lime, the more reliable the slag cement will be, and practice has proved, that the slag from a gray-iron furnace gives the best results. The slag cement which has given the results shown in the annexed table, under pressure tests, was composed of six parts of slag, from a blast-furnace producing No. 3 foundry iron, and one part of lime, of medium hydraulic properties. The above-described class of cement bears storing as long as most Portland cements, and the cheapness of its production is self-evident.

By R. P. Alger, M. P. Tixier, D. R. Tanguy

In the combination induction-electrical log used at present in the field, the induction logging tool is appropriate for the investigation of moderately invaded formations. A new induction sonde with a radius and investigation about twice as large has been developed recently for the case of deep invasion. It has very nearly the same vertical resolution as the present sonde so that thin beds are defined as accurately as before. The characteristics of the new tool are described, the corresponding interpretation charts are given and field examples are discussed. The design of the sonic logging tool has been modified to improve the calibration and the reliability. The fact that porosity can be accurately recorded by means of the sonic log has prompted new interpretation procedures for saturation estimation, wherein the data concerning the various permeable beds in a given well are correlated. One approach consists of plotting transit time vs true resistivity, with an appropriate scale. With this approach, saturations can be estimated conveniently even in cases where formation water resistivity is not well known. In another approach, a comparison is made of the values of the formation waters computed from the re- sistivity and sonic logs. Using the concept of continuity, this procedure makes possible a quick determination of zones of saturation in shaly sands and/or in case of appreciable variations of formation salinities with depth. It has been found that the comparison of porosity from the sonic log with the apparent porosity computed from a short-investigation resistivity log may reveal, in many cases, the presence of residual oil and thus detect potentially productive formations; this procedure is valuable when the true formation resistivity and the resistivity of the formation water are in doubt. INTRODUCTION During the past year, the efficiency of log interpretation has been vastly improved. The improvements have largely resulted from the introduction of a deep-investigation induction device and from the application of new interpretation techniques that utilize sonic vs resistivity readings. Since the new interpretation techniques depend, in part, upon good values of true formation resistivity, the new induction log will be discussed under Part I. The sonic interpretation techniques will be studied under Part 11. Early in 1959, the 6FF40 induction equipment was introduced in the field. This device was designed for a better approach to true formation resistivities in deeply invaded zones. The greatly improved radial investigation of the 6FF40 equipment has been achieved without sacrificing vertical resolution. The first combination induction-electrical log, the 5FF40, was introduced as a standard tool in 1956 for the logging of wells drilled with fresh muds. The tool has received wide industry acceptance in the United States. 'The 5FF40 induction log has a radial investigation sufficient to overcome average depths of mud filtrate invasion. At 5d invasion, for example, the 5FF40 induction log will read about 1.4 R. in a water sand where R., = 10R,. At 10d invasion, such induction log would read 2.45 R. in the same water sand. In either case, the effects of invasion would not be sufficiently great to cause a water sand to be mistaken for a shale-free oil- or gas-producing zone. Some formations, however, invade deeply — in excess of 10d. Such water zones could be mistaken for oil-or gas-saturated sands unless the porosity balance' can clearly make the distinction. It is for these deeply invaded formations that the 6FF40 was developed. CHARACTERISTICS OF THE 6FF40 Radial Investigation Characteristics To describe the comparative responses of the 5FF40

By E. F. Kingsbury

One of the well-known and important uses of the noble or precious metals has been for electrical contacts. In fact, the elements of this group, comprising gold, silver and the six platinum metals, have been called "noble" because of their ability to resist heat and corrosion, properties which were early found to be of importance in good contacts. Of these, platinum and its iridium alloys, palladium, and silver, have been extensively used. In recent years the electrical industry has been responsible for about one-eighth of the total consiimption of the platinum metals in the United States. Information is not available as to what part of this one-eighth has been used for contacts but it is certainly an important one. Substitutes for Platinum During the past 20 years a considerable number of contact alloys have been proposed in the patent literature largely with the aim of more or less replacing platinum by gold or silver, and some of these alloys are finding increasingly widespread application as the reliability of their performance becomes established. From even the partial data at hand it appears that this substitution by the less precious metals has replaced at least one-half of the amount of platinum that the entire electrical industry would otherwise be consuming at the present time. While the subject is of considerable interest to the metallurgist it is even more so to the electrical engineer and it is surprising that the literature records no systematic and comprehensive investigation along these lines. Practically all the scattered work on noble metal contacts has centered around platinum as a standard and has had two more or less distinct aims. One of these has been the replacement of platinum by cheaper matcrials, and the other, the dcveloprnent of materials of superior electrical performance. The merit of a contact metal has been assumed to be dependent on its infusibility, and since thc more precious metals have the higher melting points there has naturally been an association of the two qualities in the sense that the more precious were

By M. P. Tixier

In several log interpretation methods, water saturation is evaluated by use of the ratio of the readings of a short spacing resistivity device and a long spacing resistivity device plus information on the mud filtrate and connate water resistivities, which is often derivable from the SP. These methods are valid over a certain range of conditions, usually specified by moderate invasion and R.. greater than R,. Since these methods involve no explicit evaluation of the formation factor, F, the saturafiion so found may be used in the standard Archie equation to derive a computed formation factor which, as a check, may be compared with the formation factor known from other independent measurements. Discrepancies in the two values of formation factor generally indicate that the ratio method is not within its range of applicability. If the computed F is too low, the corresponding S, is too low; if the computed F is too great, the corresponding S, is too great. By means of this "porosity balance" check, and some knowledge of the probable conditions of invasion, the interpretation can often be improved. The poqosity balance check is discussed for the cases of the Induction-Electrical Log Interpretation Method, the Rocky Mountain Method, and the R.,/R, Method. For the first method, a discussion is also given for the case of shaly sands. INTRODUCTION Al though many logging tools are available to obtain a practical value of porosity, the only logging method reflecting the water saturation is electrical logging. Water saturation is commonly computed from electrical log data by means of the basic Archie formula, Where F is the formation factor, R, the resistivity of the formation water, and R, the true resistivity of the formation. One difficulty in this procedure is a possible error in the evaluation of R, from the logs— for example, in the case of very deep invasion. There is no possibility, inherent in the procedure, which would provide a check of the accuracy of R, and, hence, of the computed saturation. Other approaches have been proposed which are based essentially on the ratio of the readings made with two dzerent devices: one device, with a short radius of investigation, provides a value close to the resistivity of the invaded zone; the second device, with a long radius of investigation, gives a value close to the true resistivity of the bed. By these methods, the determination of saturation is made independently of any direct knowledge of the formation factor. Therefore, if the formation factor (or porosity) is known from other sources (core analysis, neutron log, MicroLog . . .), it becomes possible to take advantage of this additional information in order to check the value of saturation derived from these latter methods. Inasmuch as this proposed check procedure relies on a knowledge of the value of formation factor or porosity, it has been deemed appropriate to call this verification a porosity balance. The purpose of this paper is to discuss the principle and the application of the porosity balance, depending on the method of interpretation being used. Basis of the Method There are in existence three interpretation methods by means of which the water saturation can bc found without a knowledge of the formation factor (or porosity). As just stated, they are based on the fact that the ratio of the two resistivities (one measuring the resistivity close to the hole, the other measuring the resistivity far from the hole) is independent of the formation factor since it affects the two resistivity measurements in the same way. Chronologically, the resistivity ratio methods are: (1) the Rocky Mountain Method, which has been used since 1943 and was explained in a paper in 1949;' (2) the R../R, Method or Shaly Sand Method, presented in 1954, which permits the calculation of the water saturation in both clean and shaly sands2; and (3) the Induction-Electrical Log Method, presented in 1957 for the interpretation of the induction-electrical log combination." Each method requires certain definite conditions in order for the interpretation to be reliable. The porosity balance method will let us know whether these conditions are fulfilled. In clean formations we know that the fundamental Archie relation (Eq. 1) must always be verified. This relation can be rewritten:

By William E. Souder

Behavioral accidents are a particular type of accident. They are caused by inappropriate individual behaviors and faulty reactions. Catastrophe theory is a means for mathematically modeling the dynamic processes that underlie behavioral accidents. Based on a comprehensive data base of mining accidents, a computerized catastrophe model has been developed by the Bureau of Mines. This model systematically links individual psychological, group behavioral, and mine environmental variables with other accident causing factors. It answers several longstanding questions about why some normally safe behaving persons may spontaneously engage in unsafe acts that have high risks of serious injury. Field tests with the model indicate that it has three important uses: It can be used as an effective training aid for increasing employee safety consciousness; it can be used as a management laboratory for testing decision alternatives and policies; and it can be used to help design the most effective work teams.

By Charles D. Litton

Experiments were conducted to determine the optical scattering properties of diesel particulate matter (DPM) and various combustion aerosols from both flaming and smoldering combustion sources at discrete angles of 151 and 301 in the forward direction and at a light source wave length of 635 nm using a simple light scattering module. In addition to the scattering data, simultaneous measurements were made of the total aerosol mass concentration; light extinction at an average wave length of 546nm; and the response of a common bipolar ion chamber typical of residential smoke detectors modified to allow the aerosols to flow through the chamber. The results of these experiments indicate, for DPM and combustion aerosols, the intensities per unit mass concentration depend not only upon whether the aerosolis DPM or combustion aerosol but also upon the type of combustion aerosol. The results also indicate that the ion chamber responses are greatest for DPM, followed by the response to flaming combustion aerosols (FCA)and lowest for smoldering combustion aerosols (SCA). For light scattering, the greatest intensities are found for SCA, followed by the intensities from FCA, and lowest for DPM. This report describes the experiments, their results, and the use of these results to developed sign criteria for early warning fire sensor capable of the rapid and reliable detection of fires in atmospheres that may or may not be contaminated by the products produced from diesel engines. © 2002 Elsevier Science Ltd. All rights reserved.