Search Documents

By Willard R. Jillson

Within the last ten years Kentucky has become celebrated for its low-sulfur bituminous coals. Prior to this time, many investigators had discnvered the abundance of this coal but the fact was unknown to the general public until the last extension of Kentucky's mountain railroads' was completed. The coals of. Kentucky are, broadly, separated into two geographic units—the eastern and the western coal fields. On a basis of their sulfur content these coal fields fall into three units, which are indicated on the accompanying map, by numbers, 1, 2, and 3. The area marked 1, the very southeastern part of the state, shows the lowest sulfur, all coals being averaged. In this field the maximum is 1.04 per cent. and the minimum as low as 0.68 per cent.—both averages and not individual analyses, which would of course show much greater variation. Included in this group are the following counties: Lawrence, 0.87 per cent.; Martin, 0.75 per cent.; Johnson, 0.73 per cent.; Magoffin, 0.87 per cent.; Floyd, 0.88 per cent.; Pike, 0.68 per cent.; Knott, 1.04 per cent.; Perry, 0.77 per cent.; Letcher, 0.80 per cent.; Leslie, 0.70 per cent.; Harlan, 0.79 per cent.; Knox, 0.86 per cent.; and Bell, 0.92 per cent. The total average sulfur content for these thirteen counties is 0.82 per cent. Bordering the restricted low-sulfur field of eastern Kentucky is the No. 2 belt (slightly darker), 'which is commonly known as the western border of the eastern coal field. Here the maximum sulfur is 2.91 per cent. and the minimum 1.05 per cent. Again these are averages and not individual analyses. The counties in this group are: Greenup, 2.60 per • cent.; Boyd, 1.67 per cent.; Carter, 1.07 per cent.; Morgan, 1.40 per cent.; Wolfe, 1.83 per cent.; Lee, 2.93 per cent.; Breathitt, 1.85 per cent.; Owsley, 1.09 per cent.; Jackson, 1.06 per cent.; Rockcastle, 2.30 per cent.; Clay, 1.09 per cent.; Laurel, 1.08 per cent.; Pulaski, 2.24 per cent.; Whitley, 1.05 per cent.; McCreary, 1.50 per cent.; and Wayne. The sulfur percentage average for this area including fifteen important counties is 1.65 per cent. In the western Kentucky field, the great increase in the percentage of sulfur is at once apparent. Ohio County with 3.78 per cent. is the highest and Hancock with 2.87 per cent. is the lowest. The percentages by counties is as follows: Hancock, 2.87 per cent.; Henderson, 3.22 per cent.; McLean, 3.27 per cent.; Muhlenberg, 3.31 per cent.; Hopkins,

Jan 1, 1920

By P. H. Bohart

THe answer to the title question will be found by considering the ultimate influence of the petroleum engineers on industry and by considering the tools with which petroleum engineers must be equipped. Some 2000 to 3000 petroleum engineers are being absorbed by the industry yearly, indicating the magnitude of influence which they may exert on the industry. The petroleum industry is a highly technical one and it is inevitable that management be entrusted more and more to executives with technical training. The industry, therefore, expects the petroleum engineering schools to provide a broad basic training and a social and economic orientation which goes beyond merely supplying the minimum technical qualifications necessary to discharge responsibilities for a particular job. Because of the many branches of industry, emphasis must be placed on the Student's being trained in fundamentals, as well as grounded in other engineering branches, such as mechanical, civil and electrical engineering. The petroleum engineering graduate must know how to outline and organize the material for an engineering report and should be thoroughly grounded in English, both written and spoken. Even to the end that a five-year course may be required, they should receive instructions in auxiliary cultural and social science courses, industrial relations, and public speaking. Training in the method of approach to an engineering problem is something in which industry is deeply interested. A standard of ethics, a desire to arrive at the truth, and leadership training should be stressed. Introduction The basis of an an to any ques- tion concerning the petroleum engineering schools' responsibility to the industry will be found by considering the ultimate influence of the petroleum engineers on the industry and by considering also the tools with which petroleum engineers must be equipped in order to achieve success and happiness in the role in which their training, individual talents and choice of a profession ultimately will cast them. The problem must always be considered as it relates to a large group and not as it relates to a particular individual possessing strong natural talents and destined, by virtue of them, to pursue a certain course. The question implies that petroleum engineering schools have a responsibility, and indeed they do, not only to the industry but also to the students and even to the public. Because the responsibility is concentrated in a relatively few schools whose graduates are absorbed largely by the one industry, the responsibility is more apparent and the effect of failure to discharge the responsibility fully is more serious than the effect of similar failure on the part of schools supplying a smaller portion of the industry's intake of engineering graduates. This fact is especially true with respect to an industry such as the petroleum industry whose ratio of technical to nontechnical employees is relatively high and is increasing.

Jan 1, 1948

By T. P. Meloy

The behavior of slime-sized particles in a flotation cell is neither well documented nor well understood. In general, slimes (or 'fines' as they will be alternatively called in this chapter) have been either ignored, eliminated, or treated empirically. Recently, however, some good experimental and theoretical work has been done that leads qualitatively to a clearer understanding of the problems involved. Fines are small particles that, owing to their smallness, behave differently from larger particles of the same material. The smaller the particle, the larger is its surface-to-mass ratio. The difference in the surface-to-mass ratio for a micron-sized particle and a centimeter-sized particle is 10,000 to 1. For our purposes we shall arbitrarily define any particle less than 10 µ 10-3 cm) as being a fine. This is generally considered the upper limit of colloidal-sized particles. Recently, the field of colloidal chemistry has successfully used the concept of the double layer to explain dispersion and flocculation of fine particles in a static system. Flotation, however, is a dynamic system due both to the turbulence of the liquid and to the coursing of the bubbles through the cell. Further more, a flotation system is a heterogeneous system in which there are a number of different types of particles which interact with one another. A further difficulty is that no two flotation systems are the same for minute quantities of a variety of common chemicals will alter the behavior of the system. These reasons make reproducible experimental results that can ex- plain fine-particle behavior difficult to obtain. In spite of these experimental difficulties, millions of tons of ore containing a large amount of fines are successfully floated every year. For example, at the Morning mill of American Smelting and Refining Co., a galena concentrate that was 42% -6 µ and a sphalerite concentrate that was 33% -9 µ were floated leaving a tailing that was 27% - 9 µ.l Ores that contain fines can generally be divided into two classes. In the first are those in which the valuable mineral tends to slime, such as galena and hematite. In the other class, one or more of the gangue minerals forms the slime; a typical example is barite. Though the problems due to the presence of fines are often the same in handling the two classes of ore, the method of handling them is different because of economic considerations. Nonmetallic mineral flotation is particularly affected by slimes, by either slime-coating effects or excessive reagent consumption. Though a given gangue mineral may not float under the conditions prevailing in the flotation cell, it may still be able to adsorb the collector. I£ it also slimes, a significant

Jan 1, 1962

By P. H. Bohart

THe answer to the title question will be found by considering the ultimate influence of the petroleum engineers on industry and by considering the tools with which petroleum engineers must be equipped. Some 2000 to 3000 petroleum engineers are being absorbed by the industry yearly, indicating the magnitude of influence which they may exert on the industry. The petroleum industry is a highly technical one and it is inevitable that management be entrusted more and more to executives with technical training. The industry, therefore, expects the petroleum engineering schools to provide a broad basic training and a social and economic orientation which goes beyond merely supplying the minimum technical qualifications necessary to discharge responsibilities for a particular job. Because of the many branches of industry, emphasis must be placed on the Student's being trained in fundamentals, as well as grounded in other engineering branches, such as mechanical, civil and electrical engineering. The petroleum engineering graduate must know how to outline and organize the material for an engineering report and should be thoroughly grounded in English, both written and spoken. Even to the end that a five-year course may be required, they should receive instructions in auxiliary cultural and social science courses, industrial relations, and public speaking. Training in the method of approach to an engineering problem is something in which industry is deeply interested. A standard of ethics, a desire to arrive at the truth, and leadership training should be stressed. Introduction The basis of an an to any ques- tion concerning the petroleum engineering schools' responsibility to the industry will be found by considering the ultimate influence of the petroleum engineers on the industry and by considering also the tools with which petroleum engineers must be equipped in order to achieve success and happiness in the role in which their training, individual talents and choice of a profession ultimately will cast them. The problem must always be considered as it relates to a large group and not as it relates to a particular individual possessing strong natural talents and destined, by virtue of them, to pursue a certain course. The question implies that petroleum engineering schools have a responsibility, and indeed they do, not only to the industry but also to the students and even to the public. Because the responsibility is concentrated in a relatively few schools whose graduates are absorbed largely by the one industry, the responsibility is more apparent and the effect of failure to discharge the responsibility fully is more serious than the effect of similar failure on the part of schools supplying a smaller portion of the industry's intake of engineering graduates. This fact is especially true with respect to an industry such as the petroleum industry whose ratio of technical to nontechnical employees is relatively high and is increasing.

Jan 1, 1948

By B. M. Larsen

IN spite of the great amount of operating experience and of studies directed toward elucidation of the interrelation of the several processes going on in the blast furnace, the present picture of its operation is based largely on indirect evidence and on speculation. The reason for this unsatisfactory situation is the extreme difficulty of making direct observations and significant measurements while the furnace is operating, and of translating laboratory experience to the conditions obtaining in the smelting process. Yet there are questions that can be investigated with some promise of a definitive answer; and it is the purpose of the present paper to set forth some of these questions, as a basis for discussion, which may, it is hoped, lead to active work on some of them. In order to provide a starting point for this discussion, some views held tentatively by the writer are outlined, along with a certain amount of evidence gathered from operating records and interpreted on the basis of the ideas expressed. These opinions may prove not to be pertinent, and may well not be in accord with the opinions of others; they have been included to encourage others to do likewise and on purpose to provoke informal discussion. In this way it should be possible to transcend the limitations of formal papers; to bring forth suggestions of new hypotheses or new methods of investigation; and to promote better understanding between investigators and operators in this difficult field. The statements made express the present views of one physical chemist who is not too familiar with all the difficulties of maintaining these enormous reaction vessels in reasonably smooth and economical operation; who, however, is persuaded that some of these difficulties would be less formidable if surer knowledge of the several interrelated and concurrent processes were available and made use of. For these difficulties arise largely, it is believed, from variation in conditions which the operator must meet as best be can, and should lessen as these conditions can be brought under control so that they are either made substantially constant or the variations are neutralized by compensating methods. A. VARIATION IN BURDEN AND IN MIXING OF CHARGE To the scientific investigator, who aims to control the conditions of experiment and to limit the number of variable factors, the most perplexing aspect of the blast-furnace process is that nothing seems to remain constant for long. Variations of burden seem to be unnecessarily frequent; and, even if the charge were constant from hour to hour, its constituents usually are not mixed as intimately as they might well be. It is understandable that separate layers of coke and of ore, also separate layers according to size, promote permeability and the

Jan 1, 1947

By G. R. St. Pierre, M. J. Pool, R. Speiser

The activities of chromium in solid Cr-Ti solulions contaitning from 10 to 90 at, pct Cr were measured over the temperatutre range 1250" to 1380°C. The Knzudsen effusion technique with direct weighing of the Knudsen cell was used. This was possible because the vapor pressure 0f chromium is about 600 times greater than the capor pressure of titanium so that the amount 0f titanium effusing could be neg-Lected. the activities of titanium were determined by g-raphical integration of the Gibbs-Duhertz equatiotz. The solutio~ls sllcdied exhibited positive deviations from Raoult's law, but they approached ideality as the temperature was raised from 125 to 138t C. The integral free energy of mixing was calclated from the activity values. AF has a maximum value of - 1.66 kcal a1 1390°C and a composition oj 55 at. pcl Cv. The free energy of formation of 'TiCr, is —0.76' lzcal a1 1250°C. The heats and entropies of mixing usere determined using the second law and thus tile mlues could be in error by a considerable amount. The Cv-Ti system is not a regular solution. ThE increasing emphasis on application of various models to alloy formation and solution behavior has led to an intensified interest in the thermochemistry of alloys. In a majority of the cases where different models have been applied the alloy system used was a simple one. However, the amount of data available on more complicated systems is small and in order to fully understand the behavior of metallic solutions more experimental work has to be carried out on these systems. The Cr-Ti system is reasonably complex and even though this system was being investigated as early as 1937 by Kroll,' no previous work has been done on the thermodynamic properties of either the liquid or solid state. This research was undertaken to determine the thermodynamic properties of solid Cr-Ti alloys. EXPERIMENTAL PROCEDURE The experiments on the Cr-Ti system were carried out in a water-cooled metal vacuum system. This system is composed of a brass vaporization chamber wound with copper tubing for cooling purposes, a brass liquid-nitrogen trap, an oil diffusion pump, and a Duo-Seal forepump. The pressure in the system is measured with an ionization gage. The interior of the vaporization chamber is shown schematically in Fig. 1. The Knudsen cell is supported by the molybdenum tripod, and a magnetically operated shutter is provided to protect the Pyrex sighting glass against deposition of metallic vapors in the intervals between temperature measurements. Heating is accomplished by means of a 20-kw 450-kc Westinghouse oscillator. The susceptor in this case is 10-mil tantalum seamless tubing which also forms the outside of the sample container. The remainder of the sample container consists Of an im~ervims re-crystallized alumina crucible and a 10-mil-thick tantalum lid. Two holes are drilled in the lid. The diameters of these holes were measured with the calibrated stage of a "Tukon" microhardness tester. The thicknesses were measured with a micrometer. A Leeds and Northrup disappear ing-filament optical pyrometer was used to measure the temperature. This is sighted through the smaller of the two holes in the lid and corrections for the emissivity of the sample were not necessary since the ratio bf cell length to hole diameter was large enough to insure a hohlraum. The pyrometer was calibrated at the ~~t~~~l Bureau of Standards, and the correction for the Pyrex sighting window was determined in the laboratory. The titanium used in this investigation was obtained from Mallory-Sharon and the chromium from The Union Carbide Metals Co. The reported purities were

Jan 1, 1968

By G. H. Kesler, C. E. Dryden, J. H. Oxley

Rates of heat- and mass-transfer from rods to recirculating air were determined within a one-quarter-scale model of a metals deposition bulb. The dependence of local and averaged rates of transfer upon outlet geometry, number of rods, position upon a rod, and airflow rate was established for flow patterns created by radial and tangential air inlets. The results are given a qualitative interpretation in terms of rates and distribution of metal deposition in a prototype deposition bulb. IT is well known that a number of metals can be prepared from volatile compounds containing them by causing dissociation of their compounds to occur at heated surfaces. Examples of such metals include silicon, titanium, zirconium, and aluminum, which can be prepared from their halides. While the ability to prepare high-purity metals by this technique is of importance in itself, the rate of deposition of metal also is important in considering the commercial potential of such a process. If the over-all process of deposition is broken down into separate steps of transport of materials to and from the heated deposition surface and establishment of chemical equilibrium at the surface, as was described in a recent paper,1 then since surface temperatures are usually high and surface equilibrium is attained very rapidly, the rate-controlling step can be considered as either that of transport of reactants to the surface or of transport of products other than deposited metal away from the surface. These two transport steps are related through the reaction stoichiometry. It can be shown1 that the rate and distribution of metal in such a transport-controlled process are determined by the over-all and local mass-transfer coefficients in conjunction with the equilibrium conversions at the hot surfaces. The rate relationship is: w = [KA(y a-y s)] B.E where the bracketed group is the rate of arrival of the reactant at the surface of area A as determined by the mole fraction driving force (y a - y s) and the convective transport constant K. The factor B is the weight fraction of metal in the reactant, and E is the equilibrium fractional conversion to metal. Numerical values of y, and y, are fixed by the vapor composition and the equilibrium compositim of reactant at the surface; and E can be determined by thermodynamic calculations or by experiment. However, values for the convective coefficient K are available generally only for systems of simple geometry, such as cross-flow past a rod or sphere or parallel flow past a flat plate or through a tube. Often it is necessary to work with systems of more complex geometry and in which the vapor flow may not be uniform or exactly parallel or perpendicular to the surface. It is possible, of course, to establish working correlations for these more complex systems by direct experiment; but the experiments usually are costly in both time and money. To circumvent this objection and to aid in arriving at an understanding of the relationships and interactions of the process variables, it is helpful to construct a model of the flow system in which air, water, or other convenient fluids can be substituted for the process vapors and with which mass-transfer measurements or analogous he at-transfer measurements can be made upon a simulated deposition surface. By this means, a large number of measurements can be made in a reasonable time, and the correlation of these results can be applied to the prototype system. Examples of the application of these techniques are to be found in the literature (e.g., Ref. 2). It is the purpose of this paper to illustrate these techniques by description of a model study which was done at Battelle Memorial Institute. The work to be described was of a preliminary nature; the investigation subsequently was extended to develop quantitative relationships among the significant dimensionless transport, geometric, and flow parameters. This more detailed study will be the subject of a future paper.

Jan 1, 1962

By C. M. Wayman, R. Bullough

The selection rules fir twinning in fcc, bcc, and hcp lattices as established in Part I are used to predict the positions, relative to the matrix (untwinned) reciprocal lattice, of reflections due to double diffraction from the twins. These doubly diffracted reciprocal lattice points, which often coincide with ordinary matrix reciprocal lattice points, are deduced from kinematical considerations. 1. INTRODUCTION. A Bragg diffracted beam under certain conditions may be rediffracted to produce a third beam i .e., the incident beam, the first diffracted beam, and the rediffracted beam). In such a case an initial incident beam is diffracted by a set of planes; the diffracted beam may then be rediffracted by a second set of crystallographically different planes. The initially diffracted beam acts as an incident beam for the second set of planes, thus resulting in a twice-diffracted beam. This twice-diffracted beam may be regarded as a singly diffracted beam from a set of lattice planes whose structure factor is zero. In order to find the reciprocal lattice point corresponding to the doubly diffracted beam, it is customary1,2 to combine vectorially the component reciprocal lattice vectors which head to the extra reflections. For example, hcp metals may exhibit the forbidden (0001) reflections.* As Fig. 1 shows, the (0001) reciprocal lattice point results from: (1010) + (1011) = (0001) Similarly, combination of the beams (1011) and (1010) gives rise to: (1011)- (1010) = (0001) As Pashley and stowell2 have pointed out, even if in fact the double diffraction points lie on the Ewald sphere, they may not actually represent observed doubly diffracted beams because the initially diffracted primary beam, i.e., the source, must also correspond to a reciprocal lattice point which lies on the sphere of reflection; however, the reciprocal lattice point corresponding to the planes causing the rediffraction need not lie on the sphere. A little consideration of the Ewald construction in Fig. 2 makes this argument clear. In Fig. 2 the incident beam is represented by the wave vector ko. This beam is initially diffracted by the lattice planes (hkl) which are represented by the reciprocal lattice vector r(hk1) The reciprocal lattice point (hkl) of course lies on the Ewald sphere in order that the Bragg condition be satisfied. The initially diffracted beam is represented by the wave vector kl. If now kl is rediffracted by a second set of planes (h'k'l') resulting in k2 the lattice point (h"k"l") must lie on the Ewald sphere if the double diffracted beam k2 is to be observed. It is clear that (h"k"l") is the vector sum of (hkl) and (h'k'l'), and that (h'k'l') need not in fact lie on the Ewald sphere. In the usual1,2 analysis the origin (000) of the reciprocal lattice is shifted to the point (hkl). However the translation (000)— (hkl) automatically shifts the point (h'k'l') to the Ewald sphere, i.e., (h'k'l')— (h"k"1"). It is to be seen that (hkl), the source reflection, must lie on the Ewald sphere; if it does not the translation (000) — (hkl) places the "new" origin off the sphere, not satisfying the Ewald construction and hence the Bragg condition. The point (h"k"1") of course must lie on the sphere if it corresponds to an observed doubly diffracted beam. It is to be noted that crystals having the hcp structure, as just discussed, are intrinsically capable (i.e., without lattice defects such as twins) of double diffraction, as contrasted to fcc and bcc crystals which are extrinsically capable (requiring lattice defects or reoriented lattice regions). The "allowed" bcc reciprocal lattice points (which comprise an fcc lattice) are (110), (2001, (211), (220), (310), (2221, (321), (4001, and so forth. The combination

Jan 1, 1967

By J. B. Clark

The oil industry has long recognized the need for increasing well productivity. To meet this need, a process is being developed whereby the producing formation permeability is increased by hydraulically fracturing the formation. The "Hydrafrac" process, as it is now being used, consists of two steps: (1) injecting a viscous liquid containing a granular material, such as sand for a propping agent, under high hydraulic pressure to fracture the formation; (2) causing the viscous liquid to change from a high to a low viscosity so that it may be readily displaced from the formation. To date the process has been used in 32 jobs on 23 wells in 7 fields, resulting in a sustained increase in production in 11 wells. INTRODUCTION Need For Process Although explosives, acidizing, and other methods have long been used, there still exists a need for artificial means of improving the productive ability of oil and gas wells, particularly for wells which produce from formations which do not react readily with acids. This paper discusses the development of a hydraulic fracturing process, "Hydrafrac", which shows distinct promise of increasing production rates from wells producing from any type of formation. The method is also considered applicable to gas and water injection wells, wells used for solution mining of salts and, with some modification, to water wells and sulphur wells. Requirements of Process In considering such a possible process, it appeared that certain requirements must be met. Some of these are as follows: A. The hydraulic fluid selected must be sufficiently viscous that it can be injected into the well at pressure high enough to cause fracturing. B. The hydraulic fluid should carry in suspension a propping agent, such as sand, so that once a fracture is formed, it will be prevented from closing off and the fracture created will remain to serve as a flow channel for oil and gas. C. The fluid should be an oily one rather than a water-base fluid, because the latter would be harmful to many formations. D. After the fracture is made, it is essential that the fracturing fluid be thin enough to flow hack out of the well and not stay in place and plug the crack which it has formed. E. Sufficient pump capacity must be available to inject the fluid faster than it will leak away into the porous rock formation. F. In many instances, formation packers must be used to confine the fracture to the desired level, and to obtain the advantages of multiple fracturing. Development of Process As a necessary step in the development of this process, it was deemed advisable to determine if the Hydrafrac fluids were actually fracturing the formation or whether these special fluids were merely leaking away into the surrounding formation. To determine this, a shallow well, 15 feet deep, was drilled into a hard sandstone. Casing was set, the plug drilled, and the well deepened in the conventional manner. A fracturing fluid dyed a bright red was used to break down the formation. Sand mixed with distinctively colored solids was injected into the well with the fracturing fluid to prop open any fracture made in the formation. A simulated gel breaker solution dyed a bright blue was then pumped into the well to determine if the gel breaker would follow the first solution. The results are shown in Figure 1. It was noted that a fracture was formed about the well bore, that the propping agent was transported back into the break, and that the breaker solution did actually follow the fracturing gel out into the fracture. While it is realized that this shallow well test is probably not exactly equivalent to a deep test, the results were interpreted as being a definite indication of what happens down the hole during a Hydrafrac job. Of interest in this connection is an investigation reported by S. T. Yuster and J. C. Calhoun, Jr.' This study, re~orted after the Hydrafrac work was under way, presents some excellent field data supporting the theory of fracturing a formation with hydraulic pressure. METHOD Steps of Hydrafrcu: Process Figure 2 shows a simplified cross-sectional view of a well treated by one version of the process. The first step, formation breakdown, is done with a viscous fluid, usually consisting of an oil such as crude oil or gasoline, to which has been added a bodying agent. Due to availability and price, war-surplus Napalm has been used in the majority of experiments to date. Napalm is the soap which was used in the war to make "jellied gasoline". The next step consists of breaking down the viscosity of the gel by injecting a gel-breaker solution and then after several hours, putting the well back on production. Figure 3 shows diagram-matically, a typical field hookup. The oil or gasoline is unloaded into the 10 bbl. tank shown on the left rear of the truck. This base fluid is picked up by the mixing pump and pumped through the jet mixer, where the granular soap is added. Next it goes into a small mixing tub, from which the high-pressure pump takes suction. The solution is then pumped into the well. The breaker solution is then taken from an extra tank and is displaced into the well immediately following the gel. When required, additional trucks may

Jan 1, 1949

By J. B. Clark

The oil industry has long recognized the need for increasing well productivity. To meet this need, a process is being developed whereby the producing formation permeability is increased by hydraulically fracturing the formation. The "Hydrafrac" process, as it is now being used, consists of two steps: (1) injecting a viscous liquid containing a granular material, such as sand for a propping agent, under high hydraulic pressure to fracture the formation; (2) causing the viscous liquid to change from a high to a low viscosity so that it may be readily displaced from the formation. To date the process has been used in 32 jobs on 23 wells in 7 fields, resulting in a sustained increase in production in 11 wells. INTRODUCTION Need For Process Although explosives, acidizing, and other methods have long been used, there still exists a need for artificial means of improving the productive ability of oil and gas wells, particularly for wells which produce from formations which do not react readily with acids. This paper discusses the development of a hydraulic fracturing process, "Hydrafrac", which shows distinct promise of increasing production rates from wells producing from any type of formation. The method is also considered applicable to gas and water injection wells, wells used for solution mining of salts and, with some modification, to water wells and sulphur wells. Requirements of Process In considering such a possible process, it appeared that certain requirements must be met. Some of these are as follows: A. The hydraulic fluid selected must be sufficiently viscous that it can be injected into the well at pressure high enough to cause fracturing. B. The hydraulic fluid should carry in suspension a propping agent, such as sand, so that once a fracture is formed, it will be prevented from closing off and the fracture created will remain to serve as a flow channel for oil and gas. C. The fluid should be an oily one rather than a water-base fluid, because the latter would be harmful to many formations. D. After the fracture is made, it is essential that the fracturing fluid be thin enough to flow hack out of the well and not stay in place and plug the crack which it has formed. E. Sufficient pump capacity must be available to inject the fluid faster than it will leak away into the porous rock formation. F. In many instances, formation packers must be used to confine the fracture to the desired level, and to obtain the advantages of multiple fracturing. Development of Process As a necessary step in the development of this process, it was deemed advisable to determine if the Hydrafrac fluids were actually fracturing the formation or whether these special fluids were merely leaking away into the surrounding formation. To determine this, a shallow well, 15 feet deep, was drilled into a hard sandstone. Casing was set, the plug drilled, and the well deepened in the conventional manner. A fracturing fluid dyed a bright red was used to break down the formation. Sand mixed with distinctively colored solids was injected into the well with the fracturing fluid to prop open any fracture made in the formation. A simulated gel breaker solution dyed a bright blue was then pumped into the well to determine if the gel breaker would follow the first solution. The results are shown in Figure 1. It was noted that a fracture was formed about the well bore, that the propping agent was transported back into the break, and that the breaker solution did actually follow the fracturing gel out into the fracture. While it is realized that this shallow well test is probably not exactly equivalent to a deep test, the results were interpreted as being a definite indication of what happens down the hole during a Hydrafrac job. Of interest in this connection is an investigation reported by S. T. Yuster and J. C. Calhoun, Jr.' This study, re~orted after the Hydrafrac work was under way, presents some excellent field data supporting the theory of fracturing a formation with hydraulic pressure. METHOD Steps of Hydrafrcu: Process Figure 2 shows a simplified cross-sectional view of a well treated by one version of the process. The first step, formation breakdown, is done with a viscous fluid, usually consisting of an oil such as crude oil or gasoline, to which has been added a bodying agent. Due to availability and price, war-surplus Napalm has been used in the majority of experiments to date. Napalm is the soap which was used in the war to make "jellied gasoline". The next step consists of breaking down the viscosity of the gel by injecting a gel-breaker solution and then after several hours, putting the well back on production. Figure 3 shows diagram-matically, a typical field hookup. The oil or gasoline is unloaded into the 10 bbl. tank shown on the left rear of the truck. This base fluid is picked up by the mixing pump and pumped through the jet mixer, where the granular soap is added. Next it goes into a small mixing tub, from which the high-pressure pump takes suction. The solution is then pumped into the well. The breaker solution is then taken from an extra tank and is displaced into the well immediately following the gel. When required, additional trucks may

Jan 1, 1949

By E. C. Barfield, D. G. McCarty

Two-dimensional analyses offer considerable promise in providing the basic information required to effect more precise control of petroleum reservoir performance. This paper describes a method for conducting some engineering analyses of this type using a high-qreed digital computer. The general approach is to (I) develop a representative numerical model of the reservoir and (2) employ a suitable numerical technique to solve the basic equation of flow and to develop the required engineering information for the particular case represented by the model. This technique has been used to calculate pattern performance in connection with several field projects involving water injection into oil reservoirs. This type of analysis involves sequentially (I) calculating the pres-' sure distribution in the reservoir and (2) tracking the progress of the interface between the displaced and displacing fluids in a step-wise manner to provide a depletion history of the operation. The particular analysis presented in this discussion is subject to the restriction that the mobility of the displaced and displacing fluids be equal and assumes that the fluids are incompressible and that gravity and capillary pressure do not affect the shape of the flood pattern. Complete two-dimensional .flexibility is maintained with regard to definition of reservoir rock and fluid characteristics, placement of physical restrains and boundary conditions, investigation of flow characteristics in the reservoir, and movement of the displaced and displacing fluid interface. The results of these studies indicate that the highspeed digital computer is well suited for conducting reservoir performance studies in two-space dimensions and that dependable numerical techniques are available for making such analyses. INTRODUCTION One of the problems facing petroleum engineers is that of providing more precise engineering control over reservoir operations at all stages during the depletion cycle. The most direct approach to this problem is to simulate the reservoir with some type of physical or mathematical model and use the information developed from the behavior of the prototype to predict the performance of the actual reservoir. This principle is, of course, not new, and a variety of simulation techniques have been developed to provide information on the behavior of oil and gas reservoirs. One straightforward refinement of present reservoir analysis techniques is the progress from one-dimensional reservoir models to those where the reservoir can be represented in two- and three-space dimensions. Multidimensional reservoir analyses are desired in order to show not only what is happening in the reservoir but in addition to show where in the reservoir the various phenomena of interest are occurring. Thus far, two-dimensional analyses have been confined predominantly to potentiometric and electrolytic model studies. However, when large numbers of wells are involved and where large-scale models are required, such problems become very tedious and instrumentation problems become acute." It is possible through use of numerical methods and general purpose digital computers to avoid many of the instrumentation difficulties and to relax some of the restrictions of previous techniques. The purpose of this paper is to discuss some recent experiences with the use of high-speed digital computers for predicting pattern performance in fluid injection operations and to describe the numerical procedures which have been used for making such analyses. THE NUMERICAL MODEL The problem of representing the reservoir with an appropriate numerical model is composed of at least two important considerations. The first requires that an adequate description of the reservoir mechanism be formulated in the form of the basic differential equation describing the process. The second consideration is definition of the physical characteristics of the reservoir. The latter requires that representative basic data be obtained, examined, and interpreted to provide sufficient detailed definition of the pertinent physical properties throughout the system. The Differential Equation One of the basic assumptions made for the purpose of performing the two-dimensional studies being dis-

By Wayne E. Glenn

The San Miguelito Field of Ventura County, California, lies along the coast approximately five miles northwest of the city of Ventura. The structure is one of several in the Ventura anticline, a series of anticlinal folds trending northwest-southeast through the San Miguelito Field. The Ventura Avenue Field lies to the southeast along the anticline, and Padre Canyon and Rin-con Fields lie to the northwest; the series then extends beneath the Pacific Ocean. Production is obtained from the Lower Pico and Repetto formations of Pliocene age, and oil accumulation is controlled by faulting, lithology and structure. The first well in the First Grubb Pool was completed in November, 1931, and to date 48 producing wells have been drilled on the 276 productive acres. Approximately 685 ft of net pay are contained in the eight sand bodies grossing 1,220 ft of section. Production of oil to date has been 15,-600,000 bbl. A gas injection project was instigated in April, 1940, but was discontinued in September, 1949. This paper is confined to a discussion of the development and performance history of the First Grubb Pool of the San Miguelito Field. LOCATION AND GEOLOGY The San Miguelito Oil Field is located in Ventura County about five miles northwest of the city of Ventura along the coastal highway to Santa Barbara. It is located within sections 23, 24, 25, and 30 of Township 3 North, Range 24 West, San Bernardino Base and Meridian. Surface elevation in the field ranges from sea level to plus 1,200 ft, and the terrain is extremely rugged. The structure is well defined on the surface and on the wall of the Amphitheatre Canyon where appears a so-called "text book" anticlinal section. The San Miguelito structure is an asymmetrical, elongated anticline about two miles long and three-fourths of a mile wide. The trend of the axis of the anticline is northwest-southeast. Fig. 9 illustrates the contours and major faults of this structure. Dips on the southwesterly flank average about 40 degrees, whereas the northeasterly flank ranges from 60 to 75 degrees. There are several faults striking generally east and west. These faults do not appear to affect materially the area of accumulation or performance behavior. It is believed that accumulation was affected by steeply dipping beds and extremely low permeability of the northeast flank of the structure. As shown by the typical log, Fig. 1. the San Miguelito Field is divided into two stratigraphic producing intervals. The First Grubb Pool comprises the interval from the "G" to "L" markers; the Second Grubb Pool from the "L" to "M" markers. The First Grubb Pool, ranging in thickness from 1,105 to 1,680 ft, with an average of 1,220 ft, includes approximately 685 ft of oil-bearing sands. The First Grubb zone consists of eight sand intervals separated by shales. The shale bodies separating the eight sand intervals within the First Grubb Pool are continuous throughout the major part of the structure and serve as an effective barrier. Many of the thin, irregular shale or siltstone streaks within these sand bodies are discontinuous and do not affect well performance appreciably. The productive horizons in the First Grubb Pool are siltstones and sandstones of the Lower Pico and Repetto formations of Pliocene age. It is apparent from inspection of the typical log that the "G" to "L" interval might have been separated into at least six producing zones and selectively pro. duced. Unfortunately the electric log and some of our modern tools were not available or were in their infancy dur. ing the period 1931 through 1937. Not until 1937 when eight wells had been completed were there sufficient geological data to permit definition of the sand bodies. Then the problem of reservoir control in this multiple-sand field, where all sands are open to a common

Jan 1, 1950

By F. W. Jessen, H. Kazemi

A washing technique has been developed to form a spherical cavity in massive salt. The technique is, basically, a process of controlling the fluid motion in the cavity, the concentration distribution of brine, the rate of dissolution at the walls of the cavity, and the particular settings of the wash pipe, or pipes (inlet and outlet pipe) employed. For a particular arrangement of tubing and casing, fresh water enters from the casing annulus at a position above the bottom of the tubing, and the dissolved material leaves the cavity through the tubing near the bottom of the projected hole. The initial configuration of the cavity is cylindrical, representing the drilled hole situated at the axis of the projected cavity. As washing progresses the original hole is enlarged radially (any direction perpendicular to the axis of the hole). An inert fluid (generally a hydrocarbon, natural gas, or air) which serves as a blanket is injected through the annular space between the original hole and the casing. The controlled downward motion of this blanket maintains the upper edge of the brine region that is expanding on the surface of the projected sphere. A large spherical glass container and a radial model were used to investigate the flow patterns and the concentration distribution of brine occurring during the progress of solution of the salt. Aluminum flakes and dye were injected and flow patterns were observed and recorded on motion picture film. Based on the results of these model studies and other work carried out in this laboratory, the controlled washing technique was formulated. INTRODUCTION Underground cavities leached in salt beds and salt domes have been used extensively for storing hydrocarbons. Little attention has been focused on the shape of the cavities as long as large volumes were obtained and the cavities were structurally stable. The general process involves circulating fresh water and removing brine. Two distinctly different circulating systems can be employed, each of which yields characteristically different shapes. By far the most generally accepted method of washing storage caverns in salt domes is the direct circulation system. Water is introduced through the innermost pipe of a concentric pipe system, and brine is withdrawn through the annulus. In the reverse-circulation method, water enters through the annulus and is discharged through the smallesc tubing string. The size of the operation, and hence the size of tubing and casing used, depends on such factors as cavity volume, length of hole, availability of water, disposal of brine and power requirements. Underground caverns are from 100 to 1,000 ft or more in length, and when completed vary from 40 to 100 ft in diameter, providing storage of 100,000 bbl or more. The volumes of water required per volume of salt removed vary from 6 to 7 for the reverse circulation method to 10 to 11 for the direct circulation method. Up to 80,000 BWPD are often circulated. In the direct circulation method, because water enters near the bottom of the projected cavity, solution of salt is greatest at this point, particularly in view of the fact that when the washing is begun turbulent flow undoubtedly exists in the immediate vicinity of the discharge of the inlet water. As the diameter of the hole increases, flow rates also are increased, so perhaps the condition of fluid turbulence is maintained for some time. During the later stages of leaching, when capacity of equipment may be the limiting factor in maintaining flow rates, there may still exist some turbulence adjacent to the point of influx; however, for the most part laminar flow conditions prevail. There appears, therefore, rather good evidence to indicate that the state of fluid motion is responsible for the typically shaped cavity which results from such washing — namely a cylindrical form, generally somewhat larger at the bottom than at the top. Depending on the number of so-called stages which might be employed, the walls or surfaces of the salt diminish rather uniformly from bottom to

Jan 1, 1965

By H. Matlock

Beam-columns with continuous or discontinuous transverse and angular loads and elastic restraints are represented mathematically in a manner corresponding to finite-element mechanical models. Solutions for linearly elastic cases are accomplished by direct recursion-equation elimination of unknowns, and for cases involving nonlinear loads and supports, by repeated trial-and-adjustment of complete elastic-type solutions. Five practical applications are presented by which various capabilities of the general method are illustrated. The examples include simulation of a pipeline launching and the bending and collapse of a conductor pipe in deep-water drilling operations. INTRODUCTION The diversified and complex operations being conducted in support of the production of offshore oil have produced a wide variety of new problems in structural mechanics. Many of these problems cannot be solved by conventional analytical or experimental approaches and have therefore stimulated inquiry into new methods of solution. It is the purpose of this paper to give a brief description of an analytical method for solving nonlinear beam-column problems and to illustrate the generality of the method by applying it to several problems related to offshore operations. Although the data for all of the examples are artificial, the cases discussed will be recognized as being based on actual situations. DESCRIPTION OF METHOD The analytical approach which is to be described stems from research on offshore foundations, in particular, the study of lateral-load behavior of foundation piles. The method represents a generalization of a technique which was suggested by Gleser' and extended by Focht and McClelland2 and Reese and Matlock." To maintain generality and therefore versatility in the method, a considerable variety of input data is considered. Fig. 1 illustrates the various types of loads and restraints which might be applied to a short segment of a beam. Each of the terms illustrated in the figure is shown acting in a positive sense. Capital letters are used to designate concentrated quantities whereas the corresponding lower-case symbols indicate effects which are distributed in some way along the length of the beam. The transverse loads and the elastic restraints may be divided into two classes: (1) loads Q and q and springs S and s, which act normal to the axis of the beam, and (2) couples T and t and rotational springs R and r, which act in an angular sense. In addition, there is an axial force P which may be constant or may vary along the length of the beam. The bending stiffness of the beam (the product EI of the modulus of elasticity and the moment of inertia) is designated by the symbol F. The assumptions of conventional elementary beam theory provide the basis for the present method. Among these are the neglect of shearing and axial deformations and the limiting of consideration to straight beams of symmetrical cross-section. Lateral deflections are considered to be small compared to original dimensions. A deformed element of such a beam is shown in Fig. 2a. The curvature of the beam is approximated by the second derivative of the deflection y with respect to distance x along the beam. It is assumed that the material of the beam is linearly elastic and therefore the relation of Eq. 1 between curvature and bending moment M is the ordinary expression from elementary beam theory. A generalized beam element is shown in Fig. 2b. The element has been deflected in a positive direction by an amount y and tilted through a small positive angle dy/dx. The various loads and elastic restraint reactions are detailed in the figure. Temporarily, only distributed effects are considered; provision for concentrated loads and restraints will be made subsequently. Eq. 2 in Fig. 2b is obtained by summing the forces and moments acting on

By C. E. Gregor

Current modern trends in mining show that hoisting shafts are being expected to fulfill an important ventilation function. However, where rigid guides and supporting structures are mandatory, ventilating power costs will be high, unless the transverse structures (shaft sets) are streamlined. In the complex configurations of hoisting shaft furniture patterns, either for original design purposes, or for the adjustment of an existing shaft, a useful degree of aerodynamic optimization, by reduction of overall resistance, can only be effectively evaluated by model studies in a special purpose wind tunnel, as described herein. The author has shown that, in one particular case, an annual power cost saving of more than $50,000 was indicated, mainly by the selection of buntons and dividers of improved aerodynamic profile. Mining companies are now realizing the advantages that can accrue through proper attention to the aerodynamic aspects of mine shaft design. With progressive demands of civilization for an ever greater supply of coal and metals, and with greater depths of working, higher rates of extraction, advanced degrees of mechanization, and improved standards of environmental health and work comfort, shaft design is becoming an extremely important aspect of mining operations. In addition, because the cost of sinking, equipping, maintaining, and operating a mine shaft is continually increasing, there is a greater incentive to use a shaft for the maximum duty, consistent with economy. The current natural tendency is therefore to sink multi-purpose shafts, to install as much equipment as is safely and effectively possible to operate; and then, over and above all these functions, to use the shaft to its maximum possible ventilating capacity, as a major, or auxiliary, airway. Clearly, this latter additional function can only be employed at a significant cost, since the most economic operating airways are smooth-lined and devoid of internal rigid furniture. It is true that hoisting functions can be carried out in such shafts if rope guides are used, but only under restricted conditions. Where a manway must be provided, fixed guides and supporting structures are mandatory, and ventilation costs become accentuated. It is generally recognized that the prime element of cost in ventilating a mine can be represented by the power required to drive the fan unit. It is also known that the largest component of total mine resistance is that of the downcast shaft, where this shaft is equipped with rigid guides for hoisting purposes. The pressure loss varies directly as the resistance of the shaft, in accordance with a modification of the Atkinson equation p = RQ2/3600. Furthermore, the air horsepower drawn is given by the equation AHP = (p x 0/33000. Therefore, the power required to drive the fan unit varies as the resistance of the shaft and the cube of the air quantity. Clearly, if Q is a specified minimum, then the only way to minimize the pressure loss (and the power costs) is to reduce the shaft resistance R. (See Table I for a list of symbols.) The resistance to airflow of a shaft equipped with rigid guides and supporting structural sets consists of two components which, for the sake of analytical simplicity, can be regarded as additive. 1) Frictional resistance, represented by the viscous frictional forces within the boundary layer, as air flows along the surfaces of the walls and the periphery of the vertical elements (guides, pipes, cables, ropes, etc.) of the furniture. The friction loss can be readily calculated by the Atkinson equation, after selection of a suitable value for the friction factor K from the McElroy compilation.' KxLxO xQ2 5.2 xA3 Friction loss can be minimized by lining the shaft walls (particularly if of circular or elliptical cross-section) with a smooth continuous lining of monolithic concrete.

Jan 1, 1968

By John E. Carr

A bewildering variety of continuous thickeners and clarifiers are used in mineral industries and other heavy industries. General topes arc: ? Conventional thickeners ? Thickeners with flocculating feedwells ? Solids contactors ? Tray thickeners ? Lamellas ? Conventional thickeners with Lamellas ? High flow, fluid bed thickeners These units must be adapted to varied applications. They must provide operating flexibility to meet process needs. Such factors multiply the decisions which must be made prior to the selection of a thickener. The seven types of thickeners listed above, with their variations, share design concepts which limit their application and which can be used to evaluate their performance. Some of these arc: ? Feed velocities must be limited and/or controlled. ? Natural or reagent flocculation is necessary. ? Horizontal settling area is the basis for capacity. ? Overflow velocities must be limited and. or controlled. ? A controlled underflow is necessary. Thickener selection is the application of these needs to the basic process concepts. The basic function of a thickener, in the minerals industry, is the separation of a thickener sludge containing fine solids from water. In addition to this primary separation function, others include: ? Separation of solids from water ? Separation of flocculated solids from slimes ? Storage of solids ? Control of percent solids ? Storage of water ? Process control by absorbing fluctuations in feed ? Process control by increase in head clue to the depth and percent solids of the thickened sludge. Most thickener designs, suitable to minerals applications, are based on the separation of minus 65-mesh, 2.7 specific gravity solids containing slimes. These designs are routinely adapted to coarse and or denser materials. Materials of lower specific gravity such as coal, wood, char, some chemical precipitates and organics act Hike slimes if a thickener is used. Perhaps the most conservative approach to sizing and application is the use of comparative data from similar plants or similar industries. Thickeners are-routinely used in a great number of industries from which there is good historical information available. In the minerals industry, thickeners are used in: clay sizing, coal refuse, heavy media, copper concentrates, copper tailings, gold flotation concentrates, and countercurrent decantation for cyanidization in gold and CCD in uranium. Concentrate and tailings thickeners are used in iron ore plants, and concentrate and tailings units in lead-zinc plants. In addition, there arc thickeners designed to handle lime mud; limestone concentrates and tailings; phosphoric acid; potash; and rutile concentrates. Thickeners are, in addition, widely used in gravel plants for water reclamation. This gives excellent data on clay settling in some cases. Thickeners are used extensively for slurries generated by flue dust from wet scrubbers. Thickeners in flue dust applications are operating with considerable variations, treating blast furnace dust. BOF flue dust, and wood char. Similar applications are met in the power industry with fly ash and SO2 scrubber systems. Even though field data is preferred, such data has limitations as listed: ? Flowsheets are seldom identical. ? Data is often incomplete or limited. ? Ores vary significantly. ? Operating needs are changing. For these reasons, mathematical and theoretical methods are needed for application of field data to a specific project. In addition, laboratory data for similar projects and, if possible, laboratory data on a representative sample should be used to specify the thickener. There is a great volume of thickener literature. The references for this paper, along with their references, and past AIME papers provide excellent background. Older literature can be reviewed in light of current methods. It is interesting to note that the requirements for higher flow rate designs, more extreme flocculation systems, etc., can be logically introduced by reading Gaudin's Sections on Flocculation and Thickening, which was written before field installations could illustrate all of his ideas. The problem, then, is to select hardware which will take into consideration the functions of flocculation, separation of the solids, water clarity, control of the underflow, and control of the overflow. Reliable operation and simple control are customary.

Jan 1, 1979

By R. F. Mehl, G. T. Horne

Diffusion coefficients and mobilities were determined as functions of concentration in the a phase of the Cu-Zn system. Use was mode of incremental diffusion couples to determine the Kirkendall effect at various concentrations. Darken's analysis was used to calculate the individual diffusion coefficients and mobilities from these data. The general diffusion coefficient is a single-valued function of the concentration in this system to within the limits of accuracy of the experimental methods used. The form of the various functions (diffusion coefficients and mobilities) of concentration is the same in every case: it is essentially the same as the usual D vs c curve. AS Darken' has pointed out, the knowledge that the movements of different atoms diffusing in the same lattice need not proceed at the same rate is old. Ionic lattices are the classical examples; application of the generalization to lattices with metallic bonding was much more slowly accepted. It remained for Kirkendall'' to show clearly that metallic diffusion need not be symmetrical with respect to both species. Several investigators" have confirmed the initial observation of a "Kirkendall effect:" marker movement accompanying diffusion. They have extended it to other systems and conditions. The theoretical interpretations have been many; ne of the first and simplest of these is that of Darken"' which is essentially mechanism independent. Darken has shown that the Kirkendall experiment permlts a calculation of the individual diffusion coefficients and, coupled with activity data, offers a means of measuring the mobilities of the separateatoms in a binary diffusion process. There is no published information on the effects of concentration and temperature on the individual diffusion coefficients or on the mobility values in binary solid solution; this paper provides such information for zinc and copper in a brass. For purposes of clarity, the overall D value obtained, for example, from the Boltzmann-Matano solution to Fick's law, will be designated as the "general diffusion coefficient;" and the diffusion co- efficients for the separate diffusing species will be designated as the "individual diffusion coefficients." Experimental Procedure The authors have accepted Darken's analysis and have measured mobilities on this basis. The findings of da Silva and Mehl are accepted: that the Matano interface, defined as the origin of coordinates for the Boltzmann-Matano solution to Fick's law, coincides with the original weld interface; that the "shift" of a marker is the distance between the Matano interface and the marker; that the observed porosity is too small sensibly to affect this result; and that the shift is a linear function of the Matano area, the area defined as the absolute magnitude of the area on either side of the Matano interface. Alloys: The a-phase of the Cu-Zn system was chosen for investigation. The measurements were made on incremental couples for two separate reasons: l—it has been questioned whether D, the general diffusion coefficient, is a single-valued function of concentration (a comparison of the D-values for couples with small concentration range, e.g., 10 vs 20 pct Zn, against those from a wide concentration range, e.g., 0 vs 30 pct Zn, will answer that question); and 2—on the Darken analysis, the mobility can be determined only at the concentration at which the marker is observed to lie (this is but one concentration in a single couple) but choosing a series of incremental couples in the range between 0 and 30 pct Zn provides mobility values at several concentrations, thus giving mobility as a function of concentration. Four alloys were chosen; nominally, pure copper, copper with 10 pct Zn, copper with 20 pct Zn, and copper with 30 pct Zn. This allowed six different types of diffusion couples to be prepared, promising mobility determinations at six different concentrations, but having a minimum range of 10 pct Zn. The copper was an especially prepared pure copper supplied through the courtesy of the Ameri-

Jan 1, 1956

By G. W. Thomas

A simplified mathematical model of underground conduction heating in a system of limited permeability is presented. The model applies to underground retorting of oil shale, or to reservoirs containing extremely heavy oils. We assume that heat is introduced at a constant rate into a horizontal fracture which communicates between wells. The radial temperature distribution along the fractured surface is approximated by a step-function. Heat transfer away from the fracture is assumed to be by vertical conduction, and all convection effects are neglected. The model also takes into account the possible temperature dependence of thermal conductivity. A general expression for calculating the growth of the step-function temperature distribution with time is derived. The use of this expression and solutions to the one-dimensional heat equation make it possible to construct isotherms. Expressions for calculating .oil recovery, well spacing and heat efficiency are also given. An example calculation is presented for the conduction heating process in oil shale. Finally, the effect of the heat transfer coefficient between the gas and the fracture boundaries is investigated. INTRODUCTION Thermal recovery processes of oil recovery fall into the four general areas of hot fluid injection, forward combustion, reverse combustion and conduction heating. The first three of* these processes have been rather extensively studied in the past decade from both the experimental and theoretical points of view. As a result, it is possible to make reasonable engineering predictions and analyses of these processes. Little attention has been devoted, however, to the conduction heating process other than to note its possible utility.I To a certain extent, conduction heating cannot be divorced from the other regimes cited above, insofar as these provide the source of heat energy. In the conduction heating process, heat is introduced (either by combustion in the forward or reverse mode or by hot fluid injection) into a small fraction of the total reservoir thickness. This fraction may be either a streak of high permeability or an interwell fracture. Heat penetrates by conduction into the adjacent, less permeable regions of the oil-bearing rock, where the direction of conduction is essentially perpendicular to the streak or fracture. The heated product then drains by gravity or is gas driven to production wells. Conduction heating is probably most applicable to systems containing immobile bitumens such as tar sands and oil shale deposits and perhaps to low-permeability reservoirs containing highly viscous crudes. The mechanism also acts in combination with other thermal processes where fingering or overriding of a bed occurs. It seems probable that in at least one field test, conduction heating of this type was very influential.2 In this presentation we give a first approximation to some of the quantitative aspects of conduction heating. Marx and Langenheim3 treated a similar problem where they focused attention upon the injection interval, which spanned the total reservoir thickness. In their model, conduction heat losses to the bounding media imposed a practical limit on the calculated heated area. In the present study, however, we shall confine the injection interval to a small fraction of the reservoir thickness and assume it has no heat capacity. We therefore direct our attention to regions outside the injection interval into which the conduction of heat is beneficial. In particular, we will endeavor to locate specific isotherms in the media bounding the injection interval. Furthermore, we will construct our model to allow the thermal conductivity to vary arbitrarily with temperature. Thus the model will be applicable to underground retorting of shale where variations in thermal conductivity may be important. THEORY We consider the rock matrix to have no or, at best, a very limited permeability; consequently a

Jan 1, 1965

THERE has been a lot of talk about the shortage of engineers and we have done our share of it; but recently we heard a spot radio commercial-between broadcasts-urging high school seniors to study engineering. We think this is going too far. Engineering is an exacting profession requiring the highest personal qualifications. Publicity for the purpose of interesting capable young people in choosing courses in engineering should emphasize qualification. The opportunities for a career in a great constructive profession which is the basis of our high standard of living should be properly evaluated. Although there is a need for more engineers today, nobody knows what the employment potential is going to be four years hence. The profession may do itself a grave disservice by attempting mass recruitment in the face of an unpredictable national economy. The mining and metals industries are in a special position as fewer engineers are needed in these fields. Only 10 pct of all engineers are mining and metallurgical engineers. Today's shortage is not so much in recent graduates, but in trained practicing engineers. The shortage has served the profession to good effect because it directed the talents of engineers into more productive channels. In some areas high school graduates have been put into routine jobs such as, surveying. If there is one thing the mining profession does not want, it is an oversupply of graduates which would relegate college men to years of work and discouragement in routine jobs. The public has a clear cut picture of the legal and medical professions and thus the respect knowledge imparts. But the engineer is variously thought of as a man who pulls the throttle lever on a locomotive, tends a furnace, or builds a bridge. The term engineer is too loosely used. For instance we have heard of "crowd engineers" and "exterminating engineers." In some unions an engineer designation has been established to describe a person doing routine checking or other control job. We do not disparage these positions but they are different from the work done by the man qualified by years of experience and study. Many of these so-called engineers are not engineers at all. The engineering profession has cooperated closely on a program of recruiting high school graduates for engineering study. It should start a cooperative public relations program to enlighten the public on the accomplishments and work of engineers. This type of publicity need not be of the bombastic type. It can be educational and dignified. Radio, television, and the periodicals can be used to convey information. It is time that laymen and interested young people learn to make a distinction between the true engineer and the technician. Such a program would instruct parents and high school students and thereby interest the right type of people in the engineering profession. The Centennial of Engineering sponsored by the American Society of Civil Engineers is a step in the direction of instructive public relations. This fine start should not be allowed to remain as a one pass venture. The engineering profession as a whole should follow this exposition with a sustained, well directed effort to inform the public of the engineer's work and the possibilities it holds for a constructive career.

Jan 1, 1952

By R. V. Higgins, A. J. Leighton

A fast method is needed for calculating thoroughly the performance of two-phase flobr~ in reservoir rock with complex geotrietry. The authors present such a method and .slzow its accuracy by romparison with laboratory performrrnce of a five-spot pattern. In making the performance c.alculations, the octant of the five-spot pattern was divided into four channels. The computer tinle required for the culculation was one rllinrcte. This is rorighly about minute per channel. INTRODUCTION By using an electronic computer it is now possible to make reservoir performance calculations in which fewer assumptions and more significant variables are included than was possible with a desk calculator. These performance calculations include such variables as the properties of the reservoir rock and its contained fluids, properties of the rock-fluid system and even different well-spacing patterns. Usually, the more variables that are included the higher the computer cost—and time on the computer is expensive. As a result, a short procedure is needed so the many variables will be included and accurate performance calculations will result. If a reservoir is composed of alternate layers of permeable and impervioi~s rock, and if each conducting layer has a different relative permeability, the combined computer cost for calculating each layer by the present long computer programs (which includes different well-spacing patterns and other variables) would be expensive. The performances of complex spacing patterns in a field often are studied by a potentiometric model, but in these model studies the analogy is to the flow of only a single phase. This paper presents a procedure to calculate the performance of a five-spot pattern with two phases flowing using a potentiometric model as a guide for the dimensions of channels and other related data. The time required to calculate the complete performance is about one minute on the IBM 7090 or four minutes on the IBM 704. The procedure can be extended to calculate the performance of more complex well spacings or an entire field with little modification. The procedure has been tested on the performance of the laboratory water floods of a five-spot pattern reported by Douglas, et alThe range of viscosity ratios in the laboratory experiments is from 0.083 to 754, which is extensive. The accuracy of the performance calculations is excellent as shown by the plotted points in Fig. I. The recovery as influenced by the geometry of the five-spot, taking into account the deviation from true radial character, has been predicted by Aronofsky and Ramey." The authors used a potentiometric model to determine the sweep efficiency to breakthrough of the flooding agent. Dyes, Caudle and Erickson,' by means of a laboratory model, show the influence of mobility ratio on recovery after breakthrough of the injected material. Dyes, et ai, used the X-ray shadowgraph technique and miscible phases in their studies. Craig, Geffen and Morse from the analysis of the performance of their laboratory model correlate many variables, including average saturation of a flooding agent determined by permeability relationship. Craig, et al, used X-ray shadowgraphs to determine the sweep efficiencies in sandstone models. In this report the sweep efficiency is not needed, and the important effect of continuously changing saturation gradients before and after breakthrough