Search Documents

By M. A. Grossman

The harden ability of most steels can be predicted within 10 to I5 per cent provided the complete chemical composition is known, including "incidental" elements; and provided the as-quenched grain size is known; and provided, finally, that the composition and heating temperatures for hardening are such as to result in austenite free from carbide particles. In the method proposed herein, the steel is considered as having a base hardenability due to its carbon content alone (the hardenability of a "pure steel" of the given carbon content, without any other elements), and this base hardenability is multiplied by a multiplying factor for each chemical element present. After multiplying all these together, the final product is the hardenability. Grain size may be taken into account either in the base hardenability or after the multiplication. Hardenability is stated in terms of "ideal critical diameter", namely, the diameter of bar, in inches, that will just harden all the way through (absence of un-hardened core) in an "ideal" quench, and the calculation may also be related to the Jominy test The data bring out certain features rather clearly. For example: i. It is quite useless to attempt to predict hardenability unless all elements, including "incidentals,,) are known; thus an "incidental,, chromium content of 0.20 per cent increases the hardenability by about 50 per cent. 2. The relative effectiveness of different alloys is sometimes unexpected; thus molybdenum when calculated in this way appears to be of the same order of effectiveness as manganese, rather than much more powerful, as would appear to be common experience; observe that an increase from no manaanese to 0.20 per cent Mn provides a multiplying factor of 1.67, and an increase from no molybdenum to 0.20 per cent Mo provides a factor of 1.63, an increase of over 60 per cent in each case. However, if to a steel containing 0.50 per cent Mn there is added "20 points of manganese," the factor is raised from 2.65 (for 0.50 per cent Mn) to 3.35 (for 0.70 per cent Mn), or an increase of only 26 per cent. Thus the fist small addition of an element has a much more powerful percentage effect than an equal further addition when some is already present, and in most cases the effect of molybdenum is considered in relation to a steel in which molybdenum is absent. 3. If two elements are equally effective, a greater hardenability will be obtained by using, for example, 0.5 per cent of each than by using I.o per cent of either of them alone. 4. A knowledge of the as-quenched grain size is essential for precise work, since a difference of only one grain size number (say No. 7 instead of No. 6) makes a difference of almost 10 per cent in hardenability (in these units); this, however, does not apply to certain It should be emphasized that in chromium steels (over 0.30 per cent Cr) and chrome-molybdenum and chrome-vanadium steels, un-dissolved carbides are likely to be present in the steel as quenched, and that in such cases the charts can indicate only a maximum Passible hardenability, whereas the extent of hardening actually obtained may be much less. Thus tests on a number of chrome-molybdenum steels have indicated a degree of hardening amounting to only 50 to 65 per cent of the maximum possible, and in chromium steels from full hardening (100 per cent) down to as low as 70 per cent. On the other hand, when the amount of such elements is small (Cr under 0.30 per cent, Mo up to 0.25 per cent in the

Jan 1, 1942

By William Hurst, A. F. van Everdingen

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By Nicholas J. Grant, H. C. Chang

Microscopic observations during creep tests were made on AI-20 pet Zn, 80 pet Ni-20 pet Cr, and 25 and 3S aluminum specimens. All these materials failed in an inter-crystalline manner under certain stress and temperature conditions. Particular atten-tion was given to the initiation and propagation of intercrystalline cracks in coarsegrained AI-20 pct Zn specimens. The geometry of the grains and grain boundaries in these specimens was known before the tests. On the basis of the observations and a knowledge of the interaction between grains and grain boundaries, a mechanism of intercrystalline fracture and propagation of the fracture has been proposed. METAL, and alloy fracture occurs essentially in two ways: transcrystalline and intemrystalline. In the past most of the work done in this field was centered on the determination of stress criteria for fracture. The specimens used for these studies were tested under varying conditions of stress state, temperature, and strain rate. Unfortunately, very little attention was given to the initiation and propagation of the fracture, whereas the nature of fracture was examined closely only after failure. Several extensive reviews of these studies have been published by Orowan, Smith, Barrett, and Petch.1, 4 It is well known that a material which exhibits a ductile and transcrystalline fracture can be made to exhibit a brittle and intercrystalline fracture by modifying its composition and by heat treatment. This modification may or may not result in microscopically observable precipitation in the grain boundaries." In the creep of any one alloy the transition of transcrystalline to intercrystalline fracture is a function of stress, strain rate, and temperature." The tendency for intercrystalline fracture to occur increases as the strain rate decreases and the testing temperature increases. Extensive work on high purity aluminum (99.995 pet) under creep conditions showed the importance of the interaction between grains and grain boundaries, and the different responses of grain and grain boundaries to the variables: stress, strain rate, and temperature.7, 8 It was also shown that grain boundary deformation necessarily has to be accommodated in the grains by various deformation means as. for example, fold formation. The type of fracture resulting from the creep of high purity aluminum was invariably found to be transcrystalline. On the other hand, 2S and 3s aluminum were, under certain conditions, found to fracture in an intercrystalline manner." One reason given for this was that the grains could accommodate to a lesser extent the deformation of the grain boundaries.19 It was decided to examine this line of reasoning in greater detail and to investigate the particular interactions between grains and grain boundaries that result in the incidence and progression of inter-ct,ystalline fracture. It is also the purpose of this paper to show where and how intercrystalline cracks start and how they propagate to result in fracture. Several alloys—Al-20 pct Zn, 80 pet Ni-20 pet Cr, 2S and 3s aluminum-— were used for this study, but attention was paid mainly to the behavior of an A1-20 pet Zn alloy, since its general creep behavior has been established." From the results of this study a mechanism of intercrystalline fracture is proposed. Experimental Procedure The experimental procedure was similar to that used in previous work and has been detailed elsewhere.' Two parallel flat surfaces were milled from an originally round specimen, giving a gage portion which was 1x0.09X0.17 in. Some round specimens, 1/4 in. round by 1 in., were also tested. The A1-20 pet Zn specimens. which were annealed at 1030°F for 24 hr, showed two to four grains across the width of the test bar, and most of the grains occupied the whole thickness of the specimenu. Before some of the specimens were subjected to creep, the grain arrangement was mapped out in a development drawing, Fig. 1. The initiation and

Jan 1, 1957

By A. F. van Everdingen, William Hurst

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By M. A. Grossman

The harden ability of most steels can be predicted within 10 to I5 per cent provided the complete chemical composition is known, including "incidental" elements; and provided the as-quenched grain size is known; and provided, finally, that the composition and heating temperatures for hardening are such as to result in austenite free from carbide particles. In the method proposed herein, the steel is considered as having a base hardenability due to its carbon content alone (the hardenability of a "pure steel" of the given carbon content, without any other elements), and this base hardenability is multiplied by a multiplying factor for each chemical element present. After multiplying all these together, the final product is the hardenability. Grain size may be taken into account either in the base hardenability or after the multiplication. Hardenability is stated in terms of "ideal critical diameter", namely, the diameter of bar, in inches, that will just harden all the way through (absence of un-hardened core) in an "ideal" quench, and the calculation may also be related to the Jominy test The data bring out certain features rather clearly. For example: i. It is quite useless to attempt to predict hardenability unless all elements, including "incidentals,,) are known; thus an "incidental,, chromium content of 0.20 per cent increases the hardenability by about 50 per cent. 2. The relative effectiveness of different alloys is sometimes unexpected; thus molybdenum when calculated in this way appears to be of the same order of effectiveness as manganese, rather than much more powerful, as would appear to be common experience; observe that an increase from no manaanese to 0.20 per cent Mn provides a multiplying factor of 1.67, and an increase from no molybdenum to 0.20 per cent Mo provides a factor of 1.63, an increase of over 60 per cent in each case. However, if to a steel containing 0.50 per cent Mn there is added "20 points of manganese," the factor is raised from 2.65 (for 0.50 per cent Mn) to 3.35 (for 0.70 per cent Mn), or an increase of only 26 per cent. Thus the fist small addition of an element has a much more powerful percentage effect than an equal further addition when some is already present, and in most cases the effect of molybdenum is considered in relation to a steel in which molybdenum is absent. 3. If two elements are equally effective, a greater hardenability will be obtained by using, for example, 0.5 per cent of each than by using I.o per cent of either of them alone. 4. A knowledge of the as-quenched grain size is essential for precise work, since a difference of only one grain size number (say No. 7 instead of No. 6) makes a difference of almost 10 per cent in hardenability (in these units); this, however, does not apply to certain It should be emphasized that in chromium steels (over 0.30 per cent Cr) and chrome-molybdenum and chrome-vanadium steels, un-dissolved carbides are likely to be present in the steel as quenched, and that in such cases the charts can indicate only a maximum Passible hardenability, whereas the extent of hardening actually obtained may be much less. Thus tests on a number of chrome-molybdenum steels have indicated a degree of hardening amounting to only 50 to 65 per cent of the maximum possible, and in chromium steels from full hardening (100 per cent) down to as low as 70 per cent. On the other hand, when the amount of such elements is small (Cr under 0.30 per cent, Mo up to 0.25 per cent in the

Jan 1, 1942

By K. H. Coats

This paper describes an approximate technique for handling the problem of percolation of evolved gas upwards through the oil column in computer simulation of natural depletion. This technique has been incorporated into a general model for simulating three-phase flow in one, two or three dimensions. The mathematical model for perfonning these calculations is described in detail in a separate article. 1 While vertical gas percolation occurs during natural depletion in a reservoir of any configuration, it is especially pronounced in the pinnacle reel or bioherm. The reel may have an areal extent less than one section with a thickness of up to 800 It. The upward percolation of evolved gas often limits the time step in calculations as little as one day. Using this small time step, computer expense for a single 30- to 40-year simulation on even a one-dimensional basis has exceeded 92,000. The method described here for handling the percolation allows time steps of 60 days or more (depending on reservoir size and production rate), resulting in a considerable reduction in computing expense. The method also allows calculations in which secondary gas caps build up in tight zones of the oil column below the main gas cap. The validity of the method is indicated in connection with an example pinnacle reef field. Calculated results using a 2-day time step (where the method in question is not needed and not invoked) are compared with results using the method and a 60-day time step. The comparison shows good agreement. Results from one- and two-dimensional simulations of the reef are presented along with the corresponding computing times on the CDC 6600 computer. A three-dimensional simulation was also performed and the required computing time is given. INTRODUCTION A general analysis for simulating three-dimensional, three-phase flow in reservoirs has been developed. 1 In applying this model to reservoirs undergoing natural depletion, a time-step restriction was encountered due to the flow of evolved gas upwards through the oil column and toward the gas cap. A remedy for this problem has been incorporated in the analysis and is described here. The time-step restriction is encountered only in calculations which include flow in the vertical or near vertical direction. The restriction occurs to some extent in reservoirs of any geometrical configuration, but it is especially pronounced in the pinnacle reef or bioherm where the ratio of thickness to areal extent is unusually large. This paper describes the problem, a method of handling it, incorporation of the method in the three-dimensional, three-phase model, and, finally, a test of the method's validity in an application to an example pinnacle reef reservoir. THE PROBLEM During the early stages of reservoir depletion, pressure falls below bubble point in progressively lower regions of the oil column. Continuing production results in evolution of dissoived gas throughout the oil column. This gas then percolates upwards toward the top of the reservoir. If the gas encounters a sufficiently low permeability zone in its travel upward, then it can accumulate and form a secondary gas cap. If no zones are encountered which are tight enough to hold gas against the gravity forces, then the gas travels on until it reaches the top of the sand or the main gas cap. Numerical simulation of multiphase flow in reservoirs is generally performed with time steps such that the flow of a fluid into or out of a block in one time step is a fraction (considerably less than one) of the total amount of that fluid present in the block. Use of a time step where the time step's

Jan 1, 1969

By G. R. Hodgson, M. P. Lebourg

The introduction in the field of a new type well completion called for the setting of tubing open-ended in the well before perforating the casing. This paper describes a new perforating tool of the shaped charge expendable type, small enough to pass through tubing and powerful enough to perforate the casing. The perforations are produced at an angle and are large enough for optimum production. Performance of this tool in targets and under various well conditions is described. Very often no lifting equipment will be available at the well. In some cases, pressure will exist in the well at the time of the perforating operation (workover wells or multiple trip perforation) ; therefore, new equipment and new technique for the perforating operation were devised. This paper describes this equipment and technique, particularly casing collar recording, sheave support, and cable pull-down device. this last item being necessary for wells under pressure. Results obtained with this tool in the field are given. INTRODUCTION One of the key requirements for the introduction in the field of a new method of permanent well completion was the development of a new perforating method involving the use of a perforating gun small enough to go through tubing and with enough penetrating power to perforate the casing satisfactorily. The equipment described in this paper was basically designed to go through open-ended two-in. tubing and to perforate 5 1/2 in. casing. The equipment is illustrated in accompanying figures. NEW COMPLETION METHOD With the new permanent type well completion, the production string of casing is run and cemented through the potential producing zone. Open-ended tubing is -then run into the well and used to replace the drilling fluid with oil or water. The amount of the oil or water "cushion" can be varied according to the expected bottom hole pressure. Thus, if desired, the hydrostatic head can be greater, equal to, or less than the expected pressure within the zone to be perforated. Mud having been displaced, the tubing is -then positioned so that the end is above the zone to be perforated. The Christmas tree and flow lines are installed and tested. The drilling rig can now be removed from the well since it is not necessary for the perforating operation, thus resulting in a saving in rig tme. The perforating equipment is then moved to the well and the gun lowered through the tubing. After the gun passes out of the open end of the tubing, the depth measurements are correlated with casing collars and the gun placed opposite the zone to be perforated. Upon firing, the gun body is shattered and falls to the bottom of the well. Since the lubricator has been closed before firing, it is possible to withdraw the casing collar locator and head and then open the well immediately to the tanks. DESCRIPTION OF GUN An investigation of two basic types of Perforators, namely a bullet perforator and a shaped charge gun, was undertaken towards the development of a gun to meet the aforementioned requirements. Preliminary tests on various bullet perforator designs within the diameter limitations did not give results encouraging enough to warrant further study. Somewhat better results were obtained by using conventional shaped charges; but, as in the case of the bullet types, the necessarily small overall diameter proved to be a severe limitation and resulted in insufficient penetration. A shaped charge employing a new design was developed which gave excellent penetration. These shaped charges were designed to produce a large hole directed upward at a 45" angle. In order to obtain the maximum space possible for these charges, it was decided to use an expendable carrier made of drillable material, usually aluminum, which would shatter at the time of firing, leaving only harmless fragments at the bottom of the hole. The fixed spacing is five shots per ft. wit11

Jan 1, 1952

By A. R. Marder

Fe-C alloys have been investigated at temperatures below the eutectoid transformation to determine whether the superplasticity phenomenon exists for these materials. As a result of void formation at the carbide-ferrite interfaces that led to void coalescence and premature failure, the large elongations characteristic of superplustic matem'als did not develop. It had generally been assumed that the conditions for obtaining superplasticity include: a) T/Tm =0.5, b) a strain rate sensitivity of the flow stress that results in m s 0.3, and c) fine-grain structure. Although superphsticity was not observed in the alloys studied, stress-strain rate measurements indicated that an additional condition for the existence of superplasticitv in Fe-C allovs is that the fine-wain structure be equiaxed in agveement with previous results on Zn-Al al1oys. MANY investigations have been made on the hot deformation of iron-base alloys in order to obtain an understanding of the variables that may affect high-temperature strength and ductility.'-6 However, only a few papers have studied the effect of carbon content on the high-temperature properties of iron-base alloys. Peretayat'ko and Zaikov and Wallquist and carlen4 studied the effect of carbon on strength in the y temperature range, and recently Robbins, Shepard, and sherby6 studied the influence of carbon on the torsional behavior of iron in the a and cu-plus-cemen-tite temperature ranges. In addition to studying hot deformation behavior, current investigations have also been concerned with the phenomenon of superplasticity in nonferrous materials, and, in particular, ferrous alloys. Superplasticity is identified by the increased resistance to necking that results from a high strain rate sensitivity of the flow tress. ackofenr proposed the strain rate hardening coefficient m = d log a/d log e as a measure of strain rate sensitivity and showed that m must be "-0.3 for superplasticity to exist. In addition, a fine-grained microstrcture,- a eutectic or eutectoid ornosition,- and a test temperature in the creep range (TITm > 0.5)" have also been advanced as necessary conditions for superplasticity. Although many workers""23 investigated superplasticity associated with the allotropic transformation in iron and the eutectoid transformation in steel, and recently chadler " and orrison have studied superplasticity in the a - y phase region, no one has measured the superplastic effect in the a + cementite temperature region. Therefore, the purpose of this investiga- tion was to determine the effect of carbon content, microstructure, and test temperature on the strain rate sensitivity of Fe-C alloys at temperatures below the eulectoid transformation. EXPERIMENTAL PROCEDURE Fe-C alloys with carbon contents of 0.2, 0.4, 0.6, 0.8, and 1.0 wt pct were prepared as 100 lb 54 in. ingots by vacuum-induction melting. The compositions are listed in Table I. All the ingots were soaked at 2100°F for 2 hr prior to rolling to 4 in. plate. The minimum finishing temperature for rolling was 1800°F. The plate was sheared in 10-in. lengths and tensile bars were rough-machined prior to heat treatment. The samples were heat-treated according to the schedule in Table I1 in order to obtain a fine-grained microstructure consisting of finely dispersed carbides in a matrix of small ferrite grains. All the tensile specimens were heat-treated prior to finish machining, and tensile testing was accomplished on an Instron machine in conjunction with a clamshell furnace. The buttonhead specimen configuration is shown in Fig. 1. An inert atmosphere, helium, was maintained during testing with a steel tube shroud in a bell jar arrangement. Predried helium was continually purged over the buttonhead specimens and exhausted through ceramic fiber at the bottom of the tube to keep oxidation at a minimum and to prevent

Jan 1, 1970

By Attilio Ligasacchi

The production of magnetite in copper smelting has been the subject of wide investigation in the past.1-9 It is accepted that most of the magnetite is produced in the converters and returned to the re-verberatory furnace with the converter slag. Part of this magnetite is circulated back to the converters in the matte. As a result, if proper fluxing practice is not followed, and temperatures in the converting process are not strictly controlled, an undesirable circulating load of magnetite may occur. This might ultimately lead to saturation in the reverberatory bath with consequent settling of the magnetite at the furnace bottom. This paper contains the cumulative results of a series of investigations on the formation and control of magnetite presently being conducted by the metallurgical department of the Chino Mines Division of Kennecott Copper Corp. THE METHOD USED FOR THE DETERMINATION OF MAGNETITE AND OTHER MAGNETIC MATERIALS Conventional wet chemical magnetite analyses are neither reproducible nor theoretically sound. These analyses are sensitive to time, temperature, and heating rate-all factors which tend to decrease their accuracy and reproducibility. To overcome the unreliability of the above method, H. W. Mossman1 developed an inductance bridge magnetite "meter" which completely replaces the wet chemical analysis with the advantage of quicker and reproducible results. Fig. 1 is a schematic diagram of the meter in use in the Chino metallurgical laboratory. When a pulverized specimen with magnetic properties is placed in coil A of the inductance bridge circuit, the permeability of the core is increased, causing a voltage drop across the coil. The unbalance is amplified to operate the balancing motor which adjusts the movable iron core in coils C and D to bring the bridge back to balance. An indicator arm connected to the movable iron core indicates directly the percent, or weight, of the magnetic material in the specimen. Coil B acts as a reference coil and a zero adjustment. The meter is standardized with a set of various amounts of -100 mesh pure magnetite crystals (99.8 pet Fe3O4) thoroughly mixed in a test tube with -100 mesh quartz sand. These mixtures are combined with inert transoptic powder and heated to 140° to 160°C (melting point of the powder) and subsequently cooled to room temperature, solidifying the whole mixture, in this manner eliminating any possible preferential orientation of the magnetite grains. In order to determine the effect of the particle size" of the magnetite and other magnetic materials on the meter determinations, five size fractions were obtained by wet screening of the 99.8 pet pure magnetite and several reverberatory matte and slag samples. The following size fractions were investigated: -80 + 100; -100 + 150; -150 + 180; -200 + 250; -250 + 300 and -100 mesh. The deviation of the meter readings for the various sizes of each sample never exceeded ±0.3 pet. This deviation is considered well within the accuracy needed for common smelter operation procedures. Further investigation proved also that grinding slag and matte samples in our laboratory disc grinder normally produced a powder constituted of 95 pet -200 + 300 mesh. For practical purposes a - 100 mesh size was used as standard for our studies. It was also found that for rapid determinations tamping of the sample in the test tube was sufficient and gave comparable results with those samples solidified with transoptic powder. For the purpose of this study we define the magnetite content in copper smelter slags and matte as the equilibrium amount of crystalline Fe3O4 existent in the specimen at a given reference temperature, excluding any secondary magnetite precipitated in the specimen during nonequilibrium conditions. MAGNETITE DETERMINATION IN SLAGS Since it is impossible to ascertain the status of equilibrium of a conventional slag sample, due to

Jan 1, 1965

By S. K. Mitra

A refinement of the Seeger model for intersection process is investigated which is in better agreement with experimental observations than the original. It is shown that, in single crystals, the strain hardening in Stage II is mainly due to the short-mnge interactions when intersection is the rate-controlling process. It is also demonstrated that the creep curves, as predicted by this theory, are in good agreement with the experimental observations. MOTT: Cottrell? and Seeger have developed an approximate theory for plastic deformation of single crystals over the range of variables where the strain rate is controlled by the rate of intersection of dislocations. Although it is now generally agreed that the low temperature deformation of many pure metals is controlled by the intersection mechanism, various dictates of the over-simplified Seeger model are not in good agreement with all the experimental facts. It is the purpose of this paper to reveal that by appropriate extensions of the Seeger model, particularly those suggested by Basinski: a much more reliable theory results. The refined model to be presented here will be shown to account, in a satisfactory way, for the effect of stress, temperature, and strain on the creep rate under constant stress and for the effect of temperature and strain on the flow stress in constant strain-rate tests. EXPERIMENTAL METHOD To Check the theoretical deductions, both creep and tension tests were conducted on single crystals of high-purity Al (99.994 pct Al) so oriented that both the active slip plane and the Burger's vector of the operative dislocations on that plane made angles of about 45 deg to the tension axis. Single crystals of aluminum (5/8 in. by 1/10 in. by 8 in.) were grown in a graphite mold under an inert atmosphere using the Bridgman method. The chemical composition of the aluminum is given in Table I. A common seed was used for all crystals and their orientation is shown in Fig 3. The extensometry consisted of two differential transformers with matched outputs mounted so as to measure the extension in a 2 in. gage length. The amplified difference in output of the two transformers was recorded by a potentiometer. The amplification was calibrated before each run so that one chart division of the recording potentiometer indicated a 105 shear strain on the specimen. Stresses were measured to the nearest 1 x l04 dynes per sq cm. ACTIVATION VOLUME The concept of Basinski that the activation volume is not constant but is a function of the force that aids thermal fluctuation in effecting the cutting of dislocations will be adopted in this section. The force acting due to an applied shear stress on the dislocation being intersected is where L is the mean spacing between the forest dislocations, is the Burger's vector and T is the back stress. A detailed discussion about the origin of back stress will be taken up later in this report. As shown by Basinski, the activation energy, U, that must be supplied by a thermal fluctuation in order to effect intersection, is equal to where x is the distance through which the dislocation must be translated for complete intersection, and Fm is the maximum force encountered in intersection. When the applied stress is decreased abrupt ly F also decreases and the activation energy increases correspondingly as documented in Eqs. [ 1 ] and [2]. The Seeger equation for the strain rate when the deformation is controlled by rate of intersection of where the shear strain rate (per sec) N = the number of points per unit vol at which intersection can take place

Jan 1, 1962

By J. L. Huitt, S. R. Darin

The use of a partial monolayer of propping agent to obtain a high flow capacity for a hydraulically induced fracture is discussed. From the results of laboratory work it was shown that a modified form of the Kozeny-Carman relation could be used to describe the flow in the partial monolayer propped fracture. With equations presented in the paper, the density pattern of the propping agent (number of particles per unit of fracture surface) that results in the maximum flow capacity for the fracture can be determined. The maxitnum flow capacity obtained with a partial monolayer is often an order of magnitude greater than the flow capacity obtained in greater width fractures containing multilayers of the propping agent. INTRODUCTION One of the predominant factors controlling the success of a hydraulic fracturing operation is the propping of the fracture. The trend in hydraulic fracturing recently has been an increase in the ratio of propping agent-to-fluid. The primary purpose of this increased ratio is to sustain a propped fracture of greater width. It was shown in a recent study' that for some formations low concentrations of propping agent would result in "closing" or "healing" of the fracture. In some cases, in an effort to insure sufficient propping agent concentration, the fracturing operation is designed to obtain a pack of the propping agent in the fracture. The placing of a pack of propping agent in the fracture provides a fracture of maximum width; however, the width alone does not control the flow capacity of the fracture. The flow capacity is dependent on the permeability of the fracture, as well as the fracture width. Thus, increasing the permeability to obtain larger flow capacities is as important as obtaining a greater width fracture. By increasing the concentration and particle diameter of the propping agent, fractures of sufficient flow capacities are obtainable for effective well stimulation in a number of formations. However, in some of the recovery techniques, the fracture capacity obtained by commonly used practices in fracturing is inadequate. One such recovery technique is gravity drainage via a horizontal fracture placed near the base of the productive zone. For this recovery technique, the fracture capacity desired is usually one to two orders of magnitude greater than the fracture capacity normally obtained in fractures propped with the size of sand commonly used (20-40 mesh). By the use of particles of larger diameter, the permeability of a pack2 is increased; however, even a multilayer pack results in a fracture flow capacity lower than that desired for some formations. For example, a fracture of 1/4-in. width, packed with 8-12 mesh sand of about 2,000 darcies permeability, would result in a fracture flow capacity (permeability X fracture width) of only 40,000 md-ft. If the fracture flow capacity needed is an order of magnitude greater than that cited in the example, an approach different from the formation of a pack in the fracture is needed. One approach is a partial monolayer of large-diameter propping agent. In working with this type of fracture system, it was found that the flow problems were unique to the system; therefore, the study discussed in this paper was made. THEORY AND RELATED STUDIES Previous studies, although related to that described in this paper, have not pertained directly to the flow in a fracture propped with a partial monolayer of propping agent. However, studies such as flow between parallel plates,' flow in packed columns,' fracture flow capacity' and flow in simulated fractures reduce the scope of this flow problem considerably. A fracture propped with a partial monolayer of dropping agent may be considered to be between two extremes — between an open fracture and a bed of particles packed In a conduit of fracture geometry. Thus, it seems appropriate to review the extreme cases.

By I. R. Kramer

The stress associated with the high-dislocation layer at the surface of deformed aluminum crystals was measured by progressively polishing the specimen and determining the change in the initial flow stress when the specimen was reloaded. This surface stress increased linearly with strain in the Stage II and Stage III region. The depth to which the high-dislocation layer extends was found to be 0.0025 in. The influence of the surface layer on Stage I is discussed and it is shown that Stage I ends when, on the secondary slip system, the difference between the applied stress and the opposing stress due to the surface layer is equal to the critical resolved shear stress. In previous papers'-3 it was shown that the surface exerts a pronounced effect on the plastic-flow characteristics of metals. For fcc metals it was possible to increase the extent and decrease the slope of Stages I and II and decrease the work-hardening coefficient of Stage II by removing the surface of the specimen continuously during the deformation process.2,3 More recently,L it was shown that the activation energy for plastic deformation and the activated volume were also influenced by the surface. These changes were related to the existence of a high concentration of dislocations in the region near the surface of the deformed specimen and it was concluded' that the effective stress, T, acting on a dislocation was T = Ta-Ti-Ts [1] where 7, is the applied shear stress, Tj is the internal stress, and TS is the stress associated with the surface region containing the high concentration of dislocations. From Eq. [I] it is seen that the effective stress, 7, acting on a dislocation is a function of 7,. Therefore 7, becomes important in any description of the work-hardened state. In this paper the value of 7, as a function of strain and its variation with depth are reported for aluminum single crystals. In particular, it will be shown that the end of Stage I occurs when the net stress, 7, on the secondary slip system is equal to the critical resolved shear stress. To measure Ts, single crystals of aluminum were deformed in tension to various strains, whereupon the load was removed and a known amount of metal was removed by electrolytic polishing. The differ- ence between the final flow stress before the specimen was unloaded and the initial flow stress upon reloading is designated as ?Tp and is the decrease in the stress due to the removal of the surface layer in whole or in part. The value for Tp is taken equal to -Ts. EXPERIMENT PROCEDURE The details of the tensile apparatus, method of electrolytic polishing, and crystal preparation were the same as those used previously,1"3 The single crystals were 1/8 by 1/8 in. in cross section with a 3-in. gage length. The initial purity of the aluminum was 99.997 pet. The specimens were held in vacuo at temperatures 50°C below their melting points for 16 hr and furnace-cooled. Just prior to the application of the tensile stress the specimens were electrolytically polished to remove 0.004 in. from the thickness. This latter operation was conducted after the specimen was placed in the tensile apparatus. The tensile tests were conducted at a strain rate, i, of 10"5 sec-' and a temperature of 3°C. The aluminum crystals used for the determination of 7, are designated A1-3-12 and A1-116. The orientation of the various crystals used throughout the investigation is shown in Fig. 1. The ?Tp values for specimens A1-3-12 as a function of Ax were determined mainly at four strain levels, i.e., ? = 0.0217, 0.0381, 0.0593, and 0.075, where Ax is one half of the total reduction of the thickness after polishing and ? is the re-

Jan 1, 1965

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 D. R. Uhlmann, T. P. Seward, K. A. Jackson, J. D. Hunt

microscopic ohservations on alloys of organic trzaterials show that dendrite arms can melt off under normal conditiorzs of growth. This occurs because of the interactiorz of' heat and matter fluxes during dendritic growth. The re melting occurs where the den-drite arms are joined onto the main deizclrite sterrz. Many detached crystals can be produced in this fashion. It is postulated that this remelting phenornenon occurs in metal alloys, giving rise to the equi-axed region in castings. Experiments peyformed orr castings of transparent materials support remelting as a primary source of detached crystals. Ob~ervations an equiaxed tones in metal castings are discussed in terms of this and other mechanisms. It is concluded that the partial remelting of dendrites is on important mechanisnl for producing the equiaxed regions in castings. THE freezing of a casting can be divided into three regions: a chill zone formed near the mold wall, a columnar region, and a third region, at the center of the casting, which is called the equiaxed zone. The origin of the grains that comprise the equiaxed region is the subject of this paper. We will first outline the mechanisms which have been proposed to account for the origin of the equiaxed region and discuss some recent experiments which cannot be explained by these theories. We will then present some microscopic observations of dendritic growth which have led us to conclude that the equiaxed zone in castings can arise from partial remelting of dendrites. This proposal is supported by observations on transparent castings. Finally, we will discuss several observations on metal castings which could not be explained previously, and which can be accounted for by the present proposal. Previous Mechanisms. Winegard and Chalmers3 pointed out that the molten central region of an ingot could be constitutionally supercooled, because of the solute layer at the growing interface. In this supercooled region, nucleation on particles could occur, giving rise to the equiaxed crystals. Another mechanism has also been suggested,4y5 that the nuclei originate in the chill zone then float or are carried by convection to the center of the casting, where they grow to produce equiaxed crystals. In this case, all the crystals nucleate immediately on casting, so this mechanism will be termed "big-bang" nucleation. Walker's Experiments. J. L. Walker has made observations on the solidification of nickel and Ni-Cu alloys which cannot be accounted for by either of the above mechanisms. In studying the grain size resulting from freezing of nickel at various undercoolings, walker6 found that at large undercoolings the grain size was small, due to mechanical nucleation, as discussed by walker7 and Horvay.' At small undercoolings, the grain size was large, and the growth dendritic. In Ni-Cu alloys, however, for initial undercoolings starting from just below the melting point and extending for a hundred or so degrees of undercooling, fine grains were also observed. Samples were undercooled 200 Centigrade degrees without nucleating, and were then heated 130 Centigrade degrees so that the undercooling was 70 Centigrade degrees, and nucleated. The samples had a small grain size, indicating multiple nucleation. There were no heterogeneous nuclei in the liquid which could operate at the temperature of growth. Constitutional supercooling could not have been greater than the supercooling of the initial melt, since the temperature rose as the freezing continued. walker -as also shown that a Ni-Cu alloy which can be undercooled to below its solidus will produce typical columnar and equiaxed regions when cast into a cold mold. This alloy could not reach a temperature below its solidus by constitutional supercooling. It could not, therefore, have nucleated heterogeneously. Neither the big-bang mechanism nor the Winegard-Chalmers mechanism can account for the fine grain

Jan 1, 1967

By E. W. Dewing

The nature of solid solutions of sodium in non-paphitic carbon at temperatures near 1000°C has been investigated by an electrolytic technique. The activity coefficient is found to vary strongly with the heat-treatment temperature of the carbon, and an analysis of the energy terms involved suggests that this is due to variations in the Fermi level. The diffusion coefficient for sodium is in the range implying that the sodium atoms (or ions) are not bound to specific sites on the carbon, in agreement with lack of evidence for specific compounds in the activity -composition curves. In the light of These results phenomena occurring in aluminum reduction cell linings are dis-cussed. It is concluded that sodium migrates by diffusion through the carbon lallice and not the pores, and that A14C9 is formed within the lining hy The excess of NaF invariably found in used linings is formed by this reaction and by air oxidation: There has recently been much interest in the swelling of carbon which occurs when it is made cathodic in an aluminum electrolysis cell, and there is general agreement that the phenomenon is due to the penetration of sodium atoms, derived from the reaction on the carbon lattice. (NaF and AlF3 are constituents of the fused-salt electrolyte.) The low-temperature lamellar compounds of potassium with graphite are well-known, and one compound between sodium and graphite has been reported.4 There is, however, virtually no information on the basic chemistry of the interaction of sodium with nongraphitic carbon at high temperatures, although Dell3 has made some fundamental studies of rates of penetration. At the time the present work was carried out these were not available. Rapoport and samoilenkol had measured rates of physical expansion in a dilatometer, using the fractional increase in length of a carbon specimen after electrolysis in a cryolite bath for 2 hr as an index of susceptibility to cathodic swelling. Work here with a similar apparatus showed that expansion (of a 1.5-in.-diam cylindrical specimen) ceased in less than 20 hr, and that the final expansion was much more reproducible than the value at some intermediate time. This suggested that an equilibrium was being set up, and that a thermodynamic approach was required. The results obtained are not as accurate as had been hoped, but they are the only ones of their type available and have enabled useful conclusions to be drawn. The experimental work was divided into two parts: an investigation of the activity-concentration diagram for sodium in carbon, and a rough determination of the diffusion rate. For the former, most of the experiments were made with petroleum coke even though this is not a form of carbon used in cell linings. It has the advantages, however, that it is a homogeneous material, that it is low in ash, and that its structure is reasonably well understood. Anthracites and coal cokes gave results of the same type as obtained with petroleum coke, but the reprodu-cibility from specimen to specimen was extremely poor. EXPERIMENTAL 1) The Activity-Composition Diagram for Sodium in Carbon. The petroleum coke was ground, pressed at 40,000 psi into plates 19 by 9 by 2 mm without added binder, and heat-treated. The sodium activities were established by an electrolytic technique. The carbon was made anodic in a cryolite bath with an aluminum cathode, and a controlled potential, E, applied so that the sodium activity at the anode was given by ln Na = 1n asa - EF/RT [21 where as, is the sodium activity generated by Reaction [I]. The relative activity is sufficient for many purposes and is simply given by The value of Na not at present known, but when the thermodynamics of the Na-A1 system have been worked out it can be derived from the equilibrium

Jan 1, 1963

By D. Turnbull, R. E. Cech

FISHER, Hollomon, and Turnbull have developed a theory for the nucleation of martensite. They first tested the theory on Fe-C alloys and low alloy steels. The major factor influencing nucleation of martensite was considered to be statistical composition fluctuations occurring in small regions at high temperature and frozen-in on quenching. These local regions of varying size and composition serve as nucleation centers. They become supercritical, one by one, as temperature is progressively lowered, resulting in temperature-dependent or athermal transformation. Fisher next applied nucleation theory to substi-tutional solid solution alloys. Detailed predictions were made for Fe-Ni alloys because of the availability of free energy data on this system. It was shown that composition fluctuations that were significant energy-wise did not occur, and nucleation frequencies could be calculated from average properties of the system. Nucleation was predicted to occur as time-dependent and having the functional relationship to give a C curve of nucleation frequency vs temperature. The analysis further predicted that the nucleation frequency was extremely sensitive to composition. Experimentally, it would be found that the transformation in some compositions is so slow that measurable amounts will not form in a reasonable length of time. With other compositions, only slightly different, the nucleation frequency becomes so great that the material becomes transformed while still distant in temperature from the maximum nucleation frequency. On quenching an alloy of such composition, the observed transformation kinetics would be similar to those found in Fe-C alloys. Cech and Hollomon repeated experiments of Kurdjumov n which the kinetics of transformation were similar to those predicted by Fisher for Fe-Ni alloys. The alloy studied in this investigation contained 73.3 pct Fe, 23.0 pct Ni, and 3.7 pct Mn. Fisher,' using an idealized model for the extent of transformation as a function of the number of martensite crystals per grain of parent phase, derived nucleation frequencies from the transformation curves of Cech and Hollomon. Complicating influences such as coupling effects between grains in the polycrystalline specimens were neglected. Nevertheless, excellent agreement was found between the theoretically and experimentally derived nucleation frequencies. These experiments, however, could not provide a critical test of the theory. Experimental nucleation frequencies could vary widely from those calculated, depending upon the extent of deviation from ideal partitioning and the extent of coupling effects. Further, since the compositions of material theoretically analyzed and experimentally determined were different, the free energy changes involved in the experimental work could only be estimated. Also, the effect of heterogeneities on the transformation kinetics was not considered. For these reasons, it was decided that experiments designed to test the validity of the Fisher analysis must be conducted on binary Fe-Ni alloys, which were the ones considered theoretically by Fisher. The martensite transformation in Fe-Ni alloys has been the subject of considerable study. Machlin and Cohen have shown that transformation proceeds in a manner quite unlike that in any other ferrous alloy system. They found that single crystals would transform to a large extent in a single burst. In large grain polycrystalline specimens, frequently more than one grain and sometimes the whole specimen would transform at the same instant in this manner. Results on filings indicated that different particles would undergo the burst transformation at widely different temperatures. These results support the conclusion that the transformation behavior could not be described by a single nucleation frequency as would be the case if the nucleation were homogeneous. It appeared that further work was necessary to define the factors responsible for burst-type transformation, so that the conditions could be altered to favor homogeneous nucleation of martensite if such could be accomplished. This report summarizes the results of some experiments conducted with powdered Fe-Ni alloys for this purpose and the re-

Jan 1, 1957

By P. F. Gnirk, J. B. Cheatham

Single bit-tooth penetration experiments under static load were conducted on six rocks at confining pressures of O to 5,000 psi using sharp wedge-shaped teeth with included angles ranging from 30 to 120° In general, the force-displacement curves for all rocks exhibit an increasingly nonlinear and discontinuous behavior with decreasing confining pressure. The confining pressure at which a rock exhibits a macroscopic transition from predominantly ductile to predominantly brittle behavior during penetration vmies from about 500 to 1,000 psi for the limestones to greater than 5,000 psi for dolomite. The correlation between calculated values of force per unit penetration based on plasticity theory and experimental values is quite encouraging, even at confining pressures as low as 1,000 psi A qualitative correlation between volume of fragmented rock per unit energy input for a single bit-tooth and drilling rate for microbits appears to exist over a confining pressure range of 0 to 5,000 psi. INTRODUCTION Laboratory experiments utilizing a small-scale drilling apparatus have demonstrated that penetration rates are reduced considerably as a result of increasing the confining pressure from atmospheric to a few thousand psi?-3 This undesirable situation can, in general, be attributed to a combination of decreased efficiency of chip removal at the bottom of the borehole, increased rock-failure strength, and a possible change in the mechanism of chip generation and rock fragmentation with increasing confining pressure. To more fully understand the principles underlying the last circumstance, it is the purpose of this investigation to experimentally study the mechanism of single bit-tooth penetration into dry rock at low confining pressures and, in particular, to establish the confining pressure at which the penetration mechanism may undergo a brittle to ductile transition for various rock types commonly encountered in drilling. Confining pressure as used here refers to the differential pressure between the borehole fluid pressure and the formation-pore fluid pressure. EXPERIMENTAL PROCEDURE Using an experimental apparatus previously described: a single, sharp wedge-shaped tool was forced under a "statically" applied load into an effectively semi-in finite dry rock sample subjected to a prescribed confining pressure. To prevent the invasion of the con fining-pressure fluid into the pores of the rock sample during penetration, the exposed surface of the rock was jacketed with a layer of silicon putty.* Electrical instrumentation incorporated into the apparatus yielded a graphical plot of force on the tool as a function of penetration or displacement of the tool into the rock during an experiment. During the course of the experimentation the following conditions were maintained constant: (1) pore pressure — atmospheric (i.e., the rock was dry); (2) temperature — 75F; (3) rate of loading — essentially static (approximately 0.002 in./sec); (4) bit tooth — a sharp wedge-shaped tool loaded normal to the rock surface; (5) rock surface — smooth and flat; (6) drilling fluid — hydraulic oil; and (7) maximum depth of penetration — approximateiy 0.1 in. In addition, each experiment was performed on a different rock sample so the rock surface is free of a layer of cuttings and of any previous indentation craters. The influence of the corners of a borehole was neglected, since each rock sample was cemented inta a section of aLuminum tubing to simulate a semi-infinite body. Indentation experiments were made on a sandstone, two limestones, a marble, a dolomite and a schist with a 60' bit-tooth over a range of confining pressures from atmospheric to 5,000 psi and at a constant confining pressure of 1,000 psi for a variety of bit-teeth with angles ranging from 30 to

Jan 1, 1966

By C. Gatlin, M. E. Chenevert

The effects of bedding plane orientation on the elastic constants and the yield strengths of three laminated rocks (one sandstone and two shales) and one isotropic rock (a limestone) were studied. The directional dependence of the elastic properties of these rocks was experimentally evaluated using a triaxial compression cell and auxiliary stress-strain measuring equipment. Symmetry of Poisson's ratios within the bedding plane suggested that horizontal isotropy exists, but the bedding planes do give rise to an appreciable difference between properties in the horizontal and vertical directions. For the three bedded rocks studied, Young's modulus was lower normal to bedding than along bedding. Yield strengths were determined at confining pressures from 0 to 12,000 psi in a triaxial compression cell. The rocks studied showed strength reductions as high as 40 per cent when the test specimen was oriented at 20°-30° to the bedding planes. The mechanical behavior of these rocks suggested that the rock properties of shear strength and/or coefficient of internal friction can vary with direction, depending on the particular rock tested. Tensile strengths were also measured and found to be lowest when failure occurred along bedding. This work shows that bedded formations exhibit sizable directional variations in both their elastic constants and yield strengths. It is suggested that these variations may be accounted for by using the "elastic laminate" model and the "variable coefficients" failure model. INTRODUCTION The nature of rock deformation at elevated pressures has been studied by many workers; papers by Handin,l,2 and Robinson3 illustrate the present state of knowledge. Most investigators have either chosen rocks which were as isotropic as possible (in order to avoid complications of data interpretation and analysis) or they have oriented their samples so that the effects of anisotropies (such as bedding planes) have been avoided. Of the studies performed, few were concerned specifically with mechanical anisotrop-ies. Griggs4 has presented limited data for specimens cut parallel, normal, and 45° to the bedding plane; he relates the strength anisotropy observed to the fabric (bedding) anisotropy. His tests were primarily concerned with large deformations (20 per cent strain); thus, no directional values of the elastic constants were reported. Handin has reported the results of similar experiments. Bott5has discussed rock strength anisotropies due to faults, cleavage, or bedding. He was concerned primari ly with determining the shear stress on such planes and did not mention the effect of friction. Jaeger6 later generalized Bott's work by taking friction into account and presented a limited theory for the failure of rocks having a "single plane of weakness", and also for rocks having a constant coefficient of friction, and a shear strength which varies with bedding plane orientation. Donath and Cohen,7 and Donathg have evaluated rock strengths from shale and slate specimens cut normal and parallel to bedding. A dependence of cohesive strength (7,) on the specimen orientation was also shown. Adler9 has also studied this problem and lists similar results. He assumes that all bedded rocks behave according to Jaeger's single-plane-of-weakness theory. Kalinin and Belorussov10 list results for strength tests parallel and perpendicular to bedding and use this information as a basis for hole deviation analysis. From the literature it is apparent that sedimentary rocks have been tested under widely varying conditions of stress; however, the assumption of isotropy is generally, but not always, made. Since geologic sedimentation often deposits sediments in very definite lavers.,. it seems that more svstema tic attention should be given to the possible effects of this natural bedding. Bedding, as used here, refers to visible regularities of grain size or orientation resulting from depositional processes. The mechanical strength anisotropies, which bedding may present, is not the only anisotronic property to be considered, since the elastic

Jan 1, 1966

By H. H. Rachford

Numerical solutions of immiscible flow problems in which dispersive effects of capillarity are dominated by convection require excessively fine grid spacing with attendant high computing costs. The use of coarser spacing reduces cost but often produces oscillation or undue dispersion associated with displacement fronts, A numerical formulation is proposed here which should be applicable to two-dimensional flow problems. it is in part analogous to an approach previously tested for miscible systems. The convective transport is approximated using a change of variables to yield a coordinate system moving approximately with the local characteristic velocity. The capillarity-induced dispersive terms in the differential system describing the process are approximated with respect to a fixed coordinate system by the usual implicit formulation. One-dimensional tests of the procedure yielded results in which the saturation profiles tended smoothly to the zero-capillary pressure solution as the ratio of viscous to capillary forces was successively increased in a sequence of calculations. This contrasted favorably with solutions by other numerical procedures which would require attendant grid refinements to approach the zero capillary pressure results. INTRODUCTION Numerical solution of displacement problems has until recently relied on applying methods developed primarily for transient heat-flow problems. Such problems are classified as parabolic in type, and where the heat transport is purely by diffusion their solutions are characterized by a high degree of smoothness. It is not surprising, therefore, that for approximating these solutions available finite difference methods are quite adequate. In flow problems the transport is partly by diffusion, partly by convection or flow. Although the problem remains of parabolic type because the dispersive effects of capillary forces or diffusion play some role in every displacement, at high flow rates the problem is dominated by convection, and solutions tend toward those of equations of the hyperbolic type. Solutions of hyperbolic problems are characterized by the translation of fronts, or discontinuities, that may progressively increase in sharpness. Numerical methods for treating parabolic problems become less and less satisfactory as displacement rates increase and the role of dispersion due to concentration or capillary pressure gradients becomes small relative to transport due to flow. In computation the difficulty manifests itself as an error associated with the grid size chosen. 1-6 In summary, if the heat-flow type approximations are to include the terms arising due to convection, one of several choices may be made: (1) an upstream (to the direction of flow) approximation for the convection terms may be used; (2) a centered-in-distance (CID) approximation may be used; or (3) a recently developed approximation based on the theory of oscillation matrices may be chosen.6 The last appears to have significant promise for one-dimensional flow problems; its extendibility to two or three dimensions is an open question. In either of the first two approaches, a suitably small ratio of v&/D must be maintained, where v is the velocity, & is the grid spacing and D the effective dispersivity in the direction of flow. In the first choice, the approximation of the convective part is only first-order correct and errors introduced appear as a numerically induced dispersivity of magnitude proportional to v?x. In the CID choice, the approximation can be second-order correct, but the difference formulation fails to satisfy the maximum principle unless a condition on v?x/D is met. Practically, this means that for high flow rates oscillatory solutions may result in the neighborhood of a front unless exceedingly small grid intervals are taken. While the procedure proposed by Stone and Brian4 permits a less severe limitation to be placed on this ratio, ultimately the flow rates increase relative to the dispersivity the oscillation obtains. Further, extensions of their approach to higher dimensional systems may be attended by considerable

Jan 1, 1967

By T. K. Perkins, R. N. Hoffman, O. C. Johnston

A simplified theory of viscous fingering in miscible systems has been developed. It predicts the correct functional relationship between pertinent variables and permits the calculation of the order of magnitude of fingering behavior for simple systems, such as linear or radial. The theory is based on the following four observations: (1) cross-flow takes place only near the ends of fingers; (2) the relative finger width is about 0.5; (3) fingers can be suppressed by transverse dispersion, the suppression being quantitatively described by a critical ratio of dispersion times; and (4) the widths of extending fingers are close to the minimum size finger that can grow at any stage of displacement. Fingering is studied in two-dimensional linear and diverging radial systems, both theoretically and experimentally. For linear systems, the length of the fingered region is proportional to mean displacement; the finger width is proportional to square root of mean displacement; and there is a small initial region void of fingers because of suppression by longitudinal dispersion. For the radial system, two limiting conditions are recognized. If the mean displacement is small compared with the wellbore radius, the length of the fingered region is proportional to the mean displacement. The width is proportional to the square root of mean displacement. If the mean displacement is large compared with wellbore radius, length of the fingered region is proportional to mean displacement, but the number of fingers approaches a constant value. Also, in the radial case there is a small initial region void of fingers because of longitudinal dispersion. INTRODUCTION The behavior of viscous fingers in miscible systems has been of interest to the oil industry for many years. Previous studies have clearly shown the existence and nature of fingers in small models.1,2,3,4 Engineering techniques for extrapolating to reservoir situations have been proposed. 5,6 Still, because of the lack of a real understanding of the mechanics of fingering, there remains uncertainty and disagreement as to the best way to scale or calculate fingering behavior in the reservoir. In this paper we discuss a study of the mechanics of fingering in miscible systems (i.e., why do fingers behave as they do?). A simplified theory is developed which we believe will (1) predict the correct functional relationship between pertinent variables, (2) permit us to calculate the order of magnitude of fingering behavior for simple systems such as linear and radial, and (3) give further insight into the problem of scaling or otherwise extending the results to more complicated reservoir conditions. The paper includes the following sections: (1) brief summary of the mixing behavior of miscible fluids in linear and radial systems; (2) four fundamental observations of fluid flow (under fingering conditions); (3) based on these four observations, a simplified theory of fingering in linear and radial systems is developed; and (4) the theoretical equations are compared with experimental fingering data measured in laboratory models. A REVIEW OF DIFFUSION AND DISPERSION As will be shown later, the behavior of viscous fingers in miscible systems is controlled in large degree by the mixing between the two fluids. A quantitative description of fingers will first require a quantitative description of mixing. Fortunately, much work has been done to clarify this subject; a fairly comprehensive review has been presented in a previous paper. 7 For present purposes let us (1) briefly summarize the quantitative description of diffusion and dispersion (both longitudinal and transverse) within a differential element; (2) present a simplification of the "width of mixed zone"; and (3) describe the effect of geometry on width of mixed zone. If one fluid is displaced from a porous medium by another miscible fluid (in the absence of fingers)

Jan 1, 1966