Search Documents

By R. J. Blackwell, J. R. Rayne, W. M. Terry

This paper presents results of an experimental investigation of factors that control the efficiency with which oil is displaced from porous media by a miscible fluid. The study was made to elucidate the relevant processes both on microscopic level (within individual or between neighboring pore spaces) and on macroscopic level (within a large sand body). Mixing of miscible fluids on the microscopic level was studied in sand-packed tubes. It was found that molecular diffusion is the dominant dispersion mechanism for reservoir conditions of rate, length and pore sizes. Macroscopic channeling was studied for various mobility ratios in reservoir mode1s scaled to relate viscous, gravitational, and diffusional forces. The formation of channels was due to viscous firzgering, gravity segregation and variations in permeability. With adverse mobility ratios, it was found for reservoirs of realistic widths that diffusion will not be effective in preventing the formation and growth of fingers, even in homogeneous sands. At sufficiently low rates channeling was eliminated by gravity segregation in tilted reservoirs. The dependence of recovery on mobility ratio, length-to-width ratio, flow rate and angle of dip is presented. INTRODUCTION Oil recovery by solvent flooding is finding increasing application in the field. While the process promises high recoveries from the region swept by solvent, under adverse conditions only a small fraction of the reservoir volume may be swept at the time solvent breaks through to the producing well. Further, the high cost of the solvent encourages its use only as a bank whose size must be kept at a minimum. Thus, two important questions arise: (1) what fraction of the reservoir can be swept, practically, by solvent? and (2) what is the minimum size solvent bank that can be used to carry out the displacement? The answers to these questions require knowledge of both macroscopic channeling processes and microscopic mixing processes. The studies described here were carried out to gain this knowledge. Microscopic mechanisms which cause mixing will be discussed first, because an understanding of these mechanisms is necessary for proper interpretation of the experimental work on channeling described later. INVESTIGATION OF MICROSCOPIC DISPLACEMENT PHENOMENA FOR MISCIBLE FLUIDS The questions to be resolved concerning microscopic displacement behavior are: (1) how much mixing occurs between oil and solvent in the direction of flow? and (2) does the solvent completely flush the oil from the pore spaces in the region invaded by solvent? The classic work of Sir Geoffrey Taylor' and its extension by Aris2 have shown that it is possible from theory to describe the amount of mixing in single straight capillaries when solvent displaces a fluid of equal viscosity and equal density. Aris showed that the lengths of the mixing zones, 81, corresponding to the 10 per cent and 90 per cent concentration levels of the solvent can be calculated from the formula, 81 = 3.62 \/Kt, .........(1) where 8l is the length of the mixing zone, cm; t is the time, sec; and K is the effective dispersion coefficient, cm2/sec, given by K/D = 1 + (1/48) (au/D)2......(2) where a is the radius of the capillary, cm; D is the molecular diffusivity, cm2/sec; u is the average velocity

By S. Fitzel

The Betts process for electrorefining of lead has been in operation since 1902 in the Trail Operations Metallurgical complex. The Betts process allows production of high purity lead in a single electrochemical step. Process improvements and plant expansions have taken place over the many decades of operation of this process which has proven to be robust, flexible, and capable of successfully treating a large range of anode bullion compositions. This paper reviews the history of the Betts lead electrorefining process, the innovations and quick industrial adoption of new discoveries that took place at the end of the nineteenth century. The process was first implemented in Trail Operations and integrated very well with metallurgical needs then and now.

By T. C. Poulter, B. M. Caldwell, D. E. Brink

The flash X-ray has been used more than a decade to study the configuration of the jet from a shaped charge. The high-speed, rotating-mirror smear camera has provided time-distance graphs of detonations and shock fronts in transparent materials, and the high-speed framing camera has given pictorial representations of the progress of explosive phenomena. Adequate means for measuring the paths and velocities of all parts of the shaped-charge cavity Iiner during the collapse phase have not existed heretofore. A technique enabling the single-lens framing camera to make stereoscopic photographs of the cavity liner while it is collapsing has been developed. Analysis of this photographic record gives the directions and velocities of various parts of the Iiner surface, permitting direct quantitative measurement where this has previously been impossible. Significant improvements in shnped-clrarge design are expected to result. INTRODUCTION The familiar oilwell jet perforating charge, generally referred to as a "shaped charge", is related to the Mun-roe charge first described by C. E. Munroe in 1888 and later by Neumann in 1910. It differs in construction from the Munroe charge in that its cavity is lined with some inert material (usually metal)', and it differs in performance by projecting a fast jet of dense liner material against the target instead of a stream of expanding detonation products. Fig. 1 represents a typical shaped charge in axial cross section. The charge is circularly symmetrical about its longitudinal axis. The cavity at the right is lined with metal, usually copper; and, with the exception of air, the liner is otherwise completely empty. High-explosive material is intimately in contact with the convex surface of the liner and extends in the other direction to the booster. This is a pellet of high explosive having somewhat different characteristics and serves to couple the detonator to the main charge. When the detonator is fired, a small hemispherical reaction zone of very high temperature and pressure occurs at its output end. The booster responds, adding its own latent energy to the reaction, thereby initiating the main charge. The reaction progresses throughout the main charge at a velocity of 7 to 8 mm/microsecond until all the explosive material has reacted. With a reaction pressure of about 160,000 atm, the surrounding material is moved rapidly away from the high-pressure zone. As the detonation proceeds along the cavity liner, moving from apex to edge, the liner is collapsed toward the axis of the cavity, where liner material from all parts of the concave surface converges. A zone of extraordinarily high pressure results. with only one escape path—outward along the axis. The particles from the liner, thus, become a jet of very high velocity, with a density characteristic of the material from which the liner was made. The kinetic energy of the jet, concentrated in a small area of the target, has great penetrating power. It overcomes the strength of any target, and the target material crumbles or deforms plastically to escape from the path of the jet. An unlined cavity charge, on the other hand, projects a turbulent body of hot gaseous material against a target and erodes the exposed surface. With equal quantities of explosive, the penetration effectiveness of an unlined charge is a small fraction of that of a lined charge. Given an explosive charge with a detonator at one end and an axially symmetrical metal-lined cavity at the other, one has a perforator of considerable potency. To maximize the penetrating and other desirable effects,

By O. C. Holbrook, George G. Bernard

A new theoretical treatment has been obtained for the behavior of pattern waterflood injection wells when closed in. Two cases are treated: Case I where oil and water are assumed to have the same properties, and Case 2 where they arc different. In applying the method, one plots log (p — p,) vs closed-in time, where p is well-bore pressure at any tims and p, is static pressure. The value of p. is determined by trial and error as that value which makes the plot linear at large time. A value for the permeability-thickness product can be determined from the intercept of this linear part, and a value of the skin factor from the injection pressure at time of closing in. Application of the method to data from water floods at three fields seems to give reasonable results. For the case of unit mobility ratio, it is proved that this new method should give the same value for permeability-thickness product as the conventional pressure build-up method. In addition, the new method gives correct values for static pressure, whereas the conventional method does not, often indicating negative static pressures. The new method may be used in cases where the surface pressure persists after closing in as well as in cases where it does not. INTRODUCTION It is of considerable interest and importance to be able to determine the characteristics of the reservoir in an area surrounding a water injection well. Thus, if we can determine early in the life of an injection well that there is a considerable "skin effect", remedial measures can be started before a full-scale pattern flood begins. Similarly, if it can be shown that a gradual buildup of skin effect is occurring with time, measures to free the water of plugging material can be taken. Determination of static pressure in the water-injection well may show that the water is entering a thief zone and not the desired reservoir. Finally, determination of the permeability of the sand around the injection well will allow estimation of the future relation between injection pressure and rate. It should be possible to determine average reservoir permeability, skin effect and static pressure from pressure fall-off data. However, at the time we began work on this subject, it was thought that no adequate theory on which to base such determinations' was available. According to the conventional method which considers the reservoir to be filled with one fluid of small compressibility (see Van Everdingen, Joers2, and Nowak2), shut-in pressure is plotted vs log where is injection time, and At is closed-in time. The physical significance of injection time, may well be questioned in this case, since in a reservoir completely filled with a single fluid (as required by this theory) and with input and output rates equal, the pressure behavior after an initial transient is independent of t,. Attempts by our Tulsa area to use this theory led to negative values of static pressure in most cases. Because of these limitations of the method discussed above, it was decided to attempt to develop a new theory of pressure behavior in water injection wells, one which would apply when there is a gas saturation, as is so often the case in water floods. In the following treatment the assumptions and basic equations are given first, then the method of application of the equations. A complete example is given to clarify details of application. All difficult mathematics has been placed in the appendices so that the reader can follow the text without difficulty. However, if he wishes only to apply the results without knowing the basis for them, he can learn how to do this from reading only the sections entitled "Plotting of Experimental Results" and "Example." ASSUMPTIONS AND BASIC EQUATIONS Statement of Problem It will be assumed that a horizontal layer of constant thickness contains in its pore system a mixture of oil, gas and water. While water is being injected into this pore- system through a well at constant rate, an oil bank is built up, gas being expelled from the space taken by the oil as shown in Fig. 1. The saturations within each

By L. D. Mullins, E. B. Elfrink, W. L. Wahl

In the past few years several articles and papers presenting results of solution gas-drive depletion calculations have appeared in the lit-erature. Such calculations are of interest to the oil industry, for investment decisions must often be made before much is known about a reservoir. At other times, an estimate of the possible benefits to be realized from alternate production methods is desirable, and theoretical depletion calculations can serve as a floor or reference level from which to work. In any case, an estimate of ultimate oil recovery based upon engineering data is commonly required. An engineer confronted with the problem of obtaining, for a specific reservoir system, an estimate of ultimate oil recovery by solution gas-drive depletion usually will be forced to perform the calculations himself. This is despite the quantity of data in the literature. Rarely will either experience or the literature provide results from a reservoir system similar in all important respects to the one under consideration, and calculated results are not so plentiful that satisfactory interpolation procedures can be devised. Performing the calculations, however, is a tedious, time-consuming task unless an electronic computer is available, and, in practice, time and manpower are not always available for this purpose. A quick, simple, consistent method was needed for reducing the uncer- tainty in estimated oil recovery from solution gas-drive reservoirs when only minimum information about the reservoir system is available. PROCEDURE Method of Calculation The usual requisite assumptions were made so that the material balance equation could be used to calculate data for the charts. The following assumptions were made: (1) the reservoir is homogeneous and isotropic; (2) oil recovery is due entirely to solution gas drive and neither a gas cap nor a water drive nor gravity drainage is present; (3) the initial reservoir pressure is the bubble-point pressure of the reservoir fluid; (4) initial total liquid saturation is 100 per cent of pore space; (5) interstitial water saturation remains at the initial value as the reservoir pressure declines from the bubble-point pressure to atmospheric pressure; (6) equilibrium gas saturation is 5 per cent of pore space; and (7) oil and gas saturations are uniformly distributed throughout the reservoir at all times. There are no saturation gradients due to a wellbore, nor is the geometry of the reservoir system considered. The material balance equation was written in the form of a differential equation' which was integrated to determine the change in oil saturation for an assigned pressure drop. Formal integration was not possible, so recourse was made to the Runge-Kutta method' of numerical integration. All computations were performed on IBM equipment. Numerical integration yielded the change in oil saturation within the reservoir as the pressure. declined from the bubble-point pressure to atmospheric pressure. The initial oil saturation minus the change in oil saturation yielded the oil saturation at atmospheric pressure. The oil originally in place was obtained by dividing the initial oil saturation by the initial formation volume factor (differential liberation) while the oil in place at atmospheric pressure was obtained by dividing the final oil saturation by the formation volume factor at atmospheric pressure. U1timate oil recovery, expressed as a percentage of the initial oil in place, was obtained by dividing the difference between the oil initially in place and the oil in place at atmospheric pressure by the oil initially in place and multiplying by 100. PVT Data Charts were based upon 135 solutions to the material balance equation. PVT properties of the reservoir fluids were variables in this equation and had to be known as functions of pressure. The required PVT data might have been obtained from either actual reservoir fluid systems or correlated data. However, correlated PVT data were developed and employed in the calculations for the following reasons: (1) it is doubtful that 135 sets of PVT data could have been obtained for the values of variables investigated in this study, and (2) results of recovery calculations from randomly obtained PVT data could not be correlated as well as rssults from selected PVT data. The PVT data used to develop the

By J. S. Lasater

Resu1ts of experinmental measurernents of heat capacities and thermal conductivities of some typical porous rocks are presented. Measured heat capacities agree closely with va1ues calculated front known chenlical . compositions of the rocks. On the basis of this agreement, heat capacitir.7 of fluid-satttrated rocks were calc~ililted. Therrnal conductivities of the rock samples were measured under varioits conditions of fluid saturation. From these data therrnal dif-fvsivities were calculated. A significant variation of thermal diflusivity with trnlperature is indicated. Knowledge of the thermal characteristics of fluid-bearing porous rocks has become increasingly important with the development of thermal oil-recovery processes. Reliable thermal data for petroleum reservoir rocks are not available and approximate values are generally assumed for reservoir calculations. 'The effects of temperature, pressure and fluid saturation on thermal properties of porous rocks have not been carefully studied. The purpose of the present work was to measure and report the thermal characteristics of some typical porous rocks. A group of eight sedimentary rock samples, including limestone, shale, sandstone and silt-stone, was selected for the study. Properties reported include heat capacities, thermal conductivities and thermal diffusivities. The effect of fluid saturation on these properties was also considered. Determinations of thermal properties require precise measurements which are time consuming and difficult to reduce to a routine basis. Methods have not been developed for measurement of thermal proper-ties under conditions of pressure, temperature, and fluid saturation, which may be encountered in thermal recovery processes. As one approach to these problems, some preliminary work on methods of estimating thermal properties from other more readily measurable characteristics of the rock-fluid system is reported. The samples selected for the studies were oilfield cores from within, or closely associated with oil-productive zones. Descriptions of individual samples are given in Table 1. Clay minerals present in the samples were determined by differential thermal analysis (DTA). Quartz content was estimated from DTA cooling curves. For all but the thermal conductivity tests, a representative 200-gm sample of each core was disaggregated, thoroughly extracted in diox-ane and oven dried at 100°C. The cleaned and dried sample was care- fully split into appropriate portions for the several tests. Chemical analyses of the samples were obtained to permit calculations of heat capacities by methods which are outlined later. Resuits of chenlical analyses are given in Table 2. Heat contents of the test samples were measured by a Runsen-type calorimeter. The heat content of the unknown sample is measured relative to the known heat content of platinum. The apparatus is referred to as a relative error type because of the comparative measurements, and has the advantage of canceling systematic errors. Experimental results are reproducible within 0.5 per cent. A temperature base of 298.16 K is taken for heat content measurements. Measurements were made at five temperatures (260, 440, 620, 800 and 980° F). It was necessary to preheal the sample above the maximum experimental temperature to eliminate volatile constituents. Assuming the weight loss on pre-heat, ing was mostly water, a correction to the observed heat content was made to account for the water originally present in the sample. This latter procedure was followed in the determination of calculated heat contents as will be described later. The apparatus and the method of measure-

By Catherine Johnson, Barbara Rutter, Phillip Mulligan

Shock tubes, shock tunnels, and blast tubes have been widely used by explosive researchers and shock wave investigators to understand the characteristics of shock waves. Rectangular, cylindrical, and conical are the three simple geometries commonly used in shock tube and shock tunnel testing. These geometries result in shock waves that resemble the Friedlander curve’s simple exponential decay. However, geometries that are more complex have not been studied. This research investigates how changes in shock tunnel geometry affect the pressure versus time waveform. The Large Arena Test Simulator (LATS), which is composed of four different rectangular sections connected by three expansion sections, was used to determine how the changes in geometry affect the shock wave profile. Shock wave profiles were measured inside and outside the LATS to investigate the geometry’s effect. A representative shock wave profile from each geometry and open-air explosive testing were compared visually showing the geometries effect on the shock wave profile. The LATS produced a waveform with four pressure peaks behind the initial shock front, which were found to be associated with the four rectangular sections. The resulting pressure trace was compared to pressure traces from simulated urban bombings. The number of turns in an urban environment equated to a new peak, which is caused by additional expansion sections of the shock tunnel. From this observation, complex geometry shock tunnels can be purpose built to simulate an explosion in an urban setting. The use of shock tunnels with complex geometry in simulating urban bombings can allow for numerous representations and a better understanding of profiles produced during an urban bombing or similar scenario.

By Dennis R. Dolinar, John L. Ellenberger, Gabriel S. Esterhuizen, Anthony T. Iannacchione

A survey of pillar conditions was carried out at 21 operating limestone mines that use the room-and-pillar method. The surveyed pillars were all located in rock that was classified as “Good” to “Very Good” using the RMR rock mass classification system. The average pillar width was 14.5m and the width to height ratio varied between 0.53 and 3.52. Pillar conditions were rated using a visual rating system that accounts for the effect of geological structures as well as stress related fracturing. Nine cases of pillar instability were recorded. Pillar instabilities were typically limited to individual pillars or small groups of pillars within an otherwise stable layout. Local geological structures such as inclined discontinuities or weak bedding planes contributed to instability in seven of nine cases that were observed. Elevated pillar stresses contributed to instability in the two remaining cases. All the unstable pillars observed had width-to-height ratios of less than 1.5. It is concluded that pillar instability is most likely to be caused by unfavorable geological structures in pillars with width to height ratios of less than 1.5. Stress related instability such as rib spalling becomes more prevalent when the average pillar stress approaches approximately 20 MPa.

By Edmund C. Pechin

I REGRET that I am unable to present this subject in definite form and detail. All I shall attempt at this meeting is to lay before you some curious facts, the bearings and explanations of which must be reserved for future consideration. I hope before the next meeting to have carried experiments far enough to furnish a more important contribution to the science of iron metallurgy. The Wickersham process consists in the introduction into the molten metal, as it runs from the furnace, of air from the hot blast. Mr. Samuel Wickersham is the originator of the plan, as applied in practice, although Mr. Martien had the idea a long time ago, and, I believe, went to England about it. But at that time Mr. Bessemer was coming before the public with his invention, which engrossed the attention of iron-masters, and Mr. Martien accomplished nothing. The principal difficulty in our experiments at Dunbar was to get a practicable runner or trough, through which the molten metal should pass, and in which it should be exposed to the action of the Mast. The runner we used was ten feet long, and perforated on the sides to permit the entrance of the air-currents. When we used a cast-iron and fire-clay runner, the heat was so great that before we were more than half done with a run, the lower part of the runner was melted, and we had a magnificent Fourth of July display of fireworks until the blast was shut off. We had to abandon fire clay. The runner with which we have had the best success is made of wrong- ht-iron tubes, 29 inches in diameter, attached to a wrought-iron frame, and lined with a mixture of crushed conglornerate (quartz) and a little fire clay to give it cohesion. The wrought-iron shell curves inward over the upper edges, so as to catch and throw back into the current the flying sparks and particles of metal thrown out by ebullition. With this we have been able to run 12 to 14 tons at a heat with great ease. Mr. Wickersham has since prepared this material (pounded quartz) into a shaped brick or tile, with the necessary perforations for the passage of the air, which promises to answer every requirement for steady working. As regards the improvement in the quality of the iron, there seems to be no question. A large number of practical tests were made in

By J. C. Melrose, W. R. Foster, J. G. Savins, E. R. Parish

Many differences can be imagined between gas-oil flow in which the gas is supplied at the face of the core and gas-oil flow in which the flowing gas was originally dissolved in the oil. If capillary pressure characteristics and flow requirements control gas saturation distribution, the gas would be expected to be located at preferred sites within the porous medium as determined by pore sizes. On the other hand, during solution gas drive the gas first appears as bubbles through a nuclea-tion process. Nothing in self-nucleation theory specifies at which sites the first bubbles should be formed. In all probability they will be randomly distributed throughout the porous medium. Furthermore, it is not at all certain that even at low rates of production the gas will redistribute itself after nucleation to the channels normally occupied by gas in simple gas flow. Stewart, et at, have shown that at least for some limestone samples, oil recoveries could not be predicted for all rates of production using any one set of relative gas and oil permeabilities. An important factor in controlling recoveries during solution gas drive was the rate of bubble formation, higher rates giving higher recoveries. Stewart, et al, attributed the increase in recovery to a better distribution of the gas phase in heterogeneous limestone samples than is obtained by simple external gas drive. Differences in recovery from these causes were not reported for sandstone cores. In the experiments to be reported here, oil recovery, pressure and producing GOR history were measured during solution gas drive for a 5-ft sandstone core. The results were compared with predictions from the Muskat method for computing solution gas-drive behavior using external gas-drive relative permeability. The effects of changing the rate of production and oil viscosity were studied. At high laboratory rates of average pressure decline, two observations were made which would not have been predicted by Muskat's depletion theory: (1) oil recovery increased with increasing rate of production for a given viscosity oil, and (2) oil recovery increased with increasing oil viscosity for a given high rate of production. Both of these observations are explained as consequences of diffusion control of gas saturations superimposed on the normal gas-oil flow requirements, Fur- thermore, discontinuous gas phase flow appears to be significant during solution gas drive. The laboratory tests were performed at rates of average pressure decline many times greater than the maximum possible rate of average pressure decline in an actual oil field. It is, therefore, not possible to draw any direct conclusions regarding the effect of rate on recovery for the solution gas-drive mechanism under actual field conditions. However, at the lower laboratory rates, recoveries were nearly independent of rate and could be predicted by the Muskat method, using external gas-drive relative permeability data. These results suggest that at normal oilfield rates the effect of rate on recovery for the solution gas-drive mechanism is negligible. EXPERIMENTAL PROCEDURES The core material for the pressure depletion studies was Bandera sandstone from an outcrop in the Mid-Continent. This sandstone was selected because of its low permeability (about 10 md), which would permit the development of substantial pressure gradients in the corn at moderate flow rates. The core was 5-ft long and 2-in. in diameter. Its properties are listed in Table 1. Relative gas-to-oil permeability ratios were measured by an external gas-drive method2. The results are shown for a short 2-in. core and for the 5-ft core in Fig. 8. Oils used in the pressure depletion experiments were kerosene and a highly refined white oil (standard white oil No. 3) with gas-free viscosities of 1.8 and 25 cp, respectively. The gas was a naturally occurring methane from Gough field, Inglewood, Calif. The oil viscosities, gas solubilities and formation volume factors are plotted as functions of pressure at 75°F in Figs. 1 and 2. Methane viscosities and compressibilities were obtained from the literature"'. A core mounting was required which could withstand up to 2,500 psi internal pressure. This was obtained by first encasing the core completely in a plastic resin (Scotch Cast, manufactured by Minnesota Mining & Manufacturing Co.) The plastic covered core was then inserted into a steel pipe equipped with screw caps so that the plastic coating could be pressured from

By D. W. Reed, R. L. Reed, Tracht

Laboratory experiments on the reverse combustion of tar sands in a linear adiabatic system have shown that a highly upgraded oil can be produced from an exceedingly viscous, immobile oil. The dependence on the air-injection rate of peak temperature, combustion-zone velocity, oil recovery, air-oil ratio, residual coke and oil, fuel burned and distribution of product gases is shown graphically. Eflects of initial temperature, oxygen concentration, oil saturation and heat loss are discussed. Experiments bearing on the coking properties of heavy oils are mentioned and some results exhibited. Field application of the process hinges on the existence of adequate air permeability and the rate of reaction under reservoir conditions. INTRODUCTION It has been established that oil can be recovered from underground reservoirs by means of at least two fundamentally distinct methods involving in situ combustion of a certain fraction of the oil. Characteristic of both of these known methods is the production of oil from one or more wells by means of hot gases formed when a high-temperature reaction zone is advanced through the reservoir. In both cases, the reaction zone is created by heating certain of the wells to a sufficiently high temperature prior to the introduction of air, and the zone is maintained and advanced through the reservoir by appropriate control of the air-injection rate. In the first of these methods, which is called "forward combustion",' the combustion zone advances in a direction which is generally the same as that of the air flow; whereas in the second method, "reverse combustion",' the combustion zone moves in a direction generally opposite to that of the air flow. Forward combustion, on the one hand, is an ideal combustion process in the sense that a minimum of the most undesirable fraction of the oil is consumed as fuel in the form of coke, a clean sand is left behind and generated heat is used as efficiently as possible. However, the applicability of forward combustion is limited. Since the products of combustion, vaporized oil and connate water must flow into relatively cold regions of the reservoir, there is an upper limit on the viscosity of oil which can be moved by this process in a practical and economical fashion.' On the other hand, it is characteristic of reverse combustion that the vaporized oil and water together with the products of combustion are produced through sand which is already hot and has had its mobile liquid content eliminated. This means there is no upper limit on oil viscosity; indeed, the oil may be an entirely immobile semi-solid. However, fuel for the process is an intermediate fraction of the original oil, and the most undesirable fraction remains on the sand surface as a substantial deposit of coke. Since this coked material is not burned during reverse combustion, it represents energy which is available for the production of oil but is not used for this purpose. It follows that one can expect economics to be somewhat less attractive with reverse combustion than with forward combustion. Nevertheless, it is a process which is designed for reservoirs where forward combustion is impossible and, as such, has become a subject of experimental and theoretical investigation. In this paper, only experiments made with tar sands are discussed. DESCRIPTION OF THE PROCESS We proceed, then, to consider the process of reverse combustion in greater detail. Fig. 1 illustrates a temperature profile defining a combustion zone which moves from right to left when air flows from left to right. In Zone 1, the temperature is the initial reservoir temperature, and the tar sand is as yet unaltered. This statement must be modified to the extent that physical properties of the oil may be changed by low-temperature oxidation at reservoir temperature. As air passes into Zone 2, which has been warmed by conduction, it assists in vaporization of the very light ends (if there are any), and oxidation occurs at a significant rate. In this region, there is almost no production of carbon monoxide or carbon dioxide because predominantly addition-type reactions take place with the formation of oxygenated compounds such as aldehydes and acids together with water. The hydrocarbon-enriched and slightly oxygen-depleted gas stream enters Zone 3 where

The following notes have been compiled to assist authors in the preparation of papers for presentation to the Institute and for publication in the Journal. All papers must meet the standards set by the Council of the Institute, and for this purpose all papers are referred to at least two referees appointed by Council. STANDARDS FOR ACCEPTANCE To merit consideration papers should conform to the high standards which have been established for publication over many years. Papers on research should contain matter that is new, interpretations that are novel or of new significance and conclusions that cast a fresh light on old ideas. Descriptive papers should not be a repetition of well-known practices or ideas but should incorporate developments which would be of real interest to technical men and of benefit to the mining and metallurgical Industry. In some cases a well prepared review paper can be of value, and will be considered for publication. All papers and particularly research papers no matter how technical the subject, should be written with the average reader of the Journal in mind, to ensure wide interest. The amount of textbook material included in a contribution should be the minimum essential to the argument. The length of a paper is not the criterion of its worth and it should be as brief and concise as possible, consistent with the lucid presentation of the subject. Only in very exceptional circumstances should a paper exceed 15 pages of the Journal (15 000 words, if there are no tables or diagrams). Six to ten pages is more normal. NOTE: Papers in the Journal are printed in 10 point type, which is larger than the 8 point type used on this page. For special publications Council may decide on page sizes smaller than A4 used for this Journal. The text should be typewritten, double-spaced, on one side only on A4 size paper, leaving a left-hand margin of 4 cm, and should be submitted in duplicate to facilitate the work of the referees and editors. LAYOUT AND STYLE Orthodox sequence Title and author's name, with author's degrees, titles, position. Synopsis, including a brief statement of conclusions. Index, only if paper is long and involved. Introduction. Development of the main substance. Conclusions, in more detail. Acknowledgements. References. Title: This should be as brief as possible, yet give a good idea of the subject and character of the paper. Style: Writing should conform to certain prescribed standards. The Institute is guided in its requirements by: Collins, F. H., Authors & Printers' Dictionary-Oxford University Press. Hart, H. Rules for Compositors and Readers. Humphrey Milford (familiarly known as the Oxford Rules). Fowler, H. W. & F. G. The King's English-Oxford University Press. General: A few well selected diagrams and illustrations are often more pertinent than an amorphous mass of text. Over-statement and dogmatism are jarring and have no place in technical writing. Avoid the use of the first person, be objective and do not include irrelevant or extraneous matter. Avoid unnecessary use of capitals and hyphens, while punctuation should be used sparingly and be governed by the needs of sense and diction. Sentences should be short, uninvolved and unambiguous. Paragraphs should also be short and serve to separate basic ideas into compact groups. Quotation marks should be of the 'single' type for quotations and "double" for quoted matter within quotations. Interpretations in the text should be marked off by parentheses ( ), whereas brackets [ ] are employed to enclose explanatory matter in the text. Words to be printed in italics should be underlined singly. For small capitals they are to be underlined DOUBLY and for large capitals TREBLY. If there is any problem in producing formulae accurately by typewriter they should be inserted by handwriting in ink in the copy forwarded by authors. Abbreviations and symbols are laid down in British Standard 1991. Abbreviations are the same for the singular and plural, e.g. cm for centimetre and centimetres, kg for kilogram and kilograms. Percentages are written in the text as per cent; the symbol %is restricted to tables. A full stop after an abbreviation is only used if there is likely to be confusion of meaning. Metric System: The Systeme International d'Unites (S.I.) is to be used for expressing quantities. This is a coherent system of metric units derived from six basic units (metre, kilogramme, second, ampere, kelvin, and candela), from which are derived all other units, e.g. the unit of force is the newton (N) for kilogramms metre per second per second (kg m/s'). Always use the standard metric abbreviations. Commas must not be used for separating groups of digits. For ease of reading digits should be grouped in threes counting from the decimal point towards the left and the right. Illustrations: Drawings and diagrams are to be in black India ink and should be about 18 cm wide. When submitting graphical representations avoid a fine grid if possible. Curves should be in heavy line to stand out. Lettering too should be bold as a reduction in size is often involved in the printing process. (A single column is 8.5 cm wide.) Numbering of tables should be in Roman numerals: I, 11, etc. and figures in Arabic numerals: Fig. I, Fig. 2, etc. (Always use the abbreviation for figure.) Photographs should be black and white glossy prints. As a guide to the printer the author should indicate by means of notes in the typescript where tables and figures, etc. are to appear in the text. Paragraphs: A decimal system of numbering paragraphs may be used when the paper is long and complicated and there is a need for frequent reference to other parts of the paper. Proof correction: Galley proofs are sent to authors for the correction of printers' errors and not for the purpose of making alterations and additions which may be expensive. Should an author make alterations which are considered excessive, he may be required to pay for them. Standardsymbols as laid down in British Standard 1219C shoulrl be used. SYNOPSIS It is most important that the synopsis should provide a clear outline of the contents of the paper, the results obtained and the author's conclusions. It should be written concisely and in normal rather than abbreviated English and should not exceed 250 words. While the emphasis is on brevity this should not be laboured to the extent of leaving out important matter or impairing intelligibility. Summaries simplify the task of abstractors and therefore should present a balanced and complete picture. It is preferable to use standard rather than proprietary terms. FOOTNOTES AND REFERENCES Footnotes should be used only when they are indispensable. In the typescript they should appear immediately below the line to which they refer and not at the foot of the page. References should be indicated by super-script, thus. . .1 . . .2. Do not use the word Bibliography. When authors cite publications of other societies or technical and trade journals, titles should be abbreviated in accordance with the standards adopted by this Journal. GENERAL The Council will consider the publication of technical notes taking up to three pages (maximum 3 000 words). Written contributions are invited to the discussion of all papers published in the Journal. The editors, however, are empowered by the Council to edit all contributions. Once a paper or a note has been submitted to the Institute, that document becomes the property of the Institute, which then holds the copyright when it is published. The Institute as a body is, however, not responsible for the statements made or opinions expressed in any of its publications. Reproduction from the Journal is permitted provided there is full acknowledgement of the source. These points should be borne in mind by authors who may submit their work to other organizations as well as to the Institute.

By Özgen C. Karacan

It is generally accepted that typical coalbed gases (methane and carbon dioxide) are sorbed (both adsorbed and absorbed) in the coal matrix causing it to swell and resulting in local stress and strain variations in a coalbed confined under overburden pressure. The swelling, interactions of gases within the coal matrix and the resultant changes in the permeability, sorption, gas flow mechanics in the reservoir, and stress state of the coal can impact a number of reservoir-related factors. These include effective production of coalbed methane, degasification of future mining areas by drilling horizontal and vertical degasification wells, injection of CO2 as an enhanced coalbed methane recovery technique, and concurrent CO2 sequestration. Such information can also provide an understanding of the mechanisms behind gas outbursts in underground coal mines. The spatiotemporal volumetric strains in a consolidated Pittsburgh seam coal sample were evaluated while both confining pressure and carbon dioxide (C02) pore pressure were increased to keep a constant positive effective stress on the sample. The changes intemal to the sample were evaluated by maps of density and atomic number determined by dual-energy X-ray computed tomography (X-ray CT). Early-time images, as soon as CO2 was introduced, were also used to calculate the macro porosity in the coal sample. Scanning electron microscopy (SEM) and photographic images of the polished section of the coal sample at X-ray CT image location were used to identify the microlitho types and microstructures. The CO2 sorption-associated swelling and volumetric strains in consolidated coal under constant effective stress are heterogeneous processes depending on the litho types present. In the time scale of the experiment, vitrite showed the highest degree of swelling due to dissolution of CO2, while the clay (kaolinite) and inertite region was compressed in response. The volumetric strains associated with swelling and compression were between ± 15% depending on the location. Although the effective stress on the sample was constant, it varied within the sample as a result of the intemal stresses created by gas sorption-related structural changes. SEM images and porosity calculations revealed that the kaolinite and inertite bearing layer was highly porous, which enabled the fastest CO2 uptake and the highest degree of compression. Keywords: Coal; Carbon dioxide; Sorption; Swelling; Sequestration; X-ray computed tomography; Volumetric strain; Gas outbursts

By R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

By O. Ryosuke Suzuki

A new metallurgical process via sulfide is proposed: a sulfide is produced at a high temperature from its metallic oxide using gaseous CS2, and this sulfide is electrochemically reduced to its metallic state using molten salt. This combined process via sulfide is effective to obtain high purity of metallic powder, even if the metal in its oxide is strongly combined with oxygen. For example, it is not easy to reduce stable oxide TiO2, and only Ca can remove oxygen to form α-Ti. However, a fairly large amount of oxygen remains as Ti–O solid solution. Because the solubility of S in Ti is very small, this proposal was examined experimentally both on the conversion of TiO2 to TiS2 and on the successive reduction of TiS2 to Ti. TiO2 powder was exposed to CS2 gas flow at 1073 K, and the conversion to TiS2 was confirmed. TiS2 could be reduced to Ti powder either by calciothermic reduction or electrolysis in a CaCl2 melt. By Ca reduction at 1133 K in CaCl2 melt, sulfur concentration decreased to 0.03 mass%S when the amount greater than twice the stoichiometric calcium amount is added. By electrochemical reduction at 1173 K in CaCl2–CaS melt, S concentration significantly decreased to 0.01 mass% S when four times larger amount of electric charge was supplied.

By K. H. Coats

This paper presents the development and solution of a mathematical model for aquifer water movement about bottom-water-drive reservoirs. Pressure gradients in the vertical direction due to router flow are taken into account. A vertical permeability equal to a fraction of the horizontal permeability is also included in the model. The solution is given in the form of a dimensionless pressure-drop quantity tabulated as a function of dimensionless time. This quantity can be used in given equations to compute reservoir pressure from a known water-influx rate, to predict water-influx rate (or cumulative amount) from a reservoir-pressure schedule or to predict gas reservoir pressure and pore-volume performance from a given gas-in-place schedule. The model is applied in example problems to gas-storage reservoirs, and the difference between reservoir performances predicted by the thick sand model of this paper and the horizontal, radial-flow model is shown to be appreciable. INTRODUCTION The calculation of aquifer water movement into or out of oil and gas reservoirs situated on aquifers is important in pressure maintenance studies, material-balance and well-flooding calculations. In gas storage operations, a knowledge of the water movemenr is especially important in predicting pressure and pore-volume behavior. Throughout this paper the term "pore volume" denotes volume occupied by the reservoir fluid, while the term "flow model" refers to the idealized or mathematical representation of water flow in the reservoir-aquifer system. The prediction of water movement requires selection of a flow model for the reservoir-aquifer system. A physically reasonable flow model treated in detail to date is the radial-flow model considered by van Everdingen and Hurst.1 In many cases the reservoir is situated on top of the aquifer with a continuous horizontal interface between reservoir fluid and aquifer water and with a significant depth of aquifer underlying the reservoir. In these cases, bottom-water drive will occur, and a three-dimensional model accounting for the pressure gradient and water flow in the vertical direction should be employed. This paper treats such a model in detail — from the description of the model through formulation of the governing partial differential equation to solution of the equation and preparation of tables giving dimensionless pressure drop as a function of dimensionless time. The model rigorously accounts for the practical case of a vertical permeability equal to some fraction of the horizontal permeability. The pressure-drop values can be used in given equations to predict reservoir pressure from a known water-influx rate or to predict water-influx rate (or cumulative amount) when the reservoir pressure is known. The inclusion of gravity in this analysis is actually trivial since gravity has virtually no efFect on the flow of a homogeneous, slightly compressible fluid in a fixed-boundary system subject to the boundary conditions imposed in this study. Thus, if the acceleration of gravity is set equal to zero in the following equations, the final result is unchanged. The pressure distribution is altered by inclusion of gravity in the analysis, but only by the time-constant hydrostatic head. The equations developed are applied in an example case study to predict the pressure and pore-volume behavior of a gas storage reservoir. The prediction of reservoir performance based on the bottom-water-drive model is shown to differ significantly from that based on van Everdingen and Hurst's horizontal-flow model. DESCRIPTION OF FLOW MODEL The edge-water-drive flow model treated by van Everdingen and Hurst1 is shown in Fig. la. The aquifer thickness b is small in relation to reservoir radius rb. water invades or recedes from the field at the latter's edges, and only horizontal radial flow is considered as shown in Fig. 1b. The bottom-water-drive reservoir-aquifer system treated herein is sketched in Fig. 2a and 2b. Here the aquifer

By E. J. Boome, A. C. Schmidt

Pert, and how Pert techniques can be used in modern mine management By E. J. Boome, A. C. Schmidt C.A. (S.A.) This paper covers the stages of planning required to be performed and the implementation of the plan from grass roots to the partial completion of the sinking of a mine shaft system. The various steps involved in the planning are also examined in some detail. The experience gained in this project and the conclusions reached convince the authors that Pert Planning is a valuable aid to top management.