Search Documents

By R. C. Neavel

One of the most important requirements to be met by any coal producer is delivery of size consists tailored to specific applications. Generally these specifications include a minimum as well as a maximum size limit. The customers of a major Indiana coal producer noticed that some shipments of sized coal contained a large proportion of undersize pieces. Also, in these shipments, the pieces which remained within the specified limits were cracked and, when handled, crumbled readily. The producer, certain that these shipments had been properly sized at the tipple, suspected that degradation had occurred in transit. The Indiana Geological Survey was requested to examine the coal from one of these shipments to attempt to determine the cause of breakage. Two basic facts about coal and coal seams led to the employment of petrographic techniques in this investigation. First, all coal seams consist of accumulations of diverse organic and inorganic constituents; and second, most coal breakage, either desirable or undesirable, can be attributed to the inherent weakness of some of these constituents or to poor adhesion between certain constituents. Analyses of degraded coal Sampling: Samples were secured from 2 stock piles of 2-by 4-inch coal in a dealer's yard. The coal in both piles had received the same treatment from mining to storage, although it had been mined and shipped on two different days. The size consist of one pile, however, appeared smaller than that of the other. The sample from the degraded pile is referred to as the slack sample, and the other is called the unslack sample. Size analyses: In order to determine the actual size consist, the samples were screened, with a minimum of rough handling, on square 2-inch, 1-inch, 1/2-inch, and 12-mesh screens. It is assumed that 100 percent of the 2-by 4-inch coal is retained on a 2-inch round screen at the tipple. Analyses of several loads at the tipple indicate that approximately 26 percent of the 2- by 4-inch coal, 100 percent of which was assumed held on a 2-inch round screen, passed through a 2-inch square screen. Twenty-six percent may be taken as the normal difference between round and square screens. Results of the analyses of stockpile coal are shown in Table 1. Cumulative percent versus size are plotted on a Rosin-Rammler chart as Figure 1. Both samples approximately follow normal size distribution; this indicates that breakage is not size preferential and also that the samples are representative of the stockpiles. Although both samples exhibit the same general size distribution, the median size of the slack sample is considerably smaller than that of the unslack sample. Fifty-five percent of the unslack sample was undersize (29 percent greater than the normal 26 percent difference), and 83 percent of the slack sample was undersize (57 percent greater than the normal 26 percent difference). Both piles contained a greater percentage than normal of undersize coal; this indicates that they were both affected by degradation, but to a different degree. Bedding-plane analyses: Coal, because of its laminated nature, will break preferentially along bedding planes, especially where adjacent constituents are poorly cemented to each other. Subsequent breakage perpendicular to the

Jan 1, 1961

By J. W. MacDonald

Review of equipment, performance, and shaft design plans including head-frames for sinking and forms for lining are given. Also included are comparative data for shafts of 11 x20 ft and 193/4 x 30 ft with 756 and 666-ft depths, respectively. Advantage from curtain wall installation with "slip forms" which warrant consideration for normal lining are outlined. A preliminary review is also made of "chemical grouting" to stop water inflow leakage with Halliburton equipment and service. Two coal mine shafts recently completed in Franklin County, Ill., afford experience and accomplishment to warrant review. These shafts provide access to a tract with a potential recovery of one hundred million tons. This work provided the initial sinking by other than a contractor's organization since 1918. Four additional shafts required during the interim were sunk under contract arrangement with others. Supervisory contact and inspection afforded experience to govern planning and conduct of the recent work. TEST DRILLING Core test holes were drilled within each of the two shaft locations. Four additional cable tool holes were drilled on a radius of 700 ft from the shaft bottom to determine the coal seam elevations and govern the underground layout plans for each shaft. A graphic log of the hoist shaft test drillhole is shown by Fig. 1. A very soft sandstone was found at a depth of 16 ft. Cored recovery of the sandstone proved much weaker than the overlying compacted sand secured with coring facilities for soil recovery. Gradation from the very soft sandstone to fairly hard material occurred with the next 14 ft, extending to a depth of 30 ft. The weak character within the upper portion entailed need of adequate protection when sinking as described later. A part of the sandstone provided sufficient porosity and permeability to warrant further investigation. Adjacent cable tool drilling for temporary water supply disclosed inflow was limited to 2 gpm with a depth of 40 ft. It was increased to 12 gpm during the next 50 ft of drilling and occasioned subsequent precaution when sinking. General arrangement of facilities for sinking the hoist shaft is shown by Fig. 2. An available steel frame building was used to accommodate the hoist, compressor, shop, and shower or change room. The stiff leg derrick with 90-ft length of boom permitted removal of forms from the wood shop to the storage area or delivery within the shaft. Reinforcing steel and other needs were similarly handled, in addition to headframe erection. The headframe was placed over one end of the shaft to enable passage of forms for lining through the adjacent half next to the stiff leg derrick. The steel headframe as shown by Fig. 3 was fabricated locally. Side anchored cables from the

Jan 1, 1961

By H. F. Silliman, Daniel R. Hull, Earl W. Palmer

THE brass-rolling industry has not had a great deal of experience with antimony in its product. There have been some recent excursions with antimony as a corrosion inhibitor in tubes, but in sheet brass it has been regarded with apprehension and generally avoided. This is not to say that [ ] contaminations have never occurred, and probably a great many brass mill men have had their own ideas on the effect of such contamination. Our own views were based upon certain experiences some years ago with strip castings that split lengthwise in the early stages of rolling. The nature of this splitting is shown in Fig. I, from a photograph taken in 1922 as a record of 20 bars of 65-35 common brass made with an antimonial brand of copper. These bars were cast 6 ¾ in. wide by 1 in. thick by 45 in. long. In the first pass through the rolls, 19 of the 20 split lengthwise throughout their length, while each half curved outward in a wide arc as it emerged from the rolls. [ ] In those days, if splitting was not a frequent occurrence, it was not unfamiliar. One can recall periods of short duration when large amounts of common brass would behave in this way. There was nothing imaginary about these spells. The affected metal would split lengthwise in the breaking rolls, and no maneuvering of mill practice made any difference. The condition would descend like a plague, demoralize the mill for several days, and disappear. It is doubtful if anyone knew the cause-however many explanations may have been current. In 1921 a quantity of split brass was traced to the use of fired small-arms

Jan 1, 1943

By R. W. Bartlett

DISILICIDE coatings, MeSi2, on refractory metals are usually grown by a solid-state diffusion reaction similar to the parabolic oxidation of metals. Two or more silicide compounds occur in each of the binary silicon-refractory metal systems. When an adequate source of silicon is available during the coating process, the lower silicide compounds are usually not formed.1 However, during service application of the coated refractory metal, continued heating causes further diffusion and growth of the lower silicides. A study of the growth kinetics of the lower silicides in the molybdenum and tungsten binary systems with silicon was previously reported.2 The present study is concerned with the growth kinetics of Ta5Si3 and Cb5Si3, formed as conjugate layers between TaSi2/Ta and CbSi2/Cb diffusion couples, respectively. Diffusion couples consisting of disilicide-coated tantalum and columbium sheet coupons were subjected to diffusion heat treatments, cross-sectioned, and scanned by an electron-microprobe analyzer to determine the relative X-ray intensity profiles. These profiles correspond closely to concentration profiles for the constituents, and were used to identify phases and locate phase boundaries. The analysis of a single specimen after a diffusion-anneal provides only one value for the Me5Si3 phase thickness, corresponding to the time and temperature of the diffusion treatment. Consequently, each isothermal curve showing the rate of accumulation or depletion of silicide phases required several samples with identical diffusion temperatures but different diffusion times. The disilicide coatings were prepared by pyrolizing commercially pure silane in a hydrogen carrier gas. The hydrogen was purified by passage through a palladium membrane. Commercially pure tantalum and columbium sheet metal coupons, 0.005 by $ by 2 in., were used. All of the coated specimens for each system were prepared together in a single group in order to minimize coating variations among individual specimens. During pyrolysis, the brightness temperature was maintained at 1200°C. With suitable flow rates at this temperature, only the disilicides of tantalum and columbium were formed in detectable amounts. Coated specimens were subjected to diffusion an-neals in a tube furnace with automatic proportional control within ±10°C. A H2-Ar mixture, deoxidized and dried over anhydrous magnesium perchlorate, was used to prevent oxidation during the diffusion anneal. Diffusion couples were cross-sectioned, polished and examined with an A.R.L. electron-microprobe X-ray analyzer. Both metal constituents in each specimen were simultaneously scanned using a 1-p beam spot size and either 1- or 2-µ steps. Analyses were performed using the Ka radiation of columbium and silicon, and the La, radiation of tantalum. For all specimens, an acceleration voltage of 30 kv was used, and integration was made against the specimen current. LiF and A DP monochrometers were used to isolate characteristic radiations. Three or more traverses were made on each specimen. The intensity scales were adjusted using the constituent metals and standard samples of the silicide compounds which had been verified by X-ray diffraction. The relative intensity profiles were used directly, in conjunction with the standard samples, to identify compounds and phase boundaries. The stoichiometries differ considerably between the silicide compounds in both systems, and the homogeneity range for these silicide compounds is narrow. These characteristics result in marked differences in relative intensities of the fluorescing elements—differences which permit ready identification of phases and phase boundaries. Projections of both electron back-scattered images and characteristic

Jan 1, 1967

By J. R. Kyte, L. A. Rapoport

The material balance equation for partiai or complete 1:trter-drivc reser1,oir.s been re-arranged to include a pressure irrferferencr ternr. This pressure interfernce term was ohtained from the theory established hy Morinda. A .successful progranl deissed to solve the resrclritlg c,qrration in a few hours on a nleelitr .sizeti, digital compuier. No prior knowdelge of the properties of the eqrrifer is required in the solution. 7' ri~rl calculations areatle using various assumed con!binuiions of value of the constant expretsing the qquifer properties correspoinds to the best fit ohtained between observed and calculated field history. Details of the method of calculaion are illus trated by an example. INTKCDUCTICN A method 01' calculating the material balance equation for a partial water-drive reservoir completed in an aquifer has been described in the literature." Pressure interference i'rc,m other pools iocated in the same nquiucr is not taken into account in this solution. When present, the pres- sure interference effect must be included in order to make a complete performance analysis. A correction for this effect has been developed using the theory cstahlished by Mor-tada.' Manual calculation of the material balance for water-drive reservoirs is a formidable task. Such manual ca1culation when interference eflects are includcd would be impossible to completc in any reasonable time period. However, a medium-sized stored program, digltal conlputer can he used satisfactorily. The method is described in the following sections. THECKETICAL CONSIDEKATIOhS The material balance equation for a reservoir producing by partial or total water drive, as proposed by Van Everdingen, This relationship is linear for Z/a from which thc slopc, Z, and the intercept, N, can be obtained by the least-squares method from the data points. It is not necessary lo know the aquifer properties to uktain oplimum values tor 2nd IV. Various assumptions are made tor At,,: where 11 is production time, seconds (oilfield between the equal successive time intervals into which the history is divided. t,, h It, times the number of equal time periods from i.e.. the total dimen-sionless time interval from It is noted that the dinlensionless time interval, At,,, involves all the relevant aquifer properties. Accordingly, any assunlption of is an assumption of the combined aquifer properties which pertain to the solution of the equation. Deviations of the data points from a straight line arc calculated for each assumption of and the value of A, correspond ing to the minimunl deviation is chosen as the optimum solution. This method is described in Ref. 1 and 2. Thc material balance solution is usually made with an initial assumption that the reservoir being studied is located in an infinite aquifer from which there are no other withdrawals. Deviation of the points by amounts greater than can be attributed to inaccuracies in field data may be the result of: (1) boundary effects of a limited aquifer, (2) pressure interference from othcr pools completed in the .sarne aquifer, and (3) a combination of these effects. Differentiation of these effects can be made only on the basis of geologic data. Pressure interference from other fields will result in a pressure drop at the original oil-water boundary having an interference component, The measured instantaneous pressure drops at the oil-water

By H. T. Kennedy, S. M. Avasthi

An equation developed for gaseous hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 0.73 percent when applied to 264 natural gas and condensate systems including 2,043 PVT points. Another equation developed for liquid hydrocarbon mixtures predicts molal volumes with an average absolute deviation of 1.12 percent when applied to 346 crude oil systems including 1,759 PVT points. Both equations require composition of the mixture to be expressed as mole fraction of methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide, together with the characteristics of the heptanes-plus fraction in addition to the temperature and pressure. The equations cover wide ranges of the variables involved, and their accuracy is considerably better than that of other available methods. The equations were differentiated to allow calculation of the coefficients of isothermal compressibility and isobaric thermal expansion. (In this paper the coefficient of isothermal com-pressibility and the coefficient of isobaric thermal expansion will be expressed as compressibility and thermal expansion coefficient, respectively.) Equations to calculate these quantities are presented. INTRODUCTION Calculations of reservoir performance for petroleum reservoirs require accurate knowledge of the volumetric behavior of hydrocarbon mixtures, both liquid and gaseous. Compressibilities are required in transient fluid flow problems, and thermal expansion coefficients are important in thermal methods of production. An accurate laboratory investigation of the PVT behavior of each reservoir fluid encountered would be costly and time consuming. For this reason various correlations for predicting fluid properties have been developed and recorded in recent literature. Correlations have been presented in the form of graphs, tables and equations. Since an increasing number of studies are being conducted with the aid of electronic computers, recent efforts have been directed toward development of correla-tions suitable for computer programming. Application of computers permits the use of more complex correlations which otherwise are not feasible. Moreover, methods for predicting reservoir performance, particularly those based on the compositional material balance, depend upon the capability of accurately expressing the molal volumes and other fluid properties as functions of pressure, temperature and composition. The coefficient of isothermal compressibility c is defined by

Jan 1, 1969

By Jerome F. Kuzmick

Condusions. From the work on aluminum powder, the following general con clusions can bc drawn: I. The physical properties of aluminum compacts, such as density. strength and ductility, increase with increasing pressure, up to the practicable pressure limits. 2. Good properties ill aluminum con pacts can be secured only by the use of high pressures, because it is neccssary to rupture the refractory oxide film on aluminum powder. 3. Considerable softening of the specimens occurs upon sintering. regardlers of the pressure used. General Conclusions 1. For relatively short-time sintering. molding pressure is a very important factor in the final physical properties of the sintered specimens. 2. The effect of pressure. is difierent for each metal, so that few general rules can be drawn that would apply to all metals. 3. In many cases, fur a particular metal powder and sintering temperature there is a critical pressure at which the best physical properties can be obtained. 4. The principles of cold work and in nealing and their effect on grain structure applied to massive metal also holds true for sintered metals. 5. Particle size, at least for electrolytic iron powder, does not have a. great effect on ultimate grain six. DISCUSSION [E. E, Schumacher presiding V.. V. Crane. *—I was particularly impressed with Dr. Kuzmick's work on tin IF not the illustration, particularly the Second one. a very good illustraition of hot preessing. in that

Jan 1, 1945

By G. W. Thomas

Jan 1, 1967

By P. G. Hallof

Only recently has the term Induced Polarization found its way into geophysical literature. Schlumberger first mentioned the possible use of the induced polarization effect in one of his early papers. Later, induced polarization was the subject of several university theses, such as Seigel's at the University of Toronto and Bliel's at the Michigan School of Technology. In addition, a group at the New Mexico School of Mining and Technobgy has done work on induced polarization. There is very little to be found in print on the subject of induced polarization. Several papers have been presented at meetings, the 1958 AIME and SEG* conventions, for instance, on laboratory measurements of the effect and theoretical considerations useful in the interpretation of field results. This paper is the first attempt to present actual field results which show the usefulness of induced polarization as a prospecting technique. There are several phenomena that can cause polarization when electrical current is passed through the ground. In a paper presented at the 1958 SEG* meetings Professor T. R. Madden of M. I. T. demonstrated that only two of these gave rise to appreciable effects. The first is a result of the membrane and diffusion effects in a sandstone that contains appre- ciable amounts of clay minerals. The unbonded charges associated with the clay mineral forma material that passes positive and negative ions with different mobilities. When current is passed through such a system for some period of time, polarization results. The second phenomenon, which gives rise to very much larger polarization effects, is the one used in the new prospecting technique. These effects are a result of the blocking action or polarization of metallic or electronic conductors in a medium of ionic solution conduction. This electrochemical polarization, called "the over-voltage effect" by electro-chemists, occurs whenever electrical current is passed through a volume of rock which contains metallic minerals such as base metal sulphides. This polarization at metallic interfaces in ionic solutions depends upon the chemical energies necessary to allow the ions to continue the flow of current by giving up or receiving an electron from the metallic surfaces. The energy stored in this reaction isanalogous to that stored when a condenser is charged. If a d. c. current flowing through such a system is interrupted, the polarization can be detected as a small, decaying current that flows after the applied current is discontinued. The decaying potential measured is a result of the current flow caused by charges returning to their equilibrium positions from the positions they assumed under the applied potential field. Because of the polarization, there is an excess of change at each metallic interface. There is an alternate method of measurement of the induced polarization effect. The electrical

Jan 1, 1961

By G. F. Butterworth

The metallographic structures most frequently encountered in rifle barrels, and which are illustrated by the accompanying photomicrographs, fall naturally into two groups, distinguished by the method used to, produce in the stock a physical condition having the requisite properties. The fist group consists of "rolled" barrels; that is, barrels subjected to hot working by rolling in or near the critical range. In the second group, the stock, which is smaller in diameter than the first, is upset to form the butt end and is then heat treated by giving it a quench and a draw. These barrels will be referred to as " heat-treated " barrels. The structure of the rolled barrels resembles closely that of the same steel after annealing. There is the same network of excess ferrite outlining the grain boundaries, but the grains themselves are composed of sorbite rather than pearlite. The grain size is closely related to the rolling temperature. The critical-point curves of this grade of steel, which is approximately 0.50 to 0.60 per cent. of carbon and 1.00 to 1.30 per cent, of manganese, show a single very pronounced point between 1300" and 1350" F. (704" and 732" C.). Barrels rolled within this temperature range give an exceedingly fine grain, shown in Fig. 1. In fact, the grains may be so fine that, at a low magnification, they may be confused with the heat-treated structure shown in Fig. 6. A magnification of 500 diameters, however, will always resolve a rolled barrel into the characteristic structure shown in Fig. 2. As the temperature is increased, the grains are found to be larger. If, on the other hand, the rolling temperature is below the critical range, the structure previous to rolling will not be obliterated, the only effect of rolling being to elongate the coarse sorbite grains in the direction of rolling, Fig. 3. This distortion is greatest at the muzzle. Heat-treated barrels are quenched in oil from above the critical range, which should give a martensitic structure, Fig. 4, but the presence of some troostite in the martensite is frequently noted when the quench has not been sufficiently drastic, Fig. 5. The structure brought about by the subsequent draw is not so striking, though quite as typical. It is sorbitic or sorbito-pearlitic, and under a low magnification appears almost homogeneous, as in the case of fine grains of pepper and salt well mixed together, Fig. 6. This structure is substantially the same for all

Jan 1, 1920

Jan 1, 1967

Jan 1, 1948

Jan 1, 1948

By F. Keller, E. H. Dix, L. A. Willey

The consideration of alloying elements for aluminum has led to a series of investigations of the equilibrium relations between aluminum and those alloying elements. Therefore, the aluminum end of the aluminum-antimony diagram has been studied. The general form of the aluminum-antimony equilibrium diagram has been worked out by Campbell and Mathewsi and by Tammann.2 Canlpbell and mathews described principally the beginning of solidification. They found the curve of beginning of solidification rising gradually from the melting point of pure alunlinum to a maximum at about 67 per cent. aluminum. After a slight lowering of the curve at a little greater antimony concentration the curve again rises to a maximum at 18.4 per cent. aluminum, which corresponds to the compound AISb. The curve then drops to the melting point of pure antimony. Tamrriarln explained the maximum occurring at about 67 per cent. aluminum as the result of the slow formation of the single compound AISb. The space lattice of the compound AISb, according to Owen and Preaton,³ is face-centered cubic, the antimony atoms being intermeshed with an identical aluminum lattice. Very little, if any, work has been done on aluminum alloys containing small amounts of antimony. The Present Investigation It has been the purpose of this investigation to determine the solid solubility of antimony in aluminum and also to investigate the eutectic temperature and concentration. The experimental work has been carried out in a manner similar to that described in a previous paper4 from this laboratory.

Jan 1, 1931

By P. K. Guerrero, N. Arbiter

The correlation of power consumed by a tumbling mill with the dimensions, speed, and load has been attempted by three principal methods. One of these, the torque formula, has been reviewed critically elsewhere.' This approach, while analytical, disregards the individual motion and speed of the tumbling bodies and assumes that the gross geometry of the entire load is sufficient to establish the power requirement. A more general approach, using dimensional analysis, has been applied in this laboratory by the present authors1'1, and more recently by Rose.3 This disregards the internal conditions within the mill entirely and in effect is an empirical procedure for correlating data. Finally, the complete analysis of tumbling mill dynamics as first made by Davis' attempts to develop the power requirement from individual ball and rod paths and velocities. Although Davis' result is applicable only at a single speed for a given size and load it could be useful, since his optimum speed is at the upper limits of conventional tumbling mill practice. Where power values predicted by the Davis equation were compared with actual mill requirements, they were found too low for small mills and much too high for the largest mills. This suggested that the agreement with actual intermediate range requirements might be fortuitous, and the development of the equation was examined critically. The probable source of error was found to lie in an erroneous assumption made by Davis in converting the energy of ball or rod motion to power. With the same approach, but with different assumptions, a modified equation has been obtained which predicts results that are usually consistently high. This is to be expected, since it is assumed that there is no slip between load and shell and no interference among tumbling bodies. Both these effects should account for the difference. Development of Equation: When a cylindrical shell charged with a substantial weight of rods or balls is slowly rotated about a horizontal axis, the load surface first assumes an inclination that gradually increases to a maximum. Continued rotation results in failure within the load, characterized by a sliding of the top layers down the slope and a reduction in the angle of inclination. Further rotation again increases the angle, and the cycle repeats. If the speed is increased, there is continuous load failure and build-up and a continuous stream of rods or balls descends over the upcoming mass, which moves in circular arcs. At higher speeds the up-

Jan 1, 1961

By G. Thodos, D. E. Matschke

Cansiderable time and effort frequently are expended to establish, with a degree of confidence, the PVT behavior of pure substances. In particular, a great deal of experimental information contributed by a number of investigators is available in the literature. The data of each investigator, although significant in themselves, cannot be used with a high degree of confidence until they are compared with the continuum of information presented by others. Current interest in thermodynamic and transport properties, for not only the gaseous state but the con-densed state as well, requires that PVT infarmation be presented for both regions. The density dependence on temperature and pressure, when presented as a continuum, is of considerable value in itself, especially since the available equations of state fail to properly account for the PVT behavior in the condensed state. For this reason and to incorporate PVT data recently presented in the literature, this investigation was carried out. In their studies concerned with the thermodynamic properties of methane, Matthews and Hurd in 1946 utilized the data of Kvalnes and Gaddy and Olds, Reamer, Sage and Lacey to establish the PVT behavim of methane in the superheated region and the saturated liquid densities of Keyes, Taylor and Smithl in conjunction with the Clapeyron equation to develop the densities for the saturated envelope. Current activity sponsored by API Project No. 44 under the direction of F. D. Rossini has led to a tabular presentation of compressibility factors for gaseous methane in the superheated region.la These values were obtained by using smoothed residuals, developed from experimental data and the Ben edict-Webb-Rubin equation,' to establish internally consistent compressibility factors for the real gas. This approach offers probably the most precise presentation of PVT information for gaseous methane to date. The PVT behavior of methane also can be represented accurately by equations of state such as the Benedict-Webb-Rubin equation or the Martin-Hoa equation.' These equations, although valuable for analytically representing the PVT behavior, lack precision in the compressed-liquid region. The PVT data available in the literature for methane were utilized to obtain reduced densities that are presented in Figs. 1 and 2 as functions of reduced temperature with parameters of reduced pressure. The critical constants utilized for methane are those reported by Kobe and Lynn." These critical values, T, = 191.1°K, PC = 45.8 atm and p, = 0.162 gm/cc, produce a critical compressibility factor, z, = 0.289. The PVT information, as obtained from the literature, in many instances was not presented in an immediately

By H. T. Kennedy, J. E. Little

An empirical equation for the prediction of the viscosity of several pure paraffin hydrocarbons and nitrogen is presented. It involves temperature, pressure and six constants of the material, and it applies reliably to both liquids and gases. The equation is similar in form to van der Waal's equation of state. For the paraffin hydrocarbons methane through r-hexane and nitrogen, an average absolute deviation of 1.9 percent was obtained on 1,006 data points described in the literature by 14 authors. When this equation is extended to complex, liquid hydrocarbon mixtures, a correlation was obtained with an average absolute deviation of 9.9 percent. INTRODUCTION Equations describing the flow of gas and liquid through porous media contain the viscosity coefficient of the fluid. If other pertinent variables remain constant, the volume rate of flow is inversely proportional to this coefficient. In dealing with condensate fluids and volatile oils, however, the compositional effects resulting from changing pressure materially affect the viscosity. The effect of compositional changes also may be significant in certain secondary recovery or pressure maintenance processes, notably miscible displacement or gas injection. Early attempts to describe the performance of reservoirs utilized a volumetric material balance method wherein gas and liquid in the reservoir were identified as produced gas and liquid at the surface. This method of analysis proved adequate for reservoirs at moderate temperature and pressure that contained gas with moderately low amounts of condensable materials. The volumetric material balance procedures for "black oil" reservoirs leave much to be desired when applied to condensate and volatile oil reservoirs because phase behavior and compositional changes are relatively more important in these cases. The alternative is a compositional material balance, which in turn, requires a correlation of properties of the reservoir fluid with composition. This paper supplies this correlation in regard to viscosity, for reservoir crude oils. REVIEW OF LITERATURE The literature contains many empirical equations describing the effects of composition, temperature and pressure on the viscosities of pure liquids and binary liquid mixtures. However, the applicability of a majority of these equations is limited to very low pressures and to a small number of systems. Most of them, when applied to complex hydrocarbon systems, are of little value. The lack of utility of the majority of equations results from the fact that they were developed to show the separate effect of temperature, pressure or composition on viscosity, but not to predict the viscosity as a function of all three variables. And with the few exceptions noted below, they were developed to apply to much simpler systems than hydrocarbon mixtures. Eakin1 developed an equation for the viscosity of pure light hydrocarbons applicable in the gas, liquid and dense fluid regions, and in its most general form is written as µ - µo =A(ebP - e-cpk) The pure component viscosity at atmospheric pressure and temperature T, µo, was calculated using the Sutherland equation: µo = BT3/2/T + S The constants in Eakin's equation permitted the term e-CPk to become quite small at high densities (above 0.215 gm/cc), and the general equation would simplify to

Jan 1, 1969

By Enrique Hauser

Discussion of the paper of ENRIQUE HAUSER, presented at the New York meeting, February, 1916, and printed in Bulletin No. 110, February, 1916, pp. 521 to 534. GEORGE A. BURRELL, Pittsburgh, Pa. Mr. Hauser determined the influence of inert gases on the properties of fire-damp mixtures, and one might take some issue with his work on the effect of hydrogen, on those mixtures. It is valuable information to have, but as far as the practical application of it goes, it is not of much significance, because one does not find it in mines, at least not in any appreciable amount. The same Fig. 1.-APPARATUS FUR DETECTING COMBUSTIBLE GASES IN AIR. thing can be said of ethylene. Ethylene is formed in considerable quantities as a product of heat reaction, but normally, in the coal mines, as a result of a great many analyses we have made at the Bureau of Mines, we cannot say that ethylene is a combustible gas found in the fire-damp in coal mines.

Jan 5, 1916

By V. N. Cox, C. H. Grant

Both the underwater method and the rock cratering method contribute useful information in evaluating and comparing new explosive compositions. Results indicate that metallized explosive systems which contain no high explosives can be formulated to offer more energy and to do potentially more work than other slurry explosive systems. Commercial field trials are presently being conducted with a metallized explosive that does not contain any high explosive. Although there have been only a few commercial blasts to date, the results are encouraging. The Dow Chemical Co. has been engaged in explosives research since 1953. The incentive for this research was based on the desire to find additional markets and uses for magnesium metal of which DOW is a heavy producer. The thermodynamic calculations for explosives incorporating light metals indicate energies approximately twice that of nonmetallized explosives currently being used. Early in 1960, additional patent rights were acquired from independent inventors and a development station at Virginia, Minn., was set up to conduct a commercial testing and market evaluation program with this new product. The actual manufacturing of the explosive was to take place in the boreholes. Due to the complexities involved in controlling the manufacture of the explosive in the borehole, performance was sometimes in doubt. Furthermore, experience indicated that the metallized explosive systems could not be evaluated by the standard laboratory methods then being used, and it was necessary for us to look for more suitable research methods. Tests were conducted on several known methods for explosive evaluation. The most promising were the Navy underwater method for energy evaluation and a modified Livingston crater method intended to compare explosive performance in rock. After shooting over 2000 charges of 8 to 60 lb under water and in taconites and granites with many types of explosives, a metallized explosive was developed which is now being manufactured at a pilot plant in Virginia, Minn. and is being evaluated commercially. UNDERWATER METHOD An underwater method of evaluating explosives has a number of advantages. A permanent site can be established and operated indefinitely. A large number of charges may be fired in a short length of time. The medium in which the explosive is shot is exactly the same for successive shots. Data are recorded electronically minimizing human judgment. A disadvantage of the underwater method is the lack of confinement provided by the water medium. Many explosive compositions require a greater degree of confinement to develop good energy release. Although the relative significance of the two components of energy release that are usually measured, the shock energy and the bubble energy, and their relationship to breaking rock are not clearly understood, the underwater method is a valuable research tool for screening explosive compositions and gaining an insight into the manner in which explosives release their energy. Fig. 1 shows schematically the underwater apparatus. The piezoelectric gages are calibrated tourmaline crystals which convert the shock wave to an electrical impulse which is transmitted to an oscilloscope and recorded on Polaroid film. The shock traces are picked up by scopes 1 and 2. Scope 1 is timed to pick up the full trace whereas scope 2 is timed to cover the first increment of the shock trace expanded, and therefore allows full inspection of the propagation and peak pressures. Scope 3 records the bubble trace. When the explosion occurs, a bubble of gaseous detonation products is formed. The outward momentum of the water keeps expanding the bubble even after

Jan 1, 1963

By A. B. Cook, C. J. Walker, G. B. Spencer

Although water will displace oil from a petroleum reservoir to a greater extent than gas will, there are some reservoirs in which gas rather than water should be used for pressure maintenance. This is indicated by the high-percentage oil recovery in the case history reports on the Pickton field1 in East Texas and the Raleigh field2 in Mississippi. In gas cycling, total oil recovery includes vaporized products from the immobile oil in addition to oil produced by displacement. If a large part of the immobile oil is vaporized, total oil recovery may be higher than that obtainable from pressure maintenance by water injection. However, the amount of immobile oil vaporized may range from almost 0 to 100 percent, and no simple and reliable method has been presented in the literature for calculating oil vaporization. Development of such a method was the purpose of this study. Several papers concerning the calculation of oil vaporization have been published,3 6 and each is based on the concept dealing with K values — equilibrium constants for the various components in reservoir oils. (K for a component in a system of vapor and liquid in equilibrium is the mole fraction in the vapor phase divided by the mole fraction in the liquid phase.) This appears to be the most logical approach. Yet, the problem is too great for a perfect solution, even with modem laboratories to analyze oils and high-speed digital computers to perform calculations; the number of components in reservoir oils is too large. Because of their complexity no two reservoir oils are exactly alike. Also, K values change not only with variations in pressure and temperature but also with composition of the reservoir oil. In addition, during gas cycling, the lighter hydrocarbons tend to vaporize first. Thus, the reservoir oil becomes more dense and less volatile as gas cycling continues. Furthermore, the greatest amount of oil vaporization occurs near the injection well. Therefore, a simplified method is required for calculating oil vaporization because a rigorous method is not practical. The present methods of calculating vaporization generally require knowledge of the reservoir oil composition in which the mole fractions of the lighter components through hexanes are given; the remaining oil is described as C, + (heptanes and others of greater molecular weights). Assigning one K value for the C,+ system provides the simplest solution, but this incorrectly assumes that K for C, + does not change as gas cycling continues. This assumption can cause large errors because the K value may become less than a thousandth of the original after a large amount of oil vaporization. Another way to solve the problem is to determine changing K values for C, + according to how much of the C, + has been vaporized. To determine these values a sample of the reservoir oil is injected into a pressure-volume cell, and dry gas is batched in and out of the cell. Appropriate measurements and tests are made to determine K

Jan 1, 1970