Search Documents

By P. J. Jacquard, C. Jain

Interpreting pressure measurements is one of the fundamental problems in the operation of oil reservoirs. The principal methods and means of study zuhich are currently used are reviewed. Emphasis is placed on the necessity of improving the quality of the information to turn it to the best possible account. One of the important aspects of the problem is the development of a system of automatic interpretation of the pressure measurements. The proposed method aims at this objective: starting from the relation which exists between a local perturbation of permeability and its effect on the velocity potentials in the wells, we determine by successive approximation the modifications to be made to the permeability map by fitting the least squares between the measured pressures and the pressures indicated by the mathematical model. The method has already been successfully employed in the case of heterogeneous reservoirs of radial- circular symmetry. Its application is extended here to the two-dimensional case, and is illustrated by two examples, for which the results are presented. The conditions required for the method to become fully operational are likewise analyzed. I — INTRODUCTION One of the basic problems of production and recovery from oil fields is interpreting pressure measurements made in wells in the single-fluid-phase regime. Simplifying assumptions, for example supposing the reservoir to be homogeneous, make it possible to use more or less simple methods which permit the estimation, from these pressure measurements, of a mean value of the permeability around the wells. Considering the actual heterogeneity of the majority of reservoirs, and the phenomena of interference between wells, these methods are only indicative. Much has been published on the studies of different cases of comparatively simple heterogeneity: some have discussed the radial heterogeneity around a single wel1,l or vertical heterogeneity.2 The results of these works have shown the influence of variations of permeability on the pressure, but unfortunately they are not very useful in practice because they all are concerned with the phenomena around a single well, without interference between wells, and they correspond to relatively simple and theoretical schemes. These studies of reservoirs or parts of reservoirs on a larger scale are made by the electric analyzer or digital model. We define a permeability map and a porosity map, then reproduce the production history of the wells on the wells of the model, and note the corresponding pressures. These pressures are then compared with actuality, and the discrepancies observed are analyzed. From this analysis we deduce the modifications to be made to the maps of permeability and porosity, and repeat until the coincidence between the pressures indicated by the model and the actually observed pressures is considered adequate. We have adjusted the model. Obviously, in this method of investigation, everything depends on the formula which permits the modification of the permeability to approach the adjustment. Considering the complicated nature of the problem, some oversimplified formulas are poorly adapted to reaching a good adjustment. For example, we have shown that the elementary operation of increasing the permeability not far from a producing well to compensate for too low pressures in the well can, in certain cases, express itself by an effect exactly opposite. 3 The project described here is for a system where we begin with the pressure measurements and the several known parameters of the reservoir, and develop the most probable permeability map, taking account of the measurements available and a subdivision of the reservoir fixed a priori. Ideally, this system is conceived as one in which no manual intervention would be necessary. The model simulating the reservoir is only one part of the whole system, and it is neither the longest nor the most complex part. The need to integrate the model into an aggregate of much more extensive processes has led to the adoption of the digital solution, even though a solution combining the digital and the analog is perfectly conceivable. The simulating model is in that case no more

Jan 1, 1966

By Gilbert Santoro

In this study tantalum carbide filaments of various compositions in the fcc region were prepared by heating a tantalum wire in a measured amount of hydrocarbon vapor. Such properties as tensile strength, transverse-rupture strength, microhard-ness, lattice parameter, electrical resistivity, and magnetic susceptibility were measured as a function of carbon content. Some of these properties showed a maximum or a minimum between a carbon-to-tantalum molar ratio of 0.80 and 0.85, the composition range at which the color changes from metallic gray to gold. The recent literature has shown that in the large homogeneity range for the fcc phase of the Ta-C system there is a marked variation for such properties as microhardness,' electrical resistivity, ' thermodynamic properties,z'3 lattice parameter and superconductivity. In addition to the technological interest of these variations within a single phase, important theoretical benefits may ensue from a continuous investigation of properties vs stoichiometry in tantalum carbide and other refractory defect compounds, for example the effect of the metal-metal distance and crystal structure upon bond formation and electronic structure. Whether there is a similar variation of the tensile strength is of practical importance. Of theoretical interest would be the knowledge of how the magnetic susceptibility varies with carbon composition. Both of these properties, as well as the transverse-rupture strength, were determined for tantalum carbide filaments in the fcc phase region. The lattice parameter, the electrical resistivity, and the microhardness were also investigated and compared with values appearing in the literature. EXPERIMENTAL PROCEDURE Material Preparation. Tantalum carbide specimens of various carbon content were prepared by statically carburizing high-purity tantalum filaments in measured quantities of pure hydrocarbon gas. The filaments, which were 10 mil in diameter and 20 or 40 cm long, were heated by their own resistance at temperatures from 1800" to 2200°C using alternating current in some cases and direct current in others. Tantalum wire from two sources was used. One had a reported purity of 99.96 wt pct Ta. The ingot from which this wire was drawn had been electron beam-melted as the final step in its purification. The reported purity of the other wire, which was prepared by a powder metallurgical process, was 99.88 wt pct Ta. The impurities of both tantalum sources are listed in Table I. Spec-trographic-grade toluene or research-grade propane was used as the source of carbon. The apparatus for the carburization of the filaments is shown schematically in Fig. 1. The system was evacuated to 1 x l0-' mm of Hg when using the electron beam-melted tantalum and to less than 1 p when using the less-pure tantalum. Neither the difference in the base vacuum, nor in the hydrocarbon gas, nor in the tantalum filaments resulted in any detectable difference in the properties studied in this work. However, only the carbides prepared from the higher-purity tantalum were used in the magnetic-susceptibility determinations. Once the base vacuum had been reached, the system was flushed with the hydrocarbon and re-evacuated. For any desired composition of the carbide, hydrocarbon gas was allowed to leak into a constant-volume flask to a predetermined pressure. These pressures were measured by a Cartesian diver-type mercury gage with a range from 0 to 20 mm and with a scale expanded to nine times. The hydrocarbon gas was allowed into the reaction tube and the filament was quickly raised to the reaction temperature. The slack in the filament due to thermal expansion was taken up by drawing out the movable electrode with a screw mechanism. The filaments were kept at temperature from a few minutes to 200 hr. In the majority of cases the heating time was 6 hr. The reaction tube was cooled by cascading water on its outer surface. By means of an ammeter and a voltmeter the electrical resistance at the reaction

Jan 1, 1963

By E. H. Lee, F. J. Fayers

Problems in reservoir analysis can usually he cupressed in terms of a system of nonlinear partial difjerential equations. A rnethod for .setting physically reasonable boundary conditions for these systerrls by ayplication of the theory of characteristics is discussed. A short description i.c. given of the technique for finding chnmcteristics of general first order and second orrler equations. Specific examples are quoted jrorn eqrrations occurring in theoretical reservoir studies. The princip1e. ossocicred wit11 the use of chnracteristic to set physically realistic boundary conditions me developed, and thr method is discussed in terms of boundary condition for the noncupillary, compressible, two-phase flow problen: jor the cupillary, ir~cowlpressible, two-phase flow problems for the tlrtergent flooding prohlenl; as for the compressible. noncapillary, two-phase flow pro11letn. INTRODUCTION In most problcms in reservoir analysis the physical situation is described by a system of nonlinear partial differential equations. The general solution to the set of equations usually comprises a wide range of functions ancl a particular solution is normally selectetl by applying suitable boundary conditions to the problem. In some cases it is clear from the physical situation what form the boundary conditions should take, but it is usually desirable to have an independent mathematical check of the suitability of the chosen conditions. From the mathematical standpoint, a system of partial differential equations with boundary conditions is said to be properly set if the combination determines a unique solution which is continuously dependent on the equation and the boundary conditions. Such a solution is sometimes termed "reasonable". In investigating the properties of the mathematical representation, it is useful to determine the structure of the system by the application of the theory of characteristics. This theory determines in the space of the independent variables, lines or surfaces which play a special role in thc system's behavior. The characteristics determine the possibility and nature of the propagation of discontinuities in the derivatives of the dependent variables. Along the characteristics certain differential relations hold which for some problems allow a so- lution to be obtained analytically. In the case of hyperbolic equations the characteristics can bc used as 3 basis for a numerical solution. If a complicated nonlinear system is to be solved, a study of behavior of the characteristics will indicate what boundary conditions should be chosen. By using the method of characteristics, several advances have rccently been made in solving reservoir engineering problems. Sheldon, Zondek and Cardwell' have discussed how the the method may be used grap-phically to solve the incompressible two-phase displaccmerit problem when the gravity term is included in the equations. Fayers and Perrine' have applied the method to study the equations representing detergent flooding. When compressibility and solution-gas effects are introduced into the two-phase flow problem, the equations and their characteristics become more complicated. It will be seen that in this problem the characteristics are useful for indicating a reasonable set of boundary conditions but cannot readily be used to compute a numerical solution. This problem has been solved on the IBM 701 using normal finite difference methods hy West, Garvin and Sheldon." The material to be prcscntcd in this report is divided into three sections. The first discusses the determination of the characteristics of first and second order partial differential equations; the second describes the selection of suitable boundary conditions by the application of the principle that the future cannot affect the past; the third discusses the use of characteristics for setting the boundary conditions of the two-phase solution-gas flow problem. A useful introduction to the method of characteristics has been given by Hildebrand.' A summary of the method will be given in this section of the paper. The First Order Equation The quasi-linear equation of first order is It is linear in the partial derivations of z. The properly linear equation has the form The characteristics of are curves along which the

By Juanito C. Fernandez

INTRODUCTION The Philippines is endowed with mineral resources, much of which still remains unexploited. The country's known reserves, at least for the major mineral commodities, have increased every decade since the first major discoveries more than 80 years ago. Mining, therefore, naturally evolved as one of the pioneer industries and has ever since been a major factor in economic development. The Philippines has entered a new era in mineral resources development and management. One significant feature of this era is the degree to which the government is involved with the industry. In recent years, the state has increasingly encouraged the private sector to carry out detailed exploration aided by governmental contributions through extensive geological and geophysical surveys, fundamental researches and even indirect financial assistance in some specific areas. Significant changes in mineral resources development policies were brought about with the promulgation of Presidential Decree 463 or the Mineral Resources Development Decree of 1974. The Decree provides for a more rational system of administration and disposition of mineral lands to promote mineral resource exploration and exploitation. Various forms of incentives are provided. Foreign entities are allowed to participate in development through service contracts with leaseholders or with the government in the case of reserved areas. The thrust in mineral resources development has been focused on the base metals which as envisioned would provide the foundation for sound economic development. In recent years, however, major strides have been made to promote the exploitation of the non-traditionally produced commodities especially the non-metals. Development likewise point to the growing determination of the government and the private sector to expand benefits from mining. This is evident in efforts to promote downstream processing operations and to encourage investments in exploration, exploitation and mineral processing ventures. The integrated mineral resources development plan, therefore, proposes the vertical integration of mining extraction activities with the processing requirements of the metal-based manufacturing industries. While mining activities shall continue to be carved out by the private sector, the government shall influence the industry's direction by establishing legitimate goals to maximize benefits through an integrated approach giving due consideration to environmental protection and resource needs of future generations. STATUS OF MINERAL RESOURCES DEVELOPMENT Extent of Exploration Application of modem exploration tools and concepts complemented by intense infusion of foreign technology and rapid development of indigenous expertise have enabled the government jointly with the private sector to geologically explore vast grounds and to expand the country's known mineral resources. Approximately 93% of the country's total land area of 30 million hectares have been geologically mapped on at least 1:250,000 scale. Of the 93% total coverage, roughly 37.54% or 11.07 million hectares have been done on scales of 1:10,000 to 1:50,000 as shown in Figure I. Only 3.07% of the total coverage have been surveyed in detail. Land gravity surveys have been undertaken on over 45% of the total, while some 77,300 line km. of aeromagnetic surveys have been flown over 110,000 sq. km. of land area. Processing of data is almost complete. Total magnetic intensity and IGRF-corrected maps are available. Geological interpretation utilizing Landsat MSS imagery and low altitude aerial photography augmented by ground truthing have been undertaken in several large areas with very encouraging results. Beach and offshore mineral surveys were initiated in 1977 and have covered substantial areas in northwest and southeast Luzon and northeastern Mindanao. Highly promising areas for detrital magnetite, chromite, gold and some heavy minerals including ilmenite and titanium have been identified. In the petroleum sector, two offshore oilfields are currently producing from Miocene limestone reefs in Northwestern Palawan. Another field from the same area will be on stream in mid-1982.

Jan 1, 1982

By C. J. Bechtoldt, H. C. Vacher

Quenched alloys, prepared by powder metallurgical techniques, were examined by microscopic and X-ray diffraction methods. The compositions and heat treatments were chosen so that the chromium and nickel solvuses could be determined and a range of composition and temperature, in which a eutectoid had been reported, could be explored in small intervals. The results showed no evidence of eutectoid reactions having taken place. ALLOYS containing iron, chromium, molybdenum, and nickel are extensively used in high-temperature applications. In order to provide basic information on the solid-state reactions that occur in such alloys, the National Bureau of Standards is engaged in an extensive study of the equilibrium phases in the binary, ternary, and quaternary systems involving these four metals. Attention was directed to the chromium-nickel system because a review of the literature showed conflicting data. A recent compendium1 of the constitution of binary alloys gives two diagrams for the chromium-nickel system, one a simple eutectic and the second a diagram showing a eutectoid and a eutectic. The eutectoid was the result of an allotropic transformation in a high-temperature solid-solution form of a chromium-rich phase called 0. The approximate temperature and composition of the eutectoid was 1200°C and 30 at. pct Ni. In the preliminary study, at the Bureau, no transformation was indicated in the chromium-rich phase of ternary and quaternary alloys in which chromium and nickel were major components. Similar results have been reported from time to time in the literature,2-8 In consequence of this situation, a systematic study was made of the constitution of heat-treated chromium-nickel alloys using metallographic and X-ray diffraction methods. The compositions of the alloys were selected so as to determine the boundaries of the chromium-rich and the nickel-rich phases and to explore in small intervals the ranges in composition and temperature in which the eutectoid had been reported. The results of this study are reported and discussed in the present paper. EXPERIMENTAL PROCEDURES Materials, Alloys, and Heat Treatments—The alloys were prepared by a powder metallurgy technique similar to that used in earlier work.9 Analyses of the materials used are given in Table I. Powders were mixed using a vibrator in 30-g batches of predetermined composition; 2-g pellets were made from the mixture by compacting at room temperature in a 1/2-in. diam mold under 60,000 lb per sq in. pressure. The pellets were sintered in dry hydrogen for 10 days at 1300°C. At the end of the sintering treatment the specimens were quenched in water, which caused a pale green film to form on the chromium-rich alloys and thin black scale to form on the nickel-rich alloys. The pellets after sintering were approximately 1/2 in. diam by 3/32 in. thick. The purification train9 was modified by inserting a tube of titanium turnings, heated at 700°C to increase the efficiency of the removal of oxgen and nitrogen. It was found in the earlier work that the sintering treatment lowered the carbon, oxygen, and nitrogen contents of chromium-rich alloys to 0.01, 0.01, and 0.004 wt pct, respectively, and that there was a loss of chromium. In order to ascertain the extent of this loss, the chromium contents of certain alloys were determined chemically. They included those alloys used to determine the parameter-composition curves and those of critical composition, used to determine the chromium-solvus microscopically; these analyses were made on samples after the final heat treatments. The mixed powder values for the chromium content of those alloys that were not analyzed were adjusted to conform to a smooth curve based on the chemically analyzed values for

Jan 1, 1962

By Henry Howe

THE data which are collected in Table 1 show that the ferrite of low-carbon steel and of electrolytic iron, like the network of hypo- and hyper-eutectoid carbon steel, inherits, either absolutely or relatively, the, grain size of the mother austenite from which it is born in cooling through the transformation range. This evidence is so abundant and concordant as to cast doubt on Prof. Jeffries' 1 hypothesis of what might be called reversed inheritance, by which he explains Mr. Ruder's observation2 that the grains of a "slug of pressed pure electrolytic iron powder"were coarser, perhaps 12,000,000µ2, after a 3-hr. exposure to 1000° C. than after a 3-hr. exposure to 1250-1300°, when they were about 35,000µs. His explanation implies a natural tendency of coarser austenite grains to yield finer ferrite ones than spring from finer austenite grains, so that the inheritance is reversed. What seems to me a simpler explanation is that the enormous grains which followed the 1000° heating were not inherited from the mother austenite, but were formed by the process of natal coarsening, by means of which Stead and Carpenter developed grains of this same order of magnitude in this same kind of material, on heating to 950° or to 1000°, the temperature used by Mr. Ruder, provided the exposure was not unduly prolonged, and it does not seem to have been in Mr. Ruder's case. That their coarsening was natal and not through inheritance is shown first by the fact that their grain coarseness was out of all proportion to their relatively low temperatures and short periods, and hence was very far in excess of the grain size which the mother austenite must have had, and second and more cogently by the fact that the shape of these coarsened grains was evidently determined by the temperature distribution which existed during and immediately after the Ar3 transformation, that during the passage from the gamma to the non-gamma state.

Jan 9, 1917

By W. A. Tiller, S. O’Hara

Grain refinement in stirred melts has previously been shown to arise from dendrite segmentation. The present work discusses experiments capable of distinguishing between remelting and mechanical effects which result from fluid motion. Nucleation of still un-dercooled Sn-1 pct Pb melts with a single-crystal seed shows that the grain size of the resulting ingots decreases as the undercooling at nucleation increases up to about 1Z0C, then increases again for larger undercoolings. Pure melts exhibit little refinement as a result of nucleation at any undercooling up to 12°C. The presence of fluid motion gives rise to grain multiplication in pure melts, but the magnitude of this effect is partly obscured by solid-state pain growth. The results indicate that two grain-multiplication pvocesses can operate, one which involves dendrite remelt ing, and another associated with fluid motion. Small ingot experiments indicated that an equiaxed transition can be caused by fluid motion. THE presence of forced convection during recales-cence of undercooled Sn-1.3 pct Pb melts has been shown to give rise to extensive grain refinement during solidification of small ingots of this alloy.''2 In Ref. 2, this refining was shown to result entirely from dendrite fragmentation and not from any nucleation mechanism. The purpose of the present work is to establish what processes lead to dendrite segmentation. Experiments were therefore carried out to determine the resulting structure of both pure tin and Sn-1 pct Pb melts, nucleated at various undercoolings by single-crystal seeds in the absence of forced or natural convection. By comparing the information obtained from these experiments with the results of Refs. 1 and 2, the effects of fluid motion on the ingot structure, and the details of the segmentation mechanisms, should become clear. The results of these experiments should also be helpful in understanding the solidification structures of larger ingots since natural convection, which is known to play an important role in determining ingot structure13 differs only in degree from the fluid-motion effects studied here. I) EXPERIMENTAL PROCEDURE The apparatus shown in Fig. 1 was constructed to permit controlled nucleation of stirred or unstirred supercooled melts using a seed oriented in the tin dendrite growth direction. Because the melts, in all cases, were capable of larger supercoolings than the one used in each experiment, the structure observed always resulted from a single nucleation event at the surface. The sequential procedure was as follows: Tin was melted under purified hydrogen and the oxide removed at a temperature of 450°C . Electromagnetic stirring was then used to clear the surface completely and during this time the guide tube was lowered underneath the melt surface. A pointed Czochralski-like crystal was then propagated from a specially prepared tin dendrite seed4 inside the guide tube while the temperature of the bulk melt was held just above the melting point, ensuring a perfectly clean seed which was epitaxial with the original dendrite. The dendrite and epitaxial seed, still in the guide tube, were withdrawn from the melt which, in turn, was undercooled to 2C° below the temperature at which the nucleation experiment was to be performed. The melt was held there for several minutes, both in stirred and in still-melt experiments. The temperature was next raised to the required value and the melt nucleated by the carrot seed inside the guide tube, the heater power being shut off just after nucleation. In experiments where quenching after nucleation was required, a high-density graphite crucible was used, but otherwise the procedure was as described above. The temperature gradient in the melt before nucleation was kept below 1C° per cm by carrying the experiment out just below the center of the heater shown in Fig. 1.

Jan 1, 1968

By V. A. Josendal, J. W. Wilson, B. B. Sandiford

Two radioactive tracers have been tested as a means of determining core saturation in multiphase flow studies. Cesium chloride was tried as a water-phase tracer, but complications in its use in low permeability cores resulted from sorption of cesium by the core or water-wet pads. IodoI3'ben-zene proved very satisfactory as an oil-phase tracer. The synthesis of iodobenzene from the sodium iodide as received from Oak Ridge is simple and direct. The tracer is insoluble in water and there was no evidence of sorption by any of the core materials used. use of the method to determine saturation profiles during capillary and dynamic desaturations and relative permeability measurements on oil-water and oil-gas systems is described. Comparisons of the dynamic and capillary methods of relative permeability determination were made using the tracer to check core saturation and saturation distribution. Other experiments are also described in which mobility of the oil phase at various saturations was measured by displacing labeled oil by flowing inactive oil. Similar experiments were made using water labeled with cesium 134. INTRODUCTION To be fundamentally sound. any laboratory method of determining relative permeability must meet, among others, the following requirements: (1) the core saturation between the pressure taps must be uniform, and (2) there must be uniform pressure difference between the phases in the region between the pressure taps. In practice it may eventually be shown that appreciable departure from either or both of these conditions may be permitted without introducing sensible error. but the magnitude of such error remains to be established. There is need then for methods of saturation determina- tion which will determine the saturation profile in a core as well as the overall average saturation given by the gravimetric or material balance methods usually employed. Methods which have received consideration include resistivity, X-ray absorption, gamma ray absorption, neutron diffraction, and radioactive tracers. The resistivity method has use in the determination of profile, but needs checking by independent means.' The X-ray method has been used successfully, but requires elaborate equipment and calibration.23 4 At present gamma ray absorption is practical only for very large cores. The neutron methodo does not appear adaptable to measurement of saturation variations in short cores since it is a scattering method. Russell, Morgan and Muskat' employed radiovanadium in a study of the mobility of interstitial water. Coomber and Tiratsook used radioiodine as an oil-phase tracer and measured profiles in unconsolidated sand packs. ADVANTAGES OF RADIOACTIVITY METHOD A thorough discussion of the theory of radioactive tracers and the technique of using them to study the movement of fluids in sands is given by Coomber and Tiratsoos and need not be detailed here. The method requires less equipment than the X-ray method; the counting problem is similar, but the generation and regulation of the radiation, which is a major problem in the X-ray method, is not required in the radioactivity method. In the X-ray method it is customary to add up to 20 per cent of absorber, usually an iodine compound, to the phase to be labeled. while in the radioactivity method a mere trace is sufficient. A possible disadvantage of the radioactivity method is that a disproportionate amount of the activity comes from the portion of the core nearest the counter, whereas the X-ray gives a true average saturation in the region traversed by the beam. 011 the other hand, this characteristic of the radioactivity method renders it peculiarly adaptable to study of certain types of radial variations in saturation. As with the X-ray method, a tracer may be added to either the oil or water-phase. An oil-phase tracer is more generally useful. however. since it permits work on oil and gas in the

Jan 1, 1952

By G. W. Thomas

A simplified mathematical model of underground conduction heating in a system of limited permeability is presented. The model applies to underground retorting of oil shale, or to reservoirs containing extremely heavy oils. We assume that heat is introduced at a constant rate into a horizontal fracture which communicates between wells. The radial temperature distribution along the fractured surface is approximated by a step-function. Heat transfer away from the fracture is assumed to be by vertical conduction, and all convection effects are neglected. The model also takes into account the possible temperature dependence of thermal conductivity. A general expression for calculating the growth of the step-function temperature distribution with time is derived. The use of this expression and solutions to the one-dimensional heat equation make it possible to construct isotherms. Expressions for calculating .oil recovery, well spacing and heat efficiency are also given. An example calculation is presented for the conduction heating process in oil shale. Finally, the effect of the heat transfer coefficient between the gas and the fracture boundaries is investigated. INTRODUCTION Thermal recovery processes of oil recovery fall into the four general areas of hot fluid injection, forward combustion, reverse combustion and conduction heating. The first three of* these processes have been rather extensively studied in the past decade from both the experimental and theoretical points of view. As a result, it is possible to make reasonable engineering predictions and analyses of these processes. Little attention has been devoted, however, to the conduction heating process other than to note its possible utility.I To a certain extent, conduction heating cannot be divorced from the other regimes cited above, insofar as these provide the source of heat energy. In the conduction heating process, heat is introduced (either by combustion in the forward or reverse mode or by hot fluid injection) into a small fraction of the total reservoir thickness. This fraction may be either a streak of high permeability or an interwell fracture. Heat penetrates by conduction into the adjacent, less permeable regions of the oil-bearing rock, where the direction of conduction is essentially perpendicular to the streak or fracture. The heated product then drains by gravity or is gas driven to production wells. Conduction heating is probably most applicable to systems containing immobile bitumens such as tar sands and oil shale deposits and perhaps to low-permeability reservoirs containing highly viscous crudes. The mechanism also acts in combination with other thermal processes where fingering or overriding of a bed occurs. It seems probable that in at least one field test, conduction heating of this type was very influential.2 In this presentation we give a first approximation to some of the quantitative aspects of conduction heating. Marx and Langenheim3 treated a similar problem where they focused attention upon the injection interval, which spanned the total reservoir thickness. In their model, conduction heat losses to the bounding media imposed a practical limit on the calculated heated area. In the present study, however, we shall confine the injection interval to a small fraction of the reservoir thickness and assume it has no heat capacity. We therefore direct our attention to regions outside the injection interval into which the conduction of heat is beneficial. In particular, we will endeavor to locate specific isotherms in the media bounding the injection interval. Furthermore, we will construct our model to allow the thermal conductivity to vary arbitrarily with temperature. Thus the model will be applicable to underground retorting of shale where variations in thermal conductivity may be important. THEORY We consider the rock matrix to have no or, at best, a very limited permeability; consequently a

Jan 1, 1965

By H. S. Price, J. E. Warren

Techniques for studying the performance characteristics of heterogeneous reservoirs have been developed. The effect of permeability variation on both the steady-state and the transient flow of a single fluid has been investigated. Limited comparisons with field and laboratory data have been made. The physical models studied consist of random three-dimensional arrays of homogeneous porous blocks. The permeabilities of the individual elements are assigned according to a specific distribution function; uniform anisotropy is introduced by varying the relative dimensions of the blocks. A particular model is perturbed simply by re-arranging its elements at random. The behavior of a physical model is determined by digitally solving its numerical analogue. Based on computational experiments, subject to the restrictions implied by the assumptions that were made, the following general conclusions have been drawn. 1. The most probable behavior of a heterogeneous system approaches that of a homogeneous system with a permeability equal to the geometric mean of the individual permeabilities. 2. The effects of flow geometry and anisotropy on the most probable value of the effective permeability of a heterogeneous medium are finite but not significant. 3. The permeability determined from a pressure build-up curve lor a heterogeneous reservoir gives a reasonable value for the effective permeability of the drainage area. 4. A qualitative measure of the degree of heterogeneity and its spatial configuration are obtained from a comparative study of core analysis and pressure build-up data. It has been indicated that these conclusions are predicated on the assumption that the core analysis and the pressure build-up data represent true reservoir characteristics. Common sources of error in these data have been discussed. One of the most significant problems of reservoir engineering is that of determining the nature and the disposition of the heterogeneities that inevitably occur in petroliferous formations. It is with this problem that this paper is concerned. Extensive research, both theoretical and experimental, has provided reasonable descriptions of the physical processes that are associated with primary and secondaryrecovery mechanisms. Similar progress in the fie1d of numerical analysis and the evolution of truly high-speed digital computers have made it feasible to calculate the behavior of these mechanisms under almost any environmental conditions. It appears to be logical, then, to conclude that the mechanistic performance of any reservoir can be predicted with some confidence. Unfortunate1y, a necessary condition for the validity of this conclusion is that the formation be adequately described in a macroscopic sense. This condition cannot be satisfied ill general; therefore, a frustrating situation currently exists — it is possible to compute the effect of reservoir heterogeneity on behavior, but it is not possible to specify the heterogeneity itself. In many cases, the predicted behavior of a project is so completely dominated by irregularities in the formation's physical properties that the gratuitous assumptior, of a particular form for the variation can reduce the solution of the problem to a mere tautological exercise. While all porous media are microscopically heterogeneous, only macroscopic fluctuations need be considered since the fundamental concepts of reservoir mechanics are based on macroscopic quantities. For the simplest case in which only one physical property varies from point to point, the degree of heterogeneity that is apparent depends on three distinct factors, i.e., the nature of the property's variation (frequency distribution), its disposition (spatial distribution) and the inherent stability of the mechanism being studied. Consequently, any attempt to devise an unambiguous scale of heterogeneity must take these factors into account. Since this investigation is intended to be exploratory, it will be restricted to single-phase flow because of the simplification implied; both steady-state and unsteady-state behavior will be considered. Furthermore, to keep the number of parameters

By A. S. Michaels, M. C. Porter, A. Stancell

Previous laboratory studies have demonstrated that the injection of small quantities of reverse wetting agents during water displacement can increase oil recovery from unconsolidated porous media.9 In the present investigation, an attempt has been made to determine more fully the effects of reverse wetting treatments and to clarify the mechanism by which increased oil recovery is effected. Water-oil displacements were performed in beds of 140-200 mesh silica sand. Hexylamine slugs (injected after 0.25 pore volume of water throughput), when adequate in size and concentration, were effective in promoting additional oil recovery. Their effectiveness increased with the quantity of amine injected. However, slugs of sufficient size and concentration to stimulate oil production at water flow rates of 34 ft/day did not do so at 4 ft/day. Visual studies in a glass grid micromodel have shown that the stimulation of oil production, via ayueous hexylamine, is a result of transient changes in the oil wettability of the pore walls. if the amine slug is of sufficient size and concentration to induce significant changes in the adhesion-tension, large continuous oil masses will be formed. if the superficial water velocity is high enough to result in rapid desorption of the amine, a favorable "wettobility gradient" may be established across the masses; under such conditions, high oil mobility is observed, and increased oil recovery results. INTRODUCTION It is generally agreed that the efficiency of oil displacement by water in porous media is limited in part by capillary forces which cause the retention of isolated masses of oil — resulting in the so-called "irreducible minimum oil saturation". Recent estimates10 indicate that there are about 220 billion bbl of petroleum in United States reservoirs which are not economically recoverable with present techniques (such as water flooding). This amounts to almost five times the known recoverable reserves. It has been recognized for some time that a suitable alteration in the water-oil interfacial tension and/or the contact angle, as measured between the water-oil interface and the solid surface, should result in better displacement efficiency. Surface active agents can be used as interfacial tension depressants to accomplish this objective,1,2,7 but unfortunately, the additional oil recovery is seldom commensurate with the treatment cost. In contrast to interfacial tension depressants, the effect of contact angle alterations on water-oil displacements has received relatively little attention in the literature. It is known that the wettability affects the displacement process.3,4,11 Displacements in water-wet systems generally result in lower residual oil saturations than those in oil-wet systems. The effect of "'transient" wettability alterations concurrent with the displacement process have been investigated by Wagner, Leach and coworkers,4,11 wherein it has been demonstrated that the establishment of water-wet conditions during water flooding of oil-wet, oil-saturated porous media is accompanied by significant increase in oil displacement efficiency. Michaels and Timminss studied the effects of transient contact angle alterations resulting from chromatographic transport of reverse wetting agents through unconsolidated sand. It was demonstrated that chromatographic transport of short-chain (C4 through C8) primary aliphatic amines can improve oil recovery and that the recovery increases with the quantity of amine injected (i.e., with either the amine concentration or the volume of the slug injected). Circumstantial evidence indicated that the increased displacement efficiency resulted primarily from transient changes in wettability of the porous medium. In the Present investigation. additional information has been obtained on the effects of reverse wetting treatments and the mechanism by which increased oil recovery is accomplished. Hexylamine slugs of varying size and concentration were chromato-

Jan 1, 1965

By Frank T. Weems

The basic operating properties of the DorrClone are discussed and certain metallurgical applications which exploit these properties are presented. An effective method of controlling the consistency of the DorrClone underflow is described. THIS paper presents in a general manner some of the applications of the DorrClone (Dutch States Mines cyclone) in the metallurgical and heavy chemical fields. Technical or mathematical analyses of the design factors and of the mechanism whereby the cyclone effects its separation have been capably treated in a series of previous papers.'-' Engineering and operational aspects of a few installations and laboratory tests have been chosen to emphasize the unique classification characteristics of the DorrClone and to illustrate its utility in certain difficult situations. Its use in coal washing and heavy-media operations are not considered. Some 25 or more industrial installations and applications of the DorrClone in the metallurgical fields are operating at present, ranging from single trial units to complete batteries composed of a multiplicity of units. The data were obtained in part from some of these, and in part from operations of units at the Westport Mill, the laboratories and testing plant of the Dorr Co. at Westport, Conn. The DorrClone was developed, in principle, by the research dept. of the State Mines, Limburg, Holland, during World War 11. Its first use was to dewater a —50 micron sand used as a dense medium in coal washing. The metallurgical wet cyclone was first publicized by Driessen and Yancey and Geer.6 ince then the cyclone has been studied as a means of de-watering, heavy media concentration, and fine size classification. The DorrClone is shown in Figs. 1 and 2. pulp to be treated is pumped through the feed nozzle, A, into the feed cylinder, B. The feed nozzle is of rectangular cross-section, the outermost edge being tangent to the feed cylinder walls. Feed shims, C, may be attached to the innermost walls to adjust the open area of the feed nozzle available to entering pulp. The tangential injection of feed under pressure supports the vortex action of the pulp that has been detained in the cyclone; the centrifugal forces of this vortex throw the coarse particles to the walls of the cone, D, where they collect and pass downward and out through the adjustable apex opening, E. The fine particles contained in most of the water that entered with the feed pulp move to the inner spiral of the vortex and flow out through the vortex finder, F, contained in the orifice plate, G, which forms the top of cyclone chamber. The overflow pulp may then spray freely into a spray chamber, H, surmounting the cyclone or pass directly into a 90' elbow, I, of somewhat larger diameter than the vortex finder. Installation Where more than one cyclone is to be installed, best overall efficiency and ease of control are obtained with equal distribution of both pulp volume and weight of solids among the cyclones, and with the same pressure being maintained at each cyclone feed nozzle. Fig. 3 is a suggested arrangement for feeding two or more cyclones in parallel. Feed pressure control: Under normal operating conditions it is not unusual for some surge to exist in the feed that is to be sent to the cyclone. As will be discussed later, surges in the solids content of the pulp are usually not of great operational consequence, but surges in feed volume are a problem. Several systems have been devised for absorbing surges in feed volume, each involving a surge tank. If the surges are small and fairly rhythmical, a tank of sufficient size might be installed to dampen out

Jan 1, 1952

INTRODUCTION The meaning of cash flow has already been defined as after-tax earnings available for use by the corporation or business. It is therefore apparent that a significant part of mineral economic evaluation is concerned with correct and accurate computation of taxable income under existing tax laws, and the proper treatment of tax allowances and deductions. The various methods of depreciation, depletion, and amortization are all means whereby the investor can earn back or recover, before taxes are imposed, the investments made for purposes of doing business or for the production of income. Each of these three types of tax deduction are applied to different investment items such as property, equipment, buildings, research, etc. The cost of buildings, equipment, machinery, and transportation equipment are representative of capitalized business costs that must be recovered by depreciation over the life of the asset. The property investment costs of acquiring natural resources such as minerals, oil and gas, and standing timber are examples of investment costs that may be recovered by depletion. Several other capitalized business costs, such as research expenses and trademark expenses, may be recovered by amortization. The term "capitalized" means that the in- vestment cost is not deducted as an operating expense in the year that it is incurred, but is spread out over a period of years through depreciation, depletion, or amortization deductions according to lawful schedules. Of course, any capitalized cost can only be deducted once. The other costs of doing business which are not capitalized, such as labor and operating costs, are said to be "expensed," i.e., they are written off for tax purposes in the year they occur. DEPRECIATION Depreciation, according to US tax law, is a deduction from taxable income permitted as a "reasonable allowance for the exhaustion, wear and tear of property out of its use or employment in a business." The original cost of the property is the usual basis for determining depreciation. Since depreciation write-off is applied over the useful life of the asset, the first step is to determine the estimated useful life or its allowable depreciable life. These two lives may or may not be the same. No average life is applicable to all types of depreciable assets under different conditions or in different types of business. It is the taxpayer's responsibility to determine the useful life of the asset according to his own operating conditions and experience, within the guidelines set by the US Internal Revenue Service (IRS). The Class Life Depreciation Range (ADR) System permits the taxpayer to choose from a range of depreciable lives at 20% below or above the present guideline life. The second step required in calculating depreciation is to determine an estimated salvage value for the asset. The salvage value is an estimate, made at the time of acquisition, of the amount of money which can be realized by the sale of the asset at the end of its useful life. The four depreciation methods most widely used in the US are: (1) the straight-line method, (2) the declining balance method (double declining balance), (3) sum of the years - digits method, and (4) units of production method. Straight-Line Depreciation Method The straight-line depreciation method is the simplest for calculating depreciation. The depreciation for each year is determined by subtracting the estimated salvage value from the first cost of the asset and dividing the result by the years of life. The method gives equal annual deduction over the depreciable life. Salvage value equal to or less than 10% of the first cost may be ignored. [first cost - salvage Straight line = years of life ] Declining Balance Depreciation Method Declining balance depreciation, or double declining balance (DDB) as most frequently applied, is used for accelerating depreciation in the early years of the asset and is the fastest depreciation method allowed by the US Internal Revenue Service. The declining balance

Jan 1, 1980

By D. E. Norquist, J. E. Moeller

The results of wear on marked balls, 4, 31/2, 3, and 2 in. diam are given. All balls were forged steel of practically the same chemical analysis and hardness. The results indicate that balls in a given mill for a given length of time will have equal diameter losses, regardless of size.

Jan 1, 1950

By H. L. Overton, J. B. Lipson

The concept of sodium single-ion equivalent activity as developed by Gondouin, Tixier and Simard,' was used to determine the filtrate resistivity-activity relationships for 150 laboratory and 49 field drilling muds. With the exception of the native and phosphate high-resistivity muds and the calcium-treated muds, the filtrate resistivity-activity relationships were found to be within 20 per cent of that for pure sodium-chloride solutions. Gondouin, et al,* previously suggested a method of treating the high-resistivity muds in quantitative electrical log analysis. A new treatment is presented which uses the P-alkalinity of the filtrate as a correlating parameter in handling the calcium-treated muds in quantitative electrical log analysis. INTRODUCTION Gondouin, Tixier, and Simard in 1957 introduced new interpretive procedures for analyzing the self-potential curve of electric logging. They re-emphasized the previous work of the Schlumberger brothers, Wyllie and others, and re-defined the activity concept in terms of a number which they termed the effective resistivity. Their procedures are particuarly effective when the connate-water resistivities are less than 0.08 ohm-m or greater than 0.3 ohm-m. For gyp muds where the calcium content is known, a procedure was presented which allows calculation of the connate-water resistivity with improved accuracy. These procedures assume that one of the fluids in the system is a pure sodium-chloride solution, but they are not limited by the assumptions. For the fresh and saline connate waters, for instance, it is assumed that the drilling mud filtrate is a pure sodium-chloride solution. For the gyp mud interpretive procedure, it is assumed that the connate water is a pure sodium-chloride solution. Brown,' on the other hand, has suggested that it would be extremely fortuitous if drilling mud filtrates could all be treated electrochemically as sodium-chloride solutions, and he cites some examples where discrepancies exist. A study has been made of the deviations of sodium single-ion equivalent activities, an of drilling mud filtrates from that of pure sodium-chloride solutions. Drilling fluids in common use on the Texas and Louisiana Gulf Coast and laboratory-controlled muds were investigated. Using a cell patterned after Gondouin, et al, (Fig. I), the equivalent activities were measured. It should be remembered here that the activity measured is merely an index which uses numbers that represent the sodium single-ion equivalent activity a,., of the mud filtrate. The membrane potential developed across the cell was determined with a null-type potentiometer. The resistivities of the solutions were measured by the four-electrode method. The resistivity system was calibrated by using saturated KC1 solution as a standard. Chemical properties of the drilling fluids studied were measured with standard API-approved equipment for the particular property in question. For the activity measurements, two shale discs were chosen to serve as cationic membranes. These two discs were both 1 1/4 in. in diameter and, respectively, 1/4 and 1/8 in. in thickness. The membrane potential, Em, was determined using the shale discs as cationic selective membranes separating two standard sodium-chloride solutions of known mean activity. Then by comparing the constant Km in the following equation, with the theoretical value of — 59 mv at 75OF, the perfection of the membrane was determined. The variation of Km for the two membranes used in this experimental work, as well as one of the many unsuitable membranes, is shown in Fig. 2. The measurement of the activity of mud filtrate involves the balancing of its combined individual activities against a standard solution of known mean activity. At the point where there exists zero potential across the membrane cell, the activity effects may be assumed equal for all practical purposes. The mechanics of making a null reading of this sort requires having a large number of different mean activity solutions. A more practical method involves

By J. L. Carne

ORE control is a matter of planning and supervision based on a foreknowledge of the content and distribution of ore. The Inspiration orebody is predominately a copper-sulphide blanket, overlain by an oxidized zone of copper silicates, copper carbonates, and barren capping. Most of the ore is treated by leaching, and the optimum requirement of this method is the goal toward which the mining schedule is aimed. Both open-pit operations and underground mining are used. The problems of ore control at Inspiration are: I—maintenance of metallurgical balance, 2—conservation of oxide ore, and 3—mining out of the orebody at an average grade which parallels that of the reserve figure. Prior to the fall of 1926 all of the Inspiration division ores were treated by flotation, but the mixed ores of the Live Oak—Keystone section carried too large a proportion of copper in the form of chrysocolla for satisfactory flotation. For this type of mixed ore the leaching process now in use had been worked out, and the plant was put into operation in the fall of 1926. The process is known as the ferric sulphate leach, and the leaching solvents contain both ferric sulphate and free sulphuric acid. Sulphide copper in the form of chalcocite is leached by the ferric sulphate, and the oxidized copper minerals are dissolved by the sulphuric acid. From the beginning of operations to the present time the plant has treated some 67 thousand tons of ore having an average copper content of 1.142 pct, of which 0.605 pct has been in the form of oxidized copper minerals. The plant is currently being operated at the rate of 11,500 tons per day, 7 days per week. However, this tonnage is produced in 6 days of mine operation. The average current grade is approximately 1 pct copper, of which about 50 pct is present as oxidized or acid-soluble copper. Close control of metallurgical procedures is a vital factor in operation of the plant and is in turn dependent upon careful regulation of feed. An excess of sulphide copper requires the presence of more ferric sulphate in the leaching solvents. On the other hand, an excess of ferric sulphate in the presence of a low sulphide feed seriously reduces the efficiency of electrolytic precipitation in tank house operations. This condition must be avoided at all costs. An additional factor must be considered. It is known that the sulphide content is predominant in the ore which remains in the reserve. To maintain the necessary sulphide-oxide proportion, so that the entire reserve will mine out at both a grade and a sulphide-oxide ratio that will be permissible throughout the remaining life of the property, it is obvious that the oxide reserve must be currently preserved to as great an extent as possible. Fortunately much of the oxide reserve is found in the open-pit operation. This happy circumstance adds flexibility to the procedure and makes it possible to maintain the close regulation necessary to insure best overall results, both mining and metallurgical. The two ore streams from the pit and underground operations join at the coarse-crushing plant, and proper control must be exercised before the ore reaches this point. The open pit is currently providing about 55 pct of the ore and underground operations the balance of production. Both operations are under a single superintendent so that the closely interrelated problems may be more readily controlled. The least flexible of the two producing units is the underground mine, where the mining method is block caving. Although a block of ground may be undercut and developed, characteristics of the material determine its caving action. The ore in the chutes must be drawn. The area currently mined from underground is in a zone predominately sulphide and is generally

Jan 1, 1953

By J. L. Carne

ORE control is a matter of planning and supervision based on a foreknowledge of the content and distribution of ore. The Inspiration orebody is predominately a copper-sulphide blanket, overlain by an oxidized zone of copper silicates, copper carbonates, and barren capping. Most of the ore is treated by leaching, and the optimum requirement of this method is the goal toward which the mining schedule is aimed. Both open-pit operations and underground mining are used. The problems of ore control at Inspiration are: I—maintenance of metallurgical balance, 2—conservation of oxide ore, and 3—mining out of the orebody at an average grade which parallels that of the reserve figure. Prior to the fall of 1926 all of the Inspiration division ores were treated by flotation, but the mixed ores of the Live Oak—Keystone section carried too large a proportion of copper in the form of chrysocolla for satisfactory flotation. For this type of mixed ore the leaching process now in use had been worked out, and the plant was put into operation in the fall of 1926. The process is known as the ferric sulphate leach, and the leaching solvents contain both ferric sulphate and free sulphuric acid. Sulphide copper in the form of chalcocite is leached by the ferric sulphate, and the oxidized copper minerals are dissolved by the sulphuric acid. From the beginning of operations to the present time the plant has treated some 67 thousand tons of ore having an average copper content of 1.142 pct, of which 0.605 pct has been in the form of oxidized copper minerals. The plant is currently being operated at the rate of 11,500 tons per day, 7 days per week. However, this tonnage is produced in 6 days of mine operation. The average current grade is approximately 1 pct copper, of which about 50 pct is present as oxidized or acid-soluble copper. Close control of metallurgical procedures is a vital factor in operation of the plant and is in turn dependent upon careful regulation of feed. An excess of sulphide copper requires the presence of more ferric sulphate in the leaching solvents. On the other hand, an excess of ferric sulphate in the presence of a low sulphide feed seriously reduces the efficiency of electrolytic precipitation in tank house operations. This condition must be avoided at all costs. An additional factor must be considered. It is known that the sulphide content is predominant in the ore which remains in the reserve. To maintain the necessary sulphide-oxide proportion, so that the entire reserve will mine out at both a grade and a sulphide-oxide ratio that will be permissible throughout the remaining life of the property, it is obvious that the oxide reserve must be currently preserved to as great an extent as possible. Fortunately much of the oxide reserve is found in the open-pit operation. This happy circumstance adds flexibility to the procedure and makes it possible to maintain the close regulation necessary to insure best overall results, both mining and metallurgical. The two ore streams from the pit and underground operations join at the coarse-crushing plant, and proper control must be exercised before the ore reaches this point. The open pit is currently providing about 55 pct of the ore and underground operations the balance of production. Both operations are under a single superintendent so that the closely interrelated problems may be more readily controlled. The least flexible of the two producing units is the underground mine, where the mining method is block caving. Although a block of ground may be undercut and developed, characteristics of the material determine its caving action. The ore in the chutes must be drawn. The area currently mined from underground is in a zone predominately sulphide and is generally

Jan 1, 1953

By U. N. Bhrany, N. Arbiter

F. C. Bond (Consulting Engineer, Processing Machin-Dept., AllisChalmers Manufacturing Co., Milwaukee) This is a very comprehensive paper. It deals with 1) size distribution functions, 2) energy-particle size relationships, and 3) criteria for grinding efficiency. The method of analysis used consists of tracing the increase and decrease of successive Tyler scale screen fractions as grinding progresses. The treatment is based upon the assumption that the work done in grinding is concentrated entirely upon the coarse size fractions in the mill. It is postulated that when a feed consisting entirely of particles of one size fraction is batch ground that finer fractions increase at a constant rate as grinding progresses, until finally each in turn begins to be ground and to decrease in amount. It is shown that under these conditions the empirical Gates-Gaudin-Schuhmann log-log size distribution plot can be logically explained. It is difficult to imagine any condition of ball mill operation, or any probability analysis, which would result in no grinding contacts on the finer particles present. However, the authors' postulate must be accepted if their conclusions are to be sustained. They conclude that the equation expressing the for mation rate of each size fraction can be resolved into the G.-G.-S. size distribution equation, and that the curvilinear portion of the log-logsize distribution line represents the only particle sizes which are being ground as well as produced. The author's statement that the G.-G.-S. Eq. 4 applies only to the finer portion of the size distribution while the Rosin-Rammler Eq. 6 applies over the whole size distribution range certainly supports the view that an exponential equation of the Rosin-Rammler type corresponds more closely to actual size distribution than does the power law G.-G.-S. equation. However, the G.-G.-S. equation is used to derive Eq. 8 and Eq. 9, which state that the energy required to grind a ton through a particular mesh varies inversely as the particle size to the exponent -a , whic h is the slope of the extended G.-G.-S straight line. This is stated to apply to both batch and continuous grinding. The G.-G.-S. slope exponent thus becomes very important in all crushing and grinding analyses by this method. The authors include some data which can be used to test their conclusions, as described below: We will consider grinding a feed composed of 14 mesh (1190u) to -200 mesh (74u). The average of 17 Allis-Chalmers grindability tests on quartz gives a work index (Wi) of 12.77. The -200 mesh product should have 80 pct passing 41.6u or P = 41.6. The equivalent 80 pct passing size of the 14 mesh feed particles is considered to be the particle diameter, as is true for spherical particles. Thus F is 1190u. In the work index equation 10 wi 10 wi 127.7 127.7 P f 41.6 1190 = 16.1 kw-hrperton. In Table I the ball mill exponent for quartz is 1.0 and from Eqs. 9 and 10: w=kw-hr/ton =Constan/(do/di) = 16.1 ton so that (1190) Constant/ Constant and the Constant equals 259. The tons per horsepower-hour value is 0.0464, in agreement with Table IV. The exponents for other minerals are listed in Tabl e 11. In Table IIA, the kilowatt hours per ton required to grind these minerals between the same sizes as the quartz, both by the authors' method and the work index equation, are listed. According to the authors kw-hrperton = 259/(16.1)".

Jan 1, 1961

By R. P. Marshall

This note describes positive evidence that hydrogen in titanium alloys diffuses under the influence of a thermal gradient. The experiments confirmed the expected similarity of this system to the H-Zr system and produced qualitative information on the magnitude of the heat of transport for hydrogen in titanium alloys. This effect has practical importance in applications where titanium alloys operate in a thermal gradient under conditions such that hydrogen embrittlement may develop within the expected lifetime of the structure, for example in supersonic aircraft service. The theories of thermal-gradient diffusion may be found in recent reviews and experimentation in this field by Shewmon,1 Sawatzky,2 Markowitz,3 and Darken and Oriani.4 The basis for Soret thermal diffusion is the "heat of transport", denoted by Q*. This is an energy term for the amount of heat (energy) transferred by atomic diffusion, and can be either positive or negative signifying a driving force for atomic diffusion to the cooler or hotter region, respectively. Experimentally, Q* is determined by plotting the natural logarithm of the point-to-point concentrations after exposure to a thermal gradient vs the respective point-to-point temperatures in units of the inverse absolute temperature scale. The slope of the straight line is Q*/R, R being the gas constant. For the titanium experiments, specimens (3 by 1 by 1/16 in.) of a) commercially pure titanium, b) the 8 Al-1 Mo-1 V titanium alloy, and c) the 6 A1-4 V titanium alloy were exposed to thermal gradients of 6° to 115°C per in. for 18 days at 250° to 550°C. The temperature gradient was in the direction of the 3-in. dimension. The hydrogen concentrations existing after exposure were determined on sections of the specimens by vacuum extraction. Three sections of unexposed sheet of each alloy were analyzed as standards to allow a check for hydrogen

Jan 1, 1965

By A. Weiss, W. A. Coster

The search for ore is characterized by a great variety of uncertainties, each of which can be resolved by obtaining information at a cost. Three exploration stages may be distinguished in which information can be obtained, based on inherently different cost structure: 1) the reconnaissance or scanning stage, 2) the examination or screening stage and 3) the project stage. This exploration sequence would roughly coincide with replies to questions as to: where, how, why and how much. As complete information in each stage is too costly to acquire, we must assess the probability that adequate information will be furnished by a limited amount of effort and attempt to maximize the expected value of ore to be found for a given exploration expenditure. Several problems are discussed relative to the mathematical formulation and possible computer procedures for determining the optimum allocation of a total exploration effort. Based on the available literature possible solutions to this problem are indicated. When scanning through published literature on mineral exploration, it seems at first that there is little room for substantial and systematic improvements in the efficiency with which the problem of discovery of new ore deposits is presently being handled by large organizations active in exploration. It appears equally difficult to offer something relevantly new which could generate a break-through in exploration thinking and cause a marked improvement in the ratio of successful undertakings to the total exploration projects started. On the other hand, a number of experts in recent years have urged that the problems in prospecting be analyzed in a more skillful and sophisticated manner than in the past. Their thesis is based mainly on the premise that the cost of finding mineral deposits will increase with increasing technical difficulties. In present large-scale exploration ventures, especially in those regions where outcrops and immediately visible geologic evidence are scarce, older methods of prospecting are increasingly being replaced by more indirect and systematic geophysical and geochemical methods with more emphasis on correct statistical design of exploration patterns. It is perhaps the oil industry which has made the greatest strides in developing statistically, well justified patterns in its search for oil; possibly because oil seepages, as indicators of the presence of oil, left off at an earlier point in time than did outcrops as indicators of base metals. Among the pioneers in the field of statistics and operations research, as applied to mineral exploration, the excellent works of Allais, Slichter, Koop-man, Mickey, de Guhnin, Krumbein, Griffith, Sichel, Krige, de Wijs, and Matheron may be mentioned. It should, however, be pointed out that three subjects of basic importance to the allocation of exploration efforts have so far received only scant attention in published literature, even by the aforementioned authors. In the first place, a more precise definition of the object of the search, i.e., a minimum-size orebody, seems to be outstanding. Secondly, a decision model is not a widely-discussed subject. Thirdly, no concentrated attempt so far seems to have been made to establish probabilities of the occurrence of economic deposits on a regional scale and to arrive at a slightly more accurate knowledge of the probability laws governing the deposition of ore. It is not the purpose of this paper to create new order in the complexity of uncertainties in which our decision processes in exploration necessarily have to be conducted. Rather, we intend to highlight some of the outstanding, hitherto unmapped territory, where diligent compilation and a utilization of present and historic data may statistically improve our chances of success. We hope that the contribution of this paper will lie in bringing allocation problems into better focus and in increasing the mineral industry's realization of the magnitude of the factors involved. THE MINIMUM-SIZE OREBODY The concept of target specifications for what constitute a minimum-size orebody in exploration is a rather variable one and depends to some extent on

Jan 1, 1963