Search Documents

By L. W. Fagg

The purpose of this paper is to present in a convenient form data and examples necessary in making dynamometer card analyses; also to outline a procedure of well weighing. Many articles and papers have been written delving into the mathematical considerations relative to the shape and characteristics of dynamometer cards. However. it is recognized that there are too many unknown factors involved in such calculations to assure a workable degree of accuracy. For this reason the accepted procedure is to take dynamometer cards on wells in question rather than try to calculate the load curve. The polished rod dynamometer is now recognized as a necessary tool for measuring loads, torque, and horsepower. It is also used to determine pump action and trouble-shoot for any seemingly abnormal pumping condition. The apparently infinite variety of a = maximum load (height x scale constant) b = minimum load. Range of load is difference between maximum and minimum load speed—taken with stop watch stroke —measured at polished rod c = beginning of down stroke in direction of arrow. (End of up stroke) d = beginning of up stroke. (End of down stroke). Polished Rod Horsepower = (Area of card) x Scale const. x Stroke x Length x spm (Length of card) 33000x 12 1.35 1.35------x5300x 24 x 12 2.53 --------------------------= 2.06 23000x12 Appraximate Peak Torque: Upstroke = (2390) (12) (1) = 28,700 in. Ib Downstroke = (1 860) (12) (.866) = 19,300 in. Ib Counterbalance should be increased to make up stroke and down stroke peak torque equal. FIG. 1 dynamometer cards that can be obtained is one reason for the general lack of usage of the dynamometer as a control instrument rather than a means for making routine measurements of leads and horsepower. When it is considered that the dynamometer card is a record of the resultant of all forces acting on the polished rod at any particular instant during the pumping stroke. the problem is then one of breaking down this resultant into its various components. As a means of a quick review we shall consider the examples shown in Figs. 1 to 9: and Tables I and II and then proceed to the interpretation of variously shaped cards caused by some abnormal operating condition. In Table 11, when we were considering the factors involved in calculating the peak polished rod load, it can be seen that the factors involved greatly oversimplify the problem. Certain assumptions are made which may or may not he even close to the actual field conditions, such as the specific gravity of the fluid generally considered as one; that the crank has constant angular velocity; that the down-hole friction is zero; and that the fluid lift is from the pump. In the following examples we shall see what a variation in fluid weight and friction can do to the general shape and magnitude of the dynamometer card. (Figs. 10 to 22.) MAKING THE WELL STUDY It is obvious that it would be impractical to consider in detail all of these factors each time a well study is made, inasmuch as each well study job could conceivably be extended into a research project rather than serve the practical requirements of finding the answer to a specific problem. For this reason it is important that some objective be established previous to the time the well study is made. Load at TU = 3070 Ib Crank angle when polished rod is at position TU = ? ? = 30° Maximum counterbalance effect at polished rod = 2760 Ib Torque at TU = (Load at TU — Max. counterbalance effect-lbs) sin 0 x Length of Stroke 2 = (3070-2760) .5 x 24 = 1860 in. 15 2 Torque at TD = (Load at TD — Max. counterbalance effect-lbs) sin ? x Length of Stroke 2 = (530 - 2760) sin 330' x 24 = 13,400 in. Ib 2 Note: sin 0 from c to d on UPSTROKE will be positive value. sin 0 from c to d on DOWNSTROKE will be negative value Torque at c and d is zero because 0 is zero. FIG. 2 —APPROXIMATE METHOD FOR CALCUlATlNG TORQUE

Jan 1, 1950

By V. T. Stringfield, H. H. Cooper

ONE of the earliest investigations of ground water in Florida was made in 1513 when Ponce de Leon arrived at St. Augustine in search of the Fountain of Youth. The history of the development of the water resources of the State shows that the large artesian reservoir that underlies Florida was discovered in the latter part of the last century. Part of that history is given by L. C. Johnson' who states that the first successful artesian well in Florida was drilled at St. Augustine between 1880 and 1882. After the City of Jacksonville failed to obtain a flow of artesian water at a depth of nearly 400 ft and abandoned the drilling, R. N. Ellis, City Engineer, and L. C. Johnson, using Johnson's knowledge of the geologic structure of the artesian reservoir, estimated correctly that artesian water could be obtained at a depth of about 500 ft. Thus began the development of the large artesian system in the northeastern part of the State. In 1908, the Florida Geological Survey issued its first report on the ground water of central Florida. More recent reports give the results of investigations that have been in progress in cooperation between the Florida Geological Survey and the U. S. Geological Survey since 1930. As a result, the ground-water geology and hydrology of Florida are now so well known that ground-water problems such as confronted early investigators no longer exist. However, new problems arise and new discoveries are made as the demand for more ground water increases with the development of the State. Water-Bearing Formations Descriptions of the geologic formations and a map showing their distribution at the surface are given in a recent report by Cooke.' The geologic formations that yield the ground-water supplies in Florida represent only a small part (about 1000 ft) of the total thickness of the sedimentary rocks (more than 15,000 ft) that underlie Florida. The water-bearing rocks that yield fresh water include more than two dozen formations that range in age from Eocene to recent. The Eocene formations, which consist chiefly of limestone, are the oldest and the deepest of the formations in Florida that yield fresh water to wells. In 1944, on the basis of a study of the foraminifera, the Applins -ivided the limestone of Eocene age into six formations, which are as follows, from top to bottom: Age or group Formation Jackson Ocala limestone icAvon Park limestone Claiborne -j Tallahassee limestone aLake City limestone Wilcox Salt Mountain limestone xOldsmar limestone The Ocala limestone and the underlying formations of Claiborne age, along with some of the overlying limestone of Oligocene and Miocene age, constitute the principal source of water in Florida and southeastern Georgia and generally may be regarded as forming one artesian aquifer or water-bearing unit. This aquifer is referred to in this paper as the Floridan aquifer, a name proposed by Parker.' In part of Seminole County alone more than 1000 flowing wells yield water from the Floridan aquifer. A yield of 6500 gpm, or about 9.5 million gal per day, by natural flow from a well penetrating the aquifer at Jacksonville, Florida, was observed in 1942. A yield of about 7000 gpm by natural flow was reported for a well 1390 ft deep at St. Augustine in 1887. The largest yield reported from a pumped well penetrating the aquifer is 7500 gpm, or about 10.5 million gal per day. The aquifer is also the source of some of the large springs, such as Silver Springs, whose discharge, according to measurements made by the U. S. Geological Survey, has ranged from 526 to 1350 sec-ft, or from 340 to 872 million gal per day. The top of the Ocala limestone, as represented by contours in Fig. 1, indicates in a general way the geologic structure of the formations that comprise the Floridan aquifer. The Ocala is at or near the surface in the areas represented by shading in the northwestern part of the peninsula and also in an area in western Florida adjacent to Alabama and

Jan 1, 1952

By H. A. Hoffman

Competition from an all-flotation plant, with demonstrated economies and efficiencies, plus a change in smelting contract and introduction of improved cyclones lead to conversion from gravity-flotation. Detailed descriptions are given of equipment installed and procedures used at St. Joseph Lead's Federal mill. Future plans include further quality control by instrumentation in all aspects of the mill operation. Also planned are improvements in materials handling systems. Competitive conditions will continue to dictate improvements and changes in the mill flowsheet. The advent of all-flotation mills in the Lead Belt was introduced by the Indian Creek mill in late 1953. This modern and efficient plant then became a pattern for the other mills in this area to re-evaluate their circuits in an effort to develop flowsheets that would improve operating conditions and metallurgy. The Indian Creek mill demonstrated that all-flotation would require considerably less operating and maintenance labor than the combination of tabling and flotation that was common in the Lead Belt at that time. Two significant changes occurred in recent years to allow all-flotation to be seriously considered in other Lead Belt mills. One was a change in the smelting contract which did not require gravity concentrate. Another was the development of cyclones which provided classification of the flotation feed in a very small space. Since the Federal mill was the largest concentrating plant in the Lead Belt it was felt that the greatest savings could be made by investigating the all-flotation possibilities at this mill. Interest was stimulated by the milling and ore dressing Depts. to determine if this mill would be converted without undue cost, provided the metallurgy could be improved. Accordingly, laboratory tests were initiated to learn what metallurgical benefit could be derived. Numerous tests were run which indicated that grinding to all-flotation would improve the tailing by as much as 0.02 pct Pb. This was quite significant when multiplied by the tonnage of ore treated. Projecting the added cost of power for extra grinding and flotation, and the additional flotation reagents required, plus additional new equipment that would have to be purchased, it would still add up to a considerable saving. The all-flotation mill would reduce man power by some 30 pct, and make a metallurgical improvement. Operating costs could be reduced bv 2$ per ton. The metallurgical improvement would amount to 4$ per ton, for a total of 6$ per ton of ore milled. An estimate of the cost of the conversion was set at $250,000. With the savings as estimated this could be paid off in about two years. On paper it therefore seemed attractive enough to justify a mill test. MILL TESTING One section of the Federal mill was made available for use as a separate test circuit. Cyclones and density controllers were borrowed from the Viburnum and Indian Creek mills. Denver flotation cells were obtained from the Desloge mill that was then shut down. Two 10-cell groups of Fagergren flotation machines, consisting of eight roughers and two cleaners in each bank, were segregated from the remainder of the plant. The 9x12-ft rod mill was continued as a primary grinding mill. The two 6%x12-ft mills on this section were converted from rods to balls, and the speed increased to 22.0 rpm. Cyclones were installed to close the circuit of these ball mills. The ground pulp from each of the ball mills was fed separately to its own bank of flotation cells so that two completely separated test circuits could be run in parallel. This test circuit operated on a three-shift basis from Dec. 22, 1959 to Mar. 21, 1960. The first few weeks were occupied in developing conditions for proper metallurgy and trouble-free operations. Results were rather erratic but as the ore dressing laboratory and the mill operators became more familiar with conditions they were able to obtain expected results, which could be duplicated day after day. The last five weeks of testing indicated the best results and are tabulated in Table I. During the first two weeks a Denver unit cell was used on No. 1 flotation circuit. Even though it recovered almost one half of

Jan 1, 1962

By J. P. Hirth, A. H. Clauer

A simple graphical method is presented for the determination of climb forces on dislocations exerted by uniaxial stresses. In conjunction with standard stereographic projections, the technique is applicable for any crystal system. Climb forces on glide dislocations in bcc crystals are determined as an application. In connection with a study1 of elevated temperature creep of molybdenum single crystals, a detailed analysis of the climb force on a dislocation was necessary. In general, such forces can be calculated for an arbitrary dislocation segment from the Peach-Koehler equation2 (or one of its variants3"6) where Fk are the components of the force per unit length of dislocation line, ?i of the unit sense vector tangent to the line, ojl of the stress tensor, bl of the dislocation Burgers vector, and is the Einstein permutation operator with components ?123 = ?231 = ?321 = —?132 = — ?321 = —?—213 = unity, and all other components zero. It is understood in the suffix notation of Eq. [I.] that repeated indices are to be summed over 1, 2, and 3. For the most general coordinate system, the matrix manipulations required to solve Eq. [I] are fairly complex. However, a simplification can be achieved by specifying a coordinate system based on the dislocation in question. The pertinent system involves cartesian coordinates xi, x;, xi, with attendant unit vectors e'i; ei being parallel to the Burgers vector b, and e'3; being parallel to the climb directions (b x ?), i.e., normal to the slip plane. In this coordnate system b = (bi, 0, 0) and ? = (?'1,?'2, 0). Hence, Eq. [I] yields for the climb force per unit length F3 = ?'2s'11b'1—?'1s'12b'1 [2] The first term on the right side of Eq. [2] is the climb force on the edge component of the dislocation, while the second is the (cross-slip) force normal to the glide plane acting on the screw component. Thus, the problem of determining the climb force reduces to that of resolving the stress components —il and s'12 in the xi coordinate system. Several authors7- l0 have shown, for the glide case, that such stress resolution can be performed most readily by graphical methods. Therefore, we present a similar method for the climb case. Of the suggested methods, we follow that of Hartley and Hirth,9 which can be performed directly with a standard stereographic projection for the crystal system in question. We treat the case of simple tension, or compression, this being the most common creep test in which climb is important. A set of cartesian coordinates (x1, x2, x3) is fixed with x3 parallel to the specimen axis. This set is connected to the x'i set by the factors defined in Table I. A perspective view of the coordinate systems relating them to the glide ellipse, is given in Fig. 1, while a stereographic projection of the coordinate systems is presented in Fig. 2. The transformation matrix connecting the xi and xi coordinates is x'i =aijXj [3] where an a11 a13 aij = a21 a22 a23 a31 a23 [4] s The corresponding transformation for the stresses is i [5] The only nonzero component of okl for the simple tension case is s33,. Hence, using Eqs. [4] and 151, we find that Eq. [2] becomes F3 = b'1s33 sin2 cos ?[cos ? sin a - sin ? cos a] [6] = m1b'1s033 sin a- rn2b's33 cos a where sin a = 5: and cos a = ?i, with a the angle measured from b to ?, i.e., a polar angle in the x1, x2 plane. The same equation with a sign change applies for compression. A graphical plot of the "Schmid" factor m1, relating to the climb force on the edge component, is presented in the stereographic projection of Fig. 3. In this projection the glide plane normal points toward the top of the page and the glide direction points normal to the plane of the page. To use the graph, one

Jan 1, 1970

By W. A. Sheaffer

Electric furnace installation and tough-pitch copper-casting operation at Canadian Copper Refiners Ltd. are described. General layout, power supply and control, refractories, induction pour hearth, casting equipment, metal temperature control, and oxygen-content control are discussed. E LECTRIC furnace at Canadian Copper Refiners Ltd. was put into operation in August 1949. The installation was designed primarily to melt electrolytic copper cathodes and to produce vertically cast tough-pitch shapes. To meet emergencies during reverberatory furnace shut-downs, provision was made for casting horizontal tough-pitch wire bars. Due to the ease with which metal temperatures, melting rates, oxygen content, and casting hours can be altered to suit production demands, the electric furnace has met the requirements admirably. General Layout The plant has been described by H. S. McKnight' and by J. H. Schloen and E. M. Elkin.2 Since these papers were written, the electric furnace has been installed in a 140 ft extension of the original casting building. The extension, 264 ft in width, is divided into one 24 ft and four 60 ft bays, each a continuation of similar bays in the older building. Fig. l, a partial floor plan of the extension, shows the location of major equipment. Power Supply and Control Power for the electric furnace and its auxiliary equipment is delivered by Quebec Hydro over two 12 kv three-phase, 60 cycle pole lines. One line is in use, while the second serves as an emergency standby. Connections for the electric furnace from these lines, which also feed the refinery power house, are made to a substation on Canadian Copper Refiners Ltd. property. Connections between this substation and the furnace transformer room are in underground conduit. The are-furnace transformer, 550 v auxiliary-equipment transformer, and 12 kv switchgear are in the transformer room near the furnace. Stepdown from 12 kv to are-furnace voltages is done by a 7500 kva oil-immersed water-cooled three-phase 60 cycle transformer. This unusually high kva rating is available because the arc furnace and transformer were adapted from a steel-melting unit. Secondary voltages are available in 16 steps between 251 and 95 v (across phases). A four-position tap changer, operated from the furnace switchboard is connected to taps yielding 199, 164.5, 127, and 115 v, respectively. The transformer primary feeder is equipped with a General Electric Magne-Blast circuit breaker. The pour hearth, casting wheel, bosh conveyor, and other auxiliary equipment are fed by a 12 kv to 550 v, 750 kva three-phase 60 cycle transformer. Graphite electrodes, 14 in., with tapered nipples are held in the are-furnace electrode arms with steel wedges. The arms, fed from the transformer secondary busbars by flexible cable bundles, are actuated by winch cables. Direct-current motors operating the winches are supplied by a 250 v motor-generator set. Power input to the furnace is controlled from the furnace switchboard. The board has the usual combination of remote controls for the 12 kv breaker and transformer tap changer, meters, overload relays, automatic and manual electrode motion controls, switchgear, etc. When power is on the furnace, automatic electrode feed is used; the power draw on a given voltage tap then is governed by current-limiting rheostats. The rheostats work through a Westinghouse bal-lanced-beam control circuit which, in turn, operates reversing contactors in the winch motor circuits. Arc Furnace The arc furnace is a standard-type NT Pittsburgh 'Lectromelt with inside shell diameter 12 ft 3¾ in. By hydraulic-lift cylinders, the furnace can be tilted forward to 39 ½ ° and backward to 5" from horizontal. A roof-swing cylinder is also installed but not used. Fig. 2 shows refractory-lining details. With the lining as shown, operating life between repairs is 6 to 9 months of two-shift casting operation. Charge slot, skim door, and roof refractories then are replaced completely, and any necessary repairs made to side walls. Bottom life is longer and the original under courses are still in service, while the top course was replaced in October 1951. The launder is lined with high chrome-magnesite trough brick backed by fireclay and insulating brick. The trough is covered with fireclay brick and asbestos sheet. At intervals, openings 12x6 in. are left in the permanent covering. These are covered with

Jan 1, 1956

By M. C. Flemings, R. Mehrabian, G. E. Nereo

Analytic expressions were developed and applied in two preuious papers to predict effects of solidification variables on macrosegvegation. In this paper, experiments are reported to test the analyses quantitatively and qualitatively. Good quantitative agreement is obtained between predicted and measured compositions for the following types of segregation (all experiments were on laboratory-size ingots of Al-4.5 pct Cu alloy): (a) inuerse segregation, (b) composition distribution in ingots solidified with unidirectional heat and fluid flow, (c) simulated centerline segregation, and (d) segregation resulting fro, change in solidification cross section. In addition, qualitative agreement of prediction of theory with experiment is obtained for (a) segregation resulting from sudden change in heat extraction during solidification, (b) centerline segregation in a bidirectionally solidified ingot, and (c) under-riser positive segregation. IN two previous papers,''2 it was shown that a variety of apparently different types of macrosegregation result from the same basic mechanism. This mechanism is the flow of solute -rich liquid to feed solidification and thermal contractions. Analytic expressions were derived and examples given of their application to predict macrosegregation. Calculations were for Al-4.5 pct Cu alloy. In this paper, we describe a series of experiments designed to test results of the foregoing calculations. All experiments were on laboratory-size ingots (weighing 10 lb or less), of A1-4.5 pct Cu alloy, nominal composition. APPARATUS AND PROCEDURE A typical mold and chilling arrangement used for unidirectional solidification is shown in Fig. 1. Before casting, this assembly was placed in an electrically heated air furnace, and heated to 680" to 700°C. When the plaster reached the desired temperature (usually 1 to 2 hr after the furnace temperature), the mold was filled through a funnel extending from the outside of the furnace into the cavity. The assembly, with liquid, was held for several minutes to allow convection to subside. The coolant, either air or water (depending on the desired solidification rate), was turned on and the ingot solidified. With the four sides and top insulated, heat removal for solidification was through the bottom surface, resulting in a fully columnar structure. Design of the casting cavity was varied, as discussed below, to produce desired macrosegregation effects. However, for castings solidified with unidirectional heat flow, the top was always open to the furnace atmosphere and temperature. The bottom sur- face rested on a chill through which the heat was extracted, and all other surfaces were insulated with plaster. In one case to be discussed, the mold design was changed to obtain bidirectional solidification. In this case, two vertical chills were used (on two faces of the mold) and the bottom was insulated. In a last case discussed, that dealing with under-riser segregation, experiments were conducted in bottom-chilled sand molds at ambient temperature. Except for this one case, all experiments were conducted in the electrically heated air furnace as described above. In many ingots, thermocouples were used to obtain continuous records of solidification. In all cases, these were chromel-alumel couples with output continuously monitored on a twelve-channel recorder. Where thermocouples were employed, these were embedded in the refractory mold walls during mold-making so their heads extended to the vertical centerline. Six to nine couples were employed, placed at intervals from the chill. METAL-MEL TING PROCEDURE Aim analysis was A1-4.5 pct Cu; alloy was prepared from high-purity aluminum (99.9 pct Al) and A1-Cu master alloy (50 pctOFHC Copper, 50 pct high-purity Al). The A1-4.5 pct Cu alloy was cast in 5- to 7-lb ingots for subsequent remelting for the macrosegregation experiments. Melting for ingot casting was as follows. After preheating the melting crucible for an hour or more, the charge was melted and superheated to a temperature of 690" to 720°C. After stirring, the melt was degassed with chlorine for 10 min. After stirring again for 1 min, a gas sample, under reduced pressure, was taken. When the sample was gas-free, the crucible

Jan 1, 1969

By H. E. Cline, E. F. Koch, J. L. Walter

It has been postulated and, in a few instances shown, that some kind of dislocation structure will be present at semicoherent interfaces to accommodate small lattice mismatches. In the present study of the NiAl-Cr eutectic, regular arrays of interface dislocations are observed at the boundary between the chromium-rich rods and the NiAl-rich matrix. The networks were examined by transmission electron microscopy and selected area diffraction. The rods and the matrix have a crystallographic relationship in which all directions and planes of the two phases are parallel. The dislocation networks are cmposed of a<100> dislocations lying on the intersections of the cylinders with (100) planes. Dislocations forming hexagonal rather than square arrays are also observed at certain areas of the network. The morphology of the network is consistent with the interpretation of mismatch being accommodated by interface dislocations in the cylindrical geometry. The measured spacing between dislocations was used to calculate an apparent lattice mismatch between the phases (˜0.35 PCt)interface network energy (-140 ergs per sq cm), and network strengthening (-10,000 psi). It has been proposed by Frank and Van der Merwe1 that dislocations should be present at the boundary between two semicoherent crystallographically related phases. The role of the interface dislocations would be to reduce the internal stresses, caused by the mismatch in atomic spacing across the interface. Such dislocations have been observed at the interface between expitaxially grown films and Substrates.2-4 Interfacial dislocations have also been observed at precipitate-matrix interfaces.&apos;-&apos; Directionally solidified eutectics have been shown to have semicoherent phases1&apos; and would, therefore, be expected to have interfacial dislocations as found by Weatherly at a lamellar fault in A1-A12Cu.11 The NiAl-Cr eutectic appeared to be a promising system to examine because both phases are cubic, the lattice mismatch is small, and the phases are crys-tallographically related. Furthermore, the eutectic is easily thinned for transmission microscopy. Indeed, interfacial dislocations were observed and this report describes the nature of the dislocation networks in the boundary between the NiAl-rich ma-trix-phase and the fine chromium-rich rods.21 I) EXPERIMENTAL PROCEDURE Ingots, 3/4 in. in diam and 6 in. long were made by melting in vacuum and casting under argon using 99.9 pct pure material. The composition, in at. pct, was 33 pct Ni, 33 pct Al, and 34 pct Cr. The ingots were then placed in A1203 crucibles on a water-cooled base, melted by means of induction, and withdrawn from the hot zone at the rate of 1 in per hr under argon.* * T his material was first directionally solidified in this laboratory by E. R. Stover Slices were taken perpendicular to the growth direction of the directionally solidified ingot for metallography and for transmission electron microscopy. The electron transmission samples were thinned mechanically, then thinned electrolytically in A-2 electrolyte* *A-2 electrolyte: 62 ml perchloric acid, 700 ml ethanol, 100 ml butylcellosolve, 137 ml distilled H20. until a hole appeared in the foil. 11) EXPERIMENTAL RESULTS A) Optical Microscopy. The microstructure, viewed on a plane perpendicular to the growth direction, is shown in Fig. 1. The structure consists of cells or colonies of parallel chromium-rich rods in the NiAl matrix. The cells occur when there are impurities present12 or, in a ternary eutectic, if the composition is slightly off the eutectic composition. The axis of the chromium-rich rods is parallel to the growth direction except near the cell boundaries. Here the rods may assume angles to the growth direction; however, examination shows that the crystallographic relationship between the rod and the matrix remains the same. Fig. 1 includes cell boundaries where the rods formed at a large angle to the growth direction. The variation of rod position across the cells made it possible to Fig. 1-Structure on plane perpendicular to growth direction. Rods near cell walls are at large angle to growth direction. Magnification 315 times.

Jan 1, 1970

By J. T. Norton, A. L. Mowry

The monocarbides of the A subgroup elements in the fourth and fifth group of the periodic table in addition to being hard and refractory are of special interest in that they are isomorphous in crystalline structure. They are cubic with a sodium chloride type structure in which the metal atoms are essentially close packed in a face-centered cubic arrangement with the carbon atoms placed in the interstices between. Interstitial structures of this close packed type were first investigated systematically by Egg1 and he gave the rule for their formation, stating that the radius ratio of the nonmetal to the metal atom should not exceed the value of 0.59. The carbides of interest are those of titanium and zirconium of the fourth group and vanadium, columbium and tantalum of the fifth group. Table 1 shows the radius ratio using the Goldschmidt radii for 12 coordination for the metal atoms and the diamond radius for the carbon atom. It will be noted that while there is considerable variation in the size of the metal atom, in all cases the ratio is smaller than the limit of 0.59 placed by Hägg. It has been known for some time that these cubic carbides are soluble in one another, at least to some extent or, in other words, the metal atoms can be replaced, one by another without destroying the stability of the structure. Since the stability of these close packed interstitial substances appears to depend more upon geometry than upon the exact chemical nature of the atoms involved, it is of interest to examine the possibilities of replacement in these carbides in some detail. Hume-Rothery2 has pointed out the importance of the difference in size of solute and solvent atom as a factor in limiting the solubility in simple binary solid solutions. Largely on an empirical basis, he states that if the difference in size between solvent and solute atom is more than 14-15 pct of the solvent atom, the range of solubility is very restricted. The atom size was based on the distance of closest approach in the elements involved. While there is some question as to how one should calculate the size of the metal atom in the carbide structures, reference to Table 1 will show that zirconium is the largest and vanadium the smallest of the group and that the difference is about 15 pct. The Ti-Zr difference is about 9 pct and the others are smaller. Thus one would predict that if the size factor controls the solubility, all of the pairs except VC-ZrC would have wide or complete solubility whereas this latter pair is on the border line and might have restricted solubility. The purpose of the present investi- gation was to examine the solubility of the several pairs of carbides by heating them together until equilibrium was established and then examining the product by X rays. Previous Work Agte3 and his associates prepared various transition metal carbides and determined the melting points of binary mixtures. He concluded from the shapes of the melting point curves that there was extensive solubility in the case of the cubic carbides. Umanskii and his colleagues made an investigation of a number of pairs of the cubic carbides, using X rays and plotted lattice parameter vs. composition curves for the systems TaC-Tic, CbC-Tic, TaC-ZrC and CbC-ZrC. All pairs showed a continuous series of solid solutions. The first two pairs gave a linear relation while the latter two showed a negative deviation from Vegard&apos;s law. Kiefer and Nowotny, in a paper which became available after the present work was well advanced, investigated the binary pairs of the five cubic carbides by means of X rays. Relatively few points were obtained and results indicated that in some cases, at least, equilibrium was not reached at the temperatures used. The results indicated that solubility in the VC-ZrC system was not complete. All of the results of previous investigations indicated the desirability of a more detailed study. Materials The raw materials used were mono-carbides of titanium, zirconium, vanadium, columbium and tantalum and were the purest which could readily be obtained commercially. Spectrographic qualitative analysis showed that the CbC and TaC contained less than 1 pct

Jan 1, 1950

By R. H. Olds, B. H. Sage

The transient flow of gases in long pipe lines is a problem of industrial interest. The present discussion deals with the application of the conservation of momentum and material to the transient flow of gas in conduits which are long with respect to their diameter. Solutions of the resulting equations include consideration of variations in friction and specific volume of the fluid as functions of time and position along the conduit. The changes in state along the pipe line are established as a function of time. The methods proposed may be extended to other situations. The flow of compressible fluids under steady conditions&apos; has been considered in detail. Primary emphasis was placed on the development of simple means of taking into account the non-uniform flow&apos; realized in long conduits under such conditions. Binder2 discussed the essential differences between the flow of compressible and uoncompressible fluids. Bonilla" proposed direct methods of solution for equations describing isothermal flow in long pipes. Similar considerations limited to steady flow have been discussed by Ruth.&apos; Palsgrove 5 sugrested a graphical means of solving specific problems relating to the flow of gases in conduits of uniform section. Some of the basic concepts of the flow of a compressible viscous fluid in straight tubes have been described by Kuo.6 Problems of this type require knowledge of the equation of state of the fluid. In many cases it liar been convenient to assume that the fluid follows the relationship ascribed to a perfect gas. The treatment of unsteady flow of compressible fluids is less complete. This state of affairs result* from the fact that the differential equations relating pressure and flow rate with time and position for unsteady flow are nonlinear in character.&apos; Muskat and Wyckoff presented an excellent treatment of the situation in laminar Bow where inertia forces were not of controlling importance. Somewhat earlier Hurst9 discussed the uncteady flow of fluids in petroleum reservoirs. Helhering-ton. MarRoberts. and Huntington&apos;" established experimentally the pressure distribution and accumulative flow rate of natural gas through an unconsolidated sand. Recently Joffe11,12 conhitlcretl the specific problem of the flow of natural gas in relatively long.,. pipe lines for unsteady conditions with emphasis upon the economic feusibility of the storage of gas in such a reservoir. These discussions were based upon approximations of the partial differential equations applying to the condititons of flow. The results were descriptive of the physical situation and indirated the effertiveness of large pipe lines as cyclic storage facilities for natural gas. The advent of analog computers" has increased materially the practicality of solving nonlinear partial differential expressions. Several iterative numerical and graphical methods14,15,16 have been developed to solve nonlinear equations. However, it appears that the pair of simultaneous nonlinear partial differential equations involved in problems of this nature are not directly susceptible to solution by presently available analog computer or formal iterative numerical techniques. Nevertheless. graphical analyses afford a means of obtaining data of engineering value without special equipment. MATHEMATICAL CONSIDERATIONS Since the length of the conduit is large in comparison to the cross section in the cases to be considered, it is possible to treat the characteristics of the flow as independent of position in the cross section of the conduit. This possibility follows since the energy- connected with the velocity profile is small compared to that required to overcome the friction associated with its movement throughout the length of the line. Average values of the pressure. velocity. kinetic energy. momentum, and specific volume may he assigned to the fluid at any particular cross section. From a mathematical viewpoint this may be expressed in the following way for velocity: µ = µ(?X)..........(1) The uncertainties introduced by such assumptions usually are negligible when the ratio of length to diameter is greater than the order of 1.000. On this basis a volume element ma!. be considered to occupy the entire cross section of the conduit. Such an element is presented in Fig. 1.

Jan 1, 1951

By R. H. Olds, B. H. Sage

The transient flow of gases in long pipe lines is a problem of industrial interest. The present discussion deals with the application of the conservation of momentum and material to the transient flow of gas in conduits which are long with respect to their diameter. Solutions of the resulting equations include consideration of variations in friction and specific volume of the fluid as functions of time and position along the conduit. The changes in state along the pipe line are established as a function of time. The methods proposed may be extended to other situations. The flow of compressible fluids under steady conditions&apos; has been considered in detail. Primary emphasis was placed on the development of simple means of taking into account the non-uniform flow&apos; realized in long conduits under such conditions. Binder2 discussed the essential differences between the flow of compressible and uoncompressible fluids. Bonilla" proposed direct methods of solution for equations describing isothermal flow in long pipes. Similar considerations limited to steady flow have been discussed by Ruth.&apos; Palsgrove 5 sugrested a graphical means of solving specific problems relating to the flow of gases in conduits of uniform section. Some of the basic concepts of the flow of a compressible viscous fluid in straight tubes have been described by Kuo.6 Problems of this type require knowledge of the equation of state of the fluid. In many cases it liar been convenient to assume that the fluid follows the relationship ascribed to a perfect gas. The treatment of unsteady flow of compressible fluids is less complete. This state of affairs result* from the fact that the differential equations relating pressure and flow rate with time and position for unsteady flow are nonlinear in character.&apos; Muskat and Wyckoff presented an excellent treatment of the situation in laminar Bow where inertia forces were not of controlling importance. Somewhat earlier Hurst9 discussed the uncteady flow of fluids in petroleum reservoirs. Helhering-ton. MarRoberts. and Huntington&apos;" established experimentally the pressure distribution and accumulative flow rate of natural gas through an unconsolidated sand. Recently Joffe11,12 conhitlcretl the specific problem of the flow of natural gas in relatively long.,. pipe lines for unsteady conditions with emphasis upon the economic feusibility of the storage of gas in such a reservoir. These discussions were based upon approximations of the partial differential equations applying to the condititons of flow. The results were descriptive of the physical situation and indirated the effertiveness of large pipe lines as cyclic storage facilities for natural gas. The advent of analog computers" has increased materially the practicality of solving nonlinear partial differential expressions. Several iterative numerical and graphical methods14,15,16 have been developed to solve nonlinear equations. However, it appears that the pair of simultaneous nonlinear partial differential equations involved in problems of this nature are not directly susceptible to solution by presently available analog computer or formal iterative numerical techniques. Nevertheless. graphical analyses afford a means of obtaining data of engineering value without special equipment. MATHEMATICAL CONSIDERATIONS Since the length of the conduit is large in comparison to the cross section in the cases to be considered, it is possible to treat the characteristics of the flow as independent of position in the cross section of the conduit. This possibility follows since the energy- connected with the velocity profile is small compared to that required to overcome the friction associated with its movement throughout the length of the line. Average values of the pressure. velocity. kinetic energy. momentum, and specific volume may he assigned to the fluid at any particular cross section. From a mathematical viewpoint this may be expressed in the following way for velocity: µ = µ(?X)..........(1) The uncertainties introduced by such assumptions usually are negligible when the ratio of length to diameter is greater than the order of 1.000. On this basis a volume element ma!. be considered to occupy the entire cross section of the conduit. Such an element is presented in Fig. 1.

Jan 1, 1951

By T. A. Pollard, C. R. Sandberg, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler&apos;, and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl&apos; have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By C. R. Sandberg, T. A. Pollard, E. B. Elfrink

This paper presents an evaluation of compressibility factor data and a discussion of their application to the estimation of gas reserves. A correlation is presented which provides compressibility factors for use in both two-phase and single-phase hydrocarbon systems. Accuracies comparable to those obtained previously for single-phase systems only can be expected. A simple means of predicting the presence or absence of a liquid phase in a condensate system of known composition is illustrated. The correlation is based on 1,030 compressibility determinations from 21 hydrocarbon samples taken from eight oil fields. Of the data used, 75 per cent were from California, 15 per cent were from the Mid-Continent area, and 10 per cent were from South America. The average numerical deviation of the experimental data from this compressibility chart is 1.22 per cent. Charts and tables are included and discussed which illustrate the errors involved through the misuse or nonuse of compressibility factors in estimates of gas-in-place. INTRODUCTION The compressibility factor is a coefficient which expresses the deviation of a gas of given composition from the Perfect Gas Laws. A factor such as this is necessary in any calculation involving volumes of gaseous mixtures, and finds extensive application in the estimation of natural gas and condensate reserves. It is used in the decline curve, the volumetric, and the material balance methods of gas reserves estimation. The behavior of gases at high pressures has been investigated extensively in recent years and considerable data have been assembled on the compressibility characteristics of a number of gases. Notable among the investigators of high pressure gas relationships are Kvalnes and Caddya, Kay" Sage and Laceys, and Brown, Standing" and Katz4. Other contributors are Smith and Watson, Roland and Kaveler&apos;, and Stevens and Vance". The data on compressibility factors assembled by these investigators have been presented and correlated in numerous ways; particularly noteworthy is the method first presented by Kay5 which relates compressibility factor "Z" to pseudo-reduced temperature and Pseudo-reduced pressure. Other investigators have used this empirical relationship successfully and numerous charts covering a large number of gases have been constructed on this basis. It is the purpose of this paper to evaluate compressibility factor data and to illustrate the importance of their proper use in the estimation of gas-in-place for both gas and condensate reservoirs. EVALUATION OF COMPRESSIBILITY FACTOR DATA The compressibility chart presented by Standing, Katz, Brown and Holcomb in 1941 is probably the most recent and reliable chart for natural gases. The reliability of the Katz chart has been investigated by the writers for a fairly large number of gases (single-phase) in the pressure range from 500 to 10,000 psia and temperature range from 80 to 400°F. A summary of this comparison is presented in Table I. The numerical average deviation of the chart compressibility factors from the experimental compressibility factors is 1.05 per cent, and the algebraic deviation is -0.28 per cent. The Katz correlation appears to be quite accurate for single-phase hydrocarbon systems. It was based on single-phase systems and does not necessarily provide a means of accurately predicting the behavior of two-phase systems which might be encountered in condensate work. Data published recently by Sage and Oldsl&apos; have enlarged the knowledge of phase behavior in condensate systems by covering extensively the Compressibilities of systems in the two-phase as well as single-phase regions for a wide range of temperatures, pressures and gas-oil ratios. An evaluation of these data indicates that conditions can arise, especially in condensate systems having fairly low gas-oil ratios, in which two phases occur and the usual gaseous compressibility factors do not accurately apply. The measured PVT data of Sage and Olds covering two-phase systems have been correlated with additional such data from the Mid-Continent area and South America; these are presented in Fig. 1. This chart is the result of 1,030 compressibility* determinations from 21 hydrocarbon samples taken from eight oil and condensate fields. Of the data used, 74.75 per cent were from California, 15.45 per cent were from the Mid-Continent area, and 9.80 per cent were from South America. The chart construction is based on the original

Jan 1, 1949

By H. H. Podgurski, F. N. Davis

In silver alloys containing less than 0.2 wt pet Cu. the reaction 9 + 1/2 0, = CuO(s) was found to proceed to equilibrium between 700o and 808oC. From measurements of the equilibrium dissociation pressures of the CuO at several temperatures, the differential heat of solution and the activity coeflicients for copper in silver were calculated. These values were found to be in reasonable agreement with those calculated from data appearing in the literature. A metastable "copper oxide" with an atom ratio of oxygen to copper as high as 1.7 was formed by internal oxidation at 300°C of these same dilute Ag-Cu alloys. No anomalous behavior was noted in the temperature dependence of oxygen solubility in silver. The solubility minima between 300" and 800°C reported several years ago can be accounted for, at least in part, by reactions with trace impurities, such as copper. Cold-worked siloer exhibits an enhanced permeability to oxygen. THIS investigation was undertaken because of our interest in interactions between solute and lattice defects in metals. The reason for the choice of the O-Ag system for study is the anomaly in the solubility of oxygen in silver reported by Steacie and Johnson1 in 1926, specifically that isobars between 100 and 800 Torr show solubility minima at 400°C. It was also claimed that the copper impurity in the silver was not responsible for the minima. Recently, Eichenauer and Müller2 proposed that surface adsorption might have been responsible for this solubility anomaly, but no adsorption isotherms are available to check the pressure and temperature dependence observed at the low temperatures. Surface-tension measurements on silver made by Buttner, Funk, and udin3 suggest that a considerable fraction of a monolayer of oxygen exists on silver at 922°C and at an oxygen pressure of 150 Torr. To account for the pressure sensitivity reported by Steacie and Johnson1 at 300°C, the oxygen bound to the surface at 922°C cannot be considered relevant; the existence of a surface layer at 300°C characterized by a lower energy of binding would be required to explain the effect. All of our attempts to detect an interaction of this type with surface sites have failed. In this investigation we have not been able to associate the dissolution of oxygen in silver with a dislocation interaction. Evidently, severely cold-worked silver does not contain a sufficient number of trapping -sites for oxygen. Indeed, our work shows that oxidation of trace impurities in the silver was probably responsible for Steacie and Johnson&apos;s results. In addition, we have been able to establish the nature of the reaction with copper impurity in silver by thermodynamic considerations. EXPERIMENTAL PROCEDURE Measurements of both the internal oxidation rates and the solubility were made volumetric ally in a system equipped with a gas burette, a mercury manometer, a McLeod Gage, and a mass spectrometer. The silver and the silver-alloy samples were protected from contamination by mercury vapor from our pressure gages by suitably placed refrigerated (-78°C) traps. Silica reaction vessels were employed for experiments performed above 500°C. Spectroscopically pure gases were used in this investigation. The highest-purity silver used in the solubility measurements was 99.999 pet. By chemical analysis 2 ppm of Cu were found in this silver. The two Ag-Cu alloys used in this research were made up to 0.14 and 0.15 wt pet Cu starting with 99.999 pet Ag. An unexpected source of error was discovered in the course of the work. At temperatures near 800°C silver distilled from the hot silica reaction tubes into cooler regions of our system. Although the weight loss of silver from the sample was not important in itself, oxygen was being consumed by the silver distilled into the cooler part of the system to form a stable oxide phase from which the oxygen was not recovered. In a separate experiment conducted to determine the necessary correction, losses between 0.2 and 0.3 cu cm (stp) of 0, in 24 hr were measured at 810°C. On the assumption that the flux of silver from the vessel to form Ag2O de-

Jan 1, 1964

By F. X. Tartaron

The comparison of actual energy of comminution with theoretical surface energy presents a wide gap. On the other hand, Solid State Theory presents a viewpoint that places actual energy of breakage in understandable relationship with the theoretical strength of materials. Therefore the physical background of comminution theory favors Kick&apos;s Law. A11 discussions of comminution theory appear farfetched when comparison is made of the actual energy of crushing with the energy input called for by thermodynamic theory. Quoting Taggart:&apos; "The work required to crush halite was of the general order of 0.10 kg cm per sq cm of new surface produced (based on air percolation) whereas the theoretical energy of this new surface, based on thermodynamic calculation by a number of investigators is 1 x 10&apos;4 kg cm per sq cm or the work expended was 1000 times the new surface energy — an indicated crushing efficiency of 0.10 pct." It is frustrating to probe into comminution theory with the prospect that the energy one is measuring is 99.9 pct waste energy. Do the laws of comminution deal with energy consumed by mineral breakage or with waste heat energy? There seems to be something wrong. Solid State Theory may provide the answer. The following statement is particularly apt to the situation: "The fundamental fact about metals and indeed all materials, is that they are much weaker than they should be. By this, we mean that if the stress required to cause a given plane of atoms to slip past a neighboring plane is calculated from basic atomic and Solid State Theory, the resulting number is some 1000 times too large to be accounted for by errors in approximations or assumptions. We know now that this is due to dislocations, and we even have a fair understanding of the details of the process by which dislocations make a crystal weak. An intriguing possibility which suggests itself, results from the discovery of very minute metallic crystals which are apparently free of dislocations and which have a strength approaching the ideal value for a perfect crystal." Here we have an extraordinary reversal of the situation. The actual energy for breakage is only one-thousandth the theoretical instead of 1000 times as given by Taggart. But the theoretical energy presented by Solid State Theory is the binding energy between atoms (or ions) in a crystal lattice. This is the kind of energy envisaged by Kick&apos;s law. The surface energy or the other hand corresponds to Rittinger&apos;s law. Comparison of theoretical with actual, on the basis of surface or on the basis of Rittinger&apos;s law, seems to present an absurdity; yet Rittinger&apos;s law is confirmed by an impressive array of experimental evidence. Comparison of theoretical and actual on the basis of Solid State Theory or, in other words, on the basis of Kick&apos;s law seems reasonable; yet there is not a scrap of evidence to verify Kick&apos;s law. The situation appears paradoxical, but the developments in the next section of this paper show that the paradox is a hint to a strange turn of events. ENERGY-WEIGHT FRACTION FUNCTION charles3 has presented a final equation showing the relation of energy consumed in comminution to size of product n{-1) a-n + 1 E is energy, C is proportionality constant, a and n are constants and k is the coarsest size in the ground product (size modulus). In the course of study of a large amount of experimental evidence, Charles discovered that in the general case, a-n + 1 approximated zero. For some reason, a and n differed only by unity. But if a—n + 1 = 0 within experimental error, then it is obvious that Eq. 1 states that it requires an infinite amount of energy to grind any sample to any size. This is a case where nature gives an absurd answer to a man-made mathematical equation, hence something is wrong with the equation. An adjustment must be made to fit the equation with experience. TO reveal the nature of the adjustment, the differential equation preceding the second integration in Charles&apos; group of equations is given When Charles integrated the quantity xa-n it was almost impossible for him to realize that he was dealing with a special case of integration and that

Jan 1, 1962

By J. Lohrenz, C. R. Clark, B. G. Bray

Procedures to calculate the viscosities of in situ reservoir gases and liquids from their composition have been developed and evaluated. Given a composition expressed in methane through heptanes-plus, hydrogen sulfide, nitrogen and carbon dioxide together with the molecular weight and specific gravity of the heptanes-plus fraction, the procedures are capable of calculating the viscosity of the gas or liquid at the desired temperature and pressure. The procedure for reservoir liquids was developed using the residual viscosity concept and the theory of corresponding states, and was evaluated by comparing experimental and calculated results for 260 different reservoir oils ranging from black to highly volatile. The average absolute deviation was 16 per cent. This is the first known procedure for calculating the viscosity of reservoir liquids from their compositions as normally available, i.e., including the heptanes-plus fraction. The procedure for reservoir gases uses a sequence of previously published correlations. Evaluation of the procedure was accomplished by comparison of 300 calculated and experimental viscosities for high-pressure gar mixtures in the literature. The average absolute deviation was 4 per cent. The calculations are useful for (I) determining viscosities in compositional material balance computations and (2) predicting the viscosity decrease which occurs when gases, LPG, or carbon dioxide dissolve in reservoir oils. INTRODUCTION Methods to predict viscosities of reservoir fluids from the normally available field-measured variables have been presented. Beal,1 Standing,2 and Chew and Connally3 orrelated oil viscosities with temperature, pressure, oil gravity and gas-oil ratio. Carr, Kobayashi, and Burrows4 and Katz et al.5 have presented correlations for reservoir gas viscosities as a function of temperature, pressure and gas gravity. Lie all intensive physical properties, viscosity is completely described by the following function: Eq. 1 simply states that viscosity is a function of pressure, temperature and composition. These previous correlations1- hay be viewed as modifications of Eq. 1, wherein one assumes more simple functions may be used. The assumptions are practical, because the composition is frequently not known. Further, the assumptions are sufficiently valid so that these correlations are frequently used for reservoir engineering computations. In compositional material balance"&apos; computations, the compositions of the reservoir gases and oils are known. The calculation of the viscosities of these fluids using this composition information is required for a true and complete compositional material balance. For reservoir gases, Carr, Kobayashi and Burrows4 have presented a suitable compositional correlation. For reservoir oils, no correlation is available, and data from reservoir fluid analyses have been used7-9 for compositional material balance calculations.* From a theoretical point of view, this is entirely invalid. The reservoir fluid analysis, whether flash, differential, or other process, does not duplicate the compositions which occur during the actual reservoir depletion process, therefore the viscosities measured during reservoir fluid analysis are not those which occur in the reservoir. From a practical point of view, the "error" of using viscosities from reservoir fluid analysis is of varying and unknown significance. One can say qualitatively that the error is greatest where compositional effects are greatest, i.e., for volatile oil and gas condensate reservoirs and pressure maintenance operations. The first requirement to obtain a quantitative estimate of the significance of the error is to develop a reliable compositional correlation for the viscmities of reservoir oils. No such correlation has been available. Consistent with this requirement, the objective of this study was to develop a procedure to predict the viscosity of reservoir fluids from their compositions. Normally, the compositions of reservoir fluids are available expressed as mole fractions of hydrogen sulfide, nitrogen, carbon dioxide and the hydrocarbons methane through the heptane-plus fraction, with the average molecular weight and specific gravity of the latter. The final correlation was to use the composition in this form. While the more challenging objective of the study was the development of a correlation for the viscosities of reservoir oils, the viscosities of reservoir gases were also studied. The end result of the study was a procedure to calculate the viscosities of reservoir gases and liquids suitable

Jan 1, 1965

By R. L. Perrine

A stability theory is developed for miscible liquid-liquid displacement within a porous medium. In the usual case considered, a high-density high-viscosity "oil" is displaced downdip by a low-density low-viscosity "solvent&apos;! Perturbation methods are used to find the conditions under which the spreading mechanism changes from the stable dispersion process to unstable viscous fingering. We find that instability is conditional, that there is a dependence on the shape of a disturbance leading to a"diameter" effect and that very difficult experimental scaling problems may result. A useful consequence is the definition of a minimum "slug size" for stable miscible displacement. This should make possible optimum use of the solvent process for oil recovery. The results may apply to many situations in which one fluid displaces another of somewhat different fluid properties within a porous medium. INTRODUCTION The process of miscible liquid-liquid displacement looked very favorable when it first received serious consideration as a means to increase petroleum recovery. Some early laboratory results were interpreted to mean that only a small "slug" of solvent was needed, perhaps 2 to 3 per cent of a hydrocarbon pore volume. This "slug" could recover the oil in an entire reservoir provided the fluids remained miscible.1,2 Only a slow growth of the "mixing zone", or gradual transition from oil to solvent, should occur. This would follow naturally from mixing by a dispersion mechanism such as that described by Scheidegger.3 And, indeed, laboratory cores have shown this kind of behavior provided solvent viscosity was at least as great as that of the oi1.2,4-6 The same kind of behavior has also been observed with solvent viscosity lower than that of the oil after the distance traversed became large. Not all laboratory results have been this favorable, however.7-9 Mixing zone growth may become proportional to the distance traveled rather than the square root of this quantity. Such behavior accompanies high rates in short systems of large diameter, provided the viscosity ratio is adverse. If this more rapid spreading were to persist in a reservoir, the solvent requirement for miscible "slug" displacement could exceed 30 per cent of a pore volume. There is considerable economic importance in the difference between 3 and 30 per cent solvent. A 1ogical conclusion is that two different kinds of flow behavior are possible, with each kind leading to a different result. One kind gives efficient displacement, the other does not. In the efficient case, the solvent bank spreads out only by the mechanism we have termed dispersion. Under these conditions, displacement is as near piston-like as possible. In the second kind of flow, found only with adverse viscosity ratios, viscous fingering occurs. That is, permeability variations cause a small finger of low-viscosity solvent-rich fluid to move ahead of its average position within the mixing zone. This creates a path of low resistance to flow, and an even greater amount of solvent-rich fluid follows. Thus, the process is autocatalytic. Once started, the fingering mechanism rapidly becomes dominant. The question to be answered is this. Under what conditions will each of these two different kinds of flow occur? In particular, what conditions denote the transition from dispersion to viscous fingering as the solvent spreading mechanism? The answer may tell us whether or not the desirable, near piston-like miscible displacement process is practical. An answer to this problem can be obtained from theory by the use of perturbation methods. The procedure is as follows. We first formulate a mathematical representation of the system. Then a small disturbance (or perturbation) in solvent concentration profile is introduced and observed to see what happens. Unstable flow is indicated when a small disturbance will grow larger. This will lead to eventual viscous fingering. If on the other hand the disturbance dies out, the displacement is stable. In this latter case, with dispersion as the spreading mechanism, near piston-like displacement is possible. Thus, this paper presents the development of a stability theory for miscible liquid-liquid displace-

By John D. Price

COAL blending, in the preparation of coal before coke making, is so commonly practiced as to be almost universal. But the reasons underlying this practice, the benefits resulting from it, and the materials used in blending vary widely. This paper will outline the various phases of the subject and present information which may be correlated with work that has been done elsewhere. It will deal entirely with work done on the high-volatile coking coals of the western part of the United States, special emphasis being given to the coals of Colorado and Utah. A surveyL of the 86 coke plants in active operation in the United States during 1949 indicates that only 9 plants, or 10.5 pct of the total, charged one single rank of coal into their ovens, while the remaining 89.5 pcl made use of blending in some form. This report indicated that of these total plants 5 used straight high-volatile coal, 4 used straight medium-volatile coal, 47 used blends of high and low-vola-tile coals, 25 used blends of high, medium, and low-volatile coals, 2 used blends of high and medium-volatile coals, 3 used blends of medium and low-volatile coals. The fact that certain plants operated on a single kind of coal should not be interpreted to mean that no blending was practiced there, for invariably such plants secure their coal from more than one source and in the interest of uniformity do blend the coals as received. The general term coal blending covers two fields, the first of which is the mechanical mixing of a number of coals to secure uniformity. Often it is found necessary to secure coal for coke production from a number of different mines; these coals, though of the same general type or rank, may differ in their chemical composition or in the physical qualities they impart to coke made from them. Again, it is not unusual to find that coal from different sections of the same mine may show variations in quality. Under such conditions it may be necessary, in the interest of a uniform final product, to introduce a system of blending bins, a bedding yard, or other mechanical methods of securing a uniform mixture. Unfortunately this form of blending has received very little attention up to the present time; it has not received the consideration its value merits. The second type of blending, while also for the purpose of coke improvement, deals more particularly with the use of a blending agent differing in character from the base coal: it is this form of blending that will be discussed here. To consider only the western coals, for blending may be found necessary for other reasons with other coals, blending has been practiced experimentally or commercially under the following conditions: 1—When a single coal or mixture of coals of the same rank does not produce a satisfactory coke. For example, a high-volatile coal when used alone is likely to contract when coked so that a comparatively weak coke is formed. Or, if of very low rank, the coal may be deficient in the necessary bitumens required for good coke production. 2—When a product of some special quality is required, for example, when a plant ordinarily producing blast furnace coke must operate at slow coking rate to produce a high-grade foundry coke. Under this condition the reduced daily production of all products which accompanies slow coking time may be undesirable, and the use of some blending agent to increase the size of coke made at faster coking rates may be necessary. 3—When greater yield of coke or its coproducts is needed. Depending upon economic values of the products it may be found desirable to increase the yield of one or the other. 4—When supply of a particular coal must be used, either to protect reserves of high quality coking coal or to utilize a surplus or inferior product not otherwise usable. Many materials have been used for blending purposes, the exact agent to be used depending both upon the condition to be corrected and the nature of the base coal. No universal blending agent that can

Jan 1, 1954

By R. E. Ogilvie, E. T. Peters

The X-ray Kossel-line method has been used preaioz~sly for measuring lattice parameters to accuracies of 1 part in 100,000.5 A second application of this method is described for determining the crystallographic orientation of a randomly positioned single-crysta1 spherical volume that can he as small as 50 µ in diameter, within accuracy limits of ±1/2 deg. The theory, experimental procedure, and interpretation of Kossel-line patterns and an experimenta1 verification of the predicted orientation relationship of the cph k and fee a Cu-Si phases are presented. This orientation relationship can be described as (111)a, (00.1)k; [110]a [11.0]k. In addition, the lattice parameter of the a phase was found to he a. = 3.62154 ± 0.00014Å. The Kossel-line method when used in conjunction with electron microanalysis is shown to he capable of providing a complete chemical and structural analysis of a given crystal. THE most easily accomplished methods for determining the orientation of a single crystal are variations of the Laue X-ray diffraction method. Although certain materials can be oriented to within ± 1 deg accuracy by observation of exterior macroscopic features, such as etch pits, cleavage faces, or growth features, the Laue method is generally preferred for routine laboratory application. Both back-reflection and transmission patterns are coordinated by appropriate reference charts and are plotted in terms of reflection plane poles (normals) on a stereographic projection. Orientation is deduced by relating the pole distribution (which is fixed by the crystallographic symmetry of the specimen) to two specified external reference directions. The Laue method has several limitations. 1) Orientation can rarely be determined to better than ±l deg accuracy. Principal errors involve inaccuracy in specimen-to-film distance, measurement confidence of individual Laue spots, and stereographic plotting. 2) Because of the indeterminacy of the X-ray wave length diffracted to a given spot, it is not possible to determine supplementary crystallographic data, such as interplanar spacings and lattice parameters. 3) The method is generally limited to specimens of cubic symmetry or to specimens of high symmetry and known structure. 4) The relatively large cross-sectional area of the incident X-ray beam generally precludes the measurement of relative grain orientation in a polycrystalline material.* _____ Several of these limitations can be overcome by application of the Kossel-line method, which has been previously employed for precision lattice-parameter determinations.&apos;-= For this method, a point source of divergent monochromatic X-rays is generated within the crystal by means of an incident electron beam. The divergent X-rays which fulfill the Bragg law are diffracted by the specimen and are recorded on a film placed either in transmission or back reflection. Analysis of the film yields a direct measurement of orientation and lattice-spacing values to an accuracy of ±1/2 pct. As analyses can be obtained from spherical specimen volumes as small as 50 µ in diameter, the method provides a means for structural analysis of second-phase or impurity precipitates within a given matrix. The primary limitation of the Kossel-line method is the requirement for an electron microanalyzer or similar apparatus capable of producing a finely focused electron beam. This paper is designed to present the theory, experimental procedure, and geometrical interpretation of Kossel patterns. The experimentally determined orientation relationship between the k and a phases occurring in the Cu-Si system and a precision measurement of the a lattice parameter are presented as a practical application of the method. THEORY OF KOSSEL LINES Divergent X-ray beam photography utilizes an effective point source of characteristic X-rays which, when diffracted from a single crystal, form numerous diffraction and absorption cones that are recorded on film.7 The cones generated from a source lying within the crystal are called Kossel lines.&apos; Although the X-ray scattering from a divergent point source contained within a crystal is described in terms of Laue dynamical theory,9 the directions of the diffracted spectra can be ade-

Jan 1, 1965

By E. T. Hayes, A. H. Roberson, M. H. Davies

R. F. Dornagala and D. J. McPherson (Armour Research Foundation, Chicago)—I should like to compliment the authors for a workmanlike job in determining the partial phase diagram of a system comprised of two rnetals which are certainly not easy to work with. We are completing work at Armour Research Foundation on an Atomic Energy Commission-sponsored project for the determination of eight zirconium binary diagrams. Work on the Zr-Cr system has been completed and should be published within the next year. For our work, Westinghouse Grade 3 iodide crystal bar served as the zirconium melting stock. Johnson-Matthey, electrolytic chromium, specially treated for oxygen removal, was employed. The overall constitution of the system determined at Armour Research Foundation is in very good agreement with the present work. We found a eutectic at 18 pct Cr and 1280 °C, somewhat lower than the value reported. This temperature was confirmed by thermal analysis, incipient melting studies, and regular isothermal anneals. The eutectoid was located close to 1 pct Cr and 835°C by metallographic analysis of annealed specimens. Maximum solubility of chromium in /S zirconium was 4.5 pct at the eutectic temperature. Chromium solubility in a zirconium was less than 0.28 pct at all temperatures. We found the compound at 53 pct to melt around 1700°C, with an open maximum, but determined its crystal structure to be hexagonal close-packed (MgZn, type). The lattice parameters were in excellent agreement with those determined by Wallbaum in 1942. The diagrams are in substantial agreement, and .part of the differences are undoubtedly due to the use of different zirconium melting stock. M. K. McQuillan (Birmingham, England)—I read this paper with a great deal of interest, as it covered the same field as some work of my own.&apos; There are a number of points in the present paper on which I would like to comment. First, I should say that I, too, used zirconium prepared by magnesium reduction of the tetrachloride and electrolytically prepared chromium, and melted the alloys in a Kroll-type arc furnace. The purity of my alloys should, therefore, be comparable with the purity of those of the present authors, and any differences in our observations would not be expected to be attributable to this cause. The differences between my observations and theirs concern the presence of the eutectic, the temperature of the eutectoid, and the melting point of the compound. I would be very much interested in any further evidence the authors may have for the occurrence of the eutectic at 1380°C. During the course of my work I noted that a number of my alloys containing 60 to 90 pct Zr melted at about 1400 °C, and for a time assumed that a eutectic occurred at this temperature as described in the paper. On further investigation, however, I found that the structures of the as-melted alloys could not be made to fit in with this interpretation of the system. If a eutectic exists in this region of the system it would be expected that the as-melted alloys would show the usual type of cast structure, i.e., dendrites of the compound plus eutectic. This, however, does not occur, as may be seen from Fig. 9. The compound seen there is not dendritic in form, and the remaining material is by no means certainly eutectic. It may be argued that a compound such as ZrCrl would not form dendrites but would tend to crystallize in geometric shapes. In this case, however, I have evidence to the contrary, as on the chromium side of the compound, where a eutectic occurs at about 1545"C, the compound formed from the liquid takes on a conventional dendritic form, and the eutectic is observed in the interdendritic spaces in the usual way. There is no reason to suppose that the compound would behave differently in an alloy lying on the zirconium side of the compound composition if a eutectic existed there too.

Jan 1, 1953

By J. C. Wilhoit, J. B. Cheatham

Although an oil well is a long cylindrical hole with an irregular bottom, it appears likely that the nature of the stress concentration at the bottom of the hole can be ascertained from an analysis of the stresses around a short cylindrical cavity with rounded comers and smooth bottom. Such a cavity is studied primarily because it leads more readily to a solution to the problem by the use of stress functions. In this paper the stress distribution around a short cylindrical cavity subjected to bit loading, overburden and drilling-fluid pressures is determined by means of an analytical solution which approximately satisfies the boundary conditions of the problem. From this solution the stresses at the corner of the hole are calculated to be about 35 per cent lower than comparable results obtained by photoelastic and relaxation analyses. This difference is apparently due to the large radius of curvature at the corner of the cavity in the present analysis. Since good agreement is obtained between the results of this analysis and the stresses calculated for a similar loading on a semi-infinite elastic solid, it is concluded that the bit action in the region near the center of the hole is not appreciably affected by the presence of the sides of the hole. INTRODUCTION Much has been written concerning drilling "under down-hole conditions" and pertaining to the stress distribution in the rock at the bottom of an oil well.1-5 For example, it is known that identical rocks can be drilled more rapidly at the surface than under subsurface conditions of pressure and stress.6 Information on the behavior of rocks under loading can be obtained from triaxial test data.7-9 From such tests it is found that rocks exhibit brittle failure when the confining pressure and pore pressure are equal, but the mode of failure may change to ductile as the difference between the confining pressure and the pore pressure is increased. Brittle failure implies that there is very little permanent deformation before fracture, whereas ductile failure indicates that permanent deformation takes place before fracture. Some rocks are ductile at differential pressures of 5,000 psi , but other rocks are brittle even at differential pressures of more than 50,000 psi. Cuttings embedded in mud at the bottom of the hole may act as a plastic mass which the bit teeth must penetrate in order to attack the virgin rock below. 10 During drilling, the mud pressure acts as the confing pressure, and in many cases the difference between the mud pressure and the formation pore pressure IS sufficiently low so that the rock most likely falls in a brittle manner. Very little penetratiori by the bit teeth may be required for brittle failure of the rock. For these brittle materials, the elastic stress distribution is of practical interest in determicing the possible effects of a given loading on failure. The penetration of bit teeth into ductile or plastic rock has been analyzed previously ll,12 and will not be considered further in the present work. During actual drilling, many teeth act at various points ovel an irregular hole bottom. In the present analysis a solution is obtained for an idealized problem of determining the elastic stress distribution caused by only one tooth acting alone at the center of a smooth hole bottom. It obviously is not proposed that any driller should put a special one-tooth bit on bottom, but it is hoped that extensions of the simpler problems can eventually lead to a better understanding of more complex actual drilling phenomena. To obtain a solution by the use of stress functions, a short cylindrical cavity with rounded corners and smooth bottom is studied. The solution to this problem should give an insight into the nature of the stress concentration at the bottom of an oil well which is in reality a long cylindrical hole with an irregularr bottom. After selection of the proper curvilinear co-ordinate system, the stress functions are expressed as series of Legendre polynomials. The coefficient of each term in the series is then determined to satisfy the boundary conditions at a finite number of discrete points on the boundary of the cavity with least square error. Results are compared with bottom-hole stresses obtained by photoelastic and relaxation analyses.