Search Documents

By M. Sheffield

This paper presents a technique lor predicting the flow of oil, gas and water through a petroleum reservoir. Gravitational, viscous arid capillary lorces are considered, and all fluids are considered to be slightly compressible. Some theoretical work concerning the fluid flow in one-, two- and three-space dimensions is given along with example performance predictions in one- and two-space dimensions. INTRODUCTION Since the introduction of high speed computing equipment, one of the goals of reservoir engineering research has been to develop more accurate methods of describing fluid movement through underground reservoirs. Various mathematical methods have been developed or used by reservoir engineers to predict reservoir performance, 193-5,7,8 The work reported in this paper extends previously published work on three-phase fluid flow (1) by including a rigorous treatment of capillary forces and (2) by showing certain theoretical mathematical results proving that these equations can be approximated by certain numerical techniques and that a unique solution exists. DISCUSSION The method of predicting three-phase compressible fluid flow. in a reservoir can be summarized briefly by the following steps. 1. The reservoir, or a section of a reservoir, is characterized by a series of mesh points with varying rock and fluid properties simulated at each mesh point. 2. Three partial differential equations are written to describe the movement at any point in the reservoir of each of the three compressible fluids. All forces influencing movement are considered in the equations. 3. At each mesh point, the partial differential equations are replaced by a system of analogous difference equations. 4. A numerical technique is used to solve the resulting system of difference equations. Capillary forces have been included in two-phase flow calculations.1,4,8 The literature, however, does not contain examples of prediction techniques for three-phase flow that include capillary forces. Where capillary forces are considered, each of the three partial differential equations previously discussed has a different dependent variable, namely pressure in one of the three fluid phases. Therefore, three difference equations must be solved at each point in the reservoir. Where large systems of equations must be solved simultaneously, an engineer might question whether a unique solution to this system of equations actually exists and, if so, what numerical techniques may be used to obtain a good approximation to the solution. It is shown in the Appendix that a unique solution to the three-phase flow problem, as formulated, always exists. It is also shown that several methods may be used to obtain a good approximation to the solution. The partial differential equations and differenceequations used are shown in the Appendix. Matrix notation has been used in developing the mathematical results. Two sample problems were solved on a CDC 1604. They illustrate the type of problems that can be solved using a three-phase prediction technique. SAMPLE ONE-DIMENSIONAL RESERVOIR PERFORMANCE PREDICTION A hypothetical reservoir was studied to provide an example of a one-dimensional problem that can be solved. Of the several techniques available, the direct method of solution as shown in the Appendix was used. The reservoir section studied was a truncated, wedge-shaped section, 2,400 ft long, with a 6' dip. (A schematic is shown in Fig. 1.) This section was represented by 49 mesh points, uniformly spaced at 50-ft intervals. The upper end of the wedge was 2 ft wide, and the lower end was 6 ft wide. The section

Jan 1, 1970

POSITIONS VACANT (Under this heading will be published notes sent to the Secretary of the Institute by members or other persons.) Opening for man as head of laboratory of large brass manufacturing plant in Ohio; must also be competent to manage foundry making brass castings, following practice of malleable iron castings and insulating materials. Good opportunity for right man. No. 36. Professorship in Mining at a Western university; requires man competent to give instruction in ore dressing, in metal and coal mining, and in zinc smelting, and capable of gaining confidence of mining interests. Salary $2,200 to $3,000. No. 37. ENGINEERS AVAILABLE (Under this heading will be published notes sent to the Secretary of - the Institute by members or other persons introduced by members.) Member, graduate engineer, aged 44, with over 20 years' experience in connection with large coal and coke operations in different parts of the country, desires position as general superintendent of large coal company, or will take an interest with the management in a smaller company. No. 206.

Jan 6, 1915

By F. N. Rhines, A. H. Grobe

At elevated temperatures the oxide of silver is unstable in the air at atmospheric pressure, consequently no external oxide scale forms upon pure silver under conditions of high-temperature annealing. When small quantities of certain alloying elements are present in the silver, the formation of a thin external scale is possible1 and in addition there may form a subscale composed of the oxide of the solute element precipitated within the body of the silver. Norbury2 and Leroux and Raub3 have reported internal oxidation (subscaleiormation) in alloys of silver with 2, 7.5, and 30 per cent of copper. The presence of the subscale is believed to be responsible, at least in part, for the objectionable "fire mark" in Sterling silver.4 Several other alloys of silver, after oxidizing heat-treatments, are known to exhibit undesirable polishing characteristics that may be the result of internal oxidation. Except for the absence of an external scale of silver oxide, it is to be anticipated that silver alloys will prove to be very similar in their oxidation behavior to the alloys of copper, the oxidation characteristics of which have been studied in some detai1.5,6 The present research confirms this anticipation. The oxidation of a series of 20 dilute alloys of silver has been studied metallographically; some types of subscale not encountered among the copper alloys have been found. Instances of internal oxidation in alloys of most of the metals of the 1-6 and VIII groups of the periodic system are on record, but evidence of this type of oxidation in alloys of the metals of the intermediate groups is lacking. A number of the metals of the intermediate group, among them cadmium, lead, tin, and zinc, appear to provide the conditions essential to subscale formation; i.e., they form oxides with a relatively low negative free energy of formation, they dissolve other metals that form more stable oxides, and, presumably, oxygen will diffuse through them. In a study of 40 alloys of these white metals only a few cases of internal oxidation have been found. The probable reasons for this difference in behavior will be discussed presently. Experimental Procedure Silver.—The silver alloys employed in the oxidation studies were prepared in heats of 30 grams each from high-purity silver (99.993 per cent Ag)* and the purest avail-

Jan 1, 1942

By A. H. Grobe, F. N. Rhines

At elevated temperatures the oxide of silver is unstable in the air at atmospheric pressure, consequently no external oxide scale forms upon pure silver under conditions of high-temperature annealing. When small quantities of certain alloying elements are present in the silver, the formation of a thin external scale is possible1 and in addition there may form a subscale composed of the oxide of the solute element precipitated within the body of the silver. Norbury2 and Leroux and Raub3 have reported internal oxidation (subscaleiormation) in alloys of silver with 2, 7.5, and 30 per cent of copper. The presence of the subscale is believed to be responsible, at least in part, for the objectionable "fire mark" in Sterling silver.4 Several other alloys of silver, after oxidizing heat-treatments, are known to exhibit undesirable polishing characteristics that may be the result of internal oxidation. Except for the absence of an external scale of silver oxide, it is to be anticipated that silver alloys will prove to be very similar in their oxidation behavior to the alloys of copper, the oxidation characteristics of which have been studied in some detai1.5,6 The present research confirms this anticipation. The oxidation of a series of 20 dilute alloys of silver has been studied metallographically; some types of subscale not encountered among the copper alloys have been found. Instances of internal oxidation in alloys of most of the metals of the 1-6 and VIII groups of the periodic system are on record, but evidence of this type of oxidation in alloys of the metals of the intermediate groups is lacking. A number of the metals of the intermediate group, among them cadmium, lead, tin, and zinc, appear to provide the conditions essential to subscale formation; i.e., they form oxides with a relatively low negative free energy of formation, they dissolve other metals that form more stable oxides, and, presumably, oxygen will diffuse through them. In a study of 40 alloys of these white metals only a few cases of internal oxidation have been found. The probable reasons for this difference in behavior will be discussed presently. Experimental Procedure Silver.—The silver alloys employed in the oxidation studies were prepared in heats of 30 grams each from high-purity silver (99.993 per cent Ag)* and the purest avail-

Jan 1, 1942

By Homer Turner

THE nature and comparative features of anthracite and bituminous coals have been discussed by the writer in two previous papers.1 Although this paper is offered as a further contribution to the subject, it seems desirable at the outset to repeat some of the descriptive matter of the early works in order to show the actual connection between various constituents of both coals and the new photomicrographs and other physical and chemical properties.

Jan 1, 1929

By P. E. Baker

An experimental study of heat flow in steam flooding was carried out with steam displacing water in a plane-radial fluid-flow model. Temperature distributions in the model reservoir, overburden, and substratum were measured at frequent intervals during flooding by means of a fixed array of thermocouples. Heat losses to overburden and substratum and the division of heat between the steam and hot-water zones in the reservoir were calculated from the temperature data by performing numerical integrations over the appropriate volumes. Temperature profiles and heat distribution results are presented. Some results were (1) the fraction of injected heat that is lost to overburden and substratum, when expressed as a function of time, did not depend on injection rate; (2) a significant portion of the injected heat was contained in the hot-water zone ahead of and underlying the steam front in the reservoir. The division of heat between the steam and hot-water zones did depend on injection rate; and (3) considerable gravity override was observed. The experimental results are compared with those of several published theoretical studies of heat transfer in hot-fluid injection; approximate agreement on heat losses is found with some of the most idealized theories. A valid comparison cannot be made in regard to steam us hot-water zone; basic interpretation of the hot-water zone, in terms of a vertical steam front, is clouded by gravity override in the experiments. INTRODUCTION The experimental study reported here primarily concerns two problems of heat transfer in steam flooding that have economic and engineering importance. These are (1) the fraction of injected heat lost to the overburden and substratum, and (2) the division of the remaining heat between a zone containing steam and a zone of hot water at lower-than-boiling temperature. The volume of the steamed region, which has a direct effect on oil recovery, can be calculated readily from these results. Several valuable theoretical studies have been published on this problem, but no experimental tests of the theories have previously appeared. Tests were needed because the effects of the necessary simplifications in the theories could not be determined otherwise. The basic, untested theoretical problem is the effect of "decoupling" heat and mass transfer, and the resulting inability to include the effects of heat transfer and condensation on fluid flow. The theoretical problem solved so far has been a heat transfer problem in which fluid is present as a heat carrier and flows in a stipulated manner, either linear or plane-radial. The simplifications that have sometimes been made in regard to thermal conductivity by, for example, Marx and Langenheim,l Lauwerier2 and Thomas3 apparently do not have an important effect on heat loss calculations, although the results of the different methods disagree to a certain extent, as shown by Spillette4 and by Ramey.5 A recent paper by the author6 compared the results of experiments in hot-water flooding with several theoretical results and showed a rather poor agreement. It will be seen that the experimental heat loss results for steam flooding, however, show acceptable agreement. These comparisons are discussed in a later section. EXPERIMENTAL APPARATUS AND PROCEDURE The system studied was one in which a sand pack, initially water saturated and at a uniform temperature, was flooded with steam at a constant rate and temperature. The temperature distribution in the "reservoir" and adjacent "strata" was measured with an array of thermocouples, and heat losses were determined from these measurements. The sand pack contained no oil because the only concern of the experiments was heat flow, and it was important to have the flow pattern as simple

Jan 1, 1970

By Clarence Seyler

THE technical difficulties of cutting thin sections of coal for examination by transmitted light have hitherto restricted the investigation of the important subject of the microstructure of coal to the few possessed of the requisite skill and time. Apart from this, the method fails with anthracites and other coals rich in carbon, owing to the opacity of the material. Occasional use of the method of polishing a surface, with or without etching, has been made, but it was not until Winter1 published an account of his results that attention was seriously directed to the method. In May, 1923, a specimen of dull anthracite came into my hands; I tried Winter's method on it and the presence of megaspores and other vegetable structures became clearly visible. A preliminary notice was published by me;2 but I soon abandoned the use of Schulze's solution having found that a mixture of chromic and sulfuric acids gave far better results. In the meantime Turner and Randall, in America, had independently developed the method of flame etching, as applied to anthracite, with great success. I have spent a year studying the method of chromic etching, which gives excellent results with bituminous coals and cannels, as well as with anthracite. For anthracites, I have used the chromic method and Turner's method of flame etching, whichever was the more convenient or as a check one on the other. The polishing is done, according to metallurgical practice, by the use of a polishing machine and successively finer emery paper, followed by polishing on cloth moistened with a suspension of levigated green oxide of chromium, and a final rubbing on "selvyt" (a cotton velvet polishing cloth) moistened' with alcohol. The polished sample is examined by reflected light with a 2/3 and 1/6 in. objective. Bituminous coals usually show considerable structure, and even anthracites show it faintly in places. The sample is then etched by immersion in a boiling solution made by adding 10 c.c. of concentrated sulfuric acid to 30 c.c of a saturated solution of chromic acid; enough water is then added to dissolve

Jan 3, 1925

By Arthur Green

A GROUP of 13 cities situated in Los Angeles and Orange counties in Southern California is engaged in constructing an aqueduct to carry water from the Colorado River at a point near Parker, Arizona, to a distribution point sear Riverside, California. The aqueduct will have a total length of approximately 241 miles, and is designed to deliver 1605 cu. ft. of water per second, or a little more than one billion gallons per day. The aqueduct crosses a barren, desert waste and the project includes approximately 91 miles of 16-ft. tunnel. The large amount of tunnel driving that has been done in this country in the past few years has resulted in some remarkable developments in equipment and methods. The building of the Colorado River Aqueduct, with its 91 miles of tunnels, has called on the best engineering knowledge in the country. New and better equipment has been designed for all phases of the work, and this equipment has been coordinated and adapted to the various conditions with a high degree of efficiency. The result has been that records for tunnel driving have been set under adverse natural conditions that exceed the estimates for progress under good conditions. All of the tunnels in the main aqueduct, except the Valverde and Bernasconi at the western end, have a finished diameter, inside the concrete lining, of 16 ft. Fig. 4 shows a section of the tunnel with the concrete lining in place. Fig. 5 illustrates the finished tunnel in three typical sections. The standard tunnel grade is 3.432 ft. per mile in the 16-ft. tunnels, but the Valverde and Bernasconi tunnels have a steeper grade and a finished diameter of 15 ft. 3 in. The rough bore of the large main aqueduct tunnels averages 12.95 solid cubic yards of material per foot of advance in timbered sections and 10.76 solid cubic yards per foot in unsupported sections.

Jan 1, 1936

By M. L. Brashears

IN the last 15 years industrial use of ground water has more than doubled, and in 1951 amounted to 5 billion gallons per day. A similar sharp increase in the utilization of ground water for irrigation and public-water supply occurred in the same period. In many areas rapid increase in withdrawal from wells has taken place almost entirely unhampered by regulatory control and with little or no integration of effort. As might be expected, the chief interest in many regions has been maximum production rather than sustained perennial yield. As a result, widespread depletion of underground reservoirs and deterioration of the quality of the water stored in them has taken place in many areas, even though total pumpage in the United States is far below ultimate potential. Of even more concern is the fact that excessive withdrawal has drawn salt water into the reservoirs beneath many heavily populated centers along the Atlantic, Gulf, and Pacific coasts, causing costly abandonment of pumping plants.' Many hydrologists expect that consumption of water will rise rapidly in the near future, and some predict that industrial requirements will more than double in the next decade.2,3 Thus it appears likely that the draft on many already heavily pumped underground reservoirs will be greatly increased and the search for additional sources of usable ground water intensified in years to come. In view of this, industry as a whole will be forced more and more to recognize the potentialities and limitations of ground-water reservoirs and to utilize them more effectively to prevent costly water shortages and disruption of production. Through painful experience, some industries are already well aware of the need for effective water utilization, and have managed individually or through joint effort to check trends threatening to deplete underground reservoirs completely or to impair the quality of the water. Various remedial measures have been used to bring about successful management of local or regional ground-water resources. Of these, replen-ishment of aquifers by recharge wells or basins has played an important role in overcoming some ground water problems. Artificial recharge of underground reservoirs by water spreading has been practiced successfully in the United States for many years. In the West it has become an important method of salvaging flood run-off for irrigation of crops and maintenance of public water-supply reserves, and it is used to some extent in parts of the East. Artificial recharge by means of wells, on the other hand, is a relatively new development. Until recently it was employed in only a few areas, principally along the East coast. For the last few years, in the ever increasing search for additional water supplies, industry has had greater recourse to this method. Utilization of recharge wells to control the temperature and quality of underground water supplies is also being considered seriously. Operation of recharge wells, like water spreading, is governed largely by local conditions. It requires water relatively low in turbidity, whereas in some areas water spreading has been used successfully with water of high turbidity and silt content. However, water spreading must be employed in large areas and can be carried on effectively only where aquifers crop out at the surface. Recharge wells can be used in limited space. Recharge wells are similar to production wells except that the water flows in the opposite direction. Thus any water-bearing bed that will yield water to wells may be recharged by wells. Often, however, the water available for recharge is of a different character and temperature from that existing in the ground-water reservoir and if transmitted directly underground from a recharge well to a production well might require expensive or difficult treatment before it could be used. Fortunately the physical characteristics of reservoir beds, which control the movement and behavior of ground water, are generally not homogeneous. Moreover, the movement of ground water is very slow because of the frictional resistance of the reservoir beds. By taking full advantage of hydrologic and geologic conditions, it is therefore possible in many instances to bring about favorable changes of temperature and dilution as the water moves from the recharge wells underground to the production wells. Furthermore, if the natural quality or temperature of ground water is unfavorable for industrial purposes, recharge wells may be used to introduce water of more favorable quality or temperature into the ground-water reservoir. When water is discharged into a recharge well, the head in the well is increased. Because of this, a cone of elevation is produced on the water table or the artesian pressure surface in the area surrounding the well. The cone of elevation is similar to the cone of depression produced around a pumping well except that the apex of the cone is above the water table or artesian pressure surface. Thus if a recharge well and a production well tapping the same water-bearing bed are close together, as would be the case at many industrial plants, some of the water discharged from the recharge well would be drawn into the production well within a short time. Under such conditions it is apparent that water of unfavorable temperature and chemical characteristics should not be used for recharging. The more important ground-water reservoirs in the United States often consist of alternating layers of impermeable beds and porous material that will yield water readily to wells. Physical characteristics of individual beds in a ground-water reservoir may not persist over great distances, the impermeable layers grading into beds that will yield large quantities of water. Thus the water-yielding material in underground reservoirs, whether large or small,

Jan 1, 1954

By R. J. Wagner, F. F. Craig, H. G. Riley, J. D. Griffith

Peripheral water injection has been underway in the Sholem Alechem Fault Block "A" Unit, Stephens County, Okla., since 1955. In the engineering planning of the flood, it was recognized that maintenance of pressure in the large, dry gas caps was vitally necessary to prevent the loss of an estimated 27 million bbl of oil through saturation of the gas caps by the oil banks moving updip on the peripheral flood pattern. The development of an abundant, low-cost water supply and the rising cost of makeup-gas purchases indicated that gas-cap control by water injection was desirable from an economic standpoint, provided an eflective water barrier could be placed without premature water breakthrough into downdip oil producers. A laboratory model scaled to the existing reservoir conditions was comtructed to evaluate the feasibility of gas-cap water injection. Model studies indicated that (I) water injected at the edge of the gar-oil contact would form a continuous barrier between the oil and gas zones, and (2) water would advance uniformly downdip without segregation due to gravity effects provided a specified, minimum water-injection rate was maintained. The model test results also indicated that it was not necessary to fill the gas cap completely with water. Water injection at the edge of one of the gas caps was initiated in March, 1958, to field test gas-cap water injection. Water breakthrough occurred at a downstruc-ture producing well in April, 1959, after approximately 2.5 million bbl of water had been injected into the gas cap. This was the volume calculated to be required to form the barrier between the oil zone and gas cap. This performance indicates that prevention of oil encroachment into dry gas caps in a peripheral water flood by water injection at the gas-oil contact appears practical and that the model studies provided a valid prediction of field performance. INTRODUCTION The Sholem Alechem Fault Block "A" Sims Sand Unit was one of the fist major waterflood projects in Oklahoma to utilize a peripheral injection pattern. Design of control measures to minimize oil encroachment into the large crestal gas caps was one of the greatest single problems during engineering planning of this flood. Engineering studies indicated that saturation of the gas caps with oil displaced by the peripheral water flood could result in a loss of about 27 million bbl of oil, or over one-half the estimated waterflood recovery. This estimate was based on a residual oil saturation after water flooding of 35 per cent pore volume, compared to the negligible oil saturation initially in the gas caps. The Sholem Alechem Springer field, located in Stephens County, Okla., is an elongated anticline with structural dips varying from 400 to 1,000 ft/mile. The structure contains numerous faults, some of which are sealing. The portion of the field considered in this paper is Fault Block "A", which is the largest of the five major fault segments in the field. A map of the unit is shown on Fig. 1. The unit is bounded on the east and west by major sealing faults, on the north by a water-oil contact and to the south by a deterioration of sand development. Two domes, each having gas caps, exist in the unit and are separated by a north-south fault as shown in Fig. 1. The sands being flooded are the Sims sand members of the Springer series of early Pennsylvanian age. The Sims sands average over 100 ft in thickness in Fault Block "A" and have a porosity of approximately 19 per cent. Core permeabilities average 160 md, with interstitial water saturation of about 23 per cent pore volume. Three different zones can be correlated over most of the area. However, they are commingled in most of the producing wells. The two upper zones, designated as First and Second Sims, are the most important. Both the First and Second Sims have gas caps in the east dome, while only the First Sims has a gas cap in the west dome. These gas caps, whose limits are also shown on Fig. 1, initially occupied nearly 11 per cent of the total reservoir volume. The major concern of this paper is the control of the two gas caps in the First Sims sand. Fig. 2 is a cross section, to scale, across the east dome of this formation. Notable features are the relatively shallow dip, the deterioration of the First Sims sand on the south side of the gas cap and the water-oil contact on the northern edge. The formation dip at the gas-oil contact in this cross section is approximately 5".

By H. O. Jahns

This paper describes the application of regression analysis for obtaining a two-dimensional areal description of heterogeneous reservoirs from short-term pressure-time data such as that obtained in interference tests. The method replaces the time-consuming trial-and-error procedure commonly used to match field data on an electric analyzer or digital computer with a systematic search which is programmed for a computer. The computer program adjusts the properties of a reservoir model automatically until a least-squares fit is obtained between observed and calculated pressure data. The reservoir is simulated by a single-phase, compressible, two-dimensional model. It is divided into a number of homogeneous blocks whose transmissibility (kh/F) and storage (Fch) values are varied to obtain the least-squares fit. The reliability of these values is determined from their standard deviations and correlation coefficients. Although the method is rigorously applicable to single-phase flow only, multiphase flow can be handled provided saturation changes are small during the test. Possibly the method can also be used to obtain a reservoir description from pressure-production history, but this application is outside the scope of this work. The paper includes, in addition to a description of the numerical procedure, a discussion of some of the problems associated with the method. Rules are given to help in selecting the number of homogeneous blocks and deciding upon their arrangement. The uniqueness of a reservoir description is considered. Finally, the use of the method is illustrated by the interpretation of field data from two interference tests. INTRODUCTION Pressure data from short-term transient tests, such as single-well and interference tests, are widely used to obtain reservoir properties. These tests are usually analyzed by assuming a simple reservoir model; very often, a homogeneous one is used. As a result, analysis of the transient data from each well frequently gives different values for reservoir properties. The problem then arises to combine all these differing results into a more detailed picture of the reservoir. One technique is to simulate the reservoir with a digital computer or with an electrical analyzer and to adjust the reservoir parameters by trial and error until the simulated pressure data are in reasonable agreement with the observed pressure data for all wells. Although this method has been used for both transient tests and pressure-history data, it is time-consuming and subjective. A second technique uses regression analysis to replace the trial-and-error procedure with a systematic search that can be programmed for a digital computer. Use of regression analysis in reservoir description was proposed recently by Jacquard and Jain.1 They divided the reservoir into a number of homogeneous blocks whose properties are varied until a least-squares fit is obtained between observed and calculated pressures. However, they did not consider their technique to be operational, mainly because of ".. .(I) the lack of experience in using the method .. . notably for the improvement of convergence; and (2) limitations imposed by the insufficiency of available computers". 1 While the analysis presented in this paper applies the same general principle used by Jacquard and Jain, the specific method is significantly different. Some differences are (1) the regression problem is solved in a different way which requires less computer time in most cases; (2) a stepwise solution, in which the detail in the reservoir description is increased from step to step, is used to improve convergence; and (3) the reliability of the estimated reservoir .properties.. as measured by their standard deviation and correlation coefficient, is estimated'

Jan 1, 1967

By B. H. Caudle, W. J. Bernard

Factors which influence the success or failure of a waterflood can seldom be determined in the laboratory. For this reason pilot waterfloods are initiated in a repreventative portion of the oil reservoir in question. For a pilot flood to predict quantitatively the recovery to be expected in a field-wide waterflood operation, the pilot area must behave as though it were confined (surrounded by similar areas). In this study, laboratory fluid-flow models were used to determine the simplest pilot pattern, for particular conditions of mobility ratio and initial gas saturation, that would behave as though it were confined. Pilot patterns studied ranged in complexity from a single inverted five-spot to a grouping of nine regular five-spots. Only the innermost producing well in each pattern was studied. Model results showed that the optimum number of wells in the pattern depends upon the oil-water mobility ratio and the expected oil-bank size. Unfavorable mobility ratios will, in general, require more wells in the pilot pattern than will favorable mobility ratios. Pilot patterns in reservoirs which contain a dispersed, flowable, free gas saturation will require fewer wells than for the under-saturated case. The single inverted five-spot pattern was found to be unsatisfactory for predicting behavior of fully developed waterfloods. In particular, it is possible that, in reservoirs which contain a flowable, dispersed gas phase, the oil bank will never be observed at the producers due to the large amounts of free gas which continue to be produced with the oil. INTRODUCTION One method which has been used to predict the performance of a waterflood is the pilot flood. The pilot waterflood is a flood which involves only a small cluster of the reservoir wells and is located in a small, representative portion of the reservoir. The object is that oil produced from the pilot can, in some way, be related to the oil recovery to be expected from a field-wide expansion of the waterflood. However, these field pilot waterfloods have often been unreliable in the prediction of oil recovery in a fully developed waterflood. This unreliability has also been demonstrated in several laboratory studies of pilot floods. Some of the investigators have shown that there are situations in which the pilot flood oil production is far too optimistic with respect to the oil recovery in the fully de- veloped flood. Others4-G have shown that the pilot results can also be pessimistic, especially if the pilot waterflood is initiated in an oil reservoir which has been depleted by primary recovery and is at very low pressure. The major reason for this unreliability of pilot water-floods is the migration of fluids into or out of the pilot area. By the well-known method of images, if straight lines can be drawn to represent vertical planes of symmetry in a porous medium which contains pressure sources and sinks (injectors and producers), then these lines are invariant streamlines, or lines across which there is no potential gradient, and therefore no flow. In an actual reservoir, these lines of symmetry can never be established exactly because of reservoir inhomogeneities and irregular reservoir boundaries. However, if the reservoir is relatively large and contains wells in repetitive patterns, these lines of symmetry are commonly assumed to exist for the pattern units sufficiently far removed from the reservoir boundary. Lines of symmetry for the five-spot injection pattern are shown in Fig. 1. Each five-spot unit in this figure can be considered confined with respect to flow across its boundary. In pilot floods this is not the case. The lines of symmetry for the pilot patterns investigated in this study are shown in Fig. 2. It is obvious that the fluid within these pilots is not confined and is therefore able to migrate into or out of the pilot area. Intuitively, one can see that, if more wells are added to the pilot, the innermost unit tends to behave more and more like the confined pattern. However, there is a practical limit to the number of wells which should be placed in the pilot. This limit is usually determined by economic factors. It was the purpose of this study to use laboratory fluid-flow models to determine which of the previously mentioned pilot patterns will force the innermost producing well to behave as it would in a fully developed waterflood. Since fluid migration is influenced by initial saturation conditions and the mobility ratio, these factors were included in the study. The ultimate objective of this study was to develop data which would allow the operator to choose a pilot pattern and operating conditions that will yield a production history which can be applied directly as an estimate of the performance of each production well of the fully developed waterflood. BASIS FOR THE STUDY The basic problem of field pilot floods is the migration of the reservoir fluids into or out of the pilot area. This problem has been the subject of previously reported model studies on pilot floods. These studies have been concerned mainly with the development of arbitrary "correction factors" to be applied to the simple, unconfined pilot systems such as the single five-spot. The correction factors were intended to adjust the production history of the un-

The Library of the above-named Societies is open frond 9 A.M. to 10 P.M., except on holidays. It contains about 70,000 volumes and 90,000 pamphlets, including sets of technical periodicals and the publications of scientific and technical societies. Members of the Institute, with few exceptions, are forced to spend a portion of their time in localities isolated front sources of information. To these the Library, through Its Library Service Bureau, can render valuable service through correspondence; letters requesting information will receive especial attention. The Library is prepared to furnish references and photographic copies of articles on mining and metallurgical subjects; to determine the existence of mining neaps, and to furnish general information on the geology and mineral resources of all countries. All communications should he made as definite as possible so that the information received may be what is desired and not include collateral matter which may not be of interest. The time spent in searching for such collateral matter will be saved, and the information will be, sent more promptly and in more usable shape.

Jan 1, 1917

By H. A. Lesser, H. L. Stone, G. H. Bruce

ESSO PRODUCTION RESEARCH CO., HOUSTON, TEX. A mathematical model that represents the conduction heating of a rock formation of limited permeability is formulated. Heat is introduced by the injection of a hot condensing gas into horizontal fractures through the formation. The fractures are assumed to be equally spaced; the flow of fluid is linear. The model consists of three coupled partial differential equations and two functional relationships. Temperature histories for both fractures and formation are obtained by numerical procedures. Application of the model in predicting in situ heating of oil shale is illustrated with four sample cases. Influences of shale thermal diffusivity, fluid pressure in the fracture and fracture spacing upon heating rate are demonstrated. In all cases considered, several years were required to heat all the formation to 700F with 1,000F steam. Less than 32 per cent of the injected heat was utilized for heating the shale. To facilitate extensive calculations of this type, a set of seven dimension-less quantities are derived. With this dimensionless system, the effect of controllable parameters on economics can be explored. INTRODUCTION Vast deposits of hydrocarbons are contained in oil shale formations in the western United States and in other countries. These hydrocarbons, known as kerogen, are difficult to recover because they are immobile solids and because the rock formations in which they are found have very limited permeability. The kerogen can be converted through pyrolysis to recoverable liquids and gases if the formation temperature is elevated to the range 600 to 800F. To accomplish this heating, the formation might be artificially fractured and a hot fluid injected into the fractures. If the injection and production wells are vertical and the fractures horizontal, interwell communication will be established and the hot fluid will flow from an injection well to a production well. In this process, some of the heat in the fluid will be given up and conducted into the formation. Products of pyrolysis will flow from the heated rock matrix into the fractures and be produced, commingled with the injected fluid. The practicality of such a process is strongly influenced by the rate at which heat can be transferred from the hot fluids to the formation. An analysis of this heat conduction problem is the subject of this paper. Linear flow of a condensing gas through a set of horizontal, equally spaced fractures is considered. Temperature histories for both fractures and formation are obtained from numerical procedures. The analysis also applies to formations containing streaks of natural permeability instead of fractures. Heat flow problems similar to the one treated here have been studied by other workers. Thomas1 has solved a radial problem in which he considers the injection of a noncondensing gas into a single, horizontal fracture and approximates the temperature profile in the fracture by a step function. The solutions to other problems, which are less similar to the present problem than the one treated by Thomas, have been reviewed by Spillette.2 In addition to these mathematical studies, a field test of an in situ shale oil recovery process using hot natural gas as the injected fluid has been reported.3 The first part of this paper is concerned with the mathematical development of the problem. The differential equations, boundary conditions and functional relationships which describe the physical model are formulated and the equations are cast into finite difference form. The method of solution of the difference equations is discussed (Appendix B). The second part of this paper presents results computed for a typical oil shale formation for which steam at 1,000F is the injected fluid. Fracture spacing, horizontal and vertical thermal diffusivities

Jan 1, 1967

By A. J. Ardell, O. D. Sherby

The elevated-temperature steady-state creep behavior of polycrystalline a Zr was studied in vacuo under constant tensile stress. The experiments were conducted from 660° to 845°C over the stress range 870 to 5400 psi. The steady-state creep rate, is, of a Zr was found to be related to the applied stress, a, and temperature, T, by the equation where C is a constant. Two interpretations of the data are considered, each requiring different values of the parameters ß, Q, and . It is shown that this formula can be derived using Gilman's phenomeno logical equation relating dislocation velocity to the shear stress for the case where dislocation mobility is limited by the nucleation of kink motion. From the fundamental dependence of strain rate on dislocation velocity and the density of mobile dislocations, the factor exp(ßo) can be related to the stress dependence of the dislocation density, and Q to the energy of kink formation. The two interpretations considered are as follows: Case I—the parameters are adjusted so that the equation above describes the data over the entire range of stress and temperature, implying that creep is controlled exclusively by dislocation glide; Case IIthe parameters are adjusted to allow for the possibility that creep is diffusion-controlled at high stresses but glide-controlled at low stresses. For both cases the values of the parameters are within the spread permitted by the data. Factors favorable to the observation of the unusual stress -dependence of i, for a Zr are discussed and it is shown that the anomalously rapid rates of self-diffusion and the low modulus on the principal glide system are primarily responsible for the anomalous creep behavior. Metal-lograflhic observations are presented and discussed with reference to the overall creep behavior. PREVIOUS investigations of the creep behavior of a Zr have been conducted from room temperature to about 500C The results of these investigations are primarily of engineering significance because of the variety of conditions under which they were carried out, and the low testing temperatures employed. Little information is available concerning the fundamental nature of the elevated-temperature creep behavior of a Zr. This investigation of the creep behavior of a Zr above half the melting temperature (the melting temperature, T,, of a Zr is estimated to be about 1860"K was initiated because there are two significant anomalies of the physical properties of a Zr: a) the rate of self-diffusion in polycrystalline a Zr is anomalously rapid: b) both Young's modulus for polycrystalline a Zr and the shear modylus as-sociated with the prismatic slip system { 1010)( 1120) (the principal slip system in a Zr 7) have anomalously large temperature variations The importance of these quantities to the understanding of creep behavior has been demonstrated, and it was hoped that this investigation would provide a critical test of the predictions based upon the analysis of creep data for other pure metals. However, as the investigation progressed it became increasingly clear that the unusual creep behavior of a Zr could not be analyzed in a conventional manner. 1) EXPERIMENTAL EQUIPMENT AND PROCEDURES Two grades of zirconium, designated A and B, were used in this investigation. The analyses, as supplied by the vendors, are presented in Table I. Material A, with a nominal purity of 99.84 pct, was supplied by the Wah Chang Corp. in the form of cold-rolled sheet from which specimens were machined without further fabrication. Material B, with a nominal purity of 99.95 pct, was purchased from the Nuclear Materials and Equipment Corp. in the form of iodide crystal bar.

Jan 1, 1968

By Thomas T., Read

IT will be recalled that the first professor of metallurgy in the United States, appointed in 1855, never really gave any instruction in metallurgy and gradually turned into a professor of mineralogy. Metallurgy was included with. mining in the early curricula, two students who took the B.S. in mining engineering from the Polytechnic College of Pennsylvania in 1862 wrote their graduation theses on the metallurgy of copper and zinc, respectively. Even after Columbia School of Mines established, in 1868, a separate curriculum leading to the Ph.B. degree in metallurgy, very few students elected it,* preferring to take the mining degree, since they received in its curriculum nearly all the metallurgical instruction they would have in the metallurgical curriculum. For many years a large proportion of the men who attained prominence in metallurgical enterprises had taken mining rather than metallurgical degrees. METALLURGY Clifton Sorby's study of thin sections of rocks under the micro¬scope, beginning in 1850, and its subsequent development by Zirkel and others, had made microscopic petrography a recog¬nized department of geological research as early as 1873. . It was inevitable that the idea of studying polished metal surfaces-¬under the microscope by reflected instead of transmitted light would occur to investigators. Sorby himself did it in 1864, and Adolph Martens and Floris Osmond may -be especially named among the many European workers who, by the late eighties, had so developed the, necessary techniques that they were able to photograph polished metal surfaces at magnifications as high

Jan 1, 1941

By W. B. Mather

Arkansas' nonmetalliferous deposits are of many different types. The principal types briefly described in this paper are: clays, shales and slates; silica deposits; limestone and dolomite; barite; gypsum; celestite; nepheline syenite; greensand; phosphate; vermiculite; fuller's earth; bentonite; diamonds; sand and gravel. The future development of these deposits is expected to effect profoundly the industrial development of the state.

Jan 5, 1950

By C. H. Chou, Albert Sauveur

The following experiments were conducted in order to study the influence of the rate of cooling through and below solidification on the dendritic structure (macrostructure) and microstructure of some hypo-eutectoid steel. The steel selected was of commercial grade and had the following composition: C, 0.40 per cent.; Mn, 0.85; Si, 0.70; P, 0.099; S, 0.043. After melting in vacuum in an Arsem furnace, the various samples were allowed to solidify and to cool to room temperature as indicated in Fig. I. Samples were thus obtained representing all possible combinations of slow and fast cooling through the various ranges of temperatures indicated, namely: (I) zone of solidification or primary crystallization, (2) zone of granulation, and (3) zone of secondary crystallization. The location of the different temperature zones in the iron-carbon equilibrium diagram is shown in Fig. 2. The authors adopted Belaiew's notations. The designations of the various samples readily indicate their treatment, as will be apparent from Table 1. Sample SSS.—This sample was slowly solidified and cooled to room temperature in the crucible in which it was melted. The dendritic (macro) structure is shown in Fig. 3. The dendrites an: sharp and dis-

Jan 1, 1930

By J. Campbell

Various models are discussed for the evaluation of the negative pressures which may occur in solidifying materials which exhibit various deformation modes: elastic-plastic, Bingham, viscous, or creep flow. The inadequacy of the previously proposed elastic-plastic solution for solidifying metals is revealed by comparison with the more reliable creep results which are given graphically for aluminum, copper, nickel, and iron. The maximum tensions experienced in the liquid phase of solidifying spheres ranging in size from large castings to submicron powders are in the range from —10' to —105 atm for these metals. THERE has been much recent interest in the negative pressures associated with the volume change on solidification and in the possibility of the occurrence of cavitation. Considering the freezing of a highly supercooled liquid, an attempt to evaluate the stresses in the liquid ahead of the rapidly moving solidification front has been made by Horvay1 on a microscale and by Glicksman2 on a macroscale. In a casting of a wide freezing range alloy, the pressure differential due to viscous flow of residual liquid through the pasty zone has been discussed by Piwonka and Flemings,3 In a previous publication4 the author has attempted to estimate the negative pressure occurring in the residual liquid of a spherical casting, employing an elastic-plastic model to describe the collapse of the solidified shell under the internal tension. An earlier model assuming a rigid shell was shown to be inaccurate by many orders of magnitude. The elastic-plastic model is critically reviewed here, and other models are developed which are thought to be more closely related to metals and other materials near their melting points. The spherical geometry (Fig. 1) is chosen because the highest shrinkage pressures would be developed, although the analyses are readily adaptable to cylindrical geometry. A parallel sided casting experiences little internal tension because of the relatively easy dishing inward of the sides. (This commonly observed phenomenon has previously been attributed solely to atmospheric pressure.) Furthermore, small regions of confined liquid in a large solidified volume of a casting approximate reasonably well to spherical geometry. ELASTIC-PLASTIC MODEL The author has shown4 that as solidification proceeds the internal hydrostatic tension builds up until the elastic limit of the shell is exceeded. At this point the internal pressure is closely -2Y/3. Subsequently a plastic zone spreads from the inner surface toward the outer surface of the shell. When the whole casting is deforming plastically a rather more generalized analysis taking account of the externally applied pressure PA + 2y/b gives the internal pressure as: P = Pa + 2y/a + 2ys/b - 2 Y In(b/a) [1] The 2y/a and 2ys/b terms result from the tendency of the liquid-solid and solid-vapor interfaces to shrink, reducing their energy, and thereby helping to collapse the solid phase and compress the liquid phase. The 2y/b term would be important only for powders. The last term arises because of the plastic restraint of the solid, resisting collapse and so effectively expanding the residual liquid. From Eq. [I] it is easily shown that there is a minimum in the pressure at the radius amia= y/Y [2] which is of the order of 103K for the metals aluminum, copper, and iron, and corresponds to the minimum pressure Pmin = 2 Y[l-ln(bY/y [3] The results of a fully worked out elastic-plastic solution are given in a previous reporL4 The main criticism which may be leveled at this analysis when applied to metals at their melting points is the strong dependence of the yield stress on the strain rate. The strain rate varies with both solidification conditions (e.g., whether chill-cast or slowly cooled) and during solidification, as is indicated in the following section. Thus an appropriate choice of Y is very arbitrary. Before proceeding to a discussion of models which are strain-rate-dependent, it is necessary to evaluate the strain rate as a function of the rate of solidification. SOLIDIFICATION RATE Various empirical relations have been deduced5 for the rate of thickening of the solid shell by pour- out tests on partially solidified spheres. These, however, are unsatisfactory for our purposes since they become very inaccurate when the liquid core is very small. A theoretical approach is therefore necessary, and some solutions are set out below. Making the assumptions of constant surface temperature of the casting during freezing, no superheat, and a material freezing at a single temperature, Adams8 deduces the approximate solution: which becomes when b » a: Employing a semiempirical approach vallet6 finds:

Jan 1, 1969

J. B. MCGRAW, New York, N. Y. (written discussion*).-In Mr. Burke Baker's description of the process of the American Briquet Co., he speaks only of the attractive features, but every process which uses a compound or emulsified binder has certain undesirable features which must be considered. In making an emulsion, the success of which depends upon the introduction of the various ingredients at certain fixed temperatures and the complete intermixing of a soluble solution with an oil which has absolutely no affinity for it, considerable difficulty is found in keeping the various factors in strict accordance with the formula. To serve its purpose properly or, in other words, to be right for easy mixing with the coal, it should be of a certain predetermined density. Among the controlling factors of this ate the proportion of moisture contained and the temperature of the oil when introduced. If the latter is higher than intended, a greater proportion of moisture will be evaporated than is proper, and if the oil is not sufficiently hot the process of emulsification is retarded. These various factors can be controlled with more or less certainty where working in a laboratory or in a small experimental way, but when it is attempted on a large commercial scale, certain obstacles not otherwise apparent are likely to be encountered.

Jan 5, 1918