Search Documents

By U. F. Kocks

The flow stress in some latent slip systems of aluminum and silver crystals after various deformations in single slip was investigated by transverse compression and supplemented by experiments on overshoot of the compression axis and the stability of corner orientations. Work hardening was found to be remarkably isotropic, the ratio of the flow stress in the latent system to that in the primary system being between 1.0 and 1.3. The strain-rate dependence was not significantly larger in the latent system. This low anisotropy suggests that specific dislocation interactions are not the direct cause of work hardening. Either the obstacles are essentially impenetrable so that only their distance enters, or a fairly random distribution of dislocations of all systems is. in fact. formed during single slip. The latter seems to he excluded by the observations on the strain-mte dependence and on the amount of secondary strain. The amount of seconzdary strains occurring in single slip was determined by shape-change mensurements and by a new method connected with the latent-hardening experiments. It is definitely less than 1/300 times the strain in the primary svstem. WORK hardening is due to dislocation interactions; this is a hypothesis, and possibly the only one, on which all present work-harden ing theories agree. About a dozen specific interactions have been considered, and it is unlikely that any essential one has been missed. Different work-hardening theories differ to a large extent in the specific interactions which are deemed most important. The presence or absence of any particular interaction, and its degree, depend on the geometrical relationships between the two dislocations interacting. Let us enumerate some examples. Two dislocations may react, with a varying reaction energy, to form a third dislocation, and this reaction product may be sessile or glissile. Two dislocations may intersect and a jog may be formed on both, on one of them, or on neither; this jog may be sessile17Z or glissile and it may, if moved nonconservatively, lead to the production of vacancies or interstitials. Finally, a dislocation may bypass another one and thus interact only with its stress field. All these interactions depend on the specific relation between the Burgers vectors and line directions of the two dislocations. Any dislocation generated during plastic deformation of fcc crystals belongs to one of twelve slip systems. Taking any one dislocation, its relationship to each of the other eleven can fall, by symmetry, into one of only five classes. These relationships are enumerated in Fig. 1 and Table I. Each slip system is represented in the stereo-graphic projection by two triangles, namely those in one of which the specimen axis would have to lie if the specimen were to deform in single slip on this slip system in a given sense, in a tension test. Let us represent the primary slip system by the shaded triangles. The other slip systems are classified only with respect to their relationship to this primary system. The code letters (initials of established names* in German) are explained in Table I. Three specific interactions are also listed in cases where they may be strong: the product of a reaction, the formation of jogs in either system, and the elastic bypass stress. The last column relates to the types of experiments reported in this paper and summarizes their results. To decide experimentally between the relative importance of the various interactions one would like to produce a group of dislocations in one slip system and then measure the stress necessary to force another group of dislocations of a specific second system through them. Hopefully, one can produce the first set of dislocations by plastically deforming the specimen in single slip. Then one can try to achieve single slip in a specific, previously latent system and measure its flow stress. If it is different from the virgin flow stress, there has been "latent hardening". Some types of latent-hardening experiments will be discussed in the next section. The strain-rate and temperature dependence of the flow stress is an indication of the importance of dislocation intersections. By extending an investigation of this dependence to the latent systems, one can determine how much more important forest contributions are here. seeger3 has used the as-

Jan 1, 1964

By Craig S. Hartley

Expressions are obtained for the energy changes associated with the reaction of (a& (111) slip dislocations on intersecting (110)planes in anisotropic bcc metals. An energy criterion for assessing the likelihood of dissociation of the products of such reactions is also presented. It is found that the "burrier reactions" which form a(100) dislocations at the intersection of two active {110) slip planes are more energetically favorable in metals which exhibit a high value of Zener's anisotropy factor, A, than those which have a low value. The results are presented in a form which permits the stacking fault energy to be obtained from a measurement of the separation between par-tials in a dissociated configuration. However, until accurate calculations or measurements of the stacking fault energies involved are available, it is not possible to assess the physical importance of dissociated dislocations. In a recent paper,' the energy changes associated with several types of reactions between two slip dislocations, (a/2)(111){110), in bcc structures were calculated.* Isotropic elasticity and the approxima- tion v = -3- were employed. The purpose of this work is to present calculations of the energy changes for many of the same reactions using anisotropic elasticity. The problem of dissociation of a(100) and a(110) dislocations is also considered, and maximum fault energies for which dissociation will be energetically favorable are calculated for several bcc metals. Two general types of reactions are considered; those for which the reactant (a/2)(111) dislocations have long-range attractive forces and those for which the reverse is true. An example of the former is: (a/2)[lll] + (a/2)[lll]-a[l00] while the latter are typified by: (a/2)[lll] + (a/2)[111] -a[011] Only reactants lying in different slip planes are considered; therefore, the products must lie along (111) or (100) directions, which are the intersection of two {llO} planes. It will be assumed that the reactants and products are infinitely long parallel dislocations, since in this case the energy change associated with the reactions is a maximum.' THEORY The self-energy per unit length of a straight mixed dislocation in an anisotropic medium can be written? where b is the Burgers vector, K is an appropriate combination of the single-crystal elastic constants, and R and ro are, respectively, outer and inner cut-off radii of the elastic solution. The energy given by Eq. [I] does not account for any variation of the core energy with orientation. This could be manifested by an orientation dependence of the core radius or, equivalently, the Peierls width, of the dislocation. However, the energy contribution due to this source is expected to be small, and current models of the dislocation core are not sufficiently accurate to justify such a refinement. It has already been shown that for the isotropic case the energy contributions due to nonzero tractions across the cores of the reactants and products exactly cancel one another in the reaction.' Accordingly, it will be assumed that this contribution to the total energy change in the anisotropic case is small. In the subsequent discussion it is also assumed that the core radii of the reactant and product dislocation are the same and that, where stacking faults are formed, the faulted region is bounded by the centers of the partials. Consequently only changes in elastic energy due to the reactions will be considered. When the dislocation is parallel to either the (111) or the (100) directions, K may be written:375 K = (Ke sin2 a + Ks cos2 a) [2] where K, and Ks are the combination of elastic constants corresponding to an edge and screw dislocation lying along the same direction as the mixed dislocation, and a is the angle between the direction tangent to the dislocation line and the Burgers vector. Eq. [2] should not be confused with the isotropic approximation to the variation in energy with line Orientation.6 It should be noted that the essentially isotropic expression for K is a result of the characteristic symmetry of the (111) and (100) directions and is not, in general, valid for other dislocation directions in anisotropic cubic metals. The energy* change for a reaction in which the re- actant and product dislocations are parallel perfect dislocations can be written: where Ep and E, refer to the self-energies of the products and reactants, respectively. For dislocations parallel to (100) and (111) directions, Eq. [3] becomes:

Jan 1, 1969

By B. G. Matson, M. A. Whitfield, G. R. Dysart

This paper describes the development of u computer program to calculate changes that occur in the length of tubular goods due to temperature and pressure changes during stimulation operations. Due to the numerous variables involved and the uncertainty of all static and dynamic conditions that could exist, it becomes a staggering task for individuals charged with completions to perform the necessary mathematical calculations. The computer program permits advance calculations for several sets of conditions. INTRODUCTION In the Delaware basin of West Texas alone, 50 wells were contracted or drilled to 15,000 ft or deeper in 1965. Deep well activity is continuing in this and other areas on an expanding scale. Many of these deep wells require extensive stimulation for successful commercial production, and during these operations, pressures and temperatures are encountered that have a pronounced effect on the length of tubular goods. This length change during a large-volume, high-pressure stimulation treatment utilizing fluids considerably cooler than bottom-hole temperature can be of such a magnitude that permanent damage to casing and tubing will result unless mechanical design, pressures and fluid temperatures are evaluated and controlled. These pressure and temperature effects can be calculated. However, the process lends itself well to computer solutions because of the mathematical nature of the problem and the calculating hours involved in arriving at an answer. The engineering-hour demand becomes more severe as tapered strings are involved. On initial treatments on a given well, surface pressure and injection rate conditions are unknown, and offset well conditions have not proven to be a reliable method for making predictions. For these reasons, it has become rather standard procedure for operators to compensate for these uncertainties by placing unnecessary pressure and fluid temperature restrictions on stimulation design. On a number of occasions treating fluids have been preheated to as much as 160F as a means of compensating for thermal contmction resulting from pumping cool fluids. The maintenance of packer seals has been treated by Lubinski, Althouse and Logan',' and the problem of therma1 effects on pipe has been explored by Ramey." These works were expanded and the results made applicable to everyday oilfield terminology before submitting them to computer programming. The pressure and temperature effects on tubing movement previously mentioned occur simultaneously as fluid moves through the pipe. The pressure changes, for purposes of explanation, are categorized here as to the various effects these pressures have on a tubing string. These divisions are (1) the piston-like results of forces acting on horizontal surfaces exposed to pressure, (2) swelling or ballooning of the tubing along its entire length due to the forces of pressure acting against the tubing walls, (3) the elongation of tubing due to frictional drag and (4) corkscrewing of the pipe due to internal pressure. Thermal changes are also of great importance, as their results may be more significant than any of the pressure effects. Steel is an excellent conductor of heat and the earth is a relatively poor conductor. It has been calculated that pipe temperatures at depths of more than 20,000 ft approach within as little as 25" the temperfature of the surface fluid after pumping for 2 hours, or a drop in temperature in some treatments of more than 220F. The equations presented in this paper were developed for computer programming and simplicity of input information; therefore, numerical constants such as Young's modulus for steel (28 X 10\ si), the coefficient of thermal expansion of steel (6.9 X 10."IF) and Poisson's ratio for steel (0.3) are included with unit conversion factors. The moment of inertia of tubing cross-sectional area with respect to its diameter was changed to a constant times (D' — d') where D is outer diameter and d is inner diameter. Units in the equations are length in feet, diameter in inches, density in pounds per gallon, pressure in psi, rate in barrels per minute and time in hours. PISTON-LIKE REACTIONS A change in tubing internal dimensions and the exposure of other horizontal surfaces to different pressures on the inside and outside of the tubing result in a reaction much like a piston under pressure. Such is the case when the internal diameter changes in a combination string of pipe, when seals of a slick joint assembly are subject to pressure and in the end effects of a tubing string. The change in tubing length due to the piston effects of a slick joint packer is affected by the various diameters involved, the tubing pressure Ap,, the casing pressure ,Ap,, length of pipe L, densities of fluid in the tubing before and during pump-

By K. A. Slagle

Hydraulic analysis of the wellbore has become increasingly inzportant for designing cementing operations and selecting equipment, materials and techniques to complenzent modern well-c-ompletion practices. Non-Newtonian fluid technology has advanced beyond the point where former empirical methods of analysis adequately define the hydraulic system and fluid properties. In view of these factors, this paper describes a series of rheological calculations which have been found practical, through field usage, for assistance in selecting a cementing program. A relatively simple laboratory method using standard viscometric equipment is suggested for determination of the rheological properties of slurries, and clrrta are presented on some of the more common cementitrg conzposition.A. A criterion for divergence from laminar-flow characteristics has been proposed. Usefulness of the calculations is indicated by examples of cementing operations where they have been used. INTRODUCTION With the changing aspects of well-completion practices during the past few years, it has been increasingly important to have a relatively simple method of analyzing the flow conditions existing in the well during cementing operations. This is particularly true in view of the improved economics toward which most of the changes have been directed. Rheological characteristics of slurries used for cementing should be a major consideration in the trend toward smaller casing sizes, either single or multiple strings. Receiving increased attention is the practice advocated in 1948 by Howard and Clark' of attaining turbulent flow with the fluids circulated during a primary cementing operation. While there may still be a difference of opinion concerning this technique, most available information indicates that superior primary-cementing results are generally obtained when high displacement rates are employed. Fluid properties of the slurry to be used must be available, as well as calculation methods, to determine what flow rates should be attained and the probable consequences in terms of frictional pressure and horsepower utilization. It would certainly be inappropriate to attempt high displacement velocities if sufficient pressure might be developed to create lost circulation. Since cementing slurries are non-Newtonian fluids, it is not possible to define their rheological or fluid properties by the single factor of viscosity and then make calculations for the quantities just described. Because the shear stress-shear rate ratio is not constant: it becomes necessary to establish at least two parameters for adequate fluid-flow calculations. It is not the purpose of this paper to delve into the mathematical development of non-Newtonian technology, nor to discuss the arbitrary classification system under which a single fluid may resemble two or three different classes depending upon experimental conditions. Rather, the intention is to present a useful series of calculations based on a concept applicable to both Newtonian fluids and to the preponderance of non-Newtonian fluids encountered in the oil-producing industry. Development of this approach was begun some 32 years ago,' and has most recently been brought to fruition by Metzner and his co-workers at the U. of Deleware. Some non-Newtonian fluids encountered in the petroleum industry, other than cementing slurries, have also had the benefit of this method of analysis."' The two parameters required to define the fluid are usually denoted by the symbols n' and K' and, for the purposes of this discussion, are called "flow behavior index" and "consistency index", respectively. These two slurry properties permit calculation of the Reynolds' number and the "critical" velocity, or the velocity at which departure from laminar flow begins. EXPERIMENTAL DETERMINATIONS The two principal instruments used for rheological studies are the pipeline (capillary-tube) viscometer and the rotational viscometer. When conveniently possible, a capillary-tube viscometer (where the pressure drop and flow rate of the material can be measured) is the better method for rigorous determination of the flow behavior index and consistency index for non-Newtonian fluids. With pressure-drop data at various flow rates, it is then possible to prepare a logarithmic plot of shear rate as the abscissa-shear stress as the ordinate. For fluids which do not exhibit time-dependency, these data will usually produce a straight line. The flow behavior index n' represents the slope of this line, while the consistency index K' becomes the intercept of this line at unity shear rate in accordance with the mathematical derivations associated with this concept of rheology. Due to the difficulties anticipated in maintaining a uniform, pumpable cement slurry for the time interval required to obtain measurements from the pipe viscometer, the n' and K' data reported herein were obtained using a direct-indicating rotational viscometer (Fig. 2). The

By W. G. Pfann

In zone-melting, a small molten zone or zones traverse a long charge of alloy or impure metal. Consequences of this manner of freezing are examined with impurerespect to solute distribution in the ingot, with particular reference to purification and to prevention of segregation. Results are expressed in terms of the number, size, and direction of travel of the zones, the initial intermsofsolute distribution, and the distribution coefficient. IF a charge of binary solid-solution alloy is melted and then frozen slowly from one end, as for example in the Bridgman method of making single crystals,' coring usually occurs, with a resulting end-to-end variation in concentration. Such coring, or normal segregation, is undesirable where uniformity is an object. On the other hand, for certain systems, it can be utilized to refine a material by concentrating impurities at one end of the ingot.'. ' In the present paper a different manner of freezing will be examined with respect to the distribution of solute in the ingot. A number of procedures will be indicated which have in common the traversal of a relatively long charge of solid alloy by a small molten zone. Such methods will be denoted by the general term zone-,melting, while the process described in the preceding paragraph will be called normal freezing. It will be shown that, in contrast to normal freezing, zone-melting affords wide latitude in possible distributions of solute. Segregation can either be almost entirely eliminated or it can be enhanced so as to provide a high degree of sttparation of solute and solvent. A number of simplifying assumptions will be invoked which, while not entirely realizable in practice, nevertheless provide a suitable point of departure for more refined treatments. Moreover, our own experience with zone-melting has shown that, for certain systems at least, the analysis holds quite well. The present paper will be confined to a discussion of principles and a general description of procedures. Comparison with experiment is planned for later publication. Normal Freezing Before considering zone-melting, segregation during normal freezing will be reviewed briefly. If a cylinder of molten binary alloy is made to freeze from one end as in Fig. 1, there usually will be a segregating action which will concentrate the solute in one or the other end of the ingot. If the constitutional diagram for the system is like that of Fig. 2, then the distribution coefficient k, defined as the ratio of the concentration in the solid to that in the liquid at equilibrium, will be less than one and the solute will be concentrated in the last regions to freeze. If the solute raises the freezing point, then k will be greater than one and the solute will be concentrated in the first regions to freeze. The concentration in the solid as a function of g, the fraction which has solidified, can be expressed by the relation: C = kC0 (1-g)k-1 [I] where C, is the initial solute concentration in the melt. Eq 1 is based on the following assumptions: 1—Diffusion in the solid is negligible. 2—Diffusion in the liquid is complete (i.e., concentration in the liquid is uniform). 3—k is constant. Concentration curves representing eq 1 for k's from 0.01 to 5.0 are plotted in Fig. 3. This equation, in one form or another, has been treated by Gulliver,³ Scheuer,4 Hayes and Chipman5 for alloys and by McFee2 for NaCl crystals. It is derived in Appendix I. It should be pointed out that the k which is calculated from the phase diagram will be valid only in the ideal case for which the stated assumptions are correct. In all actual cases, the effective k will be larger than this value for solutes which lower the melting point, smaller for solutes which raise the melting point, and will probably vary during the beginning of the freezing process. For simplification it will be assumed that the ideal k is valid. Zone-Leveling Processes The processes of this part are designed to produce a uniform, or level, distribution of solute in the ingot. Single Pass: Consider a rod or charge of alloy whose cross-section is constant and whose composition, C2, is constant, although permissibly varying on a microscopic scale." Such a charge might be a rapidly frozen casting or a mixture of crushed or powdered constituents. Cause a molten zone of

Jan 1, 1953

By J. H. Swisher, E. O. Fuchs

Rate equations were derived to describe the kinetics of internal oxidation of cylinders and spheres. The derived equations for cylinders were checked experimentally by means of sub scale thickness and electrical conductivity measurements on Cu-Cr alloy wires. The properties of the internally oxidized samples were examined with conductivity applications in mind. It was possible to produce uniform dispersions of Cr2O3 in copper with an initial chromium content as high as 3 wt pct. While electrical conductivities only a few pct less than that of OFHC copper were obtained, the Cr2Os particle size and spacing were too large for effective dispersion hardening. T.HE process of internal oxidation has been used widely in basic studies of the permeability of gases in metals. In a review article, Rapp1 has discussed the principles of internal oxidation in considerable detail. From a technological standpoint, internal oxidation is often considered undesirable, since it is a means by which inclusions can be introduced into an otherwise clean material. Another important aspect of internal oxidation is its use as a means of dispersion hardening a material. Broutman and Krock2 discuss this and other methods for making dispersion hardened alloys. The only internally oxidized material known to the authors which is commercially available is a Cu-BeO alloy.3'4 This alloy is made from Cu-Be alloy powder, using a so-called Rhines pack. It has a tensile strength of 80,000 psi and retains its strength at relatively high temperatures. The objectives of the present study were to derive rate equations for the internal oxidation of cylinders and spheres, to check the derived equations for cylinders experimentally, and to examine the structure and properties of internally oxidized Cu-Cr alloys. The Cu-Cr system was chosen for this study because uniform dispersions are obtainable at high alloy contents, which is a desirable characteristic in dispersion hardened materials. RATE EQUATIONS FOR VARIOUS GEOMETRIES A number of authors5--9 have derived equations to describe internal oxidation kinetics. These derivations differ somewhat in mathematical assumptions and approximations, and all except one of the derivations deal exclusively with the internal oxidation of plates. The exception is a brief treatment of cylindrical and spherical geometries given by Meijering and Druy-vesteyn9 as a part of a comprehensive paper on the general subject of internal oxidation. These authors did not obtain rate data to check their derivations, although they did show that the hardness profile across an internally oxidized sample is directly related to the rate of interface movement. For cylindrical and spherical geometries, a quasi-steady-state approximation is needed to circumvent mathematical complications in obtaining a solution to the basic differential equations. In using this approximation, we consider the concentration gradient of dissolved oxygen in the internally oxidized zone or sub-scale to be the same as the gradient which would be present if there were no movement of the subscale interface. The steady-state approximation introduces an error of about 1 pct in computing the rate of internal oxidation of an Fe-1.0 pct Mn alloy plate, if the present method is compared to the more exact method of Wagner.7'10 The details of the derivations of the rate equations for cylinders and spheres are given in the Appendix, and only the results of these derivations are given below. The final equations obtained by Meijering and Druyvesteyn9 can be shown to be equivalent to our Eqs. [1] and [2], although the two approaches are somewhat different. Cylindrical Geometry. [2] where r1 is the outer radius of the cylinder or sphere, cm, r2 is the radius of the unreacted core, cm, see Fig. l(a), D is the diffusion coefficient of oxygen in copper, cm2 per sec, %O is the concentration of dissolved oxygen at the surface of the specimen, wt pct, %Cr is the initial chromium concentration in the alloy, wt pct, and t is the reaction time, sec. Plate Geometry. The analogous rate equation for a plate has been derived previously for internal oxidation of Fe-Al alloys.8'11 For Cu-Cr alloys, we may write the same equation as follows: [3] where r1 is the half-thickness of the plate, cm, and r2 is the distance from the mid-plane to the subscale intherate is An analysis of Eqs. [1], [2], and [3] shows that for a plate the rate is completely parabolic. The initial

Jan 1, 1970

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett', Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,'. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett', Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,'. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By K. J. Gill, P. Chiotti, J. T. Mason

Thermal, metallographic, and vapor pressure data were obtained to establish the pkase boundaries and the standard free energy, enthalpy, and entropy of formation for the compounds in the Y-Zn system. Three coinpounds with stoichiometric formulas of YZn, YZn2, and Y2Zn17 melt congruently at 1105", 1080°, and 890°C, respectively. Four compounds with stoiclziometric formulas of YZn3, YZn4, YZn5, and YZn,, undergo perztectic reactions at 905", 895", 870º, and 685ºC, respectively. Three eutec-tics exisl in this system with the .following eutectic temperatures and zinc contents in wtpct: 875ºC, 23.2 Zn; 1015ºC, 51 Zn; 865ºC, 82 Zn. The YZn, pkase undergoes an allotropic transformation. In the two phase YZn2 -YZn alloys the trans.formation gives a weak thermal arrest at 750°C, whereas in the two phase YZn2-YZn3 alloys no thermal arrest is observed and the transformation occurs over a temperature range below 750°C. At 500°C the free mzergies of formation per lnole vavy from —18,090 for YZn to —53,430 fov YZr11 and corresponding enthalpies vary from -24,050 to -92,080. The free energies and enthalpies per g atom as a function of composition show a maximum for the YZn2 phase; the 500°C values are -9580 and -13,180, vespectively. 1 HE only information found in the literature on Y-Zn alloys was the observation reported by Carlson, Schmidt. and speddingl that Y-20 wt pct Zn forms a low melting alloy. The alloy was produced by the bomb-reduction of YF3 and ZnF2 with calcium in an investigation of methods for producing yttrium metal. The solubility of yttrium in zinc has been determined by P. F. woerner2 and reported by Chiotti, Woerner, and Parry.3 In the temperature range 495" to 685°C the solubility may be represented by the relation In these equations N represents atom fraction of yttrium and T is the temperature in degrees Kelvin. The purpose of the present investigation was to establish the phase diagram for the Y-Zn system and to determine the standard free energy, enthalpy, and entropy of formation for the compounds formed. MATERIALS AND EXPERIMENTAL PROCEDURES The metals used in the preparation of alloys were Bunker Hill slab zinc, 99.99 pct pure, and Ames Laboratory yttrium sponge. Arc-melted yttrium buttons contained the following impurities in parts per million: C-129, N-12, 0-307, Fe-209, Ni-126, Mg-13, Ca < 10, F-105, and Ti < 50. Some of the alloys containing 70 wt pct or more of Zn were prepared from yttrium containing 5000 ppm Ti as a major impurity. Tantalum containers were found to be suitable for all alloys studied and were used throughout the investigation. The pure metals, total weight about 30 g, were sealed in 1 in. diam tantalum crucibles by welding on preformed tantalum covers. A 1/8 in. diam tantalum tube was welded in the base of each crucible for use as a thermocouple well. Welding was done with a heli-arc in a glove box which was initially evacuated and filled with argon. The sealed crucibles were enclosed in stainless steel jackets and heated in an oscillating furnace at temperatures up to 1150°C. Homogeneous liquid alloys were obtained within a half hr at these temperatures except for alloys containing less than 20 pct zinc. The latter alloys were held at 1000º to 1100°C for 2 to 3 days in order to obtain equilibrium. After the initial equilibrations the tantalum crucibles containing the alloys were removed from the steel containers and used directly for differential thermal analyses. Further annealing heat treatments for alloys in which peritectic reactions were involved were carried out in the thermal analyses furnace. After thermal analyses the tantalum crucibles were opened and the alloys sectioned and polished for metallographic examination. In the following discussion alloys referred to as "quenched" were obtained by quenching the sealed stainless steel jacket containing the tantalum crucible and alloy in water. The differential thermal analyses apparatus used was a modified version of the one described in an earlier paper., The graphite crucible was replaced by an inconel crucible, the nickel standard and sampie container were separated by a 1/8 in. MgO plate, no getter was used, and provisions were made to

Jan 1, 1963

By B. L. Landrum, J. Simmons, J. M. Pinson, P. B. Crawford

Patentiometric model data have been obtained to estimate the effect of vertical fractures on the areas swept after breakthrough in water flooding and miscible displacement programs such as gas cycling where the mobility is near one, The data are presented for the case of the fire-spot pattern in which the cemer well is fractured various lengths and orientations, the data indicate that for 10-acre spacing, fractures extetidirrg over 1300 ft in either directior1 from the fractured well may re.srrlt in reductions in sweep efficiencics from 72 to approximately 34 per cent. However. the area swept after break through may be quite largr and only 10 or 12 per cent 1ess than would be obtained if the reservoir were trot fractured. For the specific case when the volume of fluid injected is equivalent to 100 per cent of the pattern vol-unie, the swent area may vary from 80 to 88 per cent, depending on the lenght of the fracture. The former value is that which occurs when the break through or sweep efficiency was orrly 34 per cent and the latter figrrre of 88 per cent is that which is obtained if the reservoir were unfrac-ttm&apos;d. It is pointed out that although the sweep efficiency may he very low in vertically fractured five-spot patterrz.s, the area swept at low water-oil ratios may be only 5 to 10 per cent less than those achieved if the reservoir were unfractured. INTRODUCTION Since the initiation of commercial reservoir fracturing techniques it has been desirable to determine the effect of fractures on the areas swept after breakthrough. Most water flooding or gas cycling projects are continued for substantial periods after the brcakthrough of the injected fluid. Although the sweep efficiency serves as one criterion for rating various flooding patterns. the area swept after breakthrough for various water-oil ratios or percentage wet gas, if cycling. is of perhaps more importance than the sweep efficiency alone. Sweep efficiency data on the vertically fractured five-spot have been presented3. Previous work on the line-drive pattern has shown the effect of vertical fractures on the area swept after breakthrough for the case in which the distance between injection and producing wells divided by the distance between adjacent input wells was equivalent to 1.5 (see lief. 2). The data indicated that for the line-drive pattern it may be desirable to flood or cycle substantially perpendicular to the fractures in order to achieve the greatest recovery for the smallest volume of fluid injected. For this study the center well of a five-spot is assumed as the fractured well. All fractures were assumed to originate at this well and extend into the reservoir for various distances and orientations. All the fractures are straight and are of large permeability compared to the matrix proper. These data are presented to aid the engineer in estimating fractured five-spot pattern performance. ANALOGY The potentiometric model was used in making this study. The model used was 20 20 in. by approximately 1-in. deep. For certain portions of the study one corner of this model was considered to be an injection well and the opposite corner a production well. To simulate vertical fractures a copper sheet was soldered to the wire well and made to conform to the desired length and orientation. In other studies the same model was used except that the four corners of the model might be considered as the corner wells of a five-spot pattern and a fifth well was placed in the center of the model. The well placed in the center of the model was fractured. The total fracture length is L and the well spacing. d. The complimentary fracture angles will be obvious from Figs. 3 and 4. The data obtained on the potentio-metric model assumes the pay to be uniform and homogeneous, the mobility ratio is one, steady-state conditions exist and gravity effects arc neglected. The permeability of the fractures is very great compared to that of the matrix proper. The po-tentiometric model has been used widely both in water flooding and gas cycling projects, and may be used for miscible displacement; how-ever. it is believed that the poten-tiometric model data are more properly applicable to gas cycling than water flooding because the model as-

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.&apos; &apos; Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.&apos;" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-

By J. B. Conway

Stress-vuptuve data for several of- the refractory metals are frequently found to yield a linear relationship between the Larson-Miller parameter and the logarithm of the applied stress. In such cases linear stress-rupture isotherms result with slopes bearing a definite relationship to the temperature. It also follows that the stress to produce rupture in a certain period of time will be linear in temperature. Data for several refractory metals have been reviewed and excellent linearity is shown in this type of isochronal plot. In addition, the af ore - mentioned lineavity leads to a linear relation between the log of the stress to produce rupture in a certain time and the homologous temperature. This has been illustrated for the Group VI-A metals, tungsten and molybdenum. EXTENSIVE use has been made of the Larson-Miller&apos; parameter in the interpolation and extrapolation of stress-rupture and creep data. In those cases where this particular parametric approach is applicable a convenient and fairly straightforward procedure is made available for the correlation of experimental stress-rupture data. It is quite common to employ this parameter in the form of a master rupture plot in which the parameter, T(C + log tr), is expressed as a function of log stress. In many cases this functional relationship in log stress is linear within acceptable accuracy and hence the following relation results: where P is the parameter, C is the Larson-Miller constant, T is the absolute temperature, t~ is the rupture time, a is the stress, and a and b are constants. Examples of such a relationship are contained in the work of Green, Smith, and 01son2 dealing with high-temperature rupture behavior of molybdenum and in the work of Green&apos; dealing with the high-temperature behavior of tungsten. In addition, pugh4 has shown a similar linearity for some fairly low-temperature data for molybdenum. It can be shown that when the relationship in Eq. [I] is exhibited certain generalizations can be made concerning the form of the stress-rupture isotherms. For example, rearranging yields: For a given material (constant C) at a given temperature the first term on the right-hand side of Eq. [2] is a constant and hence this equation defines a straight line when log stress is plotted as a function of log-rupture time. This is recognized as the standard form usually employed in this type of data presentation. Such linearity then suggests the linear form of the Larson-Miller parameter. Or, in other words, the linear parametric relationship in Eq. [2] results only when the stress-rupture data are linear on a log-log plot of stress vs rupture time. Another interesting observation can be made in regard to Eq. [2]. It can be noted that the slope of the stress-rupture isotherms is given by - T/b and hence a direct calculation of the constant b is available. It also follows that since the value of b is the same for all temperatures the slopes of the various isotherms on the log-log stress-rupture plot cannot be the same. Indeed, the existence of the relationship in Eq. [2] precludes a system of parallel lines on this common stress-rupture plot. As a matter of fact it further specifies that in addition to being nonparallel the slope of these isotherms must decrease (i.e., become more negative) with increasing temperature. Such a condition is indeed found to exist in the case of the stress-rupture data reported for molybdenum.&apos; As a corollary to the above, it may be stated that stress-rupture data which do not lead to a linear log-log stress-rupture plot or whose isotherms do not exhibit a decrease in slope as the temperature increases will not yield the linear relationship of Eq. [I]. Applying Eq. [2] to two different temperatures and solving for C yields: Eq. [3] affords a simple and rapid method for calculating the Larson-Miller constant from the log-log stress-rupture plot. The slope of a given linear isotherm is measured and the value of "b" calculated based on Eq. [2] as: slope = - -Tb Then at an abscissa value of 1.0 hr (making log tr in Eq. [3] equal to zero) read the stress corresponding to rupture for two different temperatures. Substitution in [3] yields: A value of the Larson-Miller constant can thus be calculated from a few simple mathematical procedures employing data read directly from the log-log plot of the stress-rupture data. Of course, it is not to be overlooked that the above reasoning has been based on the linear relationship of Eq. [I] being applicable. However, if as mentioned above the log-log plot is

Jan 1, 1967

By I. E. Rosenzweig

GeophyGical prospecting by earth-resistivity methods is frequently applied to investigation of structural problems in geology. Fig. 1 indicates a scheme of the general arrangement used in these methods of measurements. The battery (or generator) E is connected with two "current electrodes" I1 and I2. The current of intensity J, supplied by this battery, passes through the earth between these elcc- trodes. A potential difference U will be measured between two other electrodes V1 and Vz, called "potential electrodes." According to the configuration of electrodes (Fig. I), the resistivity of a homogeneous soil results from the formula p — 2rca-j [1] where 1 C = -------------------------------------- a a 1 a a [2]

Jan 1, 1940

By J. C. Donnell, O&apos

Sinking an inclined ventilation shaft at the Never Sweat mine of The Anaconda Co. has required a unique underground hoisting arrangement because the shaft changes direction from vertical to 71.5º. Access to the project is made through the High Ore shaft to the 2800-ft level, then by crosscut to the bottom of the vertical Never Sweat shaft. The 71.5º incline shaft starts at this point, which makes an unusual operation as the 60-cu ft sinking skip coming up the incline changes into a vertical position above the 2800-ft level before going into the dump scrolls. This requires both rollers and guide shoes on the skip. Coming up the incline to the collar of the shaft the loaded skip enters into a guard before reaching the vertical section. The tendency of the surface- mounted Never Sweat hoist is to lift the front wheels of the skip off the track, so the guard keeps the skip on the rails during the transition.

Jan 11, 1965

By G. F. Ziesing, E. F. Loros, T. J. Smolik

The Golden Sunlight mine began full-scale operations in February 1983 with a conventional open-pit mine supplying 4536 t/d (5000 stpd) of milling ore. After crushing and grinding, the ore is leached with cyanide solution in large agitation tanks. Soluble gold values are removed from the sand portion of the leached ore with five stages of countercurrent washing cyclones. After thickening of the cyclone-wash overflow slurry, the gold is recovered in activated carbon with a six-stage carbon-in-pulp circuit and a five-stage upflow carbon column circuit. Pressure stripping, carbon reactivation, electrowinning, and refining complete the gold recovery process.

Jan 1, 1985

By Robert Glass Cleland

THE DEVELOPMENT of the Phelps Dodge copper and railroad interests began in one era of the company&apos;s history and continued unbroken into the next. Yet the two eras themselves were strikingly dissimilar and represented marked differences in the company&apos;s activities, corporate structure, and personnel. From its organization in 1834 almost to the end of the first decade of the twentieth century, the partnership known as Phelps Dodge was directed, managed, and controlled by the family dynasty established by Anson Greene Phelps. The last two directing partners of this three-generation-old regime assumed control of the company&apos;s affairs in 1879 when the economic revolution that followed the Civil War was well under way.&apos; They lived in an age of gigantic undertakings, gigantic enterprises, gigantic combinations, gi- .

Jan 1, 1952

Feb. 18 to. 21, 1918 Committee on Arrangements J. E. JOHNSON, JR., Chairman BRADLEY STOUGHTON, Vice-Chairman A. R. LEDOUX RALPH W. DEACON F. T. RUBIDGE L. W. FRANCIS E. M. SHIPP MRS. S J. JENNINGS J. V. N. DORR A. L. WALKER G. D. VAN ARSDALE R. M. RAYMOND Finance Committee ALBERT R. LEDOUX, Chairman WALTER H. ALDRIDGE J. H. JANEWAY GEORGE D. BARRON CHARLES F. RAND WALTER DOUGLAS BENJAMIN B. THAYER

Jan 4, 1918

By Elihu Thompson

The subject of welding by electricity has been so recently and so ably treated by Sir Frederick Brsmwell before the Institution of Civil Engineers* as to render evidently superfluous there-statement on this occasion of many facts and considerations which might be appropriately advanced, had they not been so well discussed in his paper already. The electrical arts grow rapidly in these days—a fact which is strikingly illustrated in the development of electric railways. In

Jan 1, 1891

In connection with the presentation of his paper in oral abstract, Mr. Hansell exhibited samples of briquettes, showing their porosity, and the change of surface produced by burning. Alfred H. Cowles, Sewaren, N. J.—I am familiar with this furnace, and I wish to point out that the larger it is, the greater is its efficiency, and the longer its life. In Europe, experiments have been made with small, short furnaces of this character, accommodating narrow cars; whereas Mr. Hansell says that there is in this country a furnace 180 ft. long .and 10 ft wide. The American furnaces, moreover, differ in design from those used in Europe; and we may expect from these larger furnaces not only greater efficiency, but also greater durability —the latter advantage being due to the fact that, in these large furnaces, projecting brick walls of greater thickness can be provided, to prevent the heat from attacking the sides of the cars, by means of bricks which, over-reaching these iron surfaces, furnish a protective cover to them. In other words, in these

Jan 1, 1913

By Henry M. Howe

Preparing the Metals.—Ordinary silver wire was hard-drawn to No. 23 Brown and Sharp gauge; it was cut in two, and one piece was carefully annealed in the ordinary way. The wire was prepared by Messrs. Guild and Delano, jewelers, Boston. The Ansonia Brass and Copper Co., Ansonia, Conn., kindly furnished me with some hard-drawn copper wire No. 31, A. W. G., and some of the same wire after annealing in the ordinary way.

Jan 1, 1885