Search Documents

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.

By Richard J. Mainiero, Harry C. Verakis

The Bureau of Mines has conducted research to determine whether the total elapsed delay time for blasting bituminous coal in underground mines could be safely expanded beyond the present 500-ms limitation without igniting a methane or methane-coal dust atmosphere. The results indicated that the increase of total delay from 500 to 1,000 ms had no detectable effect on safety relative to incendivity as long as permissible practices were observed in all other aspects. Research was also conducted to evaluate the safety of 18-in hole spacing for delay blasting in coal relative to misfires. For a variety of permissible explosives, misfires were observed for about 50 pct of the holes at 18-in spacing. Based on this observation, the Bureau has recommended that the new regulations prohibit hole spacings of less than 24 inches in underground coal mines. Misfires were also observed for 24-in spacing, but the number of misfires and test shots was too small to sup- port any firm conclusions.

By Michael Feroz

Severe drought in the Colorado River Basin over the past 10 years has caused water levels to drop in Lake Mead by more than 113 feet. The current lake level is at elevation 1,112 feet mean sea level (msl). Lake Mead is the primary source of water supply for greater Las Vegas. Water for Las Vegas from Lake Mead is accessed through two existing submerged intake shaft and tunnel systems. If the water levels continue to drop in the lake, the Intake No. 1 will become inoperable below water level elevation of 1,050 feet msl. To sustain the current levels of water delivery to all member agen-cies in the southern Nevada region, the Southern Nevada Water Authority (SNWA) has decided to construct a third intake at a deeper depth in the lake to ensure better quality water and enhanced reliability.

The Annual General Meeting of the Branch was held in the Witbank Golf Club on Thursday, 28th August, 1969 at 5.00 p.m. Mr C. G. Sowry (Chairman) was in the Chair. There were also present: Ten Members: Messrs J. Meintjes (Vice-Chairman), J. E. van Aswegen, A. C. Hofmeyr, I. Botha, B. F. v. d. Merwe, C. H. Coxon, R. Canny, R. C. J. Goode, A. Bain, G. E. Carlyon. Six Associate Members: Messrs B. T. Hosking, R. E. Powell, W. Weber, T. Bulloch, E. J. Jackson, P. J. Muller. Fourteen Visitors. Total Present: 31. The Chairman declared the meeting open and extended a welcome to all visitors including Mr R. C. J. Goode, President of the Institute. MINUTES OF PREVIOUS GENERAL MEETING The Minutes of the General Meeting held on the 26th June, 1969, were read and adopted. CHAIRMAN'S REPORT The Chairman presented a report outlining the activities of the Branch during the past year. Membership of the local Branch had increased from 80 to 99 during the year. DECLARATIONOF OFFICEBEARERS The Chairman announced that owing to his transfer out of the Witbank area, he was unable to accept the nomination as Chairman for a further year of office. The outgoing committee had elected Mr J. Meintjes Chairman and Mr A. Bain Vice-Chairman for the ensuing year. Mr J. Meintjes then took the chair. Mr J. Meintjes thanked Mr C. G. Sowry for his successful steering of the Branch during its period of inauguration. ELECTION OF COMMITTEE MEMBERS The Chairman then read the Scrutineers' Report and declared the following eight members elected as Committee Members for the 1969-1970 year: Messrs I. Botha, A. W. Brereton, G. E. Carlyon, B. T. Hosking, E. J. Jackson, B. D. McNeice, P. J. Muller, B. F. van der Merwe. TALK The Chairman then called upon Mr A. Bain to present his talk and film show on 'Impressions of the mining and metallurgy congress tour to Australasia'. CLOSURE After thanking the speaker for an interesting talk, the Chairman declared the Meeting closed at 8.40 p.m.

By Peter Swanson, John Rettkowski

As part of a research program to reduce injuries and fatalities in the U.S. mining workforce, NIOSH is investigating the potential use of noncontact vibration measurements to identify hazardous ground fall conditions. In the present study, the effectiveness of remote vibration sources in exciting resonant responses in damaged structures is examined. Initial experiments were performed on a damaged concrete slab in a controlled environment. Vibration response was measured with a laser Doppler vibrometer using both single-point and two-dimensional array recordings. Excitation sources included direct mechanical impacts to remote portions of the slab and a remotely located industrial demolition hammer. Surficial damage was effectively delineated with both sources through visualization of vibration response in spectral bands selected to discriminate wave propagation from large-scale slab resonance effects.

By Mahmood Alemrajabi

When LKAB processes their iron ores in northern Sweden, they generate apatite as a by-product. The apatite has been recovered by flotation and the concentrate contains about 0.4 wt% rare earth elements (REE). Both phosphate rock and REE are listed as critical raw materials for the European union. In the present study, the recovery of REE within the nitrophosphate process of fertilizer production has been studied. In the first stage, the REE are precipitated and isolated as a phosphate concentrate. For further purification and individual separation of REE by solvent extraction, removal of phosphorous from the concentrate to enhance the dissolution of REE in mild acidic solutions is essential. In this study, thermal treatment of the concentrate with sodium hydroxide at elevated temperatures followed by water leaching has been studied in order to separate the phosphorous and thereby facilitate dissolution of the REE. The results show that 95% of total phosphorous initially present in the REE phosphate concentrate is removed after thermal treatment with NaOH 1:1 mass ratio, at 400 °C.

By H. L. Stone, A. O. Garder

A numerical method of solving the partial differential equations which describe the one-dimensional displacement of oil by gas has been presented. Possible extension of the method to treat multidimensional flow is discussed, and the limitations of this extension are indicated. Using this method, it is possible to allow for the existence of a gas cap, the presence of any number of gas-injection or oil-production wells and the evolution of dissolved gas from the oil. It is also possible to allow for variation in the cross-sectional area, elevation, porosity and permeability of the reservoir. The influence of relative permeability and the force of gravity in the direction of flow upon the displacement is considered. The influence of capillary pressure upon the flow and the effect of gravity in the direction perpendicular to flow are neglected. The physical properties 01 the fluids are considered to be Junctions of pressure only, and equilibrium between contiguous phases is assumed. The numerical calculations can be readily carried out by the use of a digital computer. Several example analyses have been performed using the IBM 704 computer, and about one-third of an hour of computing time was required per case. Reservoir behavior predicted by use of this numerical method was compared to data obtained by other methods for three cases — complete pressure maintenance, dissolved-gas drive and gas-cap drive. The independent solutions to these problems were obtained by analytical solution, laboratory experiment and field data, respectively. Agreement of the numerical solution with data from these sources was good; this agreement establishes the convergence and accuracy of the numerical method. INTRODUCTION Most petroleum reservoirs can be produced by any one of several alternative programs. When a reservoir is produced by primary methods, production economics can be influenced by controlling the number and location of wells and the flow rate of each well. An even greater influence may be achieved by augmenting the recovery of oil obtainable by primary methods. This can be accomplished by injection of fluids such as water, natural or enriched gas or a bank of light liquid hydrocarbons. Selection of the most desirable operation requires a means of predicting the reservoir behavior which will result from each of the several alternative programs. The purpose of this paper is to present a mathematical method for predicting the behavior of reservoirs produced by gas-cap drive, dissolved-gas drive or pressure maintenance by gas injection. The method described herein takes cognizance of phase changes caused by a decline or an increase (due to gas injection) in reservoir pressure, of the presence of a gas cap and of the effect of gravity on the flow of gas and oil. Relative permeability relationships are used to define the flow properties of the rock. Allowance is made for variation in cross-sectional area, elevation, permeability and porosity of the reservoir. Both the influence of capillary pressure upon the flow and pressure gradients in the gas cap are neglected. Whenever a liquid phase and a gas phase are in contact, they are assumed to be in equilibrium. The physical properties of the fluids are considered to be functions of pressure only. Therefore, if the method is to be used to predict the effects of a gas-injection program, mixtures of the injected gas and formation crude should-have the same physical properties as mixtures of formation gas and crude. The equations to be presented in this paper apply only to a one-dimensional case; therefore, they neglect the influence of gravity in the direction transverse to the flow. As is well known, this gravitational influence may lead to overriding of oil by gas. Consequently, this procedure as presented is most applicable to long, thin reservoirs for which gravity overriding is not important. On the other hand, the equations presented can be generalized to treat multidimensional flow and, hence, to consider gravity overriding, if desired. A word of caution on two points is advisable here, however. First, the authors have not demonstrated the accuracy of the numerical technique for multidimensional flow. Second, and more important, capillary pressure will often be of importance in multidimensional problems. Obviously, in such cases a generalization of the

By Fred C. Turin, Kim M. Cornelius

Work in underground coal mines is repetitive and physically demanding. Workers commonly report a wide range of aches and pains. Management at one mine was concerned about increased reporting of aches and pains by roof bolting machine operators. An analysis of roof bolting tasks was conducted at this mine to identify early warning signs of cumulative trauma. Health and safety risk factors were identified and recommendations to reduce cumulative trauma exposure were provided. Recommendations addressed job procedures, equipment design, and worker awareness of risk factors.