Search Documents

By J. S. Levine, M. Prats

A homogeneous and uniform cylindrical reservoir containing oil and gas is fractured vertically on completion and is produced at a constant bottom-hole pressure. The fracture has an infinite flow capacity, is of limited lateral extent and is bounded above and below by the impermeable strata defining the vertical extent of the reservoir. Results show that such a fractured reservoir can be represented by a reservoir of circular symmetry having very nearly the same production history. The well radius of this circular reservoir is about 1/4 the fracture length and is essentially the same as that obtained previously for a single fluid of constant compressibility. At the same value of cumulative oil production, gas-oil ratios of fractured reservoirs producing at constant terzinal pressure are larger than those of reservoirs having no fractures. This leads to more inefficient use of the reservoir energy in fractured wells and results in lower reservoir pressures for the same cumulative oil production. The reduction in operating life due to fracturing a reservoir is not as great as that for a slightly compressible fluid. This diflerence can be accounted for by the lower reservoir pressure in the fractured reservoir and its adverse effect on the average mobility and compressibility of the oil. As anticipated, the reduction in operating life increases czs the reservoir permeability decreases. The type of results presented in this report can be used to determine the economic attractiveness of fracture treatments per se, to setect the initial spacing to be used in developing a field, and to compare the relative merits of fracturing available wells and infill drilling. INTRODUCTION The effect of vertical fractures on a reservoir producing either an incompressible or a compressible liquid has already been discussed in the 1iterature.l,2 Those results indicate that the production history of such a reservoir is essentially the same as that of a circular reservoir having an effective well radius of approximately one-fourth the fracture length. The present work reports on the effect of a vertical fracture on a reservoir producing two compressible fluids —oil and gas—by solution gas drive. Because of the empirical nature of the PVT and relative permeability data used to obtain the performance of such reservoirs, results can only be obtained numerically and with the aid of high-speed computers. Since reservoirs lose their radial symmetry when fractured vertically, pressure and saturation can no longer be given only in terms of distance from the well. Two coordinates (such as x and y) must now be used to describe the pressure and saturation within the reservoir, and, since we are dealing with compressible fluids, time is also a variable. Thus the solution of a vertically fractured reservoir requires finding two unknowns (pressure and saturation) in two space variables (say x and y) and in time (t). Since no means are readily and generally available for solving such problems at the present time, we have used the results of previous work1,2 to approximate the effect of a vertical fracture on a reservoir producing both oil and gas by depletion. The purpose of the present wmk, then, is to investigate the possibility of using available numerical techniques (limited at the moment to one space variable) to study the two-space-variable flow behavior resulting from a vertical fracture. Results obtained in the course of this investigation are also reported and discussed. Input and output data of the numerical methods used are given in practical units: BOPD, feet, psi, cp, and md. Results are discussed fist in terms of specific reservoir and crude properties and geometries. Later, dimensionless parameters are introduced in order to extend results to different values of some of the reservoir and fracture properties. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The reservoir has an impermeable, circular, cylindrical outer boundary of radius r,. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below' the producing layer (reservoir). It is assumed that there is no pressure drop in the fracture due to fluid flow. 1 indicates the general three-dimensional geometry of the fractured reservoir. Gravity effects and the effects of differential depletion resulting from variations in hydrostatic head (pressure) will be neglected. Thus, the flow behavior in the fractured reservoir is described by the

By Robert E. Mclvor

Classification system efficiency (CSE) is a simple metric that can be used to measure and increase ball milling circuit efficiency. It is defined as the percentage of the ball mill energy being used to grind target (oversize) particles. This paper describes CSE analysis conducted for a variety of operations with related improvements in circuit efficiency ranging from 3 to 30 percent. The achievable CSE for any given ball milling circuit is described, providing the means to determine the extent to which every plant can benefit from increasing grinding circuit CSE.

By Anthony Iannacchione, Christopher Mark

Residents of Maxwell Hill, a suburb of Beckley, West Virginia, reported from January until March of 1988 a series of earth tremors. The tremors rattled dishes, and the noises associated with tremors were sometimes loud enough to wake people at night. These events were puzzling because Beckley is in an area of low seismic activity. Maxwell Hill is located above abandoned coal mines, though no direct subsidence damage to surface structures has been documented. A research program was initiated to determine whether subsurface ground movements associated with the abandoned mines could be responsible for the seismic activity experienced at the surface. A hypothesis is proposed which states that seismic disturbances can be caused by abandoned mine failure when two conditions are met. First, some instability must exist at the mine level. The most likely cause of instability is the failure of pillars of coal that were left for support at the time of mining. Second, subsurface movements can result in significant seismic activity only when the rock above the mine contains at least one very strong and brittle layer that is prone to fail violently through rupture or fault slip movement. Analysis of the data indicates that both hypothesized conditions for mine-related seismicity appear to have been satisfied in the case of Maxwell Hill. It has therefore been concluded that the earth tremors experienced by Maxwell Hill residents probably are related to the abandoned mines. Analysis of the potential magnitude of these events indicates that they are unlikely to cause significant surface danlage. The probability that future events will disturb the residence of Maxwell Hill should be quite low based on the identified mechanism.

By J. K. Whyatt, B. G. White

A comprehensive survey of mine seismicity and rock bursting during development of two sublevels at the Lucky Friday Mine, Mullan, ID, USA, was conducted to better define rock failure mechanisms and sources of ground control hazards. Survey data included rock burst damage reports, seismic event locations and magnitudes and, for the most energetic events, first-motion information. Several subsets of this database, including large seismic events, rock bursts, and microseismic activity near the face, were analyzed. Elements of the large-event and rock burst data sets were nearly independent. That is, there was no relationship between the risk posed by a seismic event and its energy, although such a relationship is well established for the mine as a whole. All data sets showed that certain geologic features appear to control the spatial distribution of events and the spatial distribution of rock burst risk. The data also suggest that, within the scope of this study. the greatest risk of injury occurs when a pocket of heightened rock burst risk is first encountered and that this risk is controlled when miners adapt their practices to these conditions. Recognition of the role of particular geologic features in the spatial distribution of rock burst hazards provides an opportunity for anticipating, rather than only reacting to, a changing level of rock burst hazard- However, much work will he required to make good on this promise.

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 M. D. Arnold, H. J. Gonzalez, P. B. Crawford

A method is presented for estimating the effective directional permeability ratio and the direction of maximum and minimum permeabilities in anisotropic oil reservoirs. The method is based on the principle that production from a well in an anisotropic reservoir results in elliptical isopo-tentials about the well, rather than circular. Bottom-hole pressure data from three observation wells surrounding a producing well are required to apply the method. The method involves fitting field pressure data to a set of general charts of isopotentials and making a few simple calculations until a solution is found. The method is based on a steady-state equation for homogeneorrs fluid pow. In addition to the method, a brief discussion of the theory underlying it is presented. INTRODUCTION The existence of a different permeability in one direction than another in oil reservoirs has been mentioned in several papers. Hutchinson' reported laboratory tests on 10 limestone cores and pointed out that one-half of them showed significant, preferential, directional permeability ratios, the average being about 16:1. Johnson and Hughesz reported a permeability trend in the Bradford field in the northeast-southwest direction with flow being 25 to 30 per cent greater in that direction. Barfield, Jordan and Moore -eported an effective permeability ratio of 144:1 in the Spraberry. Crawford and Landrum4 showed that sweep efficiencies could often vary by a factor of two to four, and sometimes considerably more, due to variations in flooding direction and patterns in anisotropic media. These findings indicate that the poss'bility of anisotropy may be worthy of consideration in the development of an oil field. In considering this, it should first be determined if anisotropy exists. If it does, the direction of the maximum and minimum permeabilities and the ratio of their magnitudes are quantities which can be of value in planning the most efficient well-spacing patterns. Past methods of determining these quantities have included analysis of oriented cores and analysis of flooding performance of pilot injection patterns. In recent work, Elkins and Skov5 resented an analysis of the pressure behavior in the Spraberry which accounted for anisotropic permeability. This work was based on the transient pres- sure distribution in a porous and permeable medium, with the solution expressed as an exponential integral function involving rock and fluid properties. The purpose of this study is to provide a method, based on steady-state equations, of estimating the direction and relative magnitude of permeabilities in an oil reservoir from field pressure data and well locations only. The method presented is based on work by Muskat6 which shows that Laplace's equation represents the steady-state pressure distribution for homogeneous fluid flow in homogeneous, anisotropic media if the co-ordinates of the system are shrunk or expanded by replacing x with it is desirable that data be obtained early in the history of a field because knowledge of an anisotropic condition would allow new wells to be spaced in such a manner that reservoir development and subsequent secondary recovery programs could be planned more efficiently. THEORETICAL CONSIDERATIONS A brief discussion of the theoretical basis on which the graphical solution was developed is presented in this section. Muskat's two-dimensional6 olution for the pressure distribution in an homogeneous, anisotropic medium with an homogeneous fluid flowing can be algebraically manipulated to show that the isobaric lines are perfect ellipses. The ratio of the major axis to the minor axis, a/b, is related to the permeability ratio, k,/k,, as follows. alb = dk,/k,--...........(1) It can also be shown that the pressure varies linearly with the logarithm of the radial distance from the producing well. However, the gradient along any ray is a function of the orientation of that ray, and a ..xiable is present when anisotropy exists which cancels out for a radial (isotropic) system. For a system such as that described, a dimensionless pressure-drop ratio was developed which is completely independent of the actual magnitude of the pressures. This was done by arranging Muskat's solution in such a way that aIl variables cancelled out except k,/k, and well positions. However, this solution depends on having a co-ordinate system with axes coinciding with the major and minor axes of the elliptical isobars. Thus, it was necessary to introduce a co-ordinate system rotation factor. The two unknown variables are then k,/k. and 0, and the two measured dimensionless pressure-drop ratios are related to the unknown variables as follows.

By L. A. Wilson, N. H. Harrison, R. R. Clay, R. L. Reed. D. W. Reed

Experiments have been performed in a linear near-adiabatic system for the purpose of extending data on reverse and forward combustion from atmospheric pressure to 1,000 psig. Results obtained from reverse combustion appear to conform qualitatively with the existing description of the process when reaction kinetics are suitably modified to account for the increase in pressure, i.e., increasing tbe pressure decreases the peak temperature and increases the combustion zone velocity. Forward combustion appears to be a fuel dominated process wherein peak temperature and combustion zone velocity are not very sensitive to changes in pressure. The moderate effects of pressure that do exist at low flux virtually disappear at high flux providing all oxygen is consumed. With this provision. increasing the pressure decreases the frontal velocity and increases the peak temperature. Results me shown graphicaIIy which demonstrate the effects of presswe on peak temperature, rate of advance, oil recovery, air/oil ratio, carbon oxides produced and temperature distributions. INTRODUCTION Numerous field tests of oil recovery using the technique of underground combustion are now in progress. Operating pressures used are always substantially in excess of atmospheric pressure. Nevertheless, the literature contains no laboratory data pertaining to the effects of pressure except for those of Martin, Alexander and Dew.1 Unfortunately, these data appear to reflect heat losses which may have obscured some of the effects of increased pressure. The underground combustion processes are exceedingly complex, and general concepts relating to them have necessarily been described in relatively simple terms. Particularly with regard to forward combustion, the interaction between multiphase fluid flow, heat transport and various rate mechanisms is so involved that to date a rigorous treatment is not available. In this paper a linear near-adiabatic system is described which was found capable of operating at pressures up to 1,000 psig and temperatures of at least 1,100 F. Using this equipment, experiments on both forward and reverse combustion were accomplished at elevated pressures and the data presented as a function of air flux or temperature. It will be assumed that the reader is familiar with the papers of Wilson, Wygal, Reed, Gergins and Henderson on forward combustion and that of Reed, Reed and Tracht3 on reverse combustion, so that it is not necessary to repeat the details of the two processes. EQUIPMENT AND PROCEDURE It was decided to study the effects of pressure over the interval 0-1,000 psig using a near-adiabatic system. Since satisfactory operation of such a system requires very thin tubing walls, it proved necessary to insert the entire combustion tube and heater assembly into a high pressure jacket in such a way that the internal tube pressure could be counterbalanced by pressurizing the annulus. This meant that the numerous heater and thermocouple wires had to pass through gas-tight seals that required considerable effort to maintain effective. An unexpected difficulty arose in connection with the nature of the annulus pressuring gas. It turned out that the composition of this gas critically affected the functioning of the adiabatic control system. Operation of the equipment at elevated pressures was not particularly difficult. The problems that did arise were a direct consequence of the large numbers of connections and seals involved. EQUIPMENT The equipment consisted of the combustion tube which fitted within the pressure jacket. This jacket was in turn mounted on a wheeled carriage which facilitated movement and provided for rotation of the

By R. E. Cook

This paper presents results of a study to determine to what extent errors in estimated free gas saturation affect the results of static pressure calculations from build-up curves in two-phase systems. Use is made of the method of pressure build-up curve analysis developed by Miller, Dyes and Hutchinson1 , with modifications for two-phase flow advanced by Perrine2. Cornhined fluid compressibility of gas-oil mixtures is included, with special attention given to the discontinuity that occurs in fluid compressibility at the saturation pressure and its effect on static pressure calculations from build-up curves. The sensitivity of static pressure calculations to errors in gas saturation is considered over wide ranges of formation and fluid properties and reservoir pressure conditions. It is concluded from this study that for analyses of two-phase pressure build-up curves (I) the use of under-saturated crude compressibility data, even when producing GOR's are at or near solution values, will result in very significant negative errors in calculated static pressrrre, and (2) except for large values of buildup curve .slope and/or very low crude gravities, rather large errors in calculated gas saturation can he tolerated without .significant effect on calculated static pressures. The conc1usions arc qualified as referring only to the effects of gas saturation and excludc possible errors in the nuumerous other factors that enter into the ancllysis of pressure build-up curves. INTRODUCTION As part of the analysis of basic data for several major Eocene reservoirs in the Bolivar Coastal field (Western Venezuela), it was necessary to analyze a large number of pressure build-up curves for oil permeability and static reservoir pressure. From performance data it became apparent that gravity segregation is very active in these reservoirs. As a result of this performance it was recognized that the free gas saturation within the oil leg of the reservoirs would be small but would be subject to certain errors in calculation due to the inaccuracy of relative permeability data in the low gas saturation region. Since all of the various methods of build-up curve analysis were originally derived for sin- gle-phase, constant compressibility systems, and since the presence of a small but rather indeterminate gas saturation would greatly affect the reservoir fluid compressibility, a study was undertaken to investigate the effect of gas saturation on the results of two-phase pressure build-up analyses. TWO-PHASE BUILD-UP CURVE ANALYSIS The method of analysis used in this investigation is basically that originally developed by Miller, Dyes and Hutchinson1, with modifications for two-phase flow advanced by Perrine2. Miller, et al, showed that the original equation for single-phase permeability away from the well bore, held for two-phase flow except that the calculated permeability is the effective permeability to oil away from the wellbore. This finding was acknowledged by Per-rine, who also presented the equation for effective gas permeability* as Since Eqs. 1 and 2 contain no terms influenced by fluid compressibility, no distinct problem should develop in calculating permeability from two-phase build-up curves. In regard to the calculation of static reservoir pressure from build-up curves in single-phase systems, the modified Miller, et al, method as presented by Perrine is briefly as follows. 1. From build-up curve data, PVT data for the reservoir fluid and the well spacing and sand thickness, calculate the climensionless shut-in time corresponding to the highest value of actual shut-in time (t, hrs) on the straight-line portion of a semi-log build-up curve. Dimenrsionless shut-in time is related to actual shut-in time by the relation, t- 0.000264 kt 2. From a graph of dimension less shut-in time vs dimensionless pressure difference (Fig. 1 1 ) find the value of dimensionless pressure difference corresponding to the value of 7 calculated under (1) above. Dimen-

By C. Bryce

Rhenium is primarily recovered as a byproduct from molybdenum roasting circuits. However, rhenium is also present in copper heap leach solutions at very low concentrations, and this represents a potential new rhenium source. Freeport-McMoRan (FCX) has developed a new process using carbon adsorption followed by ion exchange or solvent extraction to concentrate rhenium from less than 1 ppm to 15 g/L. The product is compatible as a feed to an existing rhenium processing plant. A pilot plant operated at Freeport-McMoRan?s Sierrita Mine has demonstrated the technical and financial feasibility of this new process. Copper heap leaching dissolves more than just copper from the copper ore. Several metals can be identified in leach solutions at concentrations of financial interest, rhenium in particular. Rhenium spot prices ranged from a yearly average of $1,070/kg ($483/lb) in 2005 (Magyar, 2008) to $10,400/kg ($4,706/lb) in 2008 (Polyak, 2013). The spot price in September 2013 was approximately $3,000/kg ($1,357/lb) of rhenium (Metal Bulletin, 2013). Several facilities at Freeport-McMoRan have very low rhenium concentrations in their copper leach solutions, ranging from 0.1 to 0.7 mg/L. The large quantity of solution available, if coupled with an attractive rhenium price, means that the rhenium contained in leach solutions may represent a significant value if a method of extraction could be developed. Testing has confirmed that rhenium adsorbs onto activated carbon much more selectively than other metals in the pregnant leach solution (PLS). A 100-gpm pilot plant was constructed to verify the lab-scale results on a larger basis and establish operating parameters for a full-scale operation. This report summarizes the initial technology development and pilot plant operation.

By B. Habermann

Artificially consolidated sand models, representing one-quarter of a five-spot, have been developed and used to study factors aflecting misciblt. displacrmenr. Sweep efficiency at breakthrough, size of the mixing zone between two miscible liquids and per cent area contacted by drive after breakthrough were determined tor the high mobility ratios encountered in actual resvrvoirs. Quantitative relationships between the degree of viscous fingering and mobility ratio were obtained by measuring the length of the fluid interface. Scaled miscible-slug experiments, supported by field evidence reported in the literature, have shown that when the slug is followed by dry gas the process is less efficient than expected. In addition to low areal sweep efficiencies encountered for high mobility ratio displacements, the effectiveness of a miscible slug is greatly reduced. This is caused by an accelerated growth of the mixing zone between the driving and the displacing fluids. INTRODUCTION Considerable interest in miscible displacement as a secondary recovery process has been shown by the oil industry in recent years. In addition to extensive laboratory investigations, numerous field operations have been initiated so that 39 miscible displacement projects had been reported' in the United States by early 1959. By far the most popular approach in these projects is the use of an LPG slug followed by dry gas. Laboratory tests have indicated that with a small miscible slug, amounting to several per cent of the hydrocarbon pore volume, practically 100 per cent recovery of the oil in place can be achieved and that all of the slug material can be recovered by a subsequellt gas drive. However, these experiments were conducted on small-diameter long cores that can be considered as essentially uni-dimensional reservoirs that do not represent actual field conditions. It soon became apparent that the inherently adverse mobility ratios existing in the reservoir cannot be ignored in any miscible displacement field project. When a less viscous fiuid such as an LPG slug drives a more viscous reservoir oil, fingers develop that reduce the areal sweep efficiency. The slze of the miscible slug to be used is a subject of some controversy. It has been proposed that the size of siug needed be calculated by setting lengths of the mixing zone proportional to the square root of the distance traveled. This proposal is based on experimental evidence' and theoretical considerations"' in which fingering was not a factor. A prerequisite for successful application of the miscible displacement process to field operations is a better understanding of the following factors. L What is the effect of mobility ratio on the rate of deterioration of a miscible slug? 2. How large a slug should be used for given reservoir conditions? 3. What is the sweep efficiency? The purpose of this investigation was to provide some additional laboratory data that should help to answer these questions. More specifically, it was intended to study more quantitatively the effect of viscous fiinger-ing, mobility ratio and slug size on the efficiency of this relatively new secondary recovery process by means of experiments on thin, horizontal, porous medium models.

By R. A. Salathiel, A. D. Kaiser, P. H. Monaghan, B. E. Morgan

The productivities of many hydrocarbon-bearing for-rnntions appear to be reduced during drilling, completing and repairing operations. The reason for these othcarwise unexplainable low productivities may be a drastic reduction of formation permeability immediately surrounding the borehole. A possible cause of the reduced permeability is swelling and dispersion of clays present in porous media when contacted by fresh water or mud filtrates. Hence, the invasion of fresh water or mud filtrate into hydrocarbon-bearing formations that contain clays may cause low productivities as a result of formation damage. Present laboratory investigations have yielded results which supplement existing literature on formation damage. An expanded-type clay, montmorillonite, way found in hydrocarbon-bearing formations not previously analyzed. When sand columns and formation cores containing montmorillonite were contacted by fresh water, their oil permeabilities were reduced. This reduction in oil permeability way restored partially but not completely by a repair procedure, i.e., flooding with kerosene containing organic agents or flooding with water containing electrolytes. Hence, where montmorillonite is present in formations, prevention of formation danzage was indicated to be the only complete solution to the problem. To prevernt formation damage, results suggested use of an approximately 1 per cent solution of calcium chloride or stronger solutions of sodium or potassium chloride. In conventional rotary drilling operations the wall of the hole is exposed to drilling mud for days, weeks, and in some instances for as long as several months. During this time, drilling mud filtrate seeps slowly into formations surrounding the hole. The permeability and stability of many formations, especially those containing clayey materials, are impaired by the mud filtrate. Often interaction between the mud filtrate and the formation causes drilling and production difficulties. Serious reduction in productivity of oil and gas formations during drilling, completing and repairing wells has been noted often in recent years. Cases are known of development wells which produced only small volumes of fluid upon completion. However, after a period had elapsed, they began to produce oil or gas in commercial quantities. In other cases, results of drill-stem tests taken during the drilling of deep wells have indicated good production should be obtained from shallow formations. Later difficulty has been experienced in obtaining production from these shallow formations which had remained in contact with drilling muds for an extended time. In yet other cases wells have failed to attain the productivities expected from previous performance after the use of drilling mud as a circulating fluid during workover operations. It is evident from the field behavior cited that the productivities of certain oil- and gas-bearing formations have been damaged by exposure to drilling muds. This reduction in well productivity, called "formation damage", causes serious consequences. In development wells, expensive remedial operations are incurred. In wildcat wells, good formations may appear to be marginal or even non-productive. The importance of formation damage has been recognized by a number of investigators1-s who have considered the effects of fresh water, saline solutions and aqueous mud filtrates on the permeabilities of formations. A general conclusion from such studies is that interaction between water and formation clays reduces permeability of reservoir rock and that the permeability decreases as salinity of the water decreases. Data in these references, however, concern the permeability change with only a single fluid present. A more pertinent and direct consideration in the production of oil and gas from wells is the permeability change when two fluids, oil and water, are present. Recognizing this, Nowak and Kreuger2 and Bertness1" considered several ways in which mud filtrates and saline solutions caused reductions in the oil permeability of formation cores. Investigations discussed here were designed to supplement information in the literature on changes in oil

By G. Thodos, W. F. Stevens

The point-source function introduced by Horner' us a solution to the general unsteady-state equation for the flow of fluids through porous media has been utilized to calculate pressure profiles for adjacent wells as a function of time. By trial and error, the time at which the pressure waves meet to cause interference can then be calculated. For this calculation, the time of interference is arbitrarily defined as that time when both pressure waves exhibit a pressure decrease of 25 psi ar the same point within the formation. To illustrate the method, an example involving two adjacent wells producing at different rates has been worked out and is presented in detail. INTRODUCTION The initiation of production of a new well is usually accompanied by pressure drawdown measurements in order to establish the extent and characteristics of the formation. The interpretation of such pressure drawdown data is based on mathematical developments involving a single well in an infinite reservoir. If more than one well is present in a reservoir, the possibility of interaction between wells exists. Therefore, care must be taken in interpreting the pressure drawdown data. It becomes important to be able to determine when this interaction, or interference, becomes significant. This paper presents a method for the prediction of the approximate time of interference between two such adjacent wells. The basic differential equation describing the pressure distribution existing during unsteady-state flow of a slightly compressible fluid has been solved by Hor-ner.' for a constant production rate from a circular homogeneous reservoir of infinite extent. The solution utilizes the point-source function to define the pressure at various points within the formation as a function of time. The resulting expression is: q = production rate, STB/D ,;= formation volume factor, dimension-less µ = fluid viscosity, cp k = formation permeability. md H = formation thickness, ft C - fluid compressibility, psi' F = formation porosity, dimensionless R = distance from centerline of well, ft T = production time, hours Ei(—u) = Ei-function of U Eq. 1 may be used to calculate the pressure profiles generated during production from a single well in an infinite reservoir.

By G. W. Tracy, R. D. Carter

Numerical calculations were made to determine the behavior of reservoirs with high-pressure drawdown and wide well spacing where the initial productivity is low and the wells are completed by hydraulic fracturing. The two-phase flow equations were solved for the flow into a single well. This well was assumed to be producing from a reservoir with hydraulically created horizontal fractures (four different systems with fractures were studied). For comparison purposes, additional two-phase flow calculations were made assuming a reservoir with uniform rock properties. The two-phase flow results were also compared with the conventional calculation methods, which do not include the effect of saturation gradients resulting from a simultaneous flow of oil and gas which are normal to this type reservoir. It was found that the conventional methods predict (1) a high and too optimistic value of ultimate recovery, (2) a high producing rate and a high reservoir pressure at a given oil recovery and (3) a low trend of gas-oil ratio with oil recovery. Included in the two-phase flow calculations were provisions to control the oil production rate by an allowable rate and, also, by a gas-oil ratio penalty rule. For the systems with hydraulic fractures, the producing rate was controlled by the gas-oil ratio penalty rule for most of the life. This is in contrast to the system with uniform rock properties which went "on decline" almost immediately. An unexpected characteristic of the systems which included fractures was the early rise in producing gas-oil ratio from 730 cu ft/bbl to approximately 1,200 cu ft/bbl, followed by a "leveling off" before the normally expected gas-oil ratio rise began. Additional features which are a result of hydraulic fracturing are (I) greater ultimate recovery, (2) higher average producing rates and (3) a lower average reservoir pressure at a given oil recovery. INTRODUCTION Some oil fields discovered during the past few years are producing from certain volumetric ally controlled reservoirs (often referred to as solution or internal gas- drive reservoirs) which are characterized by high-pressure drawdown at the wells. Since the available pressure drawdown at a well is limited by the static reservoir pressure and the producing rate is controlled by the available drawdown, wells completed in this type of reservoir usually produce at a rate less than the allowable from the time of completion. Because of this, this type of reservoir is referred to as a low productivity reservoir. Economic considerations require the use of wide well spacing and well stimulation by hydraulic fracturing to make commercial wells in this type of reservoir. Performance predictions for volumetrically controlled reservoirs have been made using a combination of two standard equations. 1. The "Schilthuis" or "Muskat" type material balance equation is used to relate the average reservoir pressure and the cumulative oil recovery. 2. The results from the material balance equation and the productivity factor as described by Pirson' are used to relate the cumulative recovery with producing rate and time. The material balance equation assumed uniform pressure and liquid saturation conditions throughout a reservoir. The steady-state radial flow formula allows for a pressure gradient toward a well but assumes uniform liquid saturation. These calculation methods are adequate for application to reservoirs wherein the drawdown at the well to realize satisfactory producing rates is small compared to the total pressure. In low productivity, volumetrically controlled reservoirs, the pressure drawdown at the well is large cornpared to the total pressure. Although a precise number cannot be given for the magnitude of a large pressure drawdown, values in excess of 1,000 psi would definitely be included. For practical considerations, this usually occurs when the formation flow capacity is less than about 100 md-ft. However, this limit of formation flow capacity will vary with the well producing rate. The low pressure in the neighborhood of the well which results from a high drawdown causes evolution of large volumes of gas. This causes the gas saturation to be higher near the well than at a greater distance—-hence, a non-uniform gas saturation. Also, the relationship between the relative permeability to oil (K/K) and gas saturation is nonlinear but decreases approximately in an exponential way with increases in gas saturation. Because of this, the following chain reaction is established.

By P. B. Baxendell

The performance of a water-drive reservoir having a gas cap depends primarily on the movement of the gar-oil and oil-water contacts. The movement of the contacts during production depends in turn on fluid withdrawals and how the reservoir pressure changes as fluids are produced from the reservoir; that is, on how effectively the aquifer maintains pressure by replacing withdrawals. inasmuch as pressure changes, fluid withdrawals, contact positions, and produced gas-oil and oil-water ratios are interdependent, the analysis and prediction of the performance of a reservoir produced in a given way must take into account this interdependence throughout depletion. This paper presents an analysis which, within the limitations of the assumptions made, yields an engineering approach to predicting future performance based on reservoir pressure and production history. The most significant assumption is the method of extrapolation of future gas-oil and water-oil ratios. The extrapolation procedure smoothly increases both ratios to preselected values as the remaining oil column undergoes a specified decrease in thickness. This preselection is made on the busis of previous field experience in depletion of similar reservoirs under similar conditions. For computing future performance a vdumebic balance L combined with the differential equation defining pressure distribution in the aquifer to obtain positions of water-oil and gas-oil contacts. From these positions are extrapolated produced water-oil and gas-oil ratios. Reservoir performance can be investigated when oil production rates are dependent upon various factors including the performance of the reservoir itself. Examples of practical application of the procedure are included. INTRODUCTION To predict performance of water-drive reservoirs with gas caps and thin oil columns, it is necessary to describe the motions of the fluids within the reservoir during the entire production period to depletion. These motions depend on the pressure changes and on the withdrawal of oil, gas and water. Production of oil from the oil zone primarily causes water to move in to take its place; production of gas from the cap tends to cause oil to migrate into its place. In addition, the volumes of oil and gas remaining in the reservoir depend on changes in the pressure, since a decrease in pressure causes fluid expansion, gas liberation and oil shrinkage. A method of relating the future pressure to total withdrawals is used to describe the motions of fluids within the reservoir under conditions arising in possible modes of production. The analysis is based on two relationships and will reduce the problem to one amenable to digital computation. The first of these concerns the dependence of the pressure distribution in the entire aquifer furnishing the water drive upon the total reservoir withdrawals. This dependence is dictated by the permeability distribution and the extent of the aquifer; at present, such information is most readily obtained from the performance history by means of the resistancecapacitance reservoir analyzer. The second relation involves withdrawals, pressure in the reservoir, and movement of oil, gas, and water within the reservoir. The analysis is subject to certain simplifying assumptions that are necessary to permit solution of the problem. A comparison of the methods of this paper with those presently practiced is pertinent. One method of analysis is to use the reservoir analyzer to predict reservoir behavior based on aquifer characteristics determined from production history. Inasmuch as the total withdrawal rate depends upon the gas-oil and water-oil ratios, which depend in turn upon, among other things, the positions of the gas-oil and water-oil contacts, these ratios must be assumed in advance. In this paper these ratios are, instead, related to the computed positions of gas-oil and water-oil contacts. Thus, our method can be applied to problems in which the oil production rate is limited by produced gas-oil ratio, or to problems in which it is desired to determine the variable gas injection rate that will maintain the gas-oil or water-oil contact stationary. These problems cannot be worked satisfactorily on the analyzer. Another method presently used is based on a paper by Hurst.' Since his procedure is dependent upon using the solution of the heat flow equation, which requires constant permeability within the aquifer, and our procedure recognizes variations of permeability in

By Thomas E. Marshall, Michael Dunn, L. J. Prosser, Anthony T. Iannacchione

During the past few years, the Pittsburgh Research Laboratory of the National Institute for Occupational Safety and Health (NIOSH) examined and characterized conditions at a majority of the underground stone mines in the United States. Observations at these mines revealed a limited degree of roof monitoring beyond visual inspection and sounding. When monitors are used, they typically require the miner to measure movement at the roof. If conditions are unstable, the mine rmay be in a hazardous situation while recording data. Based on this scenario, researchers surmised that a simple, inexpensive monitor with the capability of recording data at a distance from the mine roof would be a safer way to gain the information. Additionally, more widespread use of monitors could potentially lead to better understanding of roof movement in general. A monitor to meet this need was designed, tested, and subsequently improved as experience was gained in its’ use at a number of underground stone mines. The Roof Monitoring Safety System (RMSS) can provide an initial indication of movement in roof beams. By understanding and measuring roof movement in underground mines, the potential for injuries and fatalities to mine workers from falls of ground can be reduced. Also, officials at a mine with a history of data are better prepared to make a decision on remedial actions in the event of ground falls. This paper will outline the evolution of the RMSS and how it can be used in a comprehensive proactive ground control safety program. Also included is a case history describing how the RMSS was used in an evaluation of the effectiveness of a mechanical impact scaling machine at an operating limestone mine.

By A. A. Stripling

The Cretaceous limestone wells of the Mara/Maru-caibo Dist. of Venezuela are extremely prolific producers. To maintain production on cessation of natural flow, large scale gas-lift operations were commenced involving high production rates on caring Flow. For the design of these gas-lift installations it was essential to have some knowledge of the pressure gradients involved in casing flow, so that the required injection pressures, optimum gas injection raws, etc., could be forecast. This paper presents a method of calculating these annular pressure gradients. Basically the method makes use of an energy balance equation coupled with an empirical energy loss factor derived from field data. The calculations have been put in a form suitable for "punch card-type" calculating machines. A set of gradient curves for Mara/Maracaibo conditions covering flow rates up to 12,000 B/D, and GOR's of 500 to 2,000 cu fd/bbl is presented. From these curves the flowing BHP can be predicted from surface data and the vertical pow performance under varying conditions of gas injection can be estimated. The method has been applied to La Paz and Mam gas-lift operations for over three years and has given consistently accurate results. INTRODUCTION The Cretaceous limestone fields of the Mara/Maracaibo area in Venezuela are characterized by the presence of widely distributed fissure systems. Under these conditions individual wells frequently have extremely high potentials and such wells are generally flowed on the casing-tubing annulus in order to reduce the pressure loss in vertical flow. To maintain production on cessation of natural flow, an extensive gas-lift system was planned. As a first approach to this gas-lift design problem, an attempt was made to construct a set of empirical flowing gradient curves from the available subsurface pressure data. Such a method has been used successfully for tubing-flow conditions in California by Gilbert'. Unfortunately, the empirical construction of casing-flow gradient curves presents additional difficulties, as pres- sure bombs cannot bc run in the casing-tubing annulus against high rates of flow. Thus, the Full flowing gradients cannot be defined by actual measurement and the best that can be done is to obtain spot pressures at the tubing shoe. Under these conditions a vast accumulation of BHP data is required to construct a set of gradients. To the La Paz and Mara fields the majority of the available flowing pressure measurements fall into a relatively narrow GOR range, in the neighborhood of the solution ratio. While there appeared to be reasonable prospects of constructing gradient curves for this limited range, measurements at higher GOR's, which were of major interest for gas-lift operations, were relatively few and the prospects of constructing empirical curves cor~spondingly remote. As an alternative approach, the possibilities of a mathematical analysis of vertical flow conditions were investigated. May and Laird2fV n 1934 presented an excellent analysis of the vertical-flow conditions in the Anglo-Iranian Oil Co.'s fields. Their analysis was based on an energy-balance equation and calculated the losses due to kinetic energy, friction and slippage separately. By combining these losses with the total calculated available energy, the required depth-pressure traverse could be deduced. On applying the method to La Paz conditions, good results were obtained for tubing flow at high tubing-head pressures. However, for lower wellhead pressures (300 psi and under) the discrepancy between calculated and measured values increased appreciably. It was considered that this was probably due to energy losses caused by slippage or liquid hang-up in a continuous gas phase, for which it was difficult to visualize any analytical approach. In 1952 Poettmann and Carpenter4 presented a method of pressure gradient calculation which, while based on an energy-balance equation similar to that employed by May and Laird, made no attempt to evaluate the various components making up the total energy loss resulting from irreversibilities of flowing fluids. Instead, they proposed a method of analysis which assumed that all significant losses for multiphase flow could be correlated in a form similar to that of the Fanning equation for frictional losses in single-phase flow. They then derived an empirical relationship linking measurable field data with a factor which, when applied to the standard form of the Fanning equation, would enable the energy losses to be determined. The total available

By D. R. Wieland, A. N. Barron, A. R. Hendrickson

A definite relationship has been found between the reactivity of flowing hydrochloric acid and its shear rate in a carbonate fracture. Both flow velocity and fracture width affect the acid reaction rate. Laboratory studies were conducted on acid reactivity at different flow velocities through horizontal-linear fractures, using 15 per cent hydrochloric acid at 80' F and approximately 1,100 psi. Fracture width varied from 0.02 to 0.20 in. These data provide a new basis from which the spending time and penetration of the acid can be estimated. Equations were derived expressing the relationship between injection rate, fracture width, acid concentration, time and fracture height, for linear and radial fracture systems. Because the penetration of the acid before spending is closely related to the extent of productivity increase resulting from an acidizing treatment, these data provide a valuable insight into some of the controlling factors that must be taken into consideration during treatment preplanning. INTRODUCTION Acidizing of carbonate reservoirs to improve production characteristics has been widely practiced since 1932. Originally, it was assumed that the acid uniformly penetrated natural formation pores and flow channels, enlarging them and thereby increasing their flow capacity. Little consideration was given to the reaction rate of the acid, or how far it would penetrate away from the wellbore into the formation, before spending. It has been shown recently' that, unless fractures are present in the rock, very little penetration is attained before spending, and the benefits of the acidizing treatment are largely confined to the immediate vicinity of the well-bore. Therefore, acid treatments may be classified into two categories: (1) matrix acidizing, in which the acid flows through multiple: formation pores; and (2) fracture acidizing, in which the bulk of the acid travels through fractures in the rock, whether natural or induced. When acidizing treatments are conducted at pressures of sufficient magnitude to open and extend such fractures, it is often desirable to inject a propping agent to hold the fracture open after the treating pressure has been released, thus providing a highly conductive flow channel through the rock. In some cases where acid attack produces surface irregularities, the resulting flow passages can be sufficient to provide high conductivity without use of a propping agent. During injection, the acid dissolves the carbonate rock with which it comes in contact until it is spent. Deeper penetration of the spent acid into the formation produces no appreciable benefit (if no propping agent is used) because the unetched fracture faces will rejoin when treating pressures are released and the fracture will "heal", with negligible resultant conductivity. The spending time of the acid thus becomes important in determining how far from the wellbore the improved-conductivity zone extends. The spending time of acid after injection into carbonate rock depends on the rate at which the acid reacts with the rock. This in turn is controlled by a number of factors, as previously reported. 2-4 These include temperature, pressure, acid concentration, rock composition, injection rate and the area-volume relationship between the acid solution volume and the surface area of the formation flow channels through which it penetrates. Thus, in matrix acidizing, where an extremely large area is exposed to the acid, spending is rapid. In contrast, spending time is prolonged in open fractures, where the area-volume ratio is much lower. Thus, greatest penetration of the acid before spending will be achieved in acid-fracturing treatments where fractures are held open by hydraulic pressure. Perkins and Kern5 have presented equations for fracture-width determinations which, for hard formations such as limestone or dolomite, indicate that high injection rates and viscous fluids are required to create wide fractures. Increasing the injection rate will, in itself, produce faster reaction rates. Because this type of reaction can be considered to be first-order diffusion-controlled, '-' the rate of movement of fluid past the rock surface affects the thickness of the diffusion layer, in turn affecting the reaction rate. Staudt, et al,10 report the results of tests in which a rotating marble cylinder was immersed in acid, and the weight loss at different speeds determined. Such test conditions, however, are not truly indicative of the situation in a formation during an acid-fracturing treatment. This paper reports reaction-rate data obtained under

By R. E. Schilson, F. H. Poettmann

The direct determination of the stabilized performance behavior of low capacity, slowly stabilizing gas wells is extremely time-consuming and wasteful of gas. From both field experience and theoretical considerations, a test procedure has been evolved by which the stabilized hack pressure behavior of such gas wells can he predicted without having to revert to long time flow tests. The method consists of using the isochrona1 test procedure to establish the slope of the back pressure curve, "n", and the short time variation of the performance coeficient, "C", with time. From this short time transient flow data and theoretical considerations, the value of C at large times can he established. By assuming the radius of drainage of a well to be half the distance between wells, one can calculate the stabilization time for various well spacing patterns. Once the stabilization time for a given .spacing has been determined, the value of C can be calculated and the stabilized back-pressure curve can he. establbrhed. The calculated perfortnance coefficient as a function of tinze was cotnpared to the experintentally measured values for a number of gas wel1.e. The deviation of the calc~*lated from the cxperimental res1tlt.s vary depending on the set of short time experimental points used to evabrate the parcrttzeters of the equation. The longer the tirne far the flow test data user1 in the calculations, the better was the agreement with the experimental results. The time necessary to obtain this data from well tests varies consitlerably, depending on the physical natrcre oj the rt3scrvoir under consideratiot~. INTRODUCTION For many years, the U. S. Bureau of Mines Monograph 7' has served as a guide for testing and evaluating the performance of gas wells by means of the back-pressure method. The back-pressure performance of a gas well is expressed by the following equation: Q~ CW-W.........(I) where the characteristics of the back-pressure equation are determined by C, the performance coefficient, and n, the exponent which corresponds to the slope of the straight line when Q and {P* - PN2) are plotted on logarithmic paper. Q is gas flow rate at standard conditions, and P, and P, are equalized and flowing bottom-hole pressures, respectively. Prior to the development of the back-pressure tcst, the "open flow" capacity method of testing a well was common. By this method, the wells were flowed wide open to the air and the flow rate measured. Such procedure was wasteful of gas and did not provide information on the deliverability of the gas to the pipe line. Monograph 7 Procedure The back-pressure method of testing wells was developed to overcome these shortcomings. Although much has been learned regarding the laws of the flow of gas through porous formations, the original development of the back-pressure relationship was based entirely on empirical methods. The back-pressure behavior provides the engineer with information essential in predicting the future development of a field. It permits him to calculate the deliverability of gas into a pipe line at predetermined line pressures, to design and analyze gas gathering lines, to determine the spacing and number of wells to he drilled during the development of a field to meet gas purchasers' requirements, and to solve many other technical and economic problems. As described in Monograph 7, the flow-after-flow method of back-pressure testing, when applied to fast stabilizing and usually high capacity wells, correctly characterized the behavior of the well. However, as the value of the gas at the wellhead increased, small capacity gas wells having slow rates of stabilization became economically operable. The flow-after-flow method of testing could not be used to describe the behavior of these slowly stabilizing wells. The procedure of Rawlins and Shellhardtl for establishing the back-pressure behavior of a gas well was based on the rcquirement that the data be obtained under stabilized flow conditions; that is, that C is constant and does not vary with time. C depends on the physical properties of the reservoir, the location, extent and geometry of the drainage radius, and the properties of the flowing fluid. In a highly permeable formation, only a very short period of time is required for the well to reach a stabilized condition, and, consequently, the requirements for the test procedure described in Monograph 7 are met. For a given well, n is also constant

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. L. Derick

On a regular basis, mine personnel should obtain air samples directly from active underground areas as well as remotely from sealed locations. These samples provide mine officials with data relative to the makeup of the various atmospheres. Should conditions underground change, these results can provide a baseline that can assist officials in determining the extent of change. Relating the concentrations of individual gases to time, applying gases to equations, and examining their change over time are tools that have been successfully used to determine if a heating event exists and, if so, the extent of the emergency. These atmospheric status equations have either been developed for the mining industry or adapted from other industries. They serve as tools to assist personnel in determining the condition of underground mine atmospheres. As long as the problem area continues to be ventilated, certain gases, gas ratios, and equations can be used to determine if a fire exists, what type of material could be burning, and what quantity of coal is burning. Once the mine or the involved section has been temporarily sealed, gases, ratios, and gas-derived equations can assist personnel in determining how conditions within the sealed volume are changing over time. To alleviate potential explosions when fresh air is once again introduced into the sealed volume, equations and a diagram have been developed to predict the likelihood of an explosive atmosphere. Similarly, prior to commencing recovery of a mine or section, equations exist that can assist in determining if the fire will rekindle when air is reintroduced. This paper serves as an overview to remind and/or instruct readers about gas-sampling methodologies and gas analyses to assist in determining the status of underground atmospheres.