Search Documents

By Axel Sahlin

Having no business interest in any of the various machines used for granulating ores, my remarks on the subject are prompted solely by my desire to contribute towards the determination of the best apparatus such information, as I have been able to gather,

Jan 1, 1893

By W. E. D. Stokes

THE application of aviation to mining and petroleum operations, on the basis of economy and attainment, has become a demonstrated fact. According to Dominion Government records, 30 Canadian companies engaged in air transportation with northland mining camps, carried 13,000 tons of freight in 1935, and in 1936 carried 16,000 tons. In addition, many mining companies used planes to transport prospectors with canoes, outboard motors, gasoline drums and supplies to remote locations, and one company maintained a plane to take radium-bearing ore from its mine on Great Bear Lake nearly 1000 miles south to the rail- way at Waterways, Alberta. Heavy pieces of mining machinery weighing 2 tons are being carried by Ford trimotor planes over the Andes, at an altitude of not less than 15,000 feet. At Bulolo, New Guinea, 17,500 tons of disassembled gold-dredging equipment has been flown over inaccessible terrain, while supplies and building materials for a community of 1000 whites have been transported by airplane. A 3 1/2-ton shaft was carried on one notable trip, while at another time an automobile, 2 tons of rice and a large safe went over the mountains in one plane. Fig. 1 shows an electric generator that was carried by plane to Bulolo. Many large American oil companies operated fleets of planes during 1936, ranging from fast, luxurious air liners for executives to elaborate geological prospecting units, ambulance planes, freight carriers and pipe-line repair crews. According to the U. S. Geological Survey, commercial firms photographed 16,291 square miles during 1935, while the Government compiled from aerial photographs by the Geological Survey in the United States 85,544 square miles. It was found that information could be determined in a fraction of the time required by surface methods. However, sounding a word of caution, a Commonwealth of Australia government report on aerial survey operations had this to say: It is necessary to repeat ad nauseam that aerial photographic survey is in no sense a substitute for geological investigation, though it is a very valuable adjunct thereto, and greatly facilitates and accelerates geological mapping. There are many details

Jan 1, 1937

By Carlton Beal

This paper presents useful charts for conversion of various viscosimeter units into centipoises and graphically summarizes published investigations of the viscosity of air, water and natural gas at high temperatures and pressures. Where possible, charts and correlations were constructed to cover a range of temperature (60°F. to 300°F.) and pressure (14.7 lb. per sq. in. abs. to 8000 lb. per sq. in. gauge) encountered in oil fields. Correlation charts, for the purpose of predicting crude oil viscosity and solubility behavior at oil field temperatures and pressures, were constructed from an analysis and correlation of 1332 viscosity and solubility observations from 953 crude oil samples taken from 747 oil fields. Of these fields, 501 are in the United States, including 75 in California. Of the 1332 observations, 1215 were viscosity values, including 786 of gas-free crude oil, 351 of oil saturated with gas, and 78 under-saturated with gas at pressures above the bubble point. Results show that crude oil viscosity at various reservoir conditions can be predicted with average deviations varying from 24.2 per cent for gas-free crude oil to as little as 2.7 per cent for undersaturated crude oil "above the bubble point." The solubility of crude oil may be predicted from oil gravity and saturation pressure with an average deviation of 22.0 per cent from observed values. In conclusion, it is shown that the viscosity of a gas and gas-saturated crude oil may be predicted within the accuracy of most reservoir computations. Viscosity is a measure of flow resistance; common units are in centipoises (cp.). The object of this report is to present conversion charts and simplified 'correlation charts showing the viscosity behavior of air, water, gas and crude oil and gas-saturated crude oil under oil field temperatures and pressures, which have been known to be as high as 300°F. and 8000 Ib. per sq. in. gauge under subsurface conditions. The viscosity of a liquid or gas is important because it is a measure of its fluidity, or its ability to flow through pipe lines or oil sands. Viscosity of air has become important because it is being used in secondaryl,2 recovery projects. The viscosity of natural gas and water at reservoir conditions is useful because they are almost universally associated and produced with crude oil. The viscosity of gas-saturated crude oil at reservoir pressures and tem peratures is of particular value in making estimates of oil reserves3 and rate of oil recovery from flush pools when production decline data are limited. Instruments Used to MeasuRe Viscosity There are many types of viscosimeters,4 and the instruments most commonly used to measure the atmospheric viscosity of liquids at moderate temperatures and their practical working ranges have been deScribed5,6 in detail. Fig. I and Fig. 2 are viscosity conversion charts which have

Jan 1, 1946

By Carlton Beal

This paper presents useful charts for conversion of various viscosimeter units into centipoises and graphically summarizes published investigations of the viscosity of air, water and natural gas at high temperatures and pressures. Where possible, charts and correlations were constructed to cover a range of temperature (60°F. to 300°F.) and pressure (14.7 lb. per sq. in. abs. to 8000 lb. per sq. in. gauge) encountered in oil fields. Correlation charts, for the purpose of predicting crude oil viscosity and solubility behavior at oil field temperatures and pressures, were constructed from an analysis and correlation of 1332 viscosity and solubility observations from 953 crude oil samples taken from 747 oil fields. Of these fields, 501 are in the United States, including 75 in California. Of the 1332 observations, 1215 were viscosity values, including 786 of gas-free crude oil, 351 of oil saturated with gas, and 78 under-saturated with gas at pressures above the bubble point. Results show that crude oil viscosity at various reservoir conditions can be predicted with average deviations varying from 24.2 per cent for gas-free crude oil to as little as 2.7 per cent for undersaturated crude oil "above the bubble point." The solubility of crude oil may be predicted from oil gravity and saturation pressure with an average deviation of 22.0 per cent from observed values. In conclusion, it is shown that the viscosity of a gas and gas-saturated crude oil may be predicted within the accuracy of most reservoir computations. Viscosity is a measure of flow resistance; common units are in centipoises (cp.). The object of this report is to present conversion charts and simplified 'correlation charts showing the viscosity behavior of air, water, gas and crude oil and gas-saturated crude oil under oil field temperatures and pressures, which have been known to be as high as 300°F. and 8000 Ib. per sq. in. gauge under subsurface conditions. The viscosity of a liquid or gas is important because it is a measure of its fluidity, or its ability to flow through pipe lines or oil sands. Viscosity of air has become important because it is being used in secondaryl,2 recovery projects. The viscosity of natural gas and water at reservoir conditions is useful because they are almost universally associated and produced with crude oil. The viscosity of gas-saturated crude oil at reservoir pressures and tem peratures is of particular value in making estimates of oil reserves3 and rate of oil recovery from flush pools when production decline data are limited. Instruments Used to MeasuRe Viscosity There are many types of viscosimeters,4 and the instruments most commonly used to measure the atmospheric viscosity of liquids at moderate temperatures and their practical working ranges have been deScribed5,6 in detail. Fig. I and Fig. 2 are viscosity conversion charts which have

Jan 1, 1946

By Günter Petzow, Siegfried Steeb

DURING an investigation of the U-Co-A1, U-Ir-A1, and U-Ni-A1 ternary systems, three previously unknown phases have been identified which correspond to the phases UCoA1, UIrA1, and UNiA1, and which be- long to the C22 (Fe2P) type (P 6 2m). In the following the structure determination of the new ternary phases is reported. The alloys have been prepared in an arc furnace under purified argon from 99.97 wt pct U, 99.83 wt pct Co, 99.99 wt pct Al, 99.95 wt pct Ni, and 99.97 wt pct Ir. The pellets obtained were bright and brittle; chemical analyses agreed with the nominal compositions. To obtain equilibrium, the samples were annealed in evacuated quartz capsules at 950°C for 30 hr and then quenched in water. Guinier photographs were processed and evaluated according to a method already described previously.1,2 The hkl's, measured and calculated sin2 6 values, and in the case of UCoAl only calculated and estimated intensities, are given in Tables I to 111. Lattice parameters refined according to the least-squares method and the respective mean-square variations are presented in Table IV. Lattice parameters, especially the c/a ratio, lead to the assumption that the three phases are crystallizing in the C22-type structure. For the calculation of the intensities given in Table I, the atomic positions recently reported by wyckoff3 have been accepted for the C22 (Fe2P) type. The uranium atoms were assumed to be in phosphorus sites; furthermore the scattering amplitudes of the boron atoms have been averaged in the stoichiometric formula UB2 according to atomic percentages; i.e., for these atoms a statistical distribution on the iron sites was assumed. The two smallest distances between a uranium and

Jan 1, 1967

By G. Sachs, G. Espey

Simple methods can be used for the determination of the residual stresses in thin walled tubing if the stresses consist of high tensile stresses at the one surface and high compressive stresses at the other surface. Such a stress distribution is approximated in sunk tubing and in deep-drawn shells. It is possible to cut sections from such articles in such a manner that a major part of the residual stress is relieved by a bending of the Part that is separated from the adjacent metal. The circumferential stress is evaluated from the change in diameter D, occurring on splitting a length of tubing (Fig. Ib) or from deflection resulting from the cutting out of a circumferential tongue (Fig. Ia) .6 The bending moment M released by such a flexure is. M - EI (I-I) = EI R1-R0 I - u2(R0-R1) I -V2 R0XR1 where E = modulus of elasticity, v = Poisson's ratio, I = moment of inertia of the section, R0 =. mean radius before splitting, R1 = mean radius after splitting. The release of this bending moment corresponds to the release of stresses, which, in a thin section, vary linearly from one surface to the other, and have the following value in a particular fiber a certain distance (x) from the surface: S = c d I - V2 x c x d 2 2 RI - R0____E (d ) R1 - R0 X R0 X R1 ~ i - u2 (2 x) Ro X R1 where C is the section modulus and d is the thickness of the section, the ratio M:C being d: 2 for a rectallgular section.* The maximum stress values at the surfaces are. Smax = +E X d X R1-R0 I — u2- 2 R0 X R1 For most purposes the simpler equations: EXd AD1 = ± I - u2 x dm2 and will be satisfactory, where D,, is the aver. age of the mean diameters of the tubing

Jan 1, 1942

By G. Sachs, G. Espey

Simple methods can be used for the determination of the residual stresses in thin walled tubing if the stresses consist of high tensile stresses at the one surface and high compressive stresses at the other surface. Such a stress distribution is approximated in sunk tubing and in deep-drawn shells. It is possible to cut sections from such articles in such a manner that a major part of the residual stress is relieved by a bending of the Part that is separated from the adjacent metal. The circumferential stress is evaluated from the change in diameter D, occurring on splitting a length of tubing (Fig. Ib) or from deflection resulting from the cutting out of a circumferential tongue (Fig. Ia) .6 The bending moment M released by such a flexure is. M - EI (I-I) = EI R1-R0 I - u2(R0-R1) I -V2 R0XR1 where E = modulus of elasticity, v = Poisson's ratio, I = moment of inertia of the section, R0 =. mean radius before splitting, R1 = mean radius after splitting. The release of this bending moment corresponds to the release of stresses, which, in a thin section, vary linearly from one surface to the other, and have the following value in a particular fiber a certain distance (x) from the surface: S = c d I - V2 x c x d 2 2 RI - R0____E (d ) R1 - R0 X R0 X R1 ~ i - u2 (2 x) Ro X R1 where C is the section modulus and d is the thickness of the section, the ratio M:C being d: 2 for a rectallgular section.* The maximum stress values at the surfaces are. Smax = +E X d X R1-R0 I — u2- 2 R0 X R1 For most purposes the simpler equations: EXd AD1 = ± I - u2 x dm2 and will be satisfactory, where D,, is the aver. age of the mean diameters of the tubing

Jan 1, 1942

Jan 1, 1893

By Kenneth L. Clifford, Thomas G. White

Many changes have occurred in lead-zinc milling in the U.S. since the 1970 Lead and Zinc Symposium. New mills have been built, novel flowsheets devised, and the use of process control techniques is increasing. The Missouri concentrators produced in 1976 more lead and zinc concentrate than those in any other state. The Tennessee mills produced nearly as much zinc concentrate in 1976 as Missouri and with some very highly developed flowsheets. Progress in the Western U.S. is centered around new operations in Colorado and Utah. New mills in New York and Missouri incorporate many of the newer control concepts. Most of the older concentrators have undergone flowsheet changes to accomodate different ore types or because of the availability of new and better equipment.

Jan 1, 1977

By E. D. Hummer

During the planning of the fine coal cleaning addition at the Gary, W. Va., Coal Preparation Plant of United States Steel Corp. careful consideration was directed toward the problem of minus 48 mesh refuse disposal. The plant flow sheet included wet screening the raw coal at ¼ in., pumping the minus ¼ in. raw coal to bowl desiltors for desliming, and then cleaning the deslimed raw coal on tables. The material to be removed from the circuit as minus 48 mesh refuse consisted of the minus 100 mesh desiltor overflow thickened in two 120 ft diam thickeners and minus 1 mm screened from the table refuse. Later a flotation plant was added to recover coal from the desiltor overflow, and a classifier replaced the screens to make a sharper size separation of the table refuse at 48 mesh. Today the minus 48 mesh refuse at Gary consists of thickened flotation tailing and thickened spiral classifier overflow. The transportation of large quantities of fine, wet refuse filter cake mixed with jig and table refuse over existing conveying equipment would have been neither dependable nor economical. Filtration was therefore impractical, and disposal of the fine materials was a problem. The decision was to pump the material to storage ponds.

Jan 3, 1965

By D. W. Gillmore, C. C. Wright

DEWATERING is a major problem in the preparation and utilization of fine-sized coals now being recovered in increasing amounts from colliery effluents, refuse banks, and silt ponds. Of the various methods which have, been proposed for dewatering fine coal, gravity drainage is probably the oldest and most widely used, and yet little information on this subject has been published. Consequently a laboratory investigation of gravity drainage phenomena was undertaken not only to provide basic data for predicting the drainage behavior of fine coals in cars, bunkers, basins, and silos, but also to obtain some understanding of the moisture retention properties of fine coal. Because of practical interest in the problem, a reproducible laboratory technique, which simulated plant conditions as closely as possible, was developed for the evaluation of drainage variables. The ap- paratus selected consisted of vertically-mounted lucite columns, nominally 4 to 8 ft in length and 2, 3 or 6 in. in diam; the bottom of the column was provided with a rubber stopper through which was inserted a -in. diam drainage tube; a wad of glass wool was placed on the upper part of this outlet tube to retain the coal and a plug was temporarily located at the lower end to prevent drainage during charging of the column, see Fig. 1. The most satisfactory testing procedure is as follows: 1—The fine coal, which has been thoroughly

Jan 1, 1953

By D. W. Gillmore, C. C. Wright

DEWATERING is a major problem in the preparation and utilization of fine-sized coals now being recovered in increasing amounts from colliery effluents, refuse banks, and silt ponds. Of the various methods which have, been proposed for dewatering fine coal, gravity drainage is probably the oldest and most widely used, and yet little information on this subject has been published. Consequently a laboratory investigation of gravity drainage phenomena was undertaken not only to provide basic data for predicting the drainage behavior of fine coals in cars, bunkers, basins, and silos, but also to obtain some understanding of the moisture retention properties of fine coal. Because of practical interest in the problem, a reproducible laboratory technique, which simulated plant conditions as closely as possible, was developed for the evaluation of drainage variables. The ap- paratus selected consisted of vertically-mounted lucite columns, nominally 4 to 8 ft in length and 2, 3 or 6 in. in diam; the bottom of the column was provided with a rubber stopper through which was inserted a -in. diam drainage tube; a wad of glass wool was placed on the upper part of this outlet tube to retain the coal and a plug was temporarily located at the lower end to prevent drainage during charging of the column, see Fig. 1. The most satisfactory testing procedure is as follows: 1—The fine coal, which has been thoroughly

Jan 1, 1953

By N. A. Gjostein

The formal descnption of the decay of an isolated scratch can be written in terms of an appropriate Fourier integral. With the application of certain approximations, this description leads to the second-moment analysis developed by King and Mullins. It is shown that the second-moment analysis does not adequately describe the behavior of a decaying, isolated scratch, and that it is difficult to obtain internally consistent values of the surface self-diffusion coefficient, D,,from this analysis. This inconsistency arises from contributions of higher-order moments to the decay process. Two new analytical methods, based on the Fourier integval representation of a scratch profile, have been developed for determining values of D, from the decay kinetics of isolated scratches. In the profile-matching method, a given initial profile is resolved into its Fourier components, the components are allowed to decay exponentially over an annealing time interval, At, and finally they are synthesized to form a decayed profile. The predicted profile is adjusted to fit the experimental profile by varying only the parameter, D,At; the contribution from volume diffusion is directly taken into account and the integrals involved in the Fourier analysis are evaluated by numerical integration on an electronic computer. In the kinetic method, the widths, W, of a series of predicted profiles are determined for various volume-diffusion contributions, and the results are plotted in the form W vs DSAt. For a fixed ratio of surface to volume diffusivities, R, these plots are linear; the slope of the lines increases as R decreases. A simi-lar analysis can be performed on the depths of the predicted profiles. In this case, plots of (l/df3 vs D,At are found to be nonlinear. The proper value of Ds is found by superimposing the data for the experimental widths and depths on Plots of this kind and adjusting Ds until a fit is obtained with the line having the appropriate value oj R. An analysis of the slopes of these plots for R —= (surface diffusion only) reveals that they do not give the correct value, 2276, for width tf second-moment standard profile. This finding allows one to modify the second-moment analysis in such a way as to obtain internally consistent values of Ds. ThERE are several processes, involving a flow of matter along a surface, due to gradients in curvature of the surface, that may be used to determine the surface self-diffusion coefficient, Ds. These processes are commonly known as 1) thermal grooving, 2) single and 3) multiple scratch decay, 4) blunting of field-emitter tips, and 5) sintering. In this series of papers, the first three of these processes are examined in detail in an attempt to subject the phenomenological theories governing these processes to a critical test and, furthermore, to make a critical comparison of the values of D, obtained from these techniques. Although the thermal grooving technique has been established'' as a reliable method for determining D, , this has not been the case for the process of scratch flattening. Mullins and coworkers374 have proposed phenomenological laws governing the growth of the width and the decay of the amplitude of an isolated scratch, but only in one study5 has the width-growth law been supported by experimental evidence. In paper I of this series detailed experimental evidence is presented that points up the deficiencies of the second-moment analysis; the origins of these deficiencies are also discussed. Following this development, some new analytical procedures, based on representing the scratch profile by Fourier integrals, are presented as a means of obtaining reliable values of D, from the kinetics of the scratch-flattening process. Further applications of the Fourier integral approach are reported in papers 11' and 111,~ where the method has been useful in interpreting the anisotropy of flattening rates of isolated scratches (paper 11) and in describing the multiple scratch decay process (paper 111). THEORY OF SCRATCH DECAY Following ~ullins" approach we consider an arbitrarily shaped scratch profile, Fig. 6, at some time, t. The profile is designated by a function z(x,t) which can be expressed in terms of a Fourier integral as follows:

Jan 1, 1967

By G. Birkbeck, T. C. Wells

Pearlite morphology in three commercially produced, low-carbon steels has been studied using optical and electron microscopy. A reduction in the cooling rate from 600° to 6°C per hr increased the interlamellar spacing but had little effect on pearlite content. The form of the lamellae was also assessed in terms of a new parameter, the ratio of lamellae length to interhmellar spacing. The variation of this parameter with cooling rate depended on steel composition. For each steel, it would appear that an optimum cooling rate exists at which lamellar pearlite can be produced to the greatest extent, though not necessarily the same extent in different steels. PeARLITE morphology depends principally on steel composition, prior austenite grain size, and either cooling rate or isothermal transformation temperature. Alloying with manganese,'" fine grain size,374 and slow cooling rates reduce the likelihood that a well-developed lamellar structure will be produced. Instead, imperfect structures usually consisting of partly globular carbides are formed in a ferrite matrix. Such structures have been described as granular.1,2,5 degenerate,6 or semi4-pearlite. Manganese additions also progressively lower the eutectoid temperature and the carbon concentration at the eutectoid so so that, for a given carbon content. the pearlite fraction is increased as manganese is added. While the effects of manganese are understood in a general way, quantitative studies of pearlite morphology have been carried out mainly on steels of eutectoid composition;9 little quantitative data is available for low-carbon steels. Certain anomalies exist in the behavior of steels containing -1.5 pct Mn, however. The pearlite morphology in niobium-treated steel differs from that which is generally accepted as typical of low-carbon/manganese steels in that a reasonably well developed lamellar structure is sometimes observed. In the present work, the effects of changes in cooling rate from 600° to 6°C per hr on the pearlite morphology in three commercially produced, low-carbon steels have been studied by both optical and electron microscopy. I) MATERIALS The steels were of the following types: 1) 0.1 pct C, 0.4 pct Mn steel (En 2, a composition similar to specification SAE 1010); 2) niobium-treated, 0.2 pct C, 1 1/2 pct Mn steel (BS 968:1962); 3) silicon-killed 0.2 pct C, 1 1/2 pct Mn steel (En 144, a composition similar to SAE 1320). The chemical analysis for each steel is given in Table I. To facilitate the comparison of pearlite morphology in the three steels, steels 1 and 3 (which were obtained as 1/2-in. plate and 1-in.-diam bar, respectively) were hot-rolled in the laboratory to $ in. thickness, finishing at around 850°C. The BS 968:1962 was received as 1/8-in. plate finished at 885°C and was not further worked. Specimens from each steel were heat-treated in vacuo at 950°C for 1 hr and cooled at rates of 600°, 60°, and 6°C per hr. 11) EXPERIMENTAL Pearlite contents were determined by point counting (500 points) on longitudinal sections mechanically polished and etched in 2 pct nital. Pearlite morphology is usually considered in terms of interlamellar spacing. For pearlite colonies randomly distributed in space and having a true, constant interlamellar spacing So Pellissier et al.9 determined the distribution curve for the apparent spacings of the lamellae on a random section plane. The areal fraction fs of the pearlite colonies occupied by apparent spacings up to a limiting value S is given by:

Jan 1, 1969

Jan 1, 1908

By J. S. Dodge, I. M. Krieger

The laminarturbulent transition was studied for monodisperse rigid polymer latices as functions of particle diameter and concentration at several tube diameters. Breaks in the graph of friction factor US Reynolds number were used to characterize the transition; rates of shear were sufficiently high so that Newtonian behavior at the high shear limiting viscosity could be assumed. No abnormal behavior was noted at concentrations below 30 volume percent. At higher concentrations, suppression of turbulence was observed with critical Reynolds numbers rising to as high as 5,500 at 50 volume percent. No definite effect on the transition was noted due to variation of particle size or tube diameter. The type of surfactant used to stabilize the latex was found to have an important influence on the critical Reynolds number. INTRODUCTION For some polymer solutions and solid-in-liquid suspensions, the transition from laminar behavior to turbulence in pipeline flow has been reported to take place at unusually low1-5 and in some cases at higher than normal6,7 critical Reynolds numbers. Extreme values of critical Reynolds numbers (NRe*) have been found to range from 1 to 6,000, as opposed to the values of 2,000 to 3,500 which are obtained for ordinary liquids. Observations have also been reported of drag reduction in fully developed turbulent flow, i.e., a decreased value of the friction factor in the Reynolds number-friction factor plot6*10-14 in colloidal and polymeric fluids. In some cases, turbulence suppression (unusually large values of Nre*) was noted to occur preceding the drag reduction,6 while a normal transition was noted in other cases of drag reduction.14 To relate the occurrence of turbulence suppression and drag reduction to the molecular or particulate structure of the fluid, it is desirable to characterize the structural variables with some precision. For most colloidal dispersions and polymeric fluids, this is very difficult; particle sizes, shapes and stiffnesses may vary over wide ranges in the same sample. Recently, however, a colloidal fluid has become available which consists of rigid spherical particles of highly uniform diameter dispersed in a simple Newtonian liquid. These are monodisperse polymer latices, prepared by carefully controlled emulsion polymerization reactions. The purpose of this research was to study the transition to turbulence in these well-characterized colloidal fluids as a function of concentration, particle diameter and tube diameter. Typical flow curves of latices of high solids content are shown in Fig. 1. While these curves were obtained with monodisperse latices, the general shape of the curve is also typical of a polydisperse latex. It is observed that a shear-thinning region occurs for low to moderate rates of shear. At sufficiently high shearing rates a limiting Newtonian viscosity is reached. The viscosity µ diminishes with decreasing volume fraction 4 at all shear rates. For 4 less than 0.25, latices are Newtonian at all rates of shear.15 An increase in average particle size (at a given extends the Newtonian viscosity region to lower rates of shear. A third factor which influences the

By AIME

Compiled by Mary Ann Gross, Department of Energy. Bruceton. PA. Supplemental data from Coal Age indicated by (*). Reprinted by permission from CoalAge magazine.

Jan 1, 1981

By Alfred H. Brooks

PetRoleUm seepages are known in Alaska at four localities, all on Pacific seaboard. These, named from east to west, are Yakataga, Katalla on Controller Bay, Iniskin Bay on Cook Inlet, and Cold Bay on the Alaska Peninsula. Besides these, a petroleum residue has been found near Smith Bay on the north Arctic coast of the Territory. At Katalla, Cold Bay, and Iniskin Bay there has been some drilling for oil, and in the first field several productive wells have been opened up. Alaskan petroleum, so far as its composition is known, is a refining oil, with paraffin base and low sulphur content. Most of what is known of the geology of Alaska oil fields is based on the investigations of Dr. George C. Martin, of the United States Geological Survey. A, G. Maddren, of the same service, has recently made a reconnaissance of the Yakataga oil field. The data to be here presented are chiefly taken from the following publications: G. C. Martin: The Petroleum Fields of the Pacific Coast of Alaska, with an account of the Bering River coal deposits, Bulletin No. 250, U. S. Geological Survey, pp. 9 to 27 (1905). G. C. Martin: Geology and Mineral Resources of the Controller Bay Region, Alaska, Bulletin No. 335, U. S. Geological Survey, pp. 112 to 130 (1908). G. C. Martin and F. J. Katz: A Geologic Reconnaissance of the Iliamna Region, Alaska, Bulletin No. 485, U. S. Geological Survey, pp. 126 to 130 (1912). A. G. Maddren: Mineral Deposits of the Yakataga District, Bulletin No. 592, U. S. Geological Survey, pp. 143 to 147 (1914). Katalla Field The Katalla field is marked by a series of seepages and gas springs occupying an east and west belt, about 25 miles in length and from 4 to 8 miles wide. (See Fig. 1.) This zone skirts the north shore of Con-

Jan 1, 1915

By J. L. Huitt, B. B. McGlothlin, D. K. Lowe

Published correlations for the principal aspects of hydraulic fracturing were combined into a digital computer program to facilitate the study of interrelated variables. The computer program includes individual relationships for fracture width during pumping, fracture area generated, propping agent embedment, flow capacities of propped fractures and transport of propping agents in horizontal fractures. The effects of more than 20 treatment and formation parameters on the predicted results of hydraulic fracturing treatments were studied. The effects of these parameters were determined for (I) fracture width during injection, (2) fracture width after the overburden comes to rest on the propping agents, assumed not to be crushed, (3) generated and propped fracture area, (4) location and concentration of propping agents in the fracture when injection ceases, (5) flow capacities of the various propped sections of the fracture and (6) expected increase in the well productivity. The effects of propping agent, formation and fracturing fluid parameters on well productivity are discussed. The parameters that were found to have the most pronounced effects on hydraulic fracturing treat~nents are injection rate, treatment volume, fracturing fluid coefficient, size and amount of propping agent, spearhead volume, well drainage radius and formation capacity. INTRODUCTION Many correlations have been published for predicting effects of various parameters that are considered in the design of hydraulic fracturing treatments. The Carter equation' can be used to predict generated fracture radius as a function of fracture width, fracturing fluid leakoff and other parameters. Fracture width can be determined by use of the Perkins and Kern correlation' in which the fracture width is related to the fracture radius, fluid injection rate and certain formation and fracturing fluid parameters. Wahl and Lowe et aL4 have reported methods of predicting the location of propping agents in fractures when pumping ceases. The former study is applicable to the case where the ratio of propping agent diameter to fracture width is less than 0.1. The latter is applicable when this ratio is greater than 0.1. These studies showed that the propping agent placement in horizontal-radial fractures depends principally on how the individual particles are transported in the fracture by the carrying fluid. Particle transport in fractures is determined by local fluid velocity in the fracture, fluid and particle properties, and the size of the particle relative to the fracture width. The distribution of propping agents, effective overburden pressure and formation rock strength control the propped fracture width6 by controlling the extent to which the propping agent particles embed into the fracture faces. From the distribution of propping agents and the propped fracture width, fracture flow capacities can be calculated or the various regions of the fracture. The flow capacities and the radial extent of these regions can be combined with reservoir information to predict the productivity increases for fractured wells. In all these studies, the effects of certain treatment and/ or reservoir parameters on one facet of fracturing can be predicted only if other facets which the parameters affect are fixed. For instance, fracture width and radius are interrelated; that is, to calculate the value of one, the value of the other must be known. Also, some parameters influence more than one aspect of fracturing. For example, prop-pant transport is a function of both fracture width and fluid viscosity. but fracture width is itself a function of fluid viscosity. Since these calculations are complex and the parameters interrelated, it is not possible to write an equation with which the over-all effects of treatment parameters can be solved explicitly. For these reasons, the correlations for determining the effects of the parameters which are most significant in hydraulic fracturing treatments have been incorporated into a digital computer program. COMPUTER PROGRAM The program, which was written for an IBM 7094 computer, can be used to predict results of most of the combinations and values for the treatment parameters that are ordinarily considered for fracturing treatments. A spearhead of fracturing fluid and a propping agent-carrying fluid with different fluid properties can be taken into account. Also, the total volumes and relative amounts of the spearhead and carrying fluids can be varied. Two different propping agents (as used in tail-in operations) and a wide range of formation properties and injection rates are considered. The computer program (Fig. 12) consists of several sets of calculations. First, the final flooded fracture radius and average fracture width at the cessation of pumping are calculated. This is done by simultaneous solution of the Perkins and Kern fracture width equation and the Carter equation for flooded fracture radius (equations used in the computer program appear in the Appendix). The next step is to determine local fluid velocity in the fracture as a function of time and radius. Since it is not possible to write this function in closed form expression, a table of velocity values is generated by the program and stored for subsequent use. The time span from the beginning of

By AIME AIME

THE report of the united. States Steel Corporation for the' year 1928 gives the, following ,figures of raw materials and' finished products . Raw materials, tons :

Jan 1, 1929