Search Documents

By R. G. Hawthorne

This report is concerned with fluid displacement in porous media, in those cases where viscous and gravitational forces control the displacement. Such a system would usually be found in a sand body of large physical dimensions such as an oil reservoir, although it is possible to create such a system in the laboratory. It is shown that the position of the fluid interface can be predicted by numerical calculations using a basic idea presented by Dietz. Fluid flow is considered in a vertical plane in a homogeneous, porous medium of sufficient thickness that the capillary transition zone is small in comparison with the total reservoir. A theory developed by Dietz' is used to make numerical calculations of the position of the fluid interface. The results for several conditions are compared with scaled model experiments. The results show that, for gas drive in a reservoir of steep dip, a relatively low flow rate can displace large volumes of oil before gas breakthrough. On the other hand, water injection at favorable mobility ratio and low dip may show best performance at high rates. Water tends to underride the oil and, given sufficient rime, will break through without much oil displacement. For certain conditions, which include relatively low flow rate, the interface is a straight line and its behavior is simple to calculate. At higher flow rates, the interface is unstable, and a numerical solution was programed for an automatic computer. In general, good agreement is shown between the fluid model and the computed results so long as gravitational forces have control. For a water drive at very unfavorable mobility ratio, many small water fingers appear. These viscous fingers are not controlled by the relatively small gravitational forces. When viscous fingering becomes the controlling factor, the mathematical model is oversimplified, and results do not check the fluid flow model. INTRODUCTION Present methods of reservoir analysis depend upon certain simplifying assumptions to obtain mathematical descriptions of practical use. Material-balance methods (Muskat2 or Tarner2) assume uniform fluid saturations in the entire reservoir, or in a few subdivisions of the reservoir. An unsteady-state flow calculation by West, er at considered pressure and saturation changes in flow to a well during solution gas drive, and neglected gravity effects. Results showed only a 4 per cent difference for ultimate oil recovery by the Muskat method, even though the case chosen for study was one in which unsteady-state effects should be high. The Buckley-Leverett5 method commonly assumes a one-dimensional flow system. It is applicable at high flow rates where viscous forces predominate over gravity forces. Simultaneous, parallel flow of the two fluids is assumed, and the concept of a fluid interface is not introduced. Permeabilities to each fluid for a given saturation must be known. The method is not applicable for a two-dimensional system where cross flow becomes possible. Less well known is the displacement equation derived by Dietz. This method is designed for two-dimensional flow systems and assumes a definable fluid interface within the porous medium. Dietz showed that, for a range of low flow rates, the interface would be stable. straight and at an angle of inclination which could be simply calculated. At a certain critical flow rate, the calculated interface tilt would equal the formation dip. For higher flow rates, a finger of displacing fluid would invade the displaced fluid. Dietz indicated that his method applied only to macroscopic reservoir behavior, while the Buckley-Leverett method applied to the small transition zone at the fluid interface. The examples worked out in this report are based on the fluid-displacement theory of Dietz. It is shown that the Dietz theory may be used to derive equations analogous to the Buckley-Leverett equations. In contrast to the Buckley-Leverett method, flow is considered in a plane rather than being limited to a line. Rather than a frontal advance, the movement of a fluid interface is followed. For flow rates substantially exceeding the critical rate and for high viscosity ratio, many fingers of invading fluid occur-—rather than the single finger assumed by Dietz. On the other hand, so long as some gravitational influence remains, the flow is not entirely parallel to the bedding planes as assumed by Buckley and Leverett; therefore, both methods fail to give an adequate descrip-

By Shiro Ban-ya, John Chipman

The effects of many alloying eletnents on the acticity coefficient of sulfur in liquid iron have-been studied by the equilibriutn in the reaction Sfin Fe) + Hz = HzS at 1550°C'. Results are expressed in terms of a concentration variable for a nonmetallic or for a substitutional metallic solute. Activity coefficient of sulfur, defined as increased by B, Al, C, Si, Ge, Sn, P, As, Sb, Mo, W, Co, and PI. It is decreased by Cu,Au, Ti. Zr, V, h'b, Ta, Cr, ,23n, and Ni. Irulues oj. the interaction coejficient 0i = 6 In are labulated. The same interactions are expressed also in terms of atom fraction and of weight percent. The thermodynamic properties of sulfur in liquid iron have been fairly well established. Our recent paper1 reported new determinations of activity coefficient up an atom fraction of 0.12 and equations based on a careful review of all published data. The effects of alloying elements on the activity of sulfur have been studied by several investigators,2"11 especially by Morris and Williams for silicon,' by Morris and Buehl for carbon,~ and by Sherman and Chipman for manganese. phosphorus, and al~minum.~ However, the effects of some important alloying elements remain to be determined. The purpose of this investigation was to determine the effects of many alloying elements on the activity of sulfur in Fe-S-j ternary systems. EXPERIMENTAL METHOD This study was based on experimental determination of equilibrium in the reaction at 1550°C: The same apparatus and procedure as described in our previous paper have been retained, and the same corrections were applied for dissociation of H2S. The alumina crucibles used in the experiments consisted of four separate compartments. In a given experimental run. kach of three alumina crucibles held four different samples of about 3 g each which reached equilibrium in the same gas atmosphere and temperature. The charges were made up of electrolytic iron, pure iron sulfide, and desired alloying elements, which were pure metal or master alloys made in the laboratory. The weighed samples were held for 4 to 12 hr in the prepared atmosphere at 1550°C. They were then low- ered to be cooled as quickly as possible. The quenched metal beads were crushed to avoid the errors of segregation in the ingot. The sulfur content was determined gravimetrically and alloy content by appropriate chemical analysis. CALCULATIONS The apparent equilibrium constants of Eq. [I] are expressed as follows from the corrected gas ratios and sulfur concentrations: The term K" is the observed equilibrium constant in any given ternary solution, K' is the value for the Fe-S binary solution. and the limiting value of K' in the infinitely dilute solution is designated by K, which is the true equilibrium constant in Eq. [I]. In the thermodynamic treatment of nonmetallic elements in a metallic solution, it has been suggested1' that the lattice ratio has certain advantages over other variables to express the concentration of solute. In interstitial solid solutions the lattice ratio zj is proportional to the ratio of filled interstitial sites to those which remain unfilled. The equations derived for the activity of the solute using zj as the concentration variable are found also to be applicable to liquid solutions containing nonmetallic solutes when the nonmetal is treated as if it were interstitial. For this purpose we adopt the following definitions: The quantity vj which is negative for interstitial solutes is taken as -1 for nonmetallic and +1 for metallic solutes. For purposes of calculation however ; j may be assigned a value which results in a linear relation such as shown in Figs. 1 to 15. The activity coefficient of sulfur, Qs, and equilibrium constant. K(z). are defined as follows: According to a Taylor series expansion. the logarithm of the activity coefficient of sulfur in Fe-S-j ternary system is: However, the value of 6 In */6zi remains constant through a broad range of dilute solutions and the terms of higher order are negligibly small. As a consequence, Eq. [6] is simplified as follows:

Jan 1, 1970

By A. K. Mukherjee, J. E. Dorn, J. D. Mole

Investigations were carried out on the effect of polyslip on the strain hardening of aluminum single crystals. The orientations investigated were those lor which the tensile axis was in the [001], [111], [112], and [012] directions plus another for which the Schmid angles for {111}(110) slip were 1 deg. The experimental data were analyzed on a model based on the intersection of dislocations with particular emphasis on the effect of polyslip on the activation volume for inter-section. It is shown that the rate of strain hardening inreases for those orientations wherein attractive dislocation intersections occur and that those orientations which produce the greater number of such intersections exhibit the greater strain hardening. Good correlation of the data is obtained with the concept that attractive junctions, as proposed by Saada, Play an important role in accounting for the rate of strain hardening. EXISTING concepts on the nature and cause of strain hardening in fcc metals have been deduced principally from experiments on the deformation of single crystals under single slip. The effect of crystal orientation on the shapes of the stress-strain curves for single slip have been summarized by seegerl and more recently by Clarebrough and Hargreaves.2 Tensile specimens whose axes fall near the center of the [001]-[011] line of the standard triangle of the stereographic projection exhibit the longest range of easy glide (Stage I) and the lowest rates of linear hardening (Stage 11) whereas specimens whose axes lie near fie [001]-[111] line of the standard triangle have limited or no easy-glide range and exhibit somewhat higher linear strain-hardening rates. Specimens whose axes lie near the [001] or the [ill] poles do not exhibit easy glide and have the highest rates of linear hardening. Kocks3 has shown that the highest rates of linear hardening Occur under Polyslip when the tensile axis coincides with the [111] or the [ 001] pole. Since the rate of strain hardening is sensitive to specimen orientation and the incidence of polyslip, these relationships might help to discriminate between various dislocation models for strain hardening in fcc metals. Previous attempts to analyze the effect of orientation on strain hardening,1"3 however, did not provide a unique answer to this problem. Consequently the present investigation was undertaken wherein additional data, particularly that for the effect of polyslip on the activation volume for intersection, was also determined in order to provide more complete information on the details of strain hardening. Whereas analyses of these data reveal that several recommended models for strain hardening are at variance with the facts, good correlation of the data is obtained with the concept that attractive junctions4 play an important role in accounting for the rate of strain hardening. I) EXPERIMENTAL APPROACH seegerl demonstrated that slip in fcc crystals at low temperatures is dependent on thermally activated intersection of glide dislocations with forest dislocations. This has been confirmed by tests on single crystals of aluminum by Mitra, Osborne and Dorn5 and on polycrystalline aluminum by Mitra and Darn.' Thus in accord with Seeger's theory, the shear strain rate, ?, below a critical temperature, T, is where .V is the number of points of contact per unit volume between forest dislocations and glide dislocations, A is the area swept out per successful intersection, b is the Burgers vector, v is the frequency of vibration of the segment of the glide dislocation undertaking intersection, U is the activation energy for intersection, k is Boltzmann's constant, and T is the absolute temperature. For aluminum, which has an extremely high staeking-fault energy, the constriction energy is negligibly small and therefore the activation energy decreases practically linearly with the stress according to as will be recomfirmed later, where Uo is jog energy at the absolute zero, G and Go are the shear moduli at the test temperature T and O°K, respectively, L is the spacing of the forest dislocations, t is the applied shear stress for slip, and tGo is the stress field that must be surmounted athermally.

Jan 1, 1965

By N. S. Stoloff, R. C. Ku, R. G. Davies

The stress-strain behavior of Fe-Co and Fe- V alloys containing up to 25 pct solute have been studied in the temperature range 25° to - 196°C. The microyield stress is independent of temperature for all alloys studied, while the temperature dependence of the ordinary yield stress is large and comparable to that for unalloyed iron. An alloy-softening phenomenon is observed for many alloys relative to unalloyed iron, and appears to be a consequence of removal of inter-stitials from solution and their subsequent agglomeration. The ductile to brittle transition temperature is increased markedly by both cobalt and vanadium. Transmission electron microscopy observations reveal that cross slip and cell formation are restricted with increasing solute, which can account for many aspects of yielding and fracture behavior. IRON and other metals of bcc structure slip on many planes of the (111) zone at moderate and elevated temperatures, giving rise to wavy glide. At low temperatures, however, there is a tendency for slip to be restricted to planes of only one or two crystallographic types; in the case of Fe-Si alloys for example, (110) glide predominates, giving rise to planar slip bands."' Restriction of cross slip at low temperatures in iron has been linked by Brown and Ekvall3 to increased strain hardening in the preyield microstrain region, leading to an apparently large temperature dependence of the yield stress. Support for this view has been provided by Lawley and Gaigher4 from a transmission electron microscope study of molybdenum. They concluded that dislocations at lower temperatures, by virtue of their inability to change planes readily, are unable to interact and annihilate, leading to an increase in internal stress and work-hardening capacity. Restriction of cross slip at low temperatures has also been suggested to be an important factor in determining the ductile to brittle transition of single-phase solids, including metals of bcc structure.5 It is pertinent to note that most solute elements tend to raise the ductile to brittle transition temperature of poly crystalline ferrite,6,7 although in some cases, e.g., Fe-Cr7 and Fe-Al,8 the dilute alloys may actually yield at lower stresses than unalloyed iron. It appears, therefore, that solute elements may have an effect analogous to lowering the test temperature with regard to restricting cross glide. Such an effect has indeed been observed in the ionic alloy system KCl-KBr,9 in which bromine atoms, although causing little misfit and consequently little effect on the yield stress of KC1, have a pronounced effect in restricting wavy glide and inducing brittleness. The purpose of the present investigation was to determine the effects of selected solute atoms on the yield stress and fracture behavior of iron, and to elucidate in turn how the slip character was related to these properties. Since it is known that interstitial atoms exert a strong influence on the mechanical behavior of iron, the solute elements chosen for study were vanadium, which reacts strong with interstitials,10 and cobalt, which does not form highly stable interstitial compounds." EXPERIMENTAL PROCEDURE Four-pound ingots of each of the following alloys were vacuum-melted in zirconia crucibles: Fe-1, 5, 10, 25 at. pct Co, and Fe-1, 4, 10, 20 at. pct V. An ingot of unalloyed iron from the same original stock as that used for the alloy was remelted under similar conditions to eliminate variations in interstitial content arising from differing thermal histories. Chemical analyses of the iron, cobalt, and vanadium melt stock are listed in Table I. All ingots were rolled and swaged, commencing at 850°C, to 1/4-in.-diameter bar. The bars were then annealed at 850°C for 1 hr (except that the 10 and 20 pct V alloys were heated for 2 hr) to give a uniform equiaxed grain size of about 0.05 mm. Tensile samples of 0.125 in. diameter by 0.8 in. gage length were utilized for studies of the ductile to brittle transition. A second set of tensile samples were utilized for microstrain experiments in which strain sensitivity of 1 x 10"6 was achieved by a capacitance method.3 Cylindrical compression samples 0.250 in. diameter by 0.4 in. long also were prepared to study the influence of composition on the macroscopic yield stress, and for metallographic analyses of defor-

Jan 1, 1965

The report of the Nominating Committee was presented and is given in this same Bulletin. A petition was received for new requirements for membership in the Institute and was ordered to take the statutory course. A petition was received for the establishment of a Tulsa Section of the Institute and was referred to a committee to report on the territory. Thirty-six Members, two Associates, and two Junior Members were elected. The status of two Junior Members was changed to Member.

Jan 12, 1917

By Porter W. Shimer

The solution used in the following method for phosphorns-determination is the filtrate obtained in the nitric and sulphuric acid method for the determination of silicon (Transactions, vii., 346). This filtrate, being easily and quickly obtained, and always quite free from silica, is a desirable one for phosphorus-precipitation. IN order, however, to get all the phosphorus in a precipitable form the solution must be made under conditions more strongly oxidizing than simple solution in boiling nitric acid. In a solution thus

Jan 1, 1889

By J. M. Mutmansky, T. K. Greer

This paper takes a look at the optimal velocity to be used in coal mine entries and determines a statistical relationship for optimal velocity that can be applied under current mining conditions. This relationship provides the optimal velocity for coal mine airways under variations in power cost, mining cost differential, airway friction factor, time value of money, airway lifespan, and quantity of air carried. This enables the ventilation engineer to make an easy, rapid determination of the optimal velocity for planning purposes.

Jan 1, 1985

By G. H. K. Schenck, H. W. Robinson

Private weather forecasting services available on a contract basis can significantly reduce the uncertainty associated with weather-induced changes in mineral production and demand. Such services offer the mineral producer the opportunity for tactical and strategic response to weather conditions. There are at least 50 private meteorological consulting services as well as numerous individual meteorology consultants. Fees for services range from $150-$600 per day for one man to over $20,000 per month where a complex forecasting or modification program is undertaken.

Jan 1, 1971

By Sir Lowthian Bell

WITH the exception of air and water, it is open to question whether there is any form of matter which the human race could less easily spare than iron. Short of going the length of asserting that, without this metal for an anchor or steel for the compass, the adventurous navigator could never have crossed the wide Atlantic, we may credit the locomotive' and the steamer, and hence iron, for that sequence of events which has peopled North America with the

Jan 1, 1891

By J. Weertman

The paper of Gifkins and Snowden considers the interesting but difficult problem of determining the stress dependence of secondary (steady-state) creep at low stresses. These authors have concluded that at stresses below 250 psi (2 x 107 dynes per sq cm) the secondary creep rate of lead is proportional to the stress (viscous creep) and is not proportional to the stress raised to about a fifth power. The experimental data considered by them were obtained on tests conducted at room temperature and at 50°C. The lowest stress employed was 50 psi (3.5 X 106 dynes per sq cm). The authors pointed out the main difficulty in determining the stress dependence of creep at low stresses. The creep tests must be run for very long lengths of time. However they made no estimates of when an experimental creep rate determination must be rejected because it does not represent a true steady-state or minimum creep rate. In order to be certain that a creep curve is within the steady-state region, the total creep strain should be of the order of 0.1 to 0.2. For a creep test of a year's duration this requirement implies that a secondary creep rate smaller than about 10-3 per hr cannot be measured reliably. The corresponding creep rate for a 10-year test is 10-6 per hr. The creep rates of the tests that were considered by the authors to prove the existence of viscous creep were of the order of or less than 10-6 per hr. One can conclude reasonably that this data does not prove unambiguously that large strain steady-state creep rate of lead is proportional to the stress in the stress range of 50 to 250 psi (3.5 x 106 to 2 X 107 dynes per sq cm). Another technique can be used to obtain the stress dependence at low stress levels. The creep rate is a very sensitive function of temperature. The creep rate can be increased by very large amounts merely by increasing the temperature. We carried out steady-state creep tests on lead single crystals25 at temperatures up to 320°C. We were able to obtain creep rate data down to stresses as low as 35 psi (2.5 x 10' dynes per sq cm). Our smallest creep rate was 8 x 10-5 per hr. Thus we obtained large strain, steady-state creep rates to even lower stresses than were considered by Gifkins and Snowden. No evidence was seen for viscous creep. The creep rate was proportional to the stress raised to about a 4.5 power down to the lowest stresses. Since there is no reason to believe that changing the temperature should change the stress dependence of steady-state creep, we feel that large strain viscous creep does not occur in the stress range quoted by the authors for lead single crystals or large-grain polycrystalline samples of lead. This conclusion does not imply that viscous creep may nat occur in a lower stress range or in the same stress range for fine grain material or at creep strains very much smaller than 0.1. Support by the U.S. Office of Naval Research is acknowledged. Authors' Reply R. C. Gifkins and K. U. Snowden We thank Dr. Weertman for his discussion and although, as we hope to show, we do not agree with his reservations, we do concur in stressing the importance of ensuring that creep rates are reliably obtained. Dr. Weertman appears to be content to accept n = 1 for low stresses with fine-grained material but not for single crystals. We believe our results show that the former result cannot be accepted without also accepting the latter. We will also show that the probable errors in our minimum creep rates are insufficient to alter our conclusions, that the criterion proposed by Dr. Weertman is arbitrarily restrictive and his alternative experimental approach possibly invalid. 1) A principal result of our Fig. 1(a) is that n = 1 for polycrystalline specimens at room temperature and 50°C for stresses below -250 psi. There was evidence that crystal slip and grain boundary sliding contributed approximately equaily to the overall strain in this low-stress regime. This implies that either a) grain boundary sliding controls slip within the grains or b) both grain boundary sliding and crystal slip independently occur according to mechanisms which give n = 1. Alternative a does not seem acceptable, so we were forced to consider b. This led us to reexamine work on bicrystals by Strutt and Gifkins and plot curves Fig. l(b). Previously Strutt et al. (loc. cit.), had merged these points with others using the Zener-Holloman parameter and thus, we now believe, had been led to overlook the behavior where n = 1. Curves C and D in Fig. l(b) did appear to confirm the hypotheses that n = 1 for crystal slip at these low stresses and the sliding curve F was similarly of the expected form. It was comparatively easy to find a quantitative theory to account for n = 1 for sliding and the similarity of curves C and D to curve F (all obtained from the same set of specimens) led us to feel that the single-crystal curves were valid. 2) We believe the secondary creep rates for both the polycrystalline and single-crystal specimens to be in error by factors c2. In Fig. 6 creep curves for polycrystalline specimens of lead(1) and lead(II) are reproduced as curves a and b, respectively, and curve c is for a single crystal at 100 psi. It is clear that, although the attainment of secondary creep rate takes 2 years for a and 150 days for b, thereafter the curve is linear for periods of 7 and -1 year, respectively. The single crystal has a linear portion commencing after 20 and extending to 90 days. Creep extension was measured directly using a traveling microscope reading to 0.01 mm on gage lengths marked on the specimens; the gage lengths

Jan 1, 1968

By H. R. Joesting

Geophysical search for minerals is stepped up to meet requirements of expanded defense program - Seismic and electrical surveys become firmly established for ground-water and engineering investigations.

Jan 2, 1953

By ALEC CROWELL

Mississippi has had oil production for only the past six years and natural gas production of minor magnitude since' 1926. The search for oil and gas commenced in 1903 and 1496 wells had been drilled in search of oil and gas up to Jan. 1, 1946. As the year 1945 ended, 659 wells were producing oil and gas from 23 separate reservoirs. The production for the year was 8,353,719,000 standard cubic feet of natural gas measured on a pressure base of 15.025 lb. per sq. in. abs., 15.316,009 bbl. of oil¬field water and 19,093,766 bbl. of crude oil and condensate. Production for the year 1944 was 2,636,067,000 standard cubic feet of gas, 14,979,793 bbl. of salt water and 16,420,346 bbl. of crude oil and condensate.

Jan 1, 1946

By Allen H. Rogers

THE need for coordinated effort on public problem by engineers has long been felt. Early in June there will assemble in Washington a conference composed of delegates from all the engineering organizations in the country. The call for this conference is issued by the Joint Conference , Committee of the four national societies as the completion of its work on the organization of a congress composed of delegates representing all the engineering organizations of the United States. The plan for this confederation was presented at a joint meeting of the administrative boards of the national founder societies on Jan. 23, 1920; and was unanimously approved in principle by that meeting. The work of the committee since that date has consisted in drafting a constitution of the organization to be presented at the conference for discussion and, after such amendments as may be adopted, acceptance.

Jan 1, 1920

By Sherwin F. Kelly

THE principal interest in the meeting sings arranged by the Geophysical Committee this year centered around the many papers dealing with the practical applications of geophysical methods to solving geological problems in many other parts of the world. The program Wednesday morning dealt with work in South Africa, Poland, and South America. Rudolph Krahmann's first paper gave a perspective of the geophysical investigations upon mineral deposits in South Africa. A striking point in his story was the tremendous increase in magnetic work during the last two or three years. This is largely accounted for by the intensive exploration for the extension of the Witwatersrand gold mining area. The author mentioned many surveys in which other methods were used for the search for ores of copper, lead, and zinc, and in the study of water-bearing and oil-bearing structures.

Jan 1, 1935

By E. A. Breitenbach

A computer program has been developed to afford rapid and complete quantitative log analysis for exploration and production decisions. The computation consists of automatic selection of tops and bottoms of porous intervals from the digitized data, and then point-by-point calculations within each selected interval. Nearly all log types can be analyzed. This paper presents the calculation techniques found to be appticable to machine evaluation and gives examples of their use. INTRODUCTION As the quantitative interpretation of well logs entails repetitive use of charts and equations, it is natural that digital computer programs would be written. A Universal Log Interpretation Computer Program (ULICP) has been developed to afiord rapid and complete analysis for exploration and production decisions. Application of this program enables the log analyst to: (I) apply rapid quantitative log analysis for exploration and production decisions on either a field, a well or on a single zone basis; (2) apply concepts of interpretation requiring detailed numerical analysis; (3) analyze all porous intervals on each log, rather than a few selected zones, and complete the analysis in less time than previously required. Computations are reported for every digitized point in each zone; (4) use empirical techniques applicable to a given area as an integral part of the computation; (5) experiment with the empirical coefficients and exponents in the interpretation equations to find the best possible solution; and (6) automatically plot both the original and computed data to scales analogous to the field prints. The primary utility of a log interpretation program stems from its ability to do an overwhelming amount of work with very little man-power. If a calculation procedure or thought process can be formalized to the extent that step-by-step logic can be written, a computer program can be developed that follows this logic. The major question is one of economics. A feasibility study of the costs for such a program indicates that digital processing is economical primarily as a means of increasing the productivity of the log analyst. In effect, the probability of missing productive intervals in any well, because of a lack of time to do detailed calculations, is greatly reduced. The nucleus of ULICP is programmed to compute large sections of a suite of logs by selecting zones automatically and then performing all pertinent computations on a data point by data point basis within each zone.'.' An entire suite of logs can be processed in this manner with very little manual intervention. Sufficient programming logic is available so that each log analyst can request computations pertinent to his area. These requests are made by simple additions or deletions of information on the input header cards. Hence, the log analyst is always in complete control of the computation process. The evaluation of a suite of logs requires pre-editing, digitization, computation, presentation of results and interpretation. Work by the log analyst is reqilired only in pre-editing and interpretation. Thus, he is allowed more time for comprehensive interpretation, rather than calculation. For continuity, the discussion of ULICP is organized sequentially: pre-editing, digitization, computation and presentation of results. DISCUSSION PRE-EDITING The extent of pre-editing prior to computation is dependent on the format of the original data. For analog prints, it requires inspection, correlation and editing of the logs, plus the entering of required data an special forms. For digitized data such as magnetic tape field recordings, only the special forms are necessary. The process for analog prints will be given here. The inspection and correlation process involves the selection of sections for digitization and the correlation of the traces to a common depth. To decrease the cost of digitization, traces that exhibit large variations in shale formations can be redrawn to a non-zero baseline. The next step is to enter all pertinent data on the input header cards. The header cards presently used for ULICP are presented as Figs. 1 through 6. Cards 1, 2, 3 and 4 (Fig. 1) are defined as the Main Header Cards. They describe the particular well for output identification and give basic information. Cards 5 through 15 (Figs. 2 through 6) are defined as the Block Header Cards. As such, they define the log types and the interpretation parameters for the block of data immediately following. Cards 1 through 10 are required for every computer run. Cards 11 through 14 pertain only to nuclear logs. Card 11 is required if any nuclear log is supplied, and cards 12 through 14 are required only when gamma-ray, density or neutron log data, respectively, are supplied. Card 15 (Fig. 6) must be supplied when cross-plot calculations are required. Either two- or three-component

Jan 1, 1967

By E. Kendall Cork

The traditional financial objective for a single mine company has been to operate as frugally as possible and to pay out most of the earnings as dividends. If the business is cyclical (as it is for most metals) the dividends might fluctuate quite widely. When the mine is exhausted the company disappears. This is still quite a viable strategy for a single mine company. It is not however a viable strategy for the world as a whole. The mining industry is built by mine development companies who can mobilize the people and capital to bring new mines into production. Their skills must include marketing, engineering, finance and other politics. It is very rare for a property to be brought in without the support of a major company that can provide all these services. The exceptions will usually have some other form of big brother support, for example the U.S. government uranium contracts at guaranteed generous prices. The mine development company will seek as a minimum to perpetuate itself by developing new mines in order to replace those which are running out. The more common and more ambitious objective is to grow -- that is to add to its ore reserves and current production by developing more new mines. The financial objectives for that company are very different. Obviously if all the earnings were paid out in dividends there would be nothing left to work with. The first financial policy then is to spend an appropriate amount on exploration for new properties. The next is to retain enough of the earnings to provide the capital for new projects at least sufficient for the equity. There is no magic formula as to what proportion of earnings should properly be distributed as dividends by a growth-oriented mine development company. As a rough rule of thumb distributing half or more will probably leave too little to work on and 30% or so is probably a good balance. However the circumstances differ widely from company to company. It may be useful to set an objective for the rate of growth of a company's earnings. Some have picked rates such as 15% per annum compounded. Others have set a target in real terms which might appear as 10 or 11% plus inflation. Obviously the arithmetic of compound interest is very attractive; however in practice there is much variation. Indeed current returns from existing operations swing widely with the business cycle and there is no assurance that economic new properties will be found according to someone's arbitrary time schedule. For example, Western Mining Corporation Limited in Australia explored for 30 years with little to show for it, but then found the great Australian nickel deposits and more recently the huge Roxby Downs copper. That long dry spell could not have fitted anyone's arbitrary calendar of growth and yet they would not have found such orebodies without that long period of effort. Should they have abandoned the search? Once a new property has been found or acquired there has to be a threshold rate of return on the new capital to be invested against which to evaluate the property's economics. Conventionally this seems to be 15% after tax, a number common in other heavy industries as well. In some cases it is expressed as a lower number plus allowance for inflation. Discounted cash flow analysis is a very useful tool but it does not make the decision. In the end a "go" decision depends on judgment of many factors some of which are numbers used in the DCF calculation whose credibility must be examined. It is curious how frequently investment proposals come in with the rates of return very close to 15%. The project advocates know that a number much less than 15% will not fly and that a number much more is not necessary. With much higher nominal and real interest rates of recent years, even though before tax, logic suggests that the hurdle rate should also rise. The power of compound interest is so great that 20% is very hard to achieve in any cash flow projection but 18% may be a sensible yard - stick. Once again it is remarkable how many project proposals come in with an 18% return. On the record the mining industry as a whole has not been overly restrictive in choosing its hurdle rates of return. This is shown by the abundance of metals in recent years and the failure of metal prices to keep up with inflation. All of the foregoing is standard text book stuff.

Jan 1, 1985

By D. M. Kehn

A method has been devloped for chromatographic analysis Of the vapor and liquid phases Of a a system containing methane to components having 20 or more carbon atoms. The method uses a windowed equilibrium cell in which volumetric phase behavior of the system can be observed accurately and from which small samples of gas or liquid can be withdrawn for analysis. Analyses are made using two chromatographs, one for the lighter and one for the heavier components in a sample. Combination of the two analyses yields a detailed analysis of the gas or liquid sample. The complexity of the condensate heavy ends is evident from the chromatograms of these fractions, and the predominance of the paraffin hydrocarborn serves as a useful marker in interpreting the chromatograms. The K-values obtained in this analytical method are presented for a high-pressure condensate system and predict closely the observed volumetric behavior of the system. INTRODUCTION Quantitative analysis of hydrocarbons from natural gas reservoirs is necessary for several reasons—to calculate the amount of sales gas produced, to calculate the amount of natural gasoline produced, to plan a liquid recovery system, or to calculate the potential economic value of a reservoir produced under one or more of several different conditions. Analysis of natural gas fluids produced to the surface consists of identifying and computing the mol fraction of each component of the mixture. Although methane is the predominant component, varying amounts of ethane, propane, butane, pentanes and heavier components are also present. Materials containing up to 30 carbon atoms occur in amounts which decrease with increasing molecular weight. However, the quantities of components in the 20 to 30 carbon atom range are usually so small that their importance is negligible, and they are undetect-able in natural gas by ordinary analytical methods. All the components up to those having 20 carbon atoms may sigsficantly affect phase behavior, however. Commonly, only the methane-through-pentane fraction is analyzed quantitatively for each component, while components heavier than Pentane are lumped and repored as "hexane- plus". Expensive, tedious techniques are required for analysis of this fraction. Consequently the detailed analyses needed for prediction of reservoir behavior are usually undertaken only when major gas fields are being developed. The need for complete analyses of condensate systems is apparent when it is recalled that most gas fields are produced by pressure depletion. As the pressure declines, some of the heavier hydrocarbons are lost as liquids which are in the reservoir. In many instances the amount of liquid in equilibrium with the gas phase at high pressure constitutes only 1 or 2 mol per cent of the total system. Flash calculations generally must predict the actual amount of liquid with an accuracy of a few per cent in order to be useful. This retrograde condensation has been understood for years, but accurate correlation methods to permit quantitative prediction of phase behavior in the retrograde region are not presently available. The increasing importance of natural gas has made accurate prediction of phase behavior and composition of produced natural gas streams an economic necessity. The work reported here was undertaken to provide a rapid, economical method for obtaining the vapor-liquid equilibrium information needed to predict accurately the composition of the fluids produced from a gas reservoir throughout its life. TO develop this method, a pressure cell equipped with windows was designed and built for observing the volumes of liquid and gas present at reservoir pressures and temperatures. Use was made of established chromatographic methods for rapid and detailed analysis of both phases. This paper describes the equipment and techniques developed for obtaining vapor-liquid equilibrium data, presents the results of analyses of a condensate system, and indicates the usefulness of these data in predicting hydrocarbon phase behavior. DESCRIPTION OF EQUIPMENT USED The equipment used in obtaining the required information on phase behavior and the complete analysis of hydrocarbon mixtures through C, will be described first, followed by a discussion of the operation of the equipment. It will be helpful, however, to consider first a brief outline of the technique used. A sample of separator gas and liquid is charged to the windowed cell (see Figs. 1 and 2) where volumetric equilibrium phase behavior at reservoir pressures and temperatures can be determined. Then samples of the coexisting phases are withdrawn. The methane-through-pentane fraction is analyzed with a chromatograph equipped with a hot-wire detector, and the pentane-plus fraction is analyzed with a second chromatograph equipped

Jan 1, 1965

By D. E. Thomas, C. E. Birchenall

ALTHOUGH Eoltzmann gave a mathematical solution for the diffusion equation (for planar diffusion in infinite 01. semi-infinite systems only) in 1894 allowing for variation of the diffusion coefficient with a change in concentration, it was not until 1933 that this solution was applied to an experimentally investigated metallic system. The calculation was carried out by Matano' on the data obtained by Grube and Jedele3 for the Cu-Ni system. Since that time concentration dependence of the diffusion coefficient has been demonstrated for many pairs of metals. However, the nature of this dependence has never been fully elucidated. Many investigators have suspected that these variations could be related to the thermodynamic properties of the solutions, one of the earliest explicit statements being contained in a discussion of irreversible transport processes by Onsager' in 1931. Development along these lines has been greatly retarded by the lack of reliable data on the variation of tliffusivity with concentration, the paucity of the thermodynamic data for the same systems at the same temperatures and compositions, and an incomplete understanding of the relation of the thermodynamic properties of the activated state for diffusion to the bulk thermodynamic properties. The last factor has been discussed by Fisher, Hollomon, and Turnbull.5 In many instances where data exist, it is difficult to know which are acceptable. This problem probably applies more strongly to diffusion data than to activity measurements. For instance, four sets of observers"-" have reported self-diffusion coefficients for copper. The average spread between extreme results is a factor of about four, though the individual sets of data are self-consistent to about 20 pct. Thus one or more factors are out of control, at least in these experiments, making estimates of internal error unreliable. The most reliable diffusion data in most systems have resulted from the use of welded couples with a plane interface from which layers for analysis are machined parallel to the interface after diffusion. The layers are analyzed, and the result is a graphical relation between distance and concentration, usually called the penetration curve. Given the same set of analytical data and distances and following the same procedure in computation, different observers will generally produce diffusion coefficients which vary appreciably, especially at the extremes of the concentration range. Experiments must be carefully designed so that the precision is good enough to answer a particular question unequivocally. In the first calculation of the dependence of the diffusion coefficient on concentration in the metallic solid solution Cu-Ni, Matano found that the coefficient was insensitive to concentration from 0 to 70 pct Cu, after which it rose more and more steeply to some undetermined value as pure copper was approached.' The same behavior was reported for Au-Ni, Au-Pd, and Au-Pt.* The data used were those of Grube and Jedele which were very good at the time, but are not considered particularly good by present standards. Furthermore, the method of calculation makes the ends of the diffusion coefficient-concentration curve unreliable. For better reliability, the high copper end of the curve has been checked by incremental couples, where the concentration spread is 67.7 to 100 atomic pct Cu. The implication of the curves calculated by Matano was that diffusion is very concentration sensitive in one dilute range of this completely isomorphous system and hardly at all in the other. Matano's result is confirmed. Later Wells and Mehll0 published data on diffusion in Fe-Ni at 1300°C, which represent a thorough test of the shape of the concentration dependence curve. They ran couples with the following ranges of nickel concentration: 0-25 pct, 1.9-20.1 pct, 0-20.1 pct, 20.1-41.8 pct, 0-99.4 pct, and 79.3-99.4 pct. Although the trend of the data indicates an S-shaped concentration dependence, their curve was drawn to the pattern set by Matano. Their original data have been recalculated for the 0-99.4 and 79.3-99.4 pct couples. Wells and Mehl's points and two independent recalculations from the raw data are plotted in Fig. 1. What appears to be the best curve is drawn through them. This curve shows little sensitivity to composition in both dilute ranges with a strong dependence at intermediate composi-tions.? Similar experiments on the Cu-Pd system are reported here at temperatures where solubility is unlimited. These lead to the same type of concentration dependence for the diffusion coefficients as was found upon recalculation of the data for the Fe-Ni system. Experimental Procedure Cu-Pd: The concentration dependence of the diffusion coefficient may be determined by the use of

Jan 1, 1953

By Rodney P. Smith

The diffusivity of carbon (0.1 wt pct C) in Fe-Nz alloys (0 to 100 pct Ni) has been determined for the temperature range 860° to 1100°C. As a function of nickel content, the diffusivity has a maximum near 60 pct Ni (the maximum diffusivity being about 1.3 times that in the absence of nickel); the activation energy has a maximum between 40 and 50 pct Ni and a maximum between 80 and 90 pct Ni. The difference between the minimum activation energy and that in iron is about 3000 cal pev g-atom; Do has a minimum between 40 and 50 pct Ni and a maximum between 80 and 90 pct Ni. The results cannot be rationalized by an approximate thermodynamic treatment. THE diffusivity of carbon has been determined in a number of iron alloys over a limited concentration range. It seemed desirable to investigate a system which allows an extended range of alloy composition within a single-phase region. The Fe-Ni system is ideal in this respect, in that all alloys from 100 pct Fe to 100 pct Ni are fee in a convenient temperature range.&apos; The carbon diffusivity was determined by a decar-burization method. The experimental procedure was identical with that used to determine the diffusivity of carbon in y Fe-Co alloys.2 The experimental data are given in Table I. A small correction (order of a few percent) has been made to the measured carbon loss to correct for the carbon lost from the ends of the cylinders.&apos; Since the diffusivity of carbon varies with carbon content the measured diffusivity is an average value for a carbon content between zero (surface) and that at the center of the sample at the end of the decarburization periods. In making the correction in D to 0.1 wt pct C it is assumed that the measured D corresponds to the arithmetical mean of the carbon content at the surface and at the center of the sample at the end of the decarburization period.3 Since this correction is small (<4 pct in D) and since for our decarburization times the changes in carbon content at the center of the sample was small the mean carbon content could have been taken as half the initial value. It is further assumed that the change in D with carbon content for the alloys is the same as that for the diffusion of carbon in iron. From the data of Wells, Batz, and Mehl4 and of smith5 the correction of D from the mean carbon content to 0.1 wt pct C is 0.3 (0.1 - mean wt pct C). The results for iron are given in Ref. 2. Within the experimental error log Do.l%C for each alloy is a linear function of 1/T; the constants for the equation determined by the method of least squares are given in Table I. The deviations of the experimental points from the least-squares line are of the order of 2 pct in D. A comparison of our results for the diffusivity of carbon in nickel with those of other investigations is shown in Fig. 1. The lower curve in Fig. 1 is a linear extrapolation of values calculated* from the equation of Diamond6 for the relaxation time (temperature range 100° to 500°C). The results indicate a small increase in the activation energy over the temperature range 100° to 1400°C; however, it is difficult to say whether the change in Q is real or experimental error. Certainly the change in Q is less than the variation of 5 kcal per g-atom in the diffusivity of carbon in a iron.6 The experimental data for all the alloys are plotted in Fig. 2. As a function of nickel content the diffusivity has a maximum near 60 wt pct Ni at all temperatures investigated and possibly a minimum between 80 and 90 wt pct Ni for temperatures below 1000°C. The activation energy, Q, and log Do are plotted as a function of the nickel content in Fig. 3. Due to the limited temperature range of our experiments neither Q nor Do can be determined precisely; the activation energies appear to be consistent to ±0.3 kcal per g-atom; however the deviation from the absolute values may be considerably larger, see Table II. The Do values probably have little significance. The solid line for Do in Fig. 3 represents the values required to reproduce the experimental values for D when Q has values represented by the upper solid line The diffusivity of carbon may be expressed in terms of the mobility B22, the activity coefficient r2,

Jan 1, 1967

By A. J. Lena, M. F. Hawkes

Changes in hardness, tensile properties, microstructure, electrical resistance, and X-ray diffraction effects indicate that lattice strains are necessary for the embrittlement of ferritic stainless steels when heated for relatively short times at 475°C (885°F). It is suggested that 475°C (885°F) embrittlement is due to the accelerated formation of an intermediate stage in the formation of s under the influence of these strains. FERRITIC stainless steels (low carbon alloys of iron with more than 15 pct Cr) are subject to two forms of embrittlement when heated in the temperature range of 375° to 750°C. The embrittlement which occurs after long time heating between 565" and 750°C is well understood; it is caused by the precipitation of the hard, brittle s phase. Sigma is an intermetallic compound of approximate equi-atomic composition with an extended range of formation in Fe-Cr alloys. The maximum temperature at which this form of embrittlement can occur is dependent upon chromium content; and is approximately 620°C for a 17 pct Cr steel and 730°C for a 27 pct Cr steel. The other form of embrittlement occurs after relatively short heating periods in the range of 375" to 565°C; in the higher chromium steels, hours may be sufficient as compared to months for s embrittlement. This phenomenon is not at all well understood and several controversial theories have been proposed. The rate and intensity of embrittlement increase with increasing chromium content but the maximum rate occurs at 475°C re-gardless of chromium content. As a result of this, the phenomenon has been termed 475°C (885°F) embrittlement. The effect of 475°C embrittlement on the properties of ferritic stainless steels has been thoroughly reviewed by Heger.1 The embrittlement causes a pronounced decrease in room temperature impact strength and ductility, a large increase in hardness and tensile strength, and a decrease in electrical resistivity and corrosion resistance. Microstructural changes accompanying embrittlement are minor and difficult to interpret with a general grain darkening, appearance of a lamellar precipitate, grain boundary widening, and precipitation along ferrite veins having been reported at various times. With the exception of reported line broadening, X-ray diffraction studies by conventional Debye analysis of solid samples have been of little value. BY making use of electron diffraction methods, Fisher, Dulis, and Car-roll&apos; have recently shown the existence of a chromi-um-rich, body-centered cubic phase in 27 pct Cr steels which had been aged at 482°C (900°F) for as long as four years. Two types of theories have been advanced to account for the embrittlement. The first of these requires the precipitation of a phase not inherent in the Fe-Cr system with various investigators suggesting a carbide,3 nitride,3 phosphide,4 or oxide." Theories of this type have difficulty accounting for the influence of alloying elements on the embrittlement and for the facts that a minimum chromium content is necessary for embrittlement and the intensity of embrittlement increases with increasing chromium content. The second type of theory that has been proposed relates 475°C embrittlement to s phase formation which is inherent in the Fe-Cr system. An assumption of this kind can adequately explain the influence of alloying elements, for they exert an effect on 475°C embrittlement similar to that on s phase for-mation as can be seen in Table I. The minimum chromium content is essentially the same for both phenomena and it has been shown12,13 that s is a stable phase in the embrittling temperature range. In addition, it has been reported14,15 that pure alloys embrittle to the same extent as commercial type alloys. There are, however, several factors which have prevented complete acceptance of a s phase theory. Foremost of these is that the embrittlement can be removed by reheating for short time periods above 600°C, which in the higher chromium steels is within the stable s region. No s has ever been observed after one of these curing treatments, nor has any s been found as a result of embrittlement at 475°C. In addition, the simple precipitation of s cannot explain the time-temperature relationships for reactions between 350°and 750°C. This behavior is shown schematically in Fig. 1. Newell 16 and Ried-rich and Loib4 have shown that 475°C embrittlement follows a C-type curve as illustrated, while Short-

Jan 1, 1955