Search Documents

By M. S. Azun

I do not at all agree with the basic points of the author's conclusion. The use of lognormal or normal model to respond to the attribute distribution function should be carefully questioned. If fitting a distribution function of regionalized variables is the main purpose, one may employ one of certain well-known families of distributions, such as Tukey's lambda functions. For instance, the distribution function of gold assay values studied by Krige (Krige, 1978) may fit the lognormal model well, but one should not exclude that the basic idea of using the distribution function of observed data is to have insight of the physical mechanism that gives rise to the regionalized variables. In Applied Geostatistics, there are many examples considering only lognormal or normal model for the distribution function of geostatistical data even if the observed distribution function has more than one mode. Of course, using lognormal or normal model always brings easy computation for the statistical properties of attribute under investigation. Still, I object to adopting any of those models without having any "statistical inference." As seen from the author's paper, lognormal and normal model produce the same average grade (0.746% copper) and almost the same sample variance. Therefore, which model explains the probabilistic behavior of ore deposit is always an important question that should be replied by the geostatistician using not the common sense but the statistical and probabilistic methods. For kriging procedures, such as linear kriging, one should describe either second order stationary properties or bivariate properties of regionalized variables. Second order properties, such as correlogram function or semivariogram function, are the inputs of the so-called kriging equations. Selecting one of any useable models for the observed second order properties, such as spherical model, is not easy since the models are not based on the correlations structure of the regionalized variables. How those models being used in Applied Geostatistics can be distinguished is another important problem. One may develop many fitting models to respond to the observed second order properties of regionalized variables. However, I suggest that the model should be based on the probabilistic behavior of geologic process. The author used spherical model for copper, wolfram, and silver grades from isotropic sector (Fig. 2, 4, 3b) and exponential model for silver grades obtained along raises in isotropic sector (Fig. 3a). It seems that those semivariograms given in Fig. 3b for isotropic sector and Fig. 4 may be described by the random model implying that regionalized variables do belong to renewal process (Azun, 1983). If I knew the total number of data points used in estimation, may prove that my conclusion is correct. Recalling the test for the first order sample correlogram value, the test procedure is introduced such that the independence is rejected if [Ir(1) I > Za/2 N1/2] where Za/2 is the a/2 percentile of standard normal random variable and N is the total number of samples used in estimation (Azun, 1983). For the other observed semivariogram functions I suggest the author might try to use the so-called "Markovian model" describing not only the correlation structure of regionalized variables, but also the physical mechanism that produces the regionalized variables. The Markovian model for correlogram function is, [p (h) = T 1ih , h>1 ,] The Markovian model for the other second order properties of ReV are also derived. [T veß], in the above equation, show, for example, structural change and mineralizational variation in a considered deposit, respectively (Azun, 1983). In selecting any model, it is not easy to search for the "best" response to the observed second order properties of ReV's. The Markovian model, based on a theoretical understanding of underlying mechanism, gives more information about the occurrence of regionalized variables and respond to all properties. Random model and the so-called hole effect structure can be easily defined as a special form of the Markovian model (Azun, 1983). An estimation of any second order properties of regionalized variables may be computed through (N-1) lag. However, the dependency between the regionalized variables may be deemed to have "died out" after the so-called stationarity width (range). In practice, the estimation is carried out through 25% of the total number samples available. Therefore, the large fluctuations around the zero level and variance for correlogram or covariogram function and semivariogram function, respectively, cannot be observed. Since the number of pairs involved in the estimation of higher second order properties is small, the estimation variance at those lags is large. There is no need to show higher order estimated second order properties. After using one of any kriging procedure for block estimation, the distribution function of average grades has the same mean as the drilling sample grades but

Jan 1, 1986

By H. E. Stauss, J. Hino

INTEREST in silicon has arisen again in the past decade as a result of improvements in crystal rectifiers.' Although the preparation of silicon was first reported by Berzelius in 1880, the early product was of relatively low purity, and only the need for rectifiers in World War II led to the production of a 99.9+ pct pure powder. This material in crystalline form was consolidated into massive silicon for use, and the method developed was to melt it with selected added constituents as "doping" agents. Melting techniques, therefore, are of great importance. There are two basic problems in producing silicon ingots free of doping additions; one is the prevention of spitting and the other is prevention of cracking of the ingot during freezing. The most satisfactory arrangement yet developed for producing massive silicon is to melt and freeze in a cylindrical quartz crucible surrounded by a concentric heating element and concentric radiation shields or insulation. For example, use can be made of a tubular heater with a high frequency generator as the source of power and reflecting shields of alundum cylinders. The spitting of silicon is related to gas evolution, and the gas comes from two primary causes—adsorbed gas and the reaction products of silicon and the crucible. Gas is also released from bubbles contained in the quartz crucible walls. Improved removal of adsorbed gas can be achieved by means of controlled melting and freezing. The seriousness of the problem in vacuo is reduced with an electrically operated mechanical movement of the high frequency power coil. The upper portion of the powder charge is melted first and the high frequency coil lowered until the powder is completely molten. During cooling the high frequency coil is raised slowly. These means also reduce the final nonviolent extrusion of large beads of metal through the ingot top during freezing. Better control of spitting and bead extrusion is obtained when melting is done under helium at. atmospheric pressure instead of in vacuo. The problem of reaction between silicon charge and crucible in practice is confined to the reaction between silicon and quartz. This2 apparently is: Si + SiO2 + 2SiO The part that this reaction plays in spitting has not been isolated for separate study. SiO is a volatile vapor at the melting point; of silicon and is released freely during melting in vacuo, but hardly at all in helium at atmospheric pressure. The cracking of ingots is a major difficulty in melting silicon, and its prevention requires special melting techniques or the addition of "toughening" agents such as aluminum or beryllium.' The cracking of the ingots has been explained as being the result of the expansion that occurs upon freezing; although direct observation of freezing ingots reveals visible cracks on the surface only after a red heat has been reached, suggesting that cracking is the result of differential contraction of silicon and quartz. Silicon wets quartz, and the ingot adheres tightly to the crucible. Therefore as ingot and crucible cool, the two either have to pull apart, or at least one must crack. Surprisingly, in spite of the relative thinness of the quartz and the thickness of the ingot, the ingot and the crucible both crack. Microscopic and X-ray4 studies fail to show any plastic flow other than twinning in the ingots. Slow cooling fails to prevent cracking. Another possible solution to cracking is to weaken the crucible. Use of thin-walled crucibles finally led to success with fused quartz crucibles with a wall thickness of 0.25 to 0.50 mm. With such thin-walled fused quartz crucibles consistently uniform success is secured in producing sound ingots 30 mm in diam from the purest available grade of silicon (99.9+) without the use of any type of addition. Melts are made in the size range of 50 to 100 g. Omission of a deliberately added doping agent is not sufficient to insure pure ingots. The reaction of silicon with crucibles and the resultant solution of impurities in the silicon is well-established." In this laboratory, the presence of Al, Be, and Zr has been found spectroscopically in ingots melted in contact with alumina, beryllia, and zircon. The best crucible materials reported in the literature are MgO and SiO2. Use of MgO in this laboratory has resulted in a heavy deposit of magnesium on the furnace walls, showing that a reduction of the magnesia occurred and the resulting magnesium removed from the melt by volatilization. In the case of quartz, the silica is reduced and SiO liberated to deposit on the equipment walls. There probably is real danger that oxygen is dissolved in the ingot when either magnesia or silica is used as the crucible material. Preliminary analyses by Dean Walter in his vacuum unit in this laboratory6 indicate the presence of oxygen in undoped silicon melted in quartz.

Jan 1, 1953

By P. A. Laxen, H. R. Spedden

WHEN a mineral particle is fractured, bonds between the atoms are broken. The unsatisfied forces that appear at the newly formed surface are considered to be responsible for the adsorption of ions at the mineral surface. A knowledge of the mechanism and extent of ion sorption from solution onto a mineral surface is of interest in the development of the theory of flotation.'*' Study of the adsorption of sodium from an aqueous solution oftheon quartz offers a simple approach to this complicated problem. The availability of a radioisotope as a tracer element meant that accurate data could be obtained."." Three main factors which appeared likely to affect the adsorption of sodium are: l—concentration of sodium in the solution, 2—concentration of onotherof cations in the solution, and 3—anions present in the solution. Hydrogen and hydroxyl ions are always present in an aqueous solution. By controlling the pH, the concentration of these two ions was kept constant. The variation in thesethe amount of sodium adsorbed with variation in sodium concentration was then determined under conditions standardized in regard to hydrogen ion. The effect of concentration of hydrogen ions and of other cations was also measured. A few experiments were made to get a preliminary idea on the effect of anions. The active isotope of sodium was available as sodium nitrate. Standard sodium nitrate solutions were used throughout these experiments except when the effects of other anions were studied. It was found that sodium adsorption increased with sodium-ion concentration, but less rapidly than in proportion to it. Increasing hydrogen-ion concentration, or conversely decreasing hydroxyl-ion, brings about a comparatively slight decrease in sodium-ion adsorption. Increasing the concentration of cations other than hydrogen or sodium decreases somewhat the adsorption of sodium ion. It would appear as if the kind of anion is a secondary factor in guiding the amount of sodium ion that is adsorbed. Materials and Methods Quartz The quartz was prepared as in previous work in the Robert H. Richards Mineral Engineering Laboratory' except for the refinement of using de-ionized distilled water for the final washing of the sized quartz, prior to drying." To minimize the laborious preparation of quartz, experiments were made to determine .whether the sodium-covered quartz could be washed free of sodium and re-used. The experiments were successful as indicated by lack of Na" activity on the repurified material and by its characteristic sodium adsorption. Table I gives the spectrographic analyses of the quartz used. The quartz ranged from 16 to 40 microns in size, averaging about 23 microns (microscope measurement), and had a surface of 1850 sq cm per g (lot I), 2210 (lot 11) and 2000 (lot 111) as determined by the Bloecher method." Radioactive Sodium Method of Beta Counting for Adsorbed Sodium: Na22, the radioisotope of sodium, possesses convenient properties.' It has a half-life of 3 years, thus requiring no allowance for decay during an experiment. On decay it emits a 0.575 mev ß+ radiation and a 1.30 mev r radiation. The decay scheme is illustrated in the following equation: ß+ NaR-------'8'77NeZ2 3 years The /3 radiation is sufficiently strong to penetrate an end-window type of Geiger-Mueller counting tube. This, in turn, makes it possible to use external counting, a great advantage in technique. Furthermore, it permits the assaying of solids arranged in infinite thickness, while assaying evaporated liquors on standardized planchets. The equipment used was standard and similar to that employed by Chang.R The original active material was 1 ml of solution containing 1 millicurie of Na" as nitrate. This active solution was diluted to 1000 ml. Five milliliters of this diluted active solution was found to give a quartz sample a sufficiently high activity for accurate evaluation of the sodium partition in the adsorption measurements. Also, 1 ml of final solution gave a sufficiently high count for precision on the liquor analyses. The sodium concentration of the diluted active solution was 1.2 mg per liter, so that 6 mg of sodium for 60 ml of test solution and 12 g of quartz was the minimum amount used. The active solution was stored in a Saftepak bottle. Procedure for Adsorption Tests: The method consisted of agitating 12 g of quartz with 60 ml of solution of known sodium concentration for enough time to establish equilibrium between the solution and the quartz surface. The quartz was separated as completely as possible from the solution by filtering and centrifuging. The activity on the quartz and in the equilibrium solution was measured and the partition of the sodium was calculated from the resulting data. The detailed procedure for the adsorption test is set forth in a thesis by Laxen." In brief, it included the following steps: 1—Ascertainment of linearity between concentration of Na" and activity measured. 2—Evaluation of factor to translate activity on solid of infinite thickness in terms of activity on an evaporated active film of minute thickness, on the various shelves of the counter shield. 3—Taking precautions to avoid evaporation of water during centrifuging

Jan 1, 1953

By P. A. Laxen, H. R. Spedden

WHEN a mineral particle is fractured, bonds between the atoms are broken. The unsatisfied forces that appear at the newly formed surface are considered to be responsible for the adsorption of ions at the mineral surface. A knowledge of the mechanism and extent of ion sorption from solution onto a mineral surface is of interest in the development of the theory of flotation.'*' Study of the adsorption of sodium from an aqueous solution oftheon quartz offers a simple approach to this complicated problem. The availability of a radioisotope as a tracer element meant that accurate data could be obtained."." Three main factors which appeared likely to affect the adsorption of sodium are: l—concentration of sodium in the solution, 2—concentration of onotherof cations in the solution, and 3—anions present in the solution. Hydrogen and hydroxyl ions are always present in an aqueous solution. By controlling the pH, the concentration of these two ions was kept constant. The variation in thesethe amount of sodium adsorbed with variation in sodium concentration was then determined under conditions standardized in regard to hydrogen ion. The effect of concentration of hydrogen ions and of other cations was also measured. A few experiments were made to get a preliminary idea on the effect of anions. The active isotope of sodium was available as sodium nitrate. Standard sodium nitrate solutions were used throughout these experiments except when the effects of other anions were studied. It was found that sodium adsorption increased with sodium-ion concentration, but less rapidly than in proportion to it. Increasing hydrogen-ion concentration, or conversely decreasing hydroxyl-ion, brings about a comparatively slight decrease in sodium-ion adsorption. Increasing the concentration of cations other than hydrogen or sodium decreases somewhat the adsorption of sodium ion. It would appear as if the kind of anion is a secondary factor in guiding the amount of sodium ion that is adsorbed. Materials and Methods Quartz The quartz was prepared as in previous work in the Robert H. Richards Mineral Engineering Laboratory' except for the refinement of using de-ionized distilled water for the final washing of the sized quartz, prior to drying." To minimize the laborious preparation of quartz, experiments were made to determine .whether the sodium-covered quartz could be washed free of sodium and re-used. The experiments were successful as indicated by lack of Na" activity on the repurified material and by its characteristic sodium adsorption. Table I gives the spectrographic analyses of the quartz used. The quartz ranged from 16 to 40 microns in size, averaging about 23 microns (microscope measurement), and had a surface of 1850 sq cm per g (lot I), 2210 (lot 11) and 2000 (lot 111) as determined by the Bloecher method." Radioactive Sodium Method of Beta Counting for Adsorbed Sodium: Na22, the radioisotope of sodium, possesses convenient properties.' It has a half-life of 3 years, thus requiring no allowance for decay during an experiment. On decay it emits a 0.575 mev ß+ radiation and a 1.30 mev r radiation. The decay scheme is illustrated in the following equation: ß+ NaR-------'8'77NeZ2 3 years The /3 radiation is sufficiently strong to penetrate an end-window type of Geiger-Mueller counting tube. This, in turn, makes it possible to use external counting, a great advantage in technique. Furthermore, it permits the assaying of solids arranged in infinite thickness, while assaying evaporated liquors on standardized planchets. The equipment used was standard and similar to that employed by Chang.R The original active material was 1 ml of solution containing 1 millicurie of Na" as nitrate. This active solution was diluted to 1000 ml. Five milliliters of this diluted active solution was found to give a quartz sample a sufficiently high activity for accurate evaluation of the sodium partition in the adsorption measurements. Also, 1 ml of final solution gave a sufficiently high count for precision on the liquor analyses. The sodium concentration of the diluted active solution was 1.2 mg per liter, so that 6 mg of sodium for 60 ml of test solution and 12 g of quartz was the minimum amount used. The active solution was stored in a Saftepak bottle. Procedure for Adsorption Tests: The method consisted of agitating 12 g of quartz with 60 ml of solution of known sodium concentration for enough time to establish equilibrium between the solution and the quartz surface. The quartz was separated as completely as possible from the solution by filtering and centrifuging. The activity on the quartz and in the equilibrium solution was measured and the partition of the sodium was calculated from the resulting data. The detailed procedure for the adsorption test is set forth in a thesis by Laxen." In brief, it included the following steps: 1—Ascertainment of linearity between concentration of Na" and activity measured. 2—Evaluation of factor to translate activity on solid of infinite thickness in terms of activity on an evaporated active film of minute thickness, on the various shelves of the counter shield. 3—Taking precautions to avoid evaporation of water during centrifuging

Jan 1, 1953

By L. D. Hall, D. R. Mash

The paper presents the results of experiments on the anelastic. behavior of gold, as manifested by grain boundary relaxation. Two grain boundary internal friction peaks are found for 99.9998 pct Au. It is found that the peaks are associated with primary and secondary recrystallization. However, the existence of two discrete peaks cannot be explained on the basis of grain size and shape alone. It is suggested that grain boundary stability, as determined by orientation, plays a role in the observed effects. EVIDENCE for the viscous behavior of grain boundaries in metals has been presented in recent years by several investigators, based upon studies of various anelastic effects, especially internal friction. KG1 has contributed greatly to this field, having put forward a coherent body of evidence for stress relaxation by the viscous intercrystalline flow mechanism. In this connection, he has made extensive use of pure aluminum (99.991 pct) as the test material, although he has also studied other metals and alloys, including pure iron (Puron).² Rotherham, Smith, and Greenough³ have studied the internal friction of pure tin, interpreting their results in a manner similar to that of KG. In view of the importance of such studies in shedding light upon the fundamental structure and behavior of the grain boundaries in pure metals, it appears that the use of a very pure test material which is inert to its environment should provide useful information on anelastic properties and the source of such behavior in pure metals. The present work was carried out on spectrograph-ically pure, 99.9998 pct Au, free of all impurities except for a trace of silver, estimated to be present to the extent of about 0.0002 pct. The term "pure gold" will hereafter refer to this very pure material. Gold of commercial purity, 99.98 pct, was also studied to observe the effects of small amounts of impurities. A pure gold "single crystal" specimen was also tested for comparison. The variation of the internal friction and rigidity modulus as a function of temperature was determined by means of a torsion pendulum apparatus employing extremely low stress amplitudes and a frequency of vibration of the order of 1 cycle per sec. A 12 in. length of 0.031 in. (20 gage) gold wire formed the suspension element. The apparatus was similar to that described by Ke.l The test procedure and the basic requirements to be met for obtaining useful experimental data by this method have been given elsewhere.1,2 It should be made clear that in all of the experiments to be described, the internal friction and rigidity were independent of the amplitude of torsional vibration. The semilog plot of amplitude of vibration vs ordinal number of vibration was a straight line. This was carefully verified for each internal friction measurement. The linear variation shows that the internal friction was independent of stress; i.e., that the specimens were not being cold-worked during testing. The reproducibility of the internal friction curves, which were obtained by cyclic heating and cooling, indicates that the gold was unaffected by its environment during the tests. The measure of internal friction adopted in the present study is the conventional "logarithmic decrement," defined as follows: log. dec. = l/n In A0/An [I] where n is the number of cycles or vibrations; A,, the initial amplitude of vibration; and An, the amplitude after the nth cycle. When the logarithmic decrement is small, the shear modulus, G, of the wire is proportional to the square of the frequency of vibration provided the length and radius of the wire are kept constant. A plot of frequency squared vs temperature gives the following ratio:' This expresses the fraction of the stress which has not been relaxed at a given temperature. Gr and Gv are the relaxed and unrelaxed moduli, respectively. The frequency of vibration in the polycrys-talline specimen is fp, and the frequency of vibration of a single crystal is f8. This latter quantity is obtained simply by extrapolating the linear, low temperature portion of the curve of frequency squared vs temperature for the polycrystalline specimens. The theory of viscous grain boundary stress relaxation as demonstrated by the anelastic behavior of metals has been discussed in detail by Zener4 and need not be reproduced here. Experimental Results Initial measurements of the internal friction of pure gold were carried out on specimens which had been drawn with no intermediate annealing, resulting in a material which had undergone approximately 99 pct reduction of area in final processing. Annealing was then carried out at successively higher temperatures starting at 400°F for 1 hr and proceeding in this manner to as high as 1600°F in 100°F intervals. After each annealing treatment an internal friction and rigidity vs temperature curve was obtained over the range from room temperature to the particular annealing temperature. The resulting internal friction curves did not exhibit well defined maxima (peaks), but rather several fairly flat

Jan 1, 1954

By George S. Ansell, John B. Hudson, Stephen A. Koffler

The permeability of hydrogen through the a phase of palladium has been measured by a low pressure permeation technique under conditions such that bulk diffusion was the rate-controlling process. The observed permeability is described by the equation: J = 1.80 x 1O-3 P½exp(-3745)/RT) cc(stp)/sec/cm2/en over a range of hydrogen pressure, P, from 2.9 x 10-5 m Hg to 5.0 x 10-3 cm Hg, and over a temperature range, T, from 300" to 709°K. The fact that the permeability shows a square root dependence on pressure and a reciprocal dependence on thickness was taken as evidence that bulk diffusion, rather than surface reactions, was the rate-controlling process. The permeability data were used in conjunction with the solubility data of Salmon et al., to determine the diffusivity of hydrogen through palludium as: D = 4.94 x 10-3 exp(-5745/RT) cm2/sec There was no influence of sub structural defects observed over the temperature range employed. From the permeability data obtained, coupled with grain size measurements, it was concluded that the ratio of grain boundary diffusivity to bulk diffusivity was less than 105 over the range of temperature investigated. THE diffusive mass transport of hydrogen in the a phase of palladium has been studied previously by numerous investigators.' In spite of the large amount of attention this system has received, there is not good agreement between the results obtained in different investigations.' This is due in part to the fact that the mass transport was surface-limited during some of these studies, rather than being diffusion-controlled3-5 The reason for the disagreement in other cases is not clear. These studies made use of such techniques as rate of absorption from solutions6 and gases,' electrochemical potential,8 time lag,' and permeation10,11 to determine the mass transport behavior. Of these, the gas permeability technique is the only method which allows an easy test to determine if diffusion is the rate-controlling mechanism, thus eliminating the uncertainty regarding the limiting transport processes inherent in the other techniques. The two most recent permeability studies are those of Toda10 and Davis." Toda determined the permeability of hydrogen in the a phase of palladium over the temperature range from 170" to 290°C, and over the pressure range from 36 to 630 mm Hg utilizing a steady-state gas-permeability technique. Toda's result was: J = 1.41 x 10-3 P½ exp(-3220/RT) where J = specific permeability in cc(stp)/sec/cm2/cm, P½ = square root of the inlet pressure in (cm Hg)½, R = Universal gas constant, and T = temperature in deg Kelvin. Davis11 also employed a steady-state gas-permeability technique over the temperature range from 200" to 700°C and over the pressure range from 0.02 to 760 mm Hg. His result for the permeability of hydrogen in the a phase of palladium was: J = 3.15 x 10-3 P½ exp(-A440/RT) In the range of overlapping temperature for these two investigations, the values of the specific permeability calculated from the above two equations differ by a factor of about 1.8. In the present investigation, the permeation of hydrogen in the a phase of palladium was determined over a wide temperature range, 27" to 436oC, and over the pressure range from 2.9 x 10-5 to 5.0 x 10-3 cm Hg. This temperature range overlaps that of the previous investigations of Toda and Davis, but also covers the lower temperature range which has never before been investigated. The lower pressure range used here avoided the interaction between the dissolved hydrogen atoms observed at higher hydrogen concentrations.' MATERIALS The as-received palladium specimens were cold rolled from a casting and were supplied as 5.08 cm discs of 0.508 and 0.762 mm thickness. According to J. Bishop and Company specifications, the composition of the discs was 99.95 pct Pd, the balance being Cu, Ag, Au, and Ir. In order to obtain samples of varying grain size, the as-received discs were then heat treated. Sample 1 was treated for eight minutes in a nitrogen atmosphere at 810°C and then air-cooled. Sample D was heated first to 550°C in helium. The helium atmosphere was then immediately replaced by hydrogen and the temperature was slowly raised to 1220°C and held for 22 hr. The sample was then cooled to 550°C where the hydrogen was replaced by helium and the disc was further furnace-cooled to room temperature. Sample H was given the same heat treatment, only with hydrogen substituted for helium below 550°C. The reason that hydrogen was not used throughout the entire annealing cycle for samples D and H was to prevent the distortion encountered by low-temperature cycling in hydrogen observed by Darling." After these heat treatments were completed, the

Jan 1, 1970

By Francis Collins, E. R. Brownscombe

A study has been made of the material balance-fluid flow method of estimating reserves and degree of water drive from pressure and production history data. By considering the effect of random pressure errors it is shown that in a particular example a standard deviation of three and one-half pounds in each of ten pressure survey? permits the determination of the reserves with a standard deviation of 8 per cent and the water drive with a standard deviation of 15 per cent, assuming that certain basic geologic data are correct. It is believed that this method of estimating reserves and water drive is useful and reliable in a number of cases. The method is particularly valuable when reservoir pressure data are accurate within a very few pounds, but may also be applied with less accurate pressure data if a relatively large reservoir pressure decline occurs early in the life of the field, as for example in an under-saturated oil field. INTRODUCTION A knowledge of the magnitude of reserves and degree of water drive present in any newly discovered petroleum reservoir is necessary to early application of proper production practices. A number of investigators have contributed to methods of relating reserves, degree of water drive, and production and pressure history. 1-8 Three types of problems of increasing complexity may be mentioned. If a reservoir is known to have no water drive. and if the ratio of the volume of the reservoir occupied by gas to the volume of the reservoir occupied by oil (which ratio permits fixing the overall compressibility of the reservoir) is known, then only one further extensive reservoir property remains to be determined, namely the magnitude of the reserves. A straightforward application of material balance considerations will permit this determination. The problem becomes very much more difficult if we wish to determine not only the magnitude of the reserves but also the magnitude of water drive, if any, which is present. In principle, a combination of material balance and fluid flow considerations will permit this evaluation. Finally, if neither the magnitude of reserves, the degree of water drive, nor the ratio of oil to gas present in the reservoir is known and it is desired to determine all three of these variables, the problem could in principle be solved by a fluid flow-material balance analysis which determines the overall compressibility of the reservoir at various points in its history. The change in compressibility with pressure would provide a means of determining the ratio of gas to liquid present, since the compressibilities of gas and liquid vary differently with pressure variation. However, in practice this problem is probably so difficult as to defy solution in terms of basic data precision apt to be available.' It is the purpose of this discussion to illustrate the second case, which involves the determination of two unknown variables, single phase reserves and degree of water drive, from pressure and production history and fluid property data, and to study the precision with which these unknowns can be determined in this manner in a particular case. Although an electric analyzer developed by Bruce as used in making the calculations to be described, numerical methods necessary in carrying out the process have been devised and have been applied for this purpose. Schilthuis,' for example, developed a comprehensive equation for the material balance in a reservoir. He combined this with a simplified water drive equation, assuming that the ratio of free gas to oil was fixed by geological data and that a period of constant pressure operation at constant rate of production was available to determine the constant for his water drive equation. On this basis he was able to compute the reserves and predict the future pressure history of the reservoir. Hurst developed a generalized equation permitting the calculation of the water drive by unsteady state expansion from a finite aquifer. He showed in a specific case how the water influx calculated by his equation, using basic geologic and reservoir data to fix the constants, matched the water influx required by material balance considerations. Old3 illustrated the simultaneous use of Schilthuis' material balance equation and Hurst's fluid flow equation for the determination of the magnitude of reserves and a water drive parameter from pressure and production history. He used this method to calculate the future pressure history of the reservoir under assumed operating conditions. As a basis for determining reserves, Old assumed a value for his water drive parameter and calculated a set of values for the reserves, using the initial reservoir pressure and each successive measured pressure. The sum of the absolute values of the deviations of the resulting reserve numbers from their mean value was taken as a criterion of the closeness of fit to the experimental data possible with the water drive parameter assumed. New values of the water drive parameter were then assumed and new sets of the reserves calculated until a set of reserves numbers having a minimum deviation from the average was established. The average value of- the re-

Jan 1, 1949

By D. G. Russell

A simple method has been developed with which flowing bottom-hole pressure data from two-rate flow tests in oil or gas wells can be analyzed to estimate the formation permeability, skin factor and average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized producing rate of the well is changed to another, higher or lower, rate. The new method yields the same information as a conventional pressure buildup analysis, but eliminates the need for closing in the well. The analysis of a two-rate flow test is of the same degree of difficulty and requires about the same engineering time for application as a conventional pressure buildup analysrs. The extended closed-in periods experienced with conventional buildups because of long, low-rate afterproduc-tion periods are eliminated by fIow tests. Other anomalous pressure buildup effects, such as "humping" due to weli-bore phase segregation, can be successfully eliminated with the new method. Generally, flow tests of about 24 hours' duration run with conventional pressure measurement equipment are sufficient for interpretation purposes. Field testing of the two-rate flow test method has established that it is a reliable and economical method which can be used in many instances to complement or even replace conventional pressure buildup methods. INTRODUCTION The principal method for estimating formation permeability and well damage, or skin factor, in a producing oil or gas well is the analysis of shut-in bottom-hole pressure buildup data.' This familiar method has been used quite successfully by reservoir engineers for many years. It is based on the solution of the radial flow equation for constant rate conditions, and requires that the well be closed in for a sufficient period of time to obtain a clearly defined linear portion on the plot of observed t + ?t bottomhole pressure vs log t + ?/?t(where At is shut-in A? time, and t is producing time to the instant of shut-in). From the slope of the plot and other normally obtainable data, the permeability, skin factor, and reservoir pressure at infinite shut-in time (if the reservoir were infinite) can be estimated. Over the years several drawbacks have become apparent in the use of conventional shut-in pressure buildups for determining permeability and skin factor. The conventional pressure buildup interpretation theory assumes that a well is closed in at the sand face and that no production into the well occurs after shut-in. In practice, of course, the well is closed in at the surface, and inflow into the well continues until the well fills sufficiently to transmit the effect of closing in to the formation. This adjustment period is commonly referred to as the "after-production" portion of the pressure buildup. In the tight reservoirs, long, low-rate after-production periods frequently occur, and the well must be shut in for several days or, in some instances, even weeks to obtain an interpretable buildup curve.' Obviously, such long shut-in times can cause loss in current income, both from reduced oil production and from the fact that personnel and pressure measurement equipment are occupied with a single well for too long a time. In other cases, even long shut-in periods do not seem to be of much aid in obtaining an interpretable buildup." If there is considerable phase redistribution (liquid fallout or bubble rise) after a well is shut in, then curves with no interpretable portion are often obtained during the buildup. In addition to instances in which wellbore effects cause trouble, there are also cases in which the major objection to use of the closed-in pressure buildup is simply the fact that the well must be shut in. When there is no proration and when the well has limited producing capacity. closing in the well means loss of income. From the foregoing discussion it is apparent that it ib desirable to have an alternative method of obtaining the same information as that derived from a conventional buildup without the need of closing in the well. One possible solution which has been offered for this problem is the use of a bottom-hole shut-in tool' which isolates the major portion of the flow string from the formation face during the buildup. In this paper an alternative method which is frequently successful in avoiding wellbore effects and which does not require the use of special equipment is presented. A new, simple method has been developed with which the flowing bottom-hole pressure data from flow tests in oil or gas wells can be used to estimate permeability, skin factor and the average reservoir pressure. The required pressure data are obtained by observation of the transient bottom-hole pressure behavior after the stabilized production rate of the well is changed to another, higher or lower rate. The need for closing in the well is eliminated, and pressure measurement periods of only 18 to 24 hours arc usually sufficient, even in tight reservoirs. Thus, the new

By E. R. Brownscombe, Francis Collins

A study has been made of the material balance-fluid flow method of estimating reserves and degree of water drive from pressure and production history data. By considering the effect of random pressure errors it is shown that in a particular example a standard deviation of three and one-half pounds in each of ten pressure survey? permits the determination of the reserves with a standard deviation of 8 per cent and the water drive with a standard deviation of 15 per cent, assuming that certain basic geologic data are correct. It is believed that this method of estimating reserves and water drive is useful and reliable in a number of cases. The method is particularly valuable when reservoir pressure data are accurate within a very few pounds, but may also be applied with less accurate pressure data if a relatively large reservoir pressure decline occurs early in the life of the field, as for example in an under-saturated oil field. INTRODUCTION A knowledge of the magnitude of reserves and degree of water drive present in any newly discovered petroleum reservoir is necessary to early application of proper production practices. A number of investigators have contributed to methods of relating reserves, degree of water drive, and production and pressure history. 1-8 Three types of problems of increasing complexity may be mentioned. If a reservoir is known to have no water drive. and if the ratio of the volume of the reservoir occupied by gas to the volume of the reservoir occupied by oil (which ratio permits fixing the overall compressibility of the reservoir) is known, then only one further extensive reservoir property remains to be determined, namely the magnitude of the reserves. A straightforward application of material balance considerations will permit this determination. The problem becomes very much more difficult if we wish to determine not only the magnitude of the reserves but also the magnitude of water drive, if any, which is present. In principle, a combination of material balance and fluid flow considerations will permit this evaluation. Finally, if neither the magnitude of reserves, the degree of water drive, nor the ratio of oil to gas present in the reservoir is known and it is desired to determine all three of these variables, the problem could in principle be solved by a fluid flow-material balance analysis which determines the overall compressibility of the reservoir at various points in its history. The change in compressibility with pressure would provide a means of determining the ratio of gas to liquid present, since the compressibilities of gas and liquid vary differently with pressure variation. However, in practice this problem is probably so difficult as to defy solution in terms of basic data precision apt to be available.' It is the purpose of this discussion to illustrate the second case, which involves the determination of two unknown variables, single phase reserves and degree of water drive, from pressure and production history and fluid property data, and to study the precision with which these unknowns can be determined in this manner in a particular case. Although an electric analyzer developed by Bruce as used in making the calculations to be described, numerical methods necessary in carrying out the process have been devised and have been applied for this purpose. Schilthuis,' for example, developed a comprehensive equation for the material balance in a reservoir. He combined this with a simplified water drive equation, assuming that the ratio of free gas to oil was fixed by geological data and that a period of constant pressure operation at constant rate of production was available to determine the constant for his water drive equation. On this basis he was able to compute the reserves and predict the future pressure history of the reservoir. Hurst developed a generalized equation permitting the calculation of the water drive by unsteady state expansion from a finite aquifer. He showed in a specific case how the water influx calculated by his equation, using basic geologic and reservoir data to fix the constants, matched the water influx required by material balance considerations. Old3 illustrated the simultaneous use of Schilthuis' material balance equation and Hurst's fluid flow equation for the determination of the magnitude of reserves and a water drive parameter from pressure and production history. He used this method to calculate the future pressure history of the reservoir under assumed operating conditions. As a basis for determining reserves, Old assumed a value for his water drive parameter and calculated a set of values for the reserves, using the initial reservoir pressure and each successive measured pressure. The sum of the absolute values of the deviations of the resulting reserve numbers from their mean value was taken as a criterion of the closeness of fit to the experimental data possible with the water drive parameter assumed. New values of the water drive parameter were then assumed and new sets of the reserves calculated until a set of reserves numbers having a minimum deviation from the average was established. The average value of- the re-

Jan 1, 1949

By C. R. Watts

The feasibility of producing composites containing nickel or molybdenum fibers in a beryllium matrix was inrestigated. The composties studied were jabricaled by powder mallurgical techniques. The 1-mil-diarr nickel fibers reacled completely below 900°C, converling the fibers .from nickel to Ni5Be2,. As the /LO/-pressing temperalure as raised above 1110oC, tlie nickel diffused outward from the beryllide fibers. The solid solubility of nickel in beryllium was clboul 20 wt pet at the 1100°C pressing temperature a1 the zone-fiber interface. The 1.5-mil-diam molybdenum fibers slzolred no evidence of reaction and little evidence of diffitsion after pressing at 900°C. Between 1000° and 1050°C pressing conditions, the fibers began lo react , producing 1ayers of MoBe2 and MoBe12, respectively surrounding the molybdenurn core. The struture remained the same at 1100°C with no evidenre of solid solubility of the molybdenum in the berylium or vice versa. In recent years a considerable amount of attention has been devoted to the determination of methods for improving the mechanical properties of materials through the use of fiber or whisker reinforcement. Previous work with metal matrix composites indicates that the study of the chemical compatibility of the fiber and matrix is an area requiring greater understanding. The metal-metal or ceramic-metal interface is frequently subject to chemical reactions that may result in the formation of hard brittle intermetal-lic compounds and/or low melting point eutectic compositions. The reaction products may reduce both the low-temperature and elevated-temperature strength of the composite by weakening the fiber-matrix bond, producing premature failure at the interface. It is well-known that most metal-metal systems and many metal-ceramic systems of interest for structural composites are thermodynamically unstable,'-" particularly at elevated temperatures. If, however, the rate of reaction under the conditions of fabrication is sufficiently low. composites can be fabricated that can be used efficiently for indefinite periods at low temperatures and for short periods at elevated temperatures. This paper presents the results of a series of tests to determine the compatibility of nickel and molybdenum fibers with beryllium at various hot-pressing temperatures. Nickel was selected as a candidate fiber material primarily because the relatively ductile fibers might be useful as crack arresters in applications such as ballistic impact where crack growth can result in catastrophic failure. The high density and the reactivity of nickel were primary factors detracting from its selection as a possible reinforcement. Molybdenum with a modulus of elasticity of 52 Xlo6 psi is one of the few metallic materials having a modulus higher than beryllium (42 X lo6). Its high modulus, coupled with its refractory characteristics, made molybdenum an attractive candidate for a relatively stable fiber reinforcement for beryllium. Its density, being higher than that of nickel and over five times that of beryllium: detracted from its other characteristics. EXPERIMENTAL PROCEDURE The specimens were prepared from beryllium powder with a dispersed phase of fibers by powder metallurgical techniques. P-20 grade powder, Table I, from Berylco was used as the matrix material. Short lengths of 0.001-in.-nominal-diam nickel fibers supplied by the Sigmund Cohn Corp. and 0.0015-in.-nominal-diam molvbdenum fibers obtained from the General Electric Co. were used as the dispersed phase. The composite constituents were combined under an argon atmosphere by mechanically mixing the powders and fibers. The compositions used were nominally 1 vol pct fibers. After mixing. the composites were hot-pressed into a-in.-diam pellets under an argon atmosphere at 900°, 1000". 1050". and 1100°C at a pressure of 6000 psi with no hold time at these temperatures so that a comparison could be made between the resultant microstructure and hot-pressing temperature. The billet was heated at a rate on the order of 30°C per min to the desired temperature and then cooled at a somewhat slower rate. The microstructure obtained should be considered as characteristic of the integrated time-temperature history of the sample, as well as the maximum temperature attained. Upon removal from the hot-pressing dies. the specimens were cut. mounted. and polished by standard procedures. No etchant was used in specimen preparation. Photomicrographs, electron microprobe scans, and electron back-scatter pictures were made. X-ray dif-fractometer patterns were made of several of the specimens. but only the lines for beryllium could be resolved. Specimens for optical and electron microprobe examination were selected partially for the roundness of the cross section. A round cross section was taken to indicate that the body of the fiber was approximately normal to the surface and that therefore effects due to fiber material immediately below the surface could be neglected. RESULTS AND DISCUSSION The microprobe scans indicated that nickel reacted as low as 900°C, converting the entire fiber cross section to NisBe21. Fig. l(a). There was no evidence of further reaction from the optical or the back-scatter pictures, Figs. 2(n) and 3(a).

Jan 1, 1970

By John A. E. M. Church

THE prospector's pan was the first implement used for saving gold, and its action is so effective that it has never been equalled for thorough work. Copper plates, blankets, sluices, and amalgamators of several patented designs have thrust it aside, but in spite of all the advances made there is no dispute among experienced men as to the comparative thoroughness of the work done by it and its rivals. When its rivals are doing their best, the pan can take their discarded refuse and prove their inefficiency by producing from it not only gold, but even the mercury which has been lost from the machines themselves. The trouble with the pan has always been that it was only a tool. Man was the machine that moved it. His work was high-priced, the quantity treated small, and it was impossible to use the pan for any but the richest sands. As a consequence the best gold-saving appliance known has been discarded for years and replaced by machines of inferior excellence, but which could more than make up for their limited ability in saving gold by their capacity for treating great quantities of ore. Mr. Peck is one of the old-time amalgamators of the Comstock district, and the first idea of the new machine came to him as he used to pan below the tailrace of his mill. No matter how fine the work of the mill was, the pan invariably proved that better could be done. The tailings always contained gold, silver, and mercury. This state of things is just as true to-day, after all the improvements that have been made, as it was then, and the much-vaunted method of the " Washoe process" consists in extracting 72 per cent. of the gold and silver by one operation, and then repeating this operation half a dozen times to get out about one-third of the remainder. With all the advances of twenty years in milling, there is nothing more crude in metallurgy than the treatment of tailings in many of the large mills of the West.

Jan 1, 1880

By D. Kunstelj, M. Stubicar, A. Tonejc, A. Bonefacic

A. Bonefafcic, D. Kunsfeli, M. Stubicar, and A. Tonejc The present communication is concerned with the variation of the texture in aluminum wires with die angle and temperature, at constant speed of extrusion. Experiments were carried out concurrently on refined samples, with a stated purity of 99.997 pct pure A1 and commercial samples of 99.5 pct Al. Ingots of refined aluminum samples were machined to 30-mm-diam by 150-mm-long billets. These billets were transformed to 5-mm-diam wires by drawing. The final form, suitable for examination, was obtained by extrusion through conical-face dies with a 1-mm hole diam and extrusion ratio 5:l in diam. The initial form of the commercial aluminum samples was drawn wire 5 mm in diam. These samples showed a poorly defined texture with (111) as a major and (001) as a minor component. A similar defined texture appeared in the refined aluminum samples after drawing to 5 mm diam. Conical-face dies with different angles (defined by the axis and the generating line of the cone) were used in our experiments. The values of the angles were 27, 35, 45, 57, 63, and 90 deg. The extrusion container was fitted with a heating element and controller permitting temperatures up to 600°C to be maintained within i5"C. Extrusion was performed at 250°, 300°, 350°, 400°, and 500°C at constant speed (approximately 1 mm per sec) and constant die reduction. The extrusion product was a wire 1 mm in diam and approximately 20 cm long. In order to remove the surface layer with the "conical" texture and to reduce the absorption by the X-ray examination of the samples, the extruded wire was etched to 0.22 mm in diam. Experiments were performed in the middle sections of the 20-cm-long wires. In addition to the die of l-mm hole diam, dies with a 1.5-, 0.7-, 0.6-, 0.5-, and 0.4-mm hole diam and 63-deg die angle were constructed. In our experiments we did not find in these ranges any important difference concerning the texture of the extruded wires and we continued our work solely on the 1.0-mm die. The diffracted X-rays (Cu K radiation) were recorded photographically. Diffracted intensities were measured on the (111) reflection with a microphotometer. The relative amounts of texture components were determined from the areas under the diffracted maxima. We found the texture of extruded aluminum wires to be strongly influenced not only by the temperature of extrusion and the purity of the sample but also by the form of the die. It is generally admitted that cold-drawn aluminum wires have mainly a (111) texture with a small amount of (001) component, Table I of Ref. 1. In our experiments with wires extruded in conditions represented by Fig. 1, in some cases a single (001) texture was obtained. If these wires were drawn repeatedly at room temperature, X-ray measurements revealed a duplex (001)-(111) fiber texture. Further drawings increased the (111) and decreased the (001) texture component. In Fig. 1 the percentage of material oriented with (001) parallel to the extrusion direction is represented as a function of the temperature and the die angle (a), for commercial and refined aluminum samples, respectively. From these diagrams we may draw the following conclusions. The slope of the die (a) influenced more strongly the texture at the lower rather than the higher temperatures. Again, a stronger influence was found in the case of the commercial in comparison with refined aluminum samples. In the case of the commercial aluminum samples the amount of material with (001) texture increases with increasing wire temperature in an approximately linear manner. This effect is less pronounced in the pure aluminum samples, with the exception of the die with a = 45 deg. In this case the (001) texture decreases with increasing temperature, as shown in Fig. 1. Component (001) is more pronounced in higher-purity aluminum samples. Our experiments led to the conclusion that both (001) and (111) components are essentially stable in extruded aluminum wires. As we obtained a single (001) texture starting with a sample of drawn wire in which the (001) component was very weak, our experiments revealed that the (001) component is not a remnant of the initial texture; this is in disagreement with the findings of Vandermeer and McHargue.1 We gratefully acknowledge discussions with Professor M. Paic. 1 R. A. Vandermeer and C. J. McHargue: Trans. 7MS-AME, 1964, vol. 230, p. 667.

Jan 1, 1969

By W. B. Gogarty

The reported decreases in water mobility do not seem unusual in view of non-Newtonian fluid properties. Shear stress vs shear rate diagrams have been reported for other solutions of water-soluble polymers. Some of these polymers are similar to the type mentioned by the author. Generally, the shear stress-shear rate is a non-linear function for these solutions. Data for plotting apparent viscosity vs shear rate can be obtained from this function. Apparent viscosity is defined as the ratio of shear stress to shear rate at a given shear rate. When plotted, the apparent viscosity decreases with increasing shear rate. This behavior is typical of a pseudoplastic fluid. For some water-soluble polymer solutions, the apparent viscosity decreases more than 50 times while the shear rate increases 1,000 times. Thus, viscosity of a pseudoplastic fluid only has meaning at a specified shear rate. Results of Fig. 1 could be explained in these terms. Viscosities measured in the Ostwald viscometer represent values at a given shear rate. Some average shear rate is affecting the polymer solutions while flowing through the core. This average value fixes the apparent viscosity as long as the flow rate remains constant. Viscosities measured by the two methods will be equal if shear rates are the same. The results indicate that shear rate in the core is lower (higher apparent viscosity) than in the viscometer. In the paper by Johnson, Bossler and Naumann, the relative permeability is independent of viscosity ratio. Thus, the relative permeability with respect to water flow at residual oil should be independent of the flowing phase viscosity. Polymer solutions will appear as Newtonian fluids The discussion emphasizes the nature of the "resistance factor effect" as discussed in the paper. Repeated anomalies arising in hundreds of experiments led us to the conclusion that non-Newtonian flow is not the only factor. Several of the key anomalies are as follows: 1. Measured viscosities over a range of shear rates from <1 sec-&apos; to 1,000 sec-&apos; do not account for but a minor fraction of the R observed in cores when compared in similar shear-rate ranges. 2. The slope of R vs flow rates in cores is always different from that expected from viscometer shear-rate measurements as shown in Fig. 2. in a core, the level of viscosity being fixed at a given flow rate. With these conditions, the definition of resistance factor R by Eq. 2 is simplified to Since , is constant with rate, R becomes a measure of the apparent viscosity in a core at a given flow rate. Variation in flow rate could easily account for the changes of R shown in Fig. 5. Also, this points to the fallacy of assuming R to be a unique parameter. The constant resistance factors at different flooding velocities appear to be in disagreement with the above discussion. The author furnishes Fig. 2 to support his arguments. As shown, the resistance factors remained substantially constant in the two cores over a considerable range of flooding velocities. However, in the 73-md core, the factor increases at lower rates. This behavior agrees with known characteristics of some pseudoplastic material. These materials act both as Newtonian and as non-Newtonian fluids in different regions of shear rate. Some exhibit first Newtonian, then non-Newtonian, finally, Newtonian character. Others are first non-Newtonian and then Newtonian. This latter type would explain the results with the 73-md core. The Fann-instrument results are not significant since shear rates in the core may be much different than with the viscometer. The higher resistance factor at high rates in the 150-md core is more difficult to explain. The greater resistance at increased flow rates could be attributed to what might be termed temporary bridging. As envisioned, changes in polymer configuration occur at the higher energy associated with the increased flow rate. These changes could cause less effective passage of polymer through the core. Correspondingly, increases in pressure drop will occur. These will be interpreted as higher resistance factors. 3. Most polymer solutions are non-Newtonian and many are more shear-rate sensitive than the polymers in question, yet only a very few polymers demonstrate useful R values. Gogarty&apos;s assumption that viscosities in cores and vis-cometers will be the same if measured at the same shear rate is only valid if non-Newtonian rheology is the only parameter. The experimental evidence does not validate this assumption. The anomalies observed in the equilibrium displacement experiment shown in Fig. 5 are not explained on the basis of varying flow rates since the rates were held constant. M

Jan 1, 1965

By D. A. Augenstein, L. V. Amundson

Operators of coal preparation plants use equipment such as large-capacity dump trucks and bulldozers to haul, spread, and compact refuse material to conform to federal and state regulations governing the disposal of solid waste. For example, federal regulations require that refuse be spread and compacted in layers not more than 0.6 m (2 ft) thick. If the refuse layers or piles have poor load-bearing characteristics, movement of equipment over them becomes extremely difficult. The increasing cost and limited space for the construction of settling ponds extensively used for fine refuse disposal are leading to the growing practice of dewatering the fine refuse by means of vacuum disk filters and combining the resulting filter cake with coarse refuse for waste bank disposal. However, this filter cake usually contains more than 25% moisture. Therefore, when it is combined with coarse refuse, the moisture level of the total refuse may become high enough to impair the load-bearing strength of the refuse. Precisely this problem and the attendant difficulties of moving heavy equipment arose when vacuum disk filters were installed at a preparation plant in West Virginia. Combining the moist filter cake with the coarse refuse generally results in a 12% moisture level for the total refuse at this plant. At this moisture level, the refuse can still be handled with some dif¬ficulty by the equipment. However, the problem has frequently been compounded after rains. The experimental program described in this paper tested the following methods of improving the load-bearing properties of coal refuse: moisture reduction, addition of crushed coarse refuse, addition of fly ash, addition of lime, and addition of a mixture of lime and fly ash. In terms of a balance between economic and technical considerations, the most effective method was demonstrated by laboratory and plant tests to be a 2% to 5% addition of lime. Test Program During formulation of a test program aimed at improving coal refuse stability, a number of treatment techniques were selected for evaluation of their effect on refuse load-bearing characteristics. The reduction of refuse moisture was the first technique considered, since moisture content plays a key role in the bearing strength of a bulk solid. The reasoning was that a very small reduction in refuse moisture, provided it could be accomplished with minimal effort and cost, would be sufficient for eliminating refuse disposal problems most of the time. At the same time, it was recognized that periods of rainy weather would offset moisture-reduction measures. Other techniques selected for evaluation because of their potential for improving refuse-bearing strength at reasonable cost were: addition of crushed coarse refuse, the addition of fly ash, and the addition of lime. A simple one-dimensional laboratory compaction test was designed for evaluating the effect of each of these techniques on refuse-bearing strength. Those techniques that gave best results in the laboratory would then be tested in the preparation plant in West Virginia. Laboratory Tests The laboratory test for evaluating each of the techniques involved the application of weight to a load module placed on a refuse sample in a container. The distance the load module sank with an increasing amount of applied weight was measured, and the relationship between module displacement and loading was obtained. The refuse container was large enough to minimize the influence of wall effects, and to decrease the effect of the container bottom, a test was terminated when the load module sank to about three fourths the depth of the container. Initially, an attempt was made to conduct compaction experiments with samples of total plant refuse. However, the presence of + 0.13 m (+ 5 in.) material in the refuse interfered with the mechanics of the test, and the large amount of material required for a representative sample made testing extremely tedious and time-consuming. To avoid this problem and because it was believed that only the fine components of the refuse were the major cause of the stability problem, tests were conducted with 13 mm X 0 (1/2 in. X 0) material from the total plant refuse. The equipment consisted of a 0.46 m (18 in.) cubic container, a loading module with 62 cm2 loading area, weights totaling 450 kg, and a portable concrete mixer for combining water and additives with the refuse. Refuse samples of about 110 kg dry weight were used for each test. Experimental compaction test data were evaluated by plotting the linear displacement of the load module vs. loading. During compaction tests, the vertical displacement of the load module was measured at each loading value. The plot of

Jan 1, 1980

By Sandford S. Cole, John S. Breitenstein

The recovery of over 80 pct of the vanadium values in titaniferous magnetite from Maclntyre Development,Tahawus, N. Y., was accomplished by an oxidizing roast with Na2O3-NaCI addition. Process description is given for leaching of roasted ore and precipitation of V2O5 and Cr2O8 from leach liquor. THE exploration and development of the Mac-Intyre orebody at Tahawus, N. Y., by the National Lead Co. provided a source of vanadium. Analyses of various composite sections of the drill cores of the MacIntyre orebody were made to establish whether or not the vanadium was constant throughout. Ten drill cores were sampled as 50 ft sections, crushed, and a portion magnetically concentrated. The head and concentrate were analyzed for total iron and vanadium. The results on the concentrates indicated that the vanadium is associated with the magnetite and maintains a close ratio to the iron content. The nominal ratio of 1:25:140 of V: TiO2:Fe was found to exist in the concentrates. Typical value for the vanadium in the magnetite both from laboratory concentration and mill production is 0.4 pct. The recovery of vanadium from the magnetite was investigated in 1942 to 1943. The research program encompassed both laboratory and pilot-plant work on sufficient scale to provide adequate data to establish the feasibility of a full scale plant. The recovery of vanadium from various ores has been reported in the literature and has been the subject of many patents. The literature dealing with recovery from titaniferous ore by roasting is quite limited. Roasting with alkaline sodium chloride, sodium chloride or alkaline earth chlorides, and sodium acid sulphate have been claimed in various processes as effective means.1-8 The reduction of the ore, followed by acid leaching, was another method proposed.&apos;-&apos; "he use of various pyrometallurgical processes for recovery of vanadium in the metal or in the slag has also been extensively investigated, but the results had little application to the problem."-" The separation of vanadium values from subsequent leach liquors and vanadium-bearing solution has been the subject of a considerable number of papers and patents. The most practical is by hydrolysis at a pH of 2 to 3 by acidifying a slightly alkaline solution. Data on solubility of V²O5 and V2O4 in water and in dilute sulphuric acid indicated a solubility of 10 g per liter in water.&apos;" Laboratory Results Magnetite Analysis: Adequate stock of magnetite was provided so that the laboratory and pilot-plant operation was on ore representative of the mill production. The ore was analyzed chemically and examined by petrographic methods to ascertain whether the vanadium was present in combined state or as an interstitial component between grain boundaries. No evidence was obtained which would indicate that the vanadium was in a free state as coulsonite.15 The analysis of the ore was as follows: Fe²O³, 47.4 pct; FeO, 29.1; TiO,, 10.1; V, 0.40; and Cr, 0.2. The screen analysis of the ore on the as-received basis was: -20 +30 mesh, 28.8 pct; —30 +40, 18.9; -40 +50, 9.7; -50 +60, 15.1; -60 4-100, 5.9; -100 + 200, 11.2; -200 +325, 3.7; and -325, 7.2. Roasting Conditions: The prior practice indicated that a chloridizing roast with or without an alkaline salt had been effective on other titaniferous magnetites. On this basis roasts with additions of sodium chloride, sodium carbonate and mixtures thereof were investigated varying the roasting temperature between 800" and 1100°C. Since the ore had shown no segregation or concentration of vanadium, the influence of particle size on the freeing of vanadium by the reagents during roasting was determined. The initial work was on silica trays in an electric resistance furnace with occasional rabbling of the charge. Subsequently, the roasting was carried out in a small Herreshoff furnace to establish the influence of products of combustion on the recovery of the vanadium. The laboratory tests showed that this ore required an alkaline chloridizing roast, in conjunction with a reduction in particle size to less than 200 mesh. When roasted in air at 900 °C with 5 pct NaCl and 10 pct Na2CO³, over 80 pct recovery of the vanadium was attained as a water-soluble salt. The presence of alkaline earth elements gave detrimental effects and care had to be exercised to avoid any contamination of the ore or roast product by such materials. The solubilization of vanadium under the various conditions is given in a series of curves in Figs. 1 to

Jan 1, 1952

By David Turnbull

IN isothermal recrystallization processes, new crystals generally grow into the matrix until they impinge upon other new crystals or an external surface, at constant linear rates G. Before impingement the perceptible course of growth can be described by the equation: 1 = G(t-7) C1I where, G = dl/dt, 1 is a crystal dimension measured in a constant direction, t is the time, and 7, the nucleation period, is a positive intercept on the time axis. Fig. 1 is a schematic representation of I as a function of time for a recrystallizing grain. G is dependent upon temperature, driving energy (strain or surface energy), relative grain and boundary orientations, but is generally independent of time. The frequency of nucleation, fi, (time" volume") can be defined by the equation: N = 1/fV [2] where ? is the mean nucleation period and V is the volume of the specimen that has not recrystallized. The kinetics of primary and secondary recrystallization generally can be described satisfactorily in terms of the parameters N and G only.&apos;-" After recrystallization is complete the average grain size 7 increases with time by "normal grain growth;" didt, the average rate of grain growth, is strongly time dependent and has not yet been precisely related to G for the motion of the individual grain boundaries constituting the system. It has been suggested4* " that the elementary act in grain boundary migration is closely related to the elementary act in grain boundary self-diffusion. Although the distance of atom movement in the two processes may be somewhat different, there is reason to expect that the activated states may be very similar, so that the free energy of activation for grain boundary migration should be of the same order of magnitude as for grain boundary self-diffusion. Therefore, it is desirable to develop a satisfactory basis for comparing data on self-diffusion and grain boundary migration and to make such comparisons where possible. Theory The formal theory of grain boundary migration rates is analogous to the theory for the rate of growth of crystals into supercooled liquids reviewed elsewhere 6-8. Boreliuss has shown that the latter theory describes, within the theoretical uncertainty, the growth of selenium crystals into supercooled liquid selenium. Motto and more recently Smolu-chowski" have derived expressions for the rate of boundary migration in recrystallization. The treatment to be presented is similar to Mott&apos;s excepting that the formalism of the absolute reaction rate theory will be used. The atomic mobility, M, in grain boundary migration is defined by: G = -M6p/6x where p is the chemical potential per atom and x is the coordinate measured in the direction of grain boundary movement. Let AF be the free energy difference per gram atom on the two sides of the boundary and k the thickness of the boundary. For RT>>AF the potential gradient across the boundary (6p/6x) is essentially linear and it follows that: SF/8x = - aF/N\ [4] where N is Avogadro&apos;s number. According to the Nernst-Einstein equation, M is related to a diffusion coefficient, Do, for matter transport in grain boundary migration by the equation: M = Da/kT [5] Substituting eqs 4 and 5 into eq 3 gives the basic relation between Do and G: G = DoaF/\RT [6] Do values may be calculated from experimental values of G from eq 6 and directly compared with the coefficient of self-diffusion within the crystal, DL, or the grain boundary self-diffusion coefficient D,. However, a more convenient, though equivalent, basis for comparing atomic mobility in grain boundary migration and self-diffusion is through the constants of the absolute reaction rate theory. According to this theory diffusion coefficients may be written:" D = k2(kT/h) exp [-AF,/RT] 171 aFa, the free energy of activation, is related to the measured energy of activation, Q, by the equation: AFA = Q - T aSx - RT [8] where aSa is the entropy of activation. Substituting eqs 8 and 7 into eq 6 gives: G = ek(kT/h) (aF/RT) exp [(AS,,)C/R] exp C-Qc/RTI C91 where the subscript G refers to boundary migration. The relationship between the driving free energy and the free energy of activation in boundary migration is indicated schematically in Fig. 2. Experience indicates that the variation of G with temperature can be described by: G= Go exp [- Qc/RT] [10] where Go and Qc are generally temperature independent over wide ranges of temperature. Comparison of eq 9 with eq 10 gives:

Jan 1, 1952

By J. P. Timmons, J. H. Moran, G. K. Miller, M. A. Coufleau

A prototype equipment has been designed and built for the digital recording of well logs on magnetic tape at the same time that the regular film recording is made. The format of the digital tape produced is such that it can be used directly at the input of the ZBM 704, 7090 or other models of ZBM computers which accept digital magnetic tape. This apparatus has been used for the experimental field recording of dipmeter tape logs which were subsequently computed by means of an ZBM 704 or 7090. In this paper the equipment and the digital tape are described briefly, and their application to the computer-interpretation of dipmeter data is discussed. A principal element in the interpretation of the dipmeter log is the correlation of the three microresirtivity dipmeter curves to determine the depth displacements between them. Several correlation methods for computer use are considered, with particular attention to their sensitivity to error and their consumption of computer time. The tape data were used to compute information content of the dipmeter microresistivity curves in terms of their frequency spectra. The results show that the sampling rate used in recording the digital information is quite adequate and illustrate a use of the digital tape in evaluating the characteristics of new tools. Some examples of field results are shown. It can be foreseen that, when digital tape recording becomes available for general field use, a whole new realm of possibilities will be opened up for the processing of other well logs through computations, which hitherto were not feasible because they were too laborious and time-con.sunzing. INTRODUCTION The last few years have seen a revolution in the design and production of data-processing equipment. Stored-pro-gram digital computers have progressed from a research curiosity to the basis of a major industry. There are now hundreds of such machines in daily use in the United States. With the acceptance of a technique that was, in fact, already clearly described by John von Neumann in 1945, the last decade has seen great strides in the development&apos;of components, reliability, programming systems and, most spectacularly, in the sheer number of machines built and in use. In 1957 there were enough digital computers available to the oil industry to justify the suggestion that it would be worthwhile to investigate the possibility of using these machines in processing well log data.&apos; The first result of this investigation was the appearance of what may be referred to as the input-output bottleneck. Well logs are customarily recorded on film. To get these data into a machine required then (and still does): a time-consuming semi-automatic reading of the film; conversion of the log data to digital form; and recording these digital data in some medium acceptable for computer input, such as cards, magnetic tape, or punched paper tape. However, the recording, reading, and re-recording could only result in deterioration of the data. Therefore, it was concluded that the fist step should be the development of a new, more direct recording technique supplemental to the film recording, which would provide easy access to the digital computer. There are many solutions to the problem of recording log data in an easily recoverable form. After careful consideration it was decided to adopt the boldest solution which, it was felt, was also the most elegant. It was decided to record well logs directly, in the field, on magnetic tape in such a way that this tape could be used without further modification as an input to the IBM 704 or 7090 computer. To realize practical field recording of magnetic tape logs, it became necessary to develop in a rather small package, an analog-to-digital converter, a tape recorder, and the necessary multiplexing and control circuits to allow the simultaneous recording of a multiplicity of logging signals. The magnetic tape recording was to be made simultaneously with the conventional logging operation in such a way as not to interfere with it. Along with the development of hardware, it was necessary to begin development of interpretation techniques and machine programs that would exploit the power of the digital computer. Here, again, there is a long list of possible applications. After much consideration it was decided to concentrate on the interpretation of the dipmeter log as a first application. It is the object of this paper to describe in some detail the developments sketched in the last three paragraphs.

By John C. Martin

The discovery of a new gas reservoir demands that the planning of a sottnd well-spacing program be initiated early in the development stage. It is the purpose of this discussion to illustrate by actual field examples the application of basic well-spacing principles, previously developed for oil reservoirs, to the problem of well spacing in natural gas fields. These studies are presented for the field use of geologists and engineers who are concerned with the initial planning of the proper development of the newly discovered gas reservoir. INTRODUCTION The phenomenal growth of a vigorous natural gas industry emphasizes the increasing importance of natural gas as a source of energy, fuel, and raw materials to our nation&apos;s economy. Since 1945 marketed production of natural gas for the U. S. has increased 21/2 times to a record high of 10.6 trillion cu ft during 1957. As major participants in the gas industry, we share an added interest to develop and produce our natural gas reserves with constantly improved efficiency. The subject of well spacing is vitally important to the gas industry, for the well itself plays a significant role in the development of the natural gas reservoir and in control of the recovery process. Maximum utilization of wells is an integral part of sound conservation practices. The discovery of a new gas reservoir demands that careful choice of well location and well spacing be made and that the planning of a sound well spacing program be initiated early an the development stage. With the drilling of the first development well, efforts of the geologist and engineer must be directed toward acquisition of adequate technical evidence upon which a firm recommendation for a spacing program may be based. With this technical appraisal as a foundation, operators and state regulatory agencies jointly can go far in providing a framework for sound development of gas fields to achieve a program of conservation that avoids the unnecessary well. WELL-SPACING CONCEPTS Through laboratory and field investigations of the mechanism of the recovery of oil and gas, of fluid behavior, and of effective control of reservoir and well, a crystallization of ideas regarding reservoir behavior has emerged as a well-developed technology. Associated with a better understanding of the fundamental principles underlying reservoir and well behavior has been the growth of concepts concerning the role of wells and their spacing in the development and operation of an oil or gas reservoir. In addition to serving as outlets for the withdrawal of fluids from the reservoir, wells are recognized as having two other important functions: (1) providing access to the reservoir to obtain information concerning the characteristics of the reservoir and its fluids, and (2) serving as a means by which the natural or induced recovery mechanism may be effectively controlled. Beyond a minimum number of wells required to fulfill these two functions, additional wells will not increase recovery. With particular emphasis upon well spacing in oil reservoirs, many studies of the well spacing-recovery relationship have evolved the concept that the ultimate oil recovery is essentially independent of the well spacing.&apos; These fundamental concepts are no different when regarding the role of wells and their spacing in the natural gas reservoir. They are equally applicable to the consideration of well spacing in gas reservoirs. For the gas reservoir, the problem of well spacing then revolves around the question of drainage and the degree or extent to which a well may drain gas from its surrounding reservoir environment. Theoretical and Experimental Work During the past 30 years, theoretical and experimental work carried on to study the physical principles involved in the flow of fluids through porous media has shed light upon the matter of drainage. Fundamental mathematical equations have been derived to describe the mechanism of flow of oil and gas through porous rocks. With the recent advent of high-speed digital computers, attempts have been made with mounting success to develop solutions, employing numerical techniques, to mathematical expressions that describe more rigorously the physical behavior and mechanism involved in the unsteady-state flow of compressible fluids, such as a gas, through porous rock. In 1953, Bruce, Peace-man, Rachford and Ricez published a stable numerical procedure for solving the equation for production of gas at constant rate. The results of these calculations are significant with respect to this matter of drainage, for they indicated (1) that depletion of the gas reservoir resulted in a drop in pressure at the extremity of the

By V. F. Leferrer

A new process for the treatment of silver crusts from desilverization of lead has been developed by Societe de Penarroya. The zinc contained in the crusts is recuperated by vacuum distillation at low temperature. The use of special retort crucibles is avoided. AT Noyelles-Godault the Parkes process (zinc addition to a silver-bearing lead melt) is used to desilverize lead bullion from the blast furnace. The crust which rises to the surface of the bath is removed and pressed by a well-known process in a press of the Howard type. The pressed crust, called zinc crust or silver crust, contains 10 pct Ag, 30 pct Zn, and 60 pct Pb. It is melted under a salt bed in a small, narrow, and deep kettle heated only in its upper part. A liquid concentrated alloy is obtained which is called T.A.C. in France (triple alliage concentr6) analyzing approximately 25 pct Ag, 65 pct Zn, 8 to 10 pct Pb, and 1 to 1.5 pct Cu. This T.A.C. is a real metallic alloy, completely deoxidized, and can be stored indefinitely. The process used at Noyelles-Godault consists in distilling this alloy in a tight furnace, under low vacuum and at low temperature, so that only the zinc distills and is condensed into liquid form. The residual silver-bearing lead is almost entirely free from zinc and oxidized dross. The vacuum dezincing equipment, see Figs. 1 and 2, consists of: 1) A furnace, constituted of a casing of steel plate lined with refractory bricks, heated internally by means of graphite resistors. The temperature of the furnace is measured by a thermocouple device and controlled by an automatic temperature regulator which acts upon the flow of the current passing through the resistors. Two doors, numbered 1 and 2 in Fig. 2, are provided for, one for charging the T.A.C., the other for tapping the high silver-bearing lead. 2) A condenser, and an ordinary room lined with good heat conductor refractory material. This condenser is provided with a graphite resistor which serves to bring to temperature when starting the furnace. A thermocouple records the temperature of the zinc melt and controls the temperature in the furnace, that is, the flow of current through the resistors, by means of the above mentioned regulator, so that the distillation rate in the furnace does not exceed that of the condensation. The condenser is connected to the vacuum pump by a filter. The whole equipment should, of course, be fluid-tight. The operation is carried out under a vacuum of about 10 mm Hg.

Jan 1, 1958

By C. R. Breck

THERE has been a running discussion over the past several years with respect to the life and adequacy of our natural gas reserves. Some of the experts agree on one phase of the subject at least-that eventually it will be necessary to augment deliveries of natural gas in the northern markets by the development of a more or less comparable synthetic gas made from coal. The purpose of this paper is to point out: 1. It may be later than a good many people think-the demand for addi¬tional domestic and commercial pipeline gas in at least one section of the North will probably be felt rather keenly before an adequate synthetic pipe¬line-gas project, developed on economic lines, could be completed, even if initiated in the near future. 2. The evidence is considerable that such a project would be technically sound and economically self-supporting, although the margin of profit at first might be so small and the magnitude of an economic enterprise so great -that financing would be impossible without Government aid in some form or other. The part of the country referred to is the heavily populated section in the northeast of the United States. To simplify an examination of the situa¬tion, let us take the potential marketing areas within range to the north and east of the three pipelines: 1. Tennessee Gas Transmission Co. as it enters Erie County, New York state. 2. Texas Eastern Transmission Corporation just before the tap that sup¬plies the general Philadelphia area. 3. Transcontinental Gas Pipe Line Corporation as it enters the state of New Jersey. The sources of gas for these three pipelines are in Texas and Louisiana. The ultimate maximum combined capacity of the three lines at the points stated is approximately 1,200,000,000 cu ft per 24 hr at the maximum load factor probable of attainment in normal operation. The areas supplied by the three lines north and east of the points stated, or assumed as eventual markets, are: 1. Central New York state, from Buffalo to Albany (both inclusive).

Jan 1, 1953