Search Documents

By J. J. Arps

This paper reviews the methods currently in use for estimating primary oil reserves and discusses the principles on which these methods are based. Particular emphasis is placed on how these methods change with the type of information available during the life cycle of an oil property. This paper contains various novel estimating methods and shortcuts heretofore unpublished. INTRODUCTION Estimating oil reserves is one of the most important phases of the work of a petroleum engineer since the solutions to the problems he deals with usually depend on a comparison of the estimated cost in terms of dollars, with the anticipated result in terms of barrels of oil. His recommendations to management regarding the best course of action are therefore normally based on the most favorable balance between these two. Specific engineering problems which require such a knowledge of recoverable oil reserves and a projection of future rates are: [a.] the exploitation and development of an oil reservoir; [b.] the construction of gasoline plants, pipelines and refineries; [c.] the division of ownership in unitized projects; [d.] the price to be paid in case of a sale or purchase of an oil property, and the magnitude of the loan which it will support; [e.] the proper depreciation rate for the investment in oil properties; and [f.] evaluation of the results of an exploration program. This discussion will be confined to the various methods and tools which are currently in use for estimating oil reserves to be obtained during the primary phase of an oil-producing reservoir and for a projection of the future production rates. Reserves which may be obtained by secondary recovery methods or fluid injection programs and gas and gas condensate reserves will not be discussed in this paper. Unfortunately, reliable oil reserve figures are most urgently needed during the early stages when only a minimum amount of information is available. Management's interest in the oil recovery from a property—aside from its use for accounting purposes—usually declines when the property approaches its economic limit, just at the time when the reliability of the estimates is at its best. To illustrate this general idea a chart is presented (Fig. 1) showing the three periods in the life of an imaginary oil property. Time is shown on the horizontal axis, while the cumulative production and estimated ultimate recovery are plotted vertically. No particular units are used and this schematical chart is not to scale. During the first period, before any wells are drilled on a property, any estimates will of necessity be of a very general nature, based on experience from similar pools or wells in the same area, and usually expressed in barrels per acre. This will

Jan 1, 1957

By P. H. Dudley

IN giving a brief account of the experiments in progress to inquire into some of the facts in regard to "railway resistances," recently commenced upon the Lake Shore and Michigan Southern Railway, with the dynagraph (which is essentially an instrument recording upon paper the force required to draw a train, at all speeds, marking the time in short intervals, so that the force absorbed due to any

Jan 1, 1876

By E. F. Poncelet

Forty years ago A. A. Griffith developed a theory explaining why brittle materials displayed such low tensile strengths.' He based his views on two points. First, he found himself compelled to assume that all brittle materials are replete with flaws, cracks, and other defects that act, although quite invisible, as large stress raisers. Second, he applied the "theorem of minimum potential energy," which says that the total potential energy of a system must pass from the unbroken to the broken condition by a process involving a continuous decrease in potential energy. By this means he satisfactorily accounted for the noted low strength of such solids and also for the wide spread obtained in experimental measurement of these strengths. So successful has the theory been that it is favored by some to this day. Unfortunately this theory is of limited use beyond the explanation of these two noted phenomena and it is keenly felt that a better theoretical insight into the physics of the fracturing process is needed as the volume of experimental evidence accumulates. The author proposes in the following to build on the fundamentally sound concepts of Griffith and, with the help of increased theoretical knowledge over that available to Griffith, develop a mechanism for frac-ture which will provide far greater understanding of the experimental evidence accumulated to date than is possible from the original Griffith idea. THE GRlFFlTH THEORY Very little progress indeed can be made without accepting the first postulate of Griffith which supposes all brittle solids to be full of microcracks. It would be difficult indeed to find a better mechanism for the small strength of such brittle materials, in conjunction with the fact that the energy that must be expended for comminution is by no means small. The postulate of the existence of the microcracks permits the breakup of the various bonds a few at a time by concentrating the stress at the tip of the progressing crack, while the total energy expended is the same as if they all had been ruptured simultaneously. The only flaw in the argument is that no reasonable explanation has been proposed to account for the genesis of such cracks. Indeed their very presence is in violation of the Griffith second postulate, the potential energy theorem. This theorem is straightforward for isothermal processes, and, in spite of Griffith, there is some doubt that treating the problem isothermally is legitimate. The surface energy of bodies is a free energy, not a potential energy as stated by Griffith, and the production of new surface free energies is not necessarily an isothermal process. There is ample evidence to the contrary. Generally speaking, if heating a body increases its surface area, then, by virtue of Le Chatelier's principle, any increase of that area by other means will tend to lower its surface temperature. Lord Kelvin calculated the actual cooling that resulted in drawing out a film of liquid.2 R. A. Houston calculated the surface cooling that resulted in stretching a metal wire.3 These calculations were made by applying the Carnot cycle to the process and evaluating the thermodynamics thereof. IRREVERSIBILITY OF THE FRACTURING PROCESS While Griffith was very careful not to say so, the impression gained from studying his papers is that he considered the fracturing process as reversible, that is, a succession of quasi-equilibrium states. There is ample evidence that it is not. The indication that the new surfaces produced by the propagation of a crack are cooler than the original body points to an irreversible heat flow from the interior to the new surfaces to equalize the temperatures. If the process be reversible, any crack accidentally formed should immediately close up as, in the absence of any strain energy, the potential energy would thereby be lowered. The fact that they do not, constitutes a paradox. Such paradoxes are nothing new where certain phenomena that propagate from minuscule nuclei are assumed to be reversible. Such is, for instance, the condensation of a pure saturated vapor that is suddenly chilled by adiabatic expansion. At the beginning the tiny droplets that are formed should be only a few angstroms in size, but the vapor pressure at such droplets is so high that they should evaporate at once. A similar situation arises if a saturated pure solution becomes super-saturated upon cooling; the first tiny crystal nuclei should dissolve as fast as

Jan 1, 1964

By O. G. Kiel, G. W. Swift

The concept of "a continuous succession of steady-states", which has been applied successfully by Aromfsky and Jenkins to obtain a solution for the nonlinear partial differential equation describing the transient Darcy flow of gas through porous media, is demonstrated to be equally valid for transient non-Darcy flow. A mathematical model, which numerically solves the partial differential equation, is used to check the validity of the succession of steady-states solution. Comparison of sand-face pressure histories compelted by the two methods shows excellent agreement. The utility of the succession of steady-states solution in predicting performance of gas wells rests in the fact that no special computation equipment is required. The development of the succession of steady-states solution leads also to a practical method for determining and analyzing field test data. A method for taking gas-well test data under constant-rate conditions is presented. Experimental data obtained in the field by employing the constant-rate method are presented and analyzed in accordance with the succession of steady-states solution. Analysis of data in this fashion is demonstrated to give direct "in situ" information for reservoir permeability, porosity and turbulence or non-Darcy coefficient. INTRODUCTION The economics of gas production are dependent upon the transient behavior of flow within the reservoir. For production from a finite reservoir, the transient flow behavior can be subdivided into two parts. At first, the transient caused by the movement of the pressure "wave" into the reservoir is of importance. Later in the production history, the pressure-wave movement ceases and the second transient stage of material depletion becomes controlling. For reservoirs of relatively high permeability, it can be shown that the pressure wave moves into the reservoir and stabilizes quite rapidly. In the case of relatively impermeable reservoirs, quite the opposite is true. Although it is theoretically possible to compute the production capability of a well from the properties of the reservoir as determined by static tests and core analyses, more reliable information is obtained by conducting flow tests on the well and thereby obtaining some measure of "in situ" formation properties. For gas wells, there are two basic types of tests in existence: the flow-after-flow method,' and the isochronal method.' Both of these techniques are tailored to obtain data that can be analyzed in accordance with the empirical performance equation: In addition, the isochronal method makes provision for the sluggish nature of pressure-wave movement in "tight" formations by requiring pressure build-up between flows and by stipulating that data obtained on successive flows be analyzed at equal elapsed flow times. It can be demonstrated that either test is valid for reservoirs of high permeability. Further, since it has been pointed out that the pressure wave stabilizes rapidly for reservoirs of this type, tests of relatively short duration will give stabilized information on the performance of a well. Further decline of sand-face pressure and/or production rate may be determined by employing material-balance techniques. Cullender' points out that for relatively impermeable reservoirs the flow-after-flow method gives invalid results. (See Appendix C.) If the isochronal method of testing is used, there are two alternatives: (1) the tests must be conducted for a sufficient length of time to obtain stabilized information (which may require months to accomplish); or (2) some method for extrapolating the results of short-term isochronal tests must be employed. The first alternative is impracticable because of manpower, conservation and economic considerations. Recourse to the second alternative requires some assurance regarding the reliability of the extrapolation technique. Poettmann and SchilsonJ present an empirical method for predicting stabilized performance. The present investigation was originally initiated to determine the reliability of this technique. To do this, a mathematical model was developed to simulate the Darcy and non-Darcy flow of gas through porous media. The model consisted of a finite-difference approximation of the nonlinear partial differential equation which was solved on an IBM 7090 computer. Long-term production histories were simulated by the model and compared against predictions obtained from the Poettmann-Schilson method. As the work progressed, it became apparent that a straightforward predictive equation could be developed by utilizing the concept of a succession of steady-states. As a result, the emphasis of the work was redirected to exploit the advantages of the new method.

By A. W. Schlechten, C. M. Hsiao

The apparent vapor pressures of a number of metal sulphides were determined by measuring their rate of weight loss when they were heated under vacuum. The calculated pressures are due in some instances to the volatility of the sulphide and in others to the decomposition pressure. FIE direct reduction of metallic sulphides, as in the precipitation method for lead or antimony, is not a widely used process today. However, it might be employed for the production of a number of metals if the operation were carried out under reduced pressures. Two general procedures are possible; the reducing agent could form a stable nonvolatile sulphide from which the valuable metal could be removed by volatilization or liquation. The other approach would be to choose a reducing agent which would form a sulphide that could be volatilized away from the valuable metal. A vacuum operation would serve to protect the sulphides and metals from oxidation and would increase the evaporation rates. A related process is the separation of metallic sulphides by selective volatilization either at atmospheric pressure or under vacuum. This procedure is already used commercially for the elimination of lead and cadmium as PbS and CdS from zinc calcines; it might well have other applications. To predict such processes it is necessary to know the relative volatilities of the metallic sulphides and to know whether they volatilize as sulphides or whether the predominant weight loss is due to decomposition. Such data would also be of value in studying problems such as the use of sulphides for refractories, as described by Brewer,] or in the evaluation of the theory of ore deposition from the vapor state, as proposed by Brown.' Qualitative statements concerning the volatility of metallic sulphides can be found frequently in the literature and are well summarized by Mellor, but accurate determinations of rates of evaporation or vapor pressures are very scarce. Kelly" lists vapor pressure data only for PbS; his reference being to the work of Schenck and Albers,4 who measured the vapor pressure of PbS by a dynamic method between 850" and 995°C. More recently Veselovskii5 has determined the vapor pressures of the sulphides of antimony, lead, cadmium, and zinc at temperatures ranging from 360" to 1000°C using Knudsen's method. A second Russian worker, Pogorelyi,6 determined the vapor pressures of CdS and ZnS at temperatures from 900" to 1250°C. The results obtained by these workers are shown in subsequent graphs in comparison with the data of this paper. Brewer and coworkers1 have studied some of the metallic sulphides to determine their suitability as refractories in vacuum. Veselovskii's determinations on PbS checked well with those of Schenck and Albers, and his determinations of CdS check almost exactly with Pogo-relyi. For the vapor pressure of ZnS the data of Veselovskii and that of Pogorelyi are in fair agreement, as shown in Fig. 1, but it should be pointed out that Lumsden7 has discovered that the equation given by Veselovskii for the vapor pressure of ZnS: - 14,200 log P9.495 is in error by a factor of ten when compared to the experimental data. In contrast with the scarcity of data on vapor pressures, the decomposition pressures of many metal sulphides have been determined and are summarized by Smithells.8 Theoretical Considerations The experiments described in this paper consisted of determinations of the rate of weight loss of certain metallic sulphides under very low pressures at different temperatures. These data plus observation and analyses of the products reveal the mechanism of the weight loss and also provide a basis for the calculation of the vapor pressure or apparent vapor pressures of the sulphides.

Jan 1, 1953

By F. C. Stanford

THE work of an engineer is to direct natural forces so that they bring about the results that he wishes to secure. Heretofore he has concerned himself chiefly with physical forces and inanimate objects, until he has secured a truly marvelous degree of control over them. The great problem of the present age is to secure a like degree of knowledge and control of the human factor. Progress has already been made along this line, but, it is only a beginning. Wilhelm Ostwald says that scientific knowledge enables a man to be a true prophet; but it must be remembered that a good many people do not seem to be able to maintain a high average in prophecy. In short, the average of knowledge is not equal to the average needs of the present time, but continually fails to keep pace with it. It has been aptly said that the really essential need of the period is for more brains; if we only had sufficient knowledge and sagacity we could not only solve all difficulties, but prevent them from arising. This generalization holds good for an industrial organization. The most vital need of the management is for more brains. The logical course for an engineer in charge of an industrial organization is, therefore, to direct whatever forces may be at his command toward the development, training, and guiding of all the brains in his organization. This is easier said than clone, for it involves the carrying on toward two objectives, the maximum of productive work, and the development of the individuals who perform it. This is essentially like the rebuilding and enlarging of a plant, keeping it in full operation the while. On the other hand, if the aim is merely for maximum production from an organization, without offering opportunity for development, the final outcome will necessarily be like that f a plant that is driven at full capacity without stopping for repairs or improvements; it eventually becomes so run down that either a complete rebuilding is necessary or else it has to be abandoned.

Jan 2, 1918

By A. Aly Selim, Bruce D. Hanson

A computer model capable of simulating the basic production system in a variety of underground mining methods has been developed at the U.S. Bureau of Mines Twin Cities Mining Research Center. The simulation of the basic production system at the working face, which includes drilling, blasting, mucking, and roof bolting, is described in this paper. The model is event-oriented; that is, the status of the system is updated at the occurrence of an event. All event or activity times, including failure and repair of equipment, are treated as stochastic variables. A location dimension, in addition to the time dimension, is used, allowing the simulation of systems with a number of stopes or faces. To aid the user in determining the efficiency of his system, the model provides several reports. An event or activity report gives the frequency of occurrence, total active time, and cumulative delay statistics for each event. An equipment report delineates the availability and utilization of all equipment. A production report provides information on various performance measures such as tons produced and holes drilled. The cost report gives the cost of various elements in the system. Operating data from an actual mining situation were used in validating the model. In addition to the description of the model, examples of simulation experiments which were performed on a room-and-pillar mine system are given. Use of the model to determine methods for improving the efficiency of the system is discussed.

Jan 1, 1976

By C. R. McEwen

This paper presents a technique for calculating the original amount of hydrocarbon in place in a petroleum reservoir, and for determining the constants characterizing the aquifer performance, based on pressure-production data. A method for doing this based on a least-squares line-fitting computationwus proposed by vanEverdingen,Timmerman andMcMahon in 1953. We found that their method would not work when there is error in the reservoir pressure data — even moderate error. The technique presented here appears to give reasonable answers when pressure data are uncertain to the degree expected in reservoir pressure determinations. The major change introduced in the present analysis is to limit the least-squares line-fitting to yield only one constant— the amount of hydrocarbon in place. The water-influx constant is then taken as proportional to the oil (or gas) in place. The constant of proportionality can be computed from estimates of effective compressibility and reservoir water saturation. It is also pointed out that the commonly used least-squares analysis assumes all of the uncertainty to be in the dependent variable. The material balance should be arranged so that this condition is fulfilled if correct inferences are to be made from statistical calculations. Examples are shown of the application of the new technique to gas reservoirs—both hypothetical and real—and to the oil-reservoir example of van Ever-dingen, Timmerman and McMahon. INTRODUCTION The amount of hydrocarbon originally in place in a petroleum resenroir can be estimated by means of the material-balance calculation. Simultaneous observations of pressure and amounts of produced fluids are required, together with the PVT data applicable to the reservoir fluids. If water encroachment is occurring, it is desirable to try to infer the behavior of the aquifer, as well as the original hydrocarbon in place, from the pressure-production data. This imposes additional demands on the method of calculation, and uncertainty in the data can result in large uncertainty in the answer. In addition, if the size of a gas cap is to be established, the whole problem becomes indeterminate, as pointed out by Woods and Muskat.1 Brownscombe and Collins2 simulated a gas reservoir and its aquifer on a reservoir analyzer and derived quantitative information on the effect of uncertainty in pressure and aquifer permeability on computed gas in place. Among the various techniques of estimating the performance of an aquifer, the method of van Ever-dingen and Hurst? based on compressible flow theory, seems to have been the most generally successful (see Ref. 4, for example). In this paper we shall confine ourselves to their representation of the aquifer. In 1953,' van Everdingen, Timmerman and McMahon5 introduced a statistical technique for deriving the amount of oil originally in place and the parameters which describe the aquifer. (We shall refer to this technique as the "VTM method'', as suggested by Mueller.6) Their example reservoir had no gas cap. It has been our experience that the VTM method gives a reasonable answer when the data are very accurate, but that inaccuracy (particularly in pressure) can cause the method to break down. The effect was first observed in gas reservoirs, but has since been seen in oil reservoirs also. In this paper we present another statistical method which has been successful in achieving a reasonable answer where the VTM method has failed. In the new method, one less parameter is derived from material-balance computations. It is assumed that values can be established for effective compressibility in the aquifer and reservoir water saturation independently of the material-balance calculation. The water-influx constant can then be obtained from these data and the quantity of hydrocarbon in place.

By Howard J. Smith

For several years this Institute has recorded in its Transactions the various discoveries of potash‡ in America, and the successive stages in the development of an independent domestic potash industry. This paper is intended primarily to continue the record through the year 1935. In retrospect it is not inappropriate to note that public interest aroused in 1910 resulted in the production of some 1090 tons of pure potash in 1915, as a small beginning in the domestic inorganic potash industry, and that the infant industry thus begun reaches its majority this year, following a record production in 1935 of about 600,000 tons of potash salts and exports to no less than 22 foreign countries, including six countries in Europe. This production, with an estimated 200,000 tons of pure potash, may be compared with a production in 1934 of 430,157 tons of salts with an equivalent of 144,342 tons of pure potash. Mine and Plant Improvement During 1935, the American Potash & Chemical Co., at Searles Lake, Calif., began to reap the benefits of its previously enlarged plant facilities, which increased its productive capacity to 180,000 tons of muriate annually. The United States Potash Co., at Carlsbad, N. Mex., improved its operations by installing trolley locomotives and other haulage facilities and is maintaining its development headings two years in advance of operations, with projected plans for orderly mining for years to come. The hoisting shaft at its mine is now equipped to hoist at the rate of 1,000,000 tons annually, and the refinery to produce 150,000 tons of muriate containing less than one per cent of impurities. The Potash Co. of America, at Carlsbad, N. Mex., completed the conversion of its plant from a prospecting to an operating basis late in 1935. The coarse-crushing plant, with a capacity of 150 tons an hour, formerly on the surface, was installed underground, with ore pockets

Jan 1, 1938

By Howard J. Smith

For several years this Institute has recorded in its Transactions the various discoveries of potash‡ in America, and the successive stages in the development of an independent domestic potash industry. This paper is intended primarily to continue the record through the year 1935. In retrospect it is not inappropriate to note that public interest aroused in 1910 resulted in the production of some 1090 tons of pure potash in 1915, as a small beginning in the domestic inorganic potash industry, and that the infant industry thus begun reaches its majority this year, following a record production in 1935 of about 600,000 tons of potash salts and exports to no less than 22 foreign countries, including six countries in Europe. This production, with an estimated 200,000 tons of pure potash, may be compared with a production in 1934 of 430,157 tons of salts with an equivalent of 144,342 tons of pure potash. Mine and Plant Improvement During 1935, the American Potash & Chemical Co., at Searles Lake, Calif., began to reap the benefits of its previously enlarged plant facilities, which increased its productive capacity to 180,000 tons of muriate annually. The United States Potash Co., at Carlsbad, N. Mex., improved its operations by installing trolley locomotives and other haulage facilities and is maintaining its development headings two years in advance of operations, with projected plans for orderly mining for years to come. The hoisting shaft at its mine is now equipped to hoist at the rate of 1,000,000 tons annually, and the refinery to produce 150,000 tons of muriate containing less than one per cent of impurities. The Potash Co. of America, at Carlsbad, N. Mex., completed the conversion of its plant from a prospecting to an operating basis late in 1935. The coarse-crushing plant, with a capacity of 150 tons an hour, formerly on the surface, was installed underground, with ore pockets

Jan 1, 1938

By S. M. Wu, R. L. Stone, T. D. Tiemann

The dissolution of silica from quartz in sodium hydroxide solutions was investigated by statistical methods. A second-order predicting equation was developed with 28 experiments expressing the dissolution of silica as a function of the variables quartz size, time, temperature and sodium hydroxide concentration. With this equation a comprehensive picture of the dissolution reaction was gained which included the effect on extraction of each variable and their interactions, fractional extractions, reaction rates and activation energies. A previous investigation of the extraction of silica from low grade siliceous iron ore by digestion in sodium hydroxide solutions presented the extraction as a function of time, temperature and sodium hydroxide concentration.' The present paper is an outgrowth of the previous work and develops by a mathematical analysis the functional relationship between the extraction of silica and time, temperature, sodium hydroxide concentration and particle size. Relatively pure quartz was used in this study in order that the basic quartz-sodium hydroxide reaction could be studied without the influence of other constituents that are commonly found in iron ores. The functional relationship between the dissolution

Jan 1, 1965

By F. A. Loving

This paper describes experiments which represent an effort to measure preciscly the pressure-time history in detonating condensed commercial explosives. The performance of explosives in mining has been a favorite subject for experimentalists and theorists for many decades. Much useful empirical information has been developed; but it has not yet been possible to predict, or even to account for, the behavior of explosives in mining on the basis of fundamental physical principals. A basic piece of information that would appear necessary for any rigorous analysis is the pressure-time history in the detonating condensed explosive. The experiments described here represent an effort to measure this relationship precisely in commercial explosives. THE EXPERIMENTAL METHOD The experimental method, first described by Cowan and Jacksonl has been extensively evaluated by Cowan and Balchan for relatively simple explosives. The experimental procedure involves the detonation of a layer of explosive between metal plates whose external surfaces are polished to a mirror finish (Fig. 1). Precise measurements of the deflection of

Jan 1, 1967

Jan 1, 1896

By James M. Barker

Jan 1, 1985

By Maxwell Gensamer, Vsevolod N. Krivobok

It is generally recognized that "low temperatures" have a definite effect on the properties of metals and alloys. Furthermore, it is also generally understood that the effect is decidedly adverse. The possibility that the same low-temperature treatment might be employed for the purpose of improving certain properties of the alloys has received less attention until recently, and the experimental and theoretical knowledge consequently is meager. The effect of low temperatures may be studied in two ways: by studying such properties as tensile strength, hardpess, impact resistance, etc., at low temperature, and by studying the properties of alloys at room temperature after these alloys have been subjected to a low-temperature treatment. Even in the latter case, certain properties will differ materially from the same properties obtained on identical material but without cold treatment. In other words, subjecting alloys to cold treatment results in an irreversible change in mechanical properties. Changes of this character are usually indications of the appearance of a new phase or phases, and would be likely to occur in alloys which, either by virtue of their composition or by heat treatment, are retained at room temperature in a metastable condition. Consequently, studies of the effect of low temperatures logically should be conducted from the physicochemical point of view, since the results of such studies not only would furnish metallurgists with interesting theoretical data but would clearly indicate the classes of alloys that would hold promise of improvement, and, possibly, of new applications. As the present investmigation deals with chromium and chromium-nickel alloys of iron of the stainless variety, it has been thought desirable to present the data concerning low-temperature properties of these alloys, available from the most recent and most comprehensive work of H. S. Russell.l A summary of these data is given in Table 1. The data are confined to a limited number of alloys of the type dealt with in the present

Jan 1, 1931

By V. V. Valleroy, A. J. Cornelius, B. T. Willman, G. W. Runberg, L. W. Powers

This paper reports the results of an investigation into the use of steam as a recovery agent. High oil recoveries by steam, as much as 100 per cent greater than by water flood, were demonstrated in laboratory linear cores. The principal mechanisms responsible for the additional oil have been identified as (1) thermal expansion of the oil, (2) viscosity reduction and (3) steam distillation. Although the first two also enhance recovery by hot water flood, they are effective with steam injection at lower volumes of produced fluids. Steam-zone residual oil saturation is independent of initial saturation. Recovery by steam is greater for lighter oils because they contain a greater fraction of steam-distillable components. Equations for the steam-zone advance in a radial homogeneous system are presented; these show the effect of reservoir parameters. The displacement in the hot liquid region moving ahead of the steam front can be estimated by a calculative procedure which is suggested. The method which takes account of the temperature gradient is applicable also to a regullar hot water flood. INTRODUCTION Considerable industry interest in the application of heat to petroleum reservoirs as a means of increasing oil recovery has been evident in the past few years. Three thermal methods which have attracted interest are (1) the injection of saturated steam, (2) the injection of hot water and (3) the injection of air into a formation to support combustion within the reservoir. Much has been published in recent years on underground burning, and there have been several papers on hot water injection. Although there are several early reports1,2 on steam injection, no recent experimental studies have been published on this recovery method. In our early thermal methods research, a marked improvement was found in oil recovery by steam injection compared with recovery by water injection in the same system. An investigation subsequently was undertaken to obtain quantitative information from which to evaluate steam injection as an oil recovery method. In addition, it was desired to investigate the recovery mechanism of hot water as compared with that of "cold" water and steam injection. An understanding of the mechanisms by which these processes operate would generally define favorable field conditions for their successful application. This paper reports the results of the laboratory study and suggests calculative approaches to estimate the field behavior of the steam-injection process. EXPERIMENTAL METHOD EQUIPMENT A schematic diagram of the equipment used for the displacement tests is shown in Fig. 1. Essentially, it consisted of a core with suitable means for recording temperatures, a feed system for injecting steam or cold or hot water, and a receiving train for condensing, cooling and measuring produced fluids. For the steam experiments, air-free saturated steam was piped from the boiler to the core face in insulated lines that were provided with steam traps and bleeds to remove any condensate formed in the lines due to heat losses. A 12-in. section of the line immediately ahead of

By Philip L. Collyer

A technique for evaluating the similarity of ore deposits is presented. A Q-mode cluster analysis is performed using information on tungsten mines in North America. This is a method of grouping several samples according to their similarity based on a simple matching coefficient. Tungsten mines are grouped or clustered on the basis of geologic variables. It is possible to interpret and evaluate different ore environments according to age, mineralogy, and other variables. It is well known that variable geologic conditions influence the occurrence of certain types of minerals. If one variable is responsible and identified, that variable can be searched for in exploration. If several interrelated variables are responsible for forming a favorable environment, then recognition and identification may be more difficult or completely masked. Therefore it is helpful in exporation to isolate and identify multiple interrelated variables which control the localization of mineral deposits. Cluster analysis is useful in mineral exploration where considerable amounts of data are available for analysis. Subtle relationships that might be overlooked because the data are voluminous can be isolated using this technique. Cluster analysis can be used to establish exploration prioriy in unknown areas on the basis of similarities to known metallogenetic provinces.

Jan 1, 1976

By F. J. Fayers, J. W. Sheldon

A theoretical discussion is presented of the behavior of the interface between two fluids of different physical properties when displacement is occurring along a thin tilted bed. An approximate equation of motion is deduced for the motion of the interface, suitable boundary conditions are derived and the conditions necessary for the approximation to be valid are determined. Steady-state behavior and the circumstances under which a steady-state interface can develop are discussed. Finite-difference solutions of the equation of motion for two characteristic water-drive problems are presenled. The theory of displacement in a thin bed may be regarded as being the porous-medium-equivalent of the classical "shallow-water wave" theory of hydrodynamics. INTRODUCTION The solution of moving interface problems are important in predicting the behavior of oil reservoirs under various conditions of production. Kidder4,5 has obtained analytic solutions for the fingering problem and for another particular case of the motion of the interface between two immiscible fluids of unequal density in a porous medium. In the present analysis, we analyze the problem of the motion of an interface along a slightly dipping bed whose length is long compared with its thickness. The requirement of a thin bed is necessary in order to give a simple differential equation for the motion. The two fluids are of unequal density and mobility. In the derivation of the equation of motion, it is assumed that the lighter fluid is being injected at constant rate at the high end of the porous layer, while heavy fluid is produced at the low end of the layer. The analysis applies equally well, however, to the reverse procedure in which the heavier fluid is injected at the low end of the tilted layer. It will be shown that the behavior of the system is critically dependent on which fluid has the higher mobility and that the interface can in certain circumstances achieve a simple progressive steady-state shape. The fluids may be considered to be immiscible or miscible. It is assumed, however, that the thickness of the interface region is negligible in comparison with other dimensions of the system, and that the strength of the saturation or concentration discontinuity remains constant during the motion. In conjunction with Buckley-Leverett theory and areal sweep-out calculations, the tilted interface problem is one of the fundamental problems of reservoir engineering; it defines the very important underpassing and overpassing phenomena associated with water-drives, gas-drives and secondary recovery operations. The tilted interface problem has been dealt with previously by Dietz.2 In the present paper, we endeavor to derive the pertinent mathematical results more rigorously, and in doing so we obtain modifications in both the differential equation and boundary conditions for those given by Dietz. We also attempt to show the connection between these results and some of the recent work on the stability of interfaces. Finally, we attempt to demonstrate the practical importance of the thin bed theory by showing the difference in behavior of two characteristic waterflood problems. In the favorable mobility case recovery is improved by increasing production rate, while in the unfavorable case undercutting or overcutting is minimized by reducing the displacement rate. PROPERTIES OF THE EQUATION OF MOTION Fig. 1 illustrates a porous layer of constant thickness H and tilt 8 with Fluid 1 of mobility ?l and density pl below Fluid 2 of mobility ?2 and density p2. The fluids are incompressible and flowing with volumetric flux q along the layer. Both the slope a of the interface with respect to the bed and the curvature of the interface are assumed to be small. What is meant by "small" is indicated in the derivation given in Appendix A of an approximate equation of motion for the interface. Other symbols are defined in the "Nomen-

By A. A. Berk

LAGGING technology and the slow spread of in- formation have been the chief obstacles to wide- spread participation in minimizing the industrial pollution load. These obstacles can be conquered by fact finding and applied research. There are three aspects to the problem of waterborne mineral waste: to minimize the amount being thrown away, to select a predisposal treatment, and to use the products or dispose of them.

Jan 7, 1957

By R. N. Hoffman

A new technique for determining capillary pressure curves has been developed and tested. The technique differs from previously reported centrifuge techniques in that the centrifuge is slowly accelerated from zero to the maximum desired speed rather than being held constant at particular, progressively higher speeds. An important advantage of this technique over other methods of determining capillary pressure curves is the short time required to obtain the desired amount of data over the chosen pressure range. For example, to obtain 30 data points between 1.2 and 104 psig with a 1.55-in. long, 3/4-in. diameter core using a brine-air system, 6.6 hours were required with this technique. An equally important development of this paper is an analytic method for the conversion of the data from the centrifuge experiment to capillary pressure curve data. Previously there has been only an approximate conversion available. Although the capillary pressure curves determined by this technique appear to be as accurate as those determined by other techniques, the accuracy could be improved if certain variables, not treated in this experiment, were investigated. Among these are the dynamic distortion of the centrifuge equipment and imperfect initial saturation of the cores. INTRODUCTION Pirson 1 defines capillary pressure in a porous medium as "the differential pressure that exists between two fluid phases at their interfaces when they are distributed under static equilibrium within a porous material". Capillary pressure in rocks is known to be a function of fluid saturation, among other things, and a capillary pressure curve is defined for the purposes of this paper as a plot of the capillary pressure-wetting phase saturation relationship for a particular rock sample. Several methods are used for determining capillary pressure curves for small rock cores. Prominent among these are the semi-permeable barrier, mercury injection, and, to a lesser extent, centrifuge methods. The semi-permeable barrier method is currently the most popular. It features simplicity in both execution and the mathematical conversion of the experimental data into a capillary pressure curve. The main disadvantages of the semi-permeable barrier method are the time required — as long as two months — to obtain several points for the curve, and a fairly low maximum pressure before breakthrough of the non-wetting phase into the barrier occurs, for example, about 32 psig for a brine-air system. It is for these reasons that other methods such as the centrifuge method have been introduced. High accelerations and the absence of a barrier result in quicker attainment of saturation equilibrium at a given pressure. However, the centrifuge method involves much more expensive equipment and more difficult procedures and calculations than does the barrier method. The purpose of this investigation has been to improve the equipment and procedures of the centrifuge method and to develop an analytic method for the conversion of the experimental data into a capillary pressure curve. Hassler and Brunner did the original work in the determination of capillary pressure using a centrifuge. In their work the centrifuge speed was increased in a step-wise manner, each speed being held constant until saturation equilibrium was reached in each core. Saturation equilibrium was indicated when the volume of liquid collected in the graduated pipette of the core holder remained constant. According to Hassler and Brunner, equilibrium was reached in "a few minutes to one-half hour or more". In the centrifuge method, as opposed to the barrier method, the fluid saturation of the core is not a constant throughout the length of the core, but varies with the radius of centrifugation. Also, the capillary pressure cannot be read directly but must be calculated from a knowledge of the centrifuge speed and other parameters. For these calculations, Hassler and Brunner derived the following equations: