Search Documents

By H. W. Fleming, R. R. Brooks

The Clearwater Deposit, a small occurrence of massive-sulphide mineralization enclosed in an envelope of disseminated-sulphide mineralization, was discovered as a result of an aeroelectromagnetic survey in New Brunswick in 1955. Subsequently, additional flights were made over the deposit and various ground follow-up methods were employed to ascertain the nature of their response. All methods proved to be capable of indicating the occurrence with varying degrees of precision. This paper is a presentation and a discussion of the magnetic and electromagnetic data, both airborne and ground, and the self-potential data obtained over this deposit. The attention of Canadian mining exploration companies was initially focussed on New Brunswick by the discovery, in 1953, of the Brunswick ore body. Ground electromagnetic surveys were of considerable assistance in locating and developing this deposit. In the following year, 1954, the Heath Steele base-metal deposits were discovered by employing an airborne electromagnetic device and ground electromagnetic surveys to follow-up the airborne indications. These discoveries demonstrated clearly that both airborne and ground electromagnetic techniques were applicable for locating the near-surface deposits of massive sulphides present in this region. Airborne electromagnetic equipment became available on a contract basis through Aeromagnetic Surveys Limited early in 1955 and was used extensively thereafter throughout a large portion of northern New Brunswick. It was as a result of one of these surveys with Aeromagnetic Surveys Limited equipment in the summer of 1955 that the Clearwater Deposit was detected. The airborne anomaly was located on the ground by means of an electromagnetic survey using a Swedish Slingram unit. During 1956, the region was investigated more thoroughly- in the hope that other sulphide bodies might be located. Airborne electromagnetic coverage was extended and the Clearwater Deposit itself was used as a testing ground for the various methods and types of equipment currently in use in New Brunswick. In this paper com prative data are presented for the airborne magnetic, airborne electromagnetic, ground magnetic and ground electromagnetic surveys over this deposit. The ground electromagnetic instruments used were the McPhar vertical-loop unit and the Slingram and Boliden horizontal-loop units. Self-potential surveys were carried out with a McPhar unit. From the geophysical point of view, all ground electrical methods employed indicated the presence and the extent of the massive-sulphide mineralization with a good degree of accuracy. The horizontal-loop methods however, defined the contacts more precisely as is to be expected when dealing with a flatly-dipping

Jan 1, 1961

By Norman F. Schulz

The technical excellence of separation achieved in a mineral concentration process, or any other process where two constituents of any kind are physically separated from each other, is expressed uniquely and quantitatively by the separation efficiency (E,): E. = (R - R,), where R is the percentage of the valuable constituent and R, the percentage of the waste appearing in the concentrate. The separation efficiency meets all the usual requirements of a true efficiency and is independent of the nature of the process used. It can be visualized as the percentage of the feed material that, in effect, actually undergoes complete separation while the rest of the feed is distributed randomly and unchanged into the two separation products. The problem of quantitatively assessing the results of a mineral concentration process recurs whenever one tries to decide which of two similar test results is more desirable. Although economic factors may greatly influence the eventual choice of process, they usually have very little to do with the technical excellence of a separation and are, therefore, not a consideration at this point. Concentrate grade and recovery, used simultaneously are the most widely accepted measures of performance, and if both the grade and the recovery are greater for one case than for another, the choice is simple. However, when the results of closely competitive operations are compared, it often occurs that one set of test results shows a higher grade and a lower recovery than another and the choice is no longer obvious. Laboratory investigations of a mineral concentration process applied to a given feed material can often be designed to yield successive paired values of concentrate grade and recovery over the practical range of the course of a separation process.' These results may be plotted as some combination of the cumulative values of grade, weight recovery, or unit recovery in the concentrate. The resulting correlations for two slightly different sets of conditions can then be compared qualitatively, or even quantitatively in particular cases.2 Such comparisons are not applicable, however, to individual concentration test results. The reduction of a separation test result to a single efficiency figure, preferably as a percentage of ideality, would afford the simplest and most precise means for comparing the results of competitive tests. The effective proportion of feed material that undergoes complete separation in any real concentration process, hereafter referred to as "separation efficiency," is such a figure. It is essily calculated from the usual feed and product data and can even be evaluated, in many cases, from data in publications describing various concentration processes. Mineral Concentration Processes and Their Evaluation Mineral Concentration Processes: The overall concentration process may be subdivided into a liberation process which severs the bonds between unlike mineral species, and a separation process, which sorts the resulting liberated particles into two products—a concentrate rich in mineral A and a tailing lean in mineral A. A perfect mineral concentration process is represented by the idealized scheme shown in Fig. 1 in which liberation and separation are not differentiated. The feed mixture consists of two mineral types which are completely liberated early in the process and then sorted by type. All the mineral particles of type A which respond uniquely to the selection forces acting in the ideal separator, report to the concentrate, while all the mineral particles of type non-A report to the tail-

Jan 1, 1971

By Frank Rieber

THE study of earthquakes long ago developed the fact that by studying the travel times of the various groups of waves from the same earthquake, as received on seismographs at varying distances, major discontinuities in the earth's crust could be reasonably well defined. Using this original discovery, and adapting the mathematical methods originally devised in seismology for determining the depth below the surface at which such discontinuities occur, other investigators were able to develop a new branch of seismology which promises to be extremely useful in extending the boundaries of our information on local conditions below the earth's surface. As an outstanding example of the applications of this work, we have the so-called "seismograph" exploration of the Gulf Coastal region, which has resulted in the discovery of a large number of salt domes, from many of which oil fields may be developed. The original seismological investigations, from which different major discontinuities in the earth's crust were deduced, as well as these later successful experiments in locating salt-dome structures, involve relatively simple problems. In each case we are dealing with surface materials in which elastic waves will be propagated with a fairly definite velocity. Under this material we have something quite different, in which the elastic waves will be propagated with a much higher velocity. The transition from one medium to the other is likewise relatively abrupt. For example, in the Gulf Coastal region the sedimentary formations have a velocity of 6000 or 7000 ft. per sec., while the velocity of elastic waves in the cap rock of a salt dome may exceed 15,000 ft. per sec. Under such conditions, delineation of the formation is a relatively simple matter, and, furthermore, there is very little opportunity for ambiguous results.

Jan 1, 1928

By E. S. Moore

ALTHOUGH nickel was in use in alloys long before the Christian era, the metal was not discovered until 1751, when Cronstedt recognized it in niccolite from Sweden. The Chinese apparently used a nickel alloy, known as "pahfong" thousands of years ago and Persian coins dating farther back than 200 B.C. consist of a copper-nickel alloy similar to some of those now in use in coins. The pure metal was first prepared about 1775 and the metal was refined for commercial purposes some years later. The nickel industry did not attain much prominence until the ores in New Caledonia were discovered. Exports of ore from that country began in 1875 and since that date there has been an ever-increasing world output of nickel and a great expansion in the uses of the metal except during short periods of commercial and industrial depression. Previous to 1875, for a number of years, Norway was the chief producer, although the metal was mined also in Sweden, Finland, Italy, Germany, Austria, France, China, the United States and other countries. With the dis-covery of nickel deposits in the Sudbury district, Ontario, a new factor was introduced into the nickel industry and Sudbury has since become the overwhelmingly dominant feature in the world's nickel markets. Although nickel was first discovered in Canada as early as 1848 near the north shore of Lake Huron and nickel-copper sulfides were found in 1856 by surveyors working close to the present location of the important mine at Creighton, in the Sudbury field, it was not until 188.3 that any importance was attached to this field. In that year copper sulfide attracted attention where the main line of the Canadian Pacific railroad was being projected across Canada to the Pacific coast. Mining began near this discovery in 1889 but in the meantime operations had begun in 1886 in other parts of the field and by 1890 about 100,000 tons of ore had been raised. In the Sudbury field, as in Norway and at the old Lancaster Gap mine in Pennsylvania, the ore was first mined for its copper content and difficulty in the refining of the copper led to the discovery of the important proportions of nickel in the ore. No. 3 property of the Inter-national Nickel Co., now known as the Frood, the greatest deposit in the

Jan 1, 1932

By A. B. Cook, C. J. Walker, G. B. Spencer

Although water will displace oil from a petroleum reservoir to a greater extent than gas will, there are some reservoirs in which gas rather than water should be used for pressure maintenance. This is indicated by the high-percentage oil recovery in the case history reports on the Pickton field1 in East Texas and the Raleigh field2 in Mississippi. In gas cycling, total oil recovery includes vaporized products from the immobile oil in addition to oil produced by displacement. If a large part of the immobile oil is vaporized, total oil recovery may be higher than that obtainable from pressure maintenance by water injection. However, the amount of immobile oil vaporized may range from almost 0 to 100 percent, and no simple and reliable method has been presented in the literature for calculating oil vaporization. Development of such a method was the purpose of this study. Several papers concerning the calculation of oil vaporization have been published,3 6 and each is based on the concept dealing with K values — equilibrium constants for the various components in reservoir oils. (K for a component in a system of vapor and liquid in equilibrium is the mole fraction in the vapor phase divided by the mole fraction in the liquid phase.) This appears to be the most logical approach. Yet, the problem is too great for a perfect solution, even with modem laboratories to analyze oils and high-speed digital computers to perform calculations; the number of components in reservoir oils is too large. Because of their complexity no two reservoir oils are exactly alike. Also, K values change not only with variations in pressure and temperature but also with composition of the reservoir oil. In addition, during gas cycling, the lighter hydrocarbons tend to vaporize first. Thus, the reservoir oil becomes more dense and less volatile as gas cycling continues. Furthermore, the greatest amount of oil vaporization occurs near the injection well. Therefore, a simplified method is required for calculating oil vaporization because a rigorous method is not practical. The present methods of calculating vaporization generally require knowledge of the reservoir oil composition in which the mole fractions of the lighter components through hexanes are given; the remaining oil is described as C, + (heptanes and others of greater molecular weights). Assigning one K value for the C,+ system provides the simplest solution, but this incorrectly assumes that K for C, + does not change as gas cycling continues. This assumption can cause large errors because the K value may become less than a thousandth of the original after a large amount of oil vaporization. Another way to solve the problem is to determine changing K values for C, + according to how much of the C, + has been vaporized. To determine these values a sample of the reservoir oil is injected into a pressure-volume cell, and dry gas is batched in and out of the cell. Appropriate measurements and tests are made to determine K

Jan 1, 1970

By H. R. Cooper

Beginning from the point of view that process control of ore flotation can be approached by on-line analysis of processing behavior, an equation was developed which has characteristics useful for a multi-variable feedback loop. The application of this equation to optimizing control with a data acquisition system will be discussed in Part 11. OBJECTIVE OF CONTROL The purpose of developing a feedback method for concentrator control is to provide a means for determining ore processing behavior in order that data describing the ore can be used to make appropriate changes in control conditions while the ore continues to flow through the flotation circuit. In absence of simple linear relationships between pairs of variables in the flotation process, attention is directed to the need for a generalized relationship which possibly would incorporate all pertinent variables. It would be possible to apply such an equation to the control of variables in a concentrator in real time according to operational requirements as ore gradually changes in quality and processing behavior. To meet this objective, it is necessary to define relationships describing flotation performance in terms of performance variables capable of being measured on-line, either with on-line instruments and analyzers or with manual methods lending themselves to quick determination of meaningful information. Relationships derived to meet this objective form a model for the flotation process. The evaluation of flotation parameters of the ore and the execution of control on the basis of the derived information while the ore is continuing through the process would constitute a feedback control system. Data for the following control and performance variables can be made available: 1) ore feed rate, 2) collector addition rate, 3) lime addition rate, 4) pulp density, 5) assay or multiple assays in feed, tails, or concentrate sufficient for material balancing, 6) particle size distribution in flotation feed, 7) froth height, and 8) frother addition rate. Other process variables, such as flotation air input and mechanical agitation intensity, are usually not considered as control variables and were not taken into account in this analysis. Other data, such as that describing the ore, would be useful for concentrator process control. However, information of this type (such as ore hardness, grain size of the mineral, chemical species of the mineral, and ore surface properties) being not presently measurable in the on-line sense, has the effect of causing disturbances in the system which can be detected only indirectly. The process control system being sought must recognize disturbances by their effect on measured performance followed by determining control changes to compensate for them in meeting the objectives of improved performance (perhaps higher throughput, maximum return, etc.). A suitable model for the process having capability of recognizing the effects of ore disturbances should provide the means for introducing proper compensation through control changes. This would be done by processing on-line data with the model to calculate coefficients reflecting ore characteristics. This procedure can only be applied in a typical flotation system when on-line data can be processed at faster rates than the usual rate at which the ore in flotation feed takes on significantly different characteristics. For example, material balance data can be accumulated in approximately 5-min intervals using on-line X-ray fluorescence analysis. The flotation system time constant might be 15 min including grinding. Ore changes having time constants longer than 1 hr could be accommodated using on-line feedback control by this procedure. FORM OF A FLOTATION MODEL A suitable form for mathematical expression of the flotation process is an equation involving rates, when each of the components entering into the separation can be treated independently by a rate equation. A first-order, kinetic rate equation can be written on

Jan 1, 1967

By P. K. Leung, S. R. Allada, P. M. Dranchuk, D. Quon

The alternating direction explicit procedure (ADEP) makes use of the boundary conditions to reduce multi-dimensional problems to a series of one-dimensional problems. The method, previously applied to reservoirs containing only an under-saturated oil, has now been extended to cover the case of natural gas reservoirs. l his involves solving a non-linear partial differential equation, application of the procedure is straightforward and no calcuIational problems were encountered. INTRODUCTION There has been a growing interest in formulating mathematical models of petroleum and natural gas reservoirs — models which permit the engineer to examine and evaluate the physical and economic consequences of various alternative production policies. The tremendous reduction in the cost of solving such models in recent years has made possible their use as an almost routine management tool. This reduction is the result not only of improved computer hardware but also of the development of more efficient mathematical techniques. The present paper is concerned with the application of a recently proposed numerical method (ADEP) to two-dimensional gas reservoirs. STATEMENT OF THE PROBLEM Given a two-dimensional Region R, bounded by a closed Curve C (Fig. 1) such that the behavior on Curve C is known, the differential mass balance for each fluid phase in R, neglecting gravitational effects and assuming Darcy's law for fluid transport, can be written as:

Jan 1, 1967

By C. F. Weinaug, R. W. Farley, J. F. Wolfe

A rapid, accurate method for predicting the dew points of gas condensate systems and their subsequent normal and retrograde phase behavior with pressure decline has been developed. The method predicts dewpoint pressure, volume percent liquid, condensate content of produced gas, and the liquid and gas compositions at pressures below the dewpoint pressure. Use of a modified form of the BWR equation of state, which was established by regression calculations on experimental data from several gas condensates, provides depletion behavior predictions that are significantly improved over predictions using K-values based on convergence pressures and Standing's correlation for density. INTRODUCTION Prediction of complex hydrocarbon-mixture thermodynamics is fundamental to reservoir and gas process engineering. In particular, accurate descriptions of phase behavior are needed to plan production operations for gas-condensate reservoirs. Although analysis of reservoir behavior during producing and cycling operations can be made for simple reservoir models, extensions of these studies to more complex reservoir systems are still difficult, both with regard to numerical techniques and the phase behavior of the fluids. This report presents a mathematical method for handling the complex phase relationships encountered in retrograde condensates; it represents an advance in reservoir behavior technology in that it will aid the development of closed-form mathematical representations of total production systems. A typical phase diagram for a complex hydrocarbon mixture is shown in Fig. 1. When the reservoir temperature T, is higher than the critical temperature Tc, and lower than the critical-condensation temperature Tct' the mixture is a gas condensate. Decreasing the pressure on a condensate while maintaining nearly constant temperature will result in two phases forming at the dew point. Further decrease in pressure will cause additional liquids to be formed. This phenomenon is called retrograde condensation. As pressure continues to decrease, normal phase behavior regions are encountered (liquid vaporizes as the pressure decreases) until the low-pressure dew point is encountered. When a gas-condensate reservoir is produced, gas is withdrawn and the volume of fluids in the reservoir remains nearly constant as the pressure declines. Heavier components tend to accumulate in the liquid phase as a result of retrograde condensation, while lighter components tend to concentrate in the vapor phase. Thus, as the gas is produced the composition of the gas and the condensate

Jan 1, 1970

Discussion, at the Virginia Beach Meeting, February, 1894, of the Paper of Prof. Posepny. (Trans., xxiii., 197, 587.) Including communications subsequently received. a T. A. Rickard, Denver, Colorado (communication to the Secretary) : The paper of Professor Posepny was printed so short a time

Jan 1, 1895

By H. J. de Haan, J. van Lookeren

In view of their size and favorable characteristics, Shell's heavy-oil fields on the eastern coast of Lake Maracaibo, Venezuela — known as the "Bolivar Coast" — are obvious candidates for the application of thermal recovery methods. During the last decade various methods have been tested, such as steam drive,', steam soak and in-situ combustion.3 Steam soak was discovered as a promising production method rather accidentally in 1969, during early steam drive testing in the Mene Grande Tar Sands1 When steam erupted at the surface due to breakdown of the overburden, the injection wells were back-flowed to relieve the reservoir pressure. This resulted in high oil production rates, all the more impressive because the reservoir is unproducible by primary means. It was concluded that injection of limited amounts of steam might be a very effective method for stimulation of heavy-oil wells. To obtain information that would be more widely applicable, testing of the process was continued in the eastern part of the Tia Juana field (Fig. 1). Here even more favorable results were obtained than in Mene Grande. This created sufficient confidence to prompt a large-scale project, which was considered essential for an adequate evaluation of ultimate recovery, steam requirements, mechanical problems and economic prospects. To obtain an early answer on the ultimate recovery a fairly depleted area was selected in the western part of the same field ("A" in Fig. 1). The first steam was injected in March, 1964. This paper describes the performance until the end of 1966. The production mechanism and well inflow performance are discussed in some detail. In view of the favorable results obtained so far, steam soak operations have expanded considerably since the start of this project. Presently, five other projects are under way, of which two are about the same size as the one under discussion, with further extensions being programmed. Description of the Project Reservoir Data Shell's part of the Main Tia Juana reservoir has a proven area of 8,800 acres and an average net oil sand thickness of 130 ft. It contained an initial quantity of 2,750 million bbl of stock tank oil (STOIIP), with an API gravity ranging from 10 to 15 API and a viscosity of 100 to 10,000 cp (at 122F). As a geological unit it is known locally as the Lower Lagun-illas member of the Lagunillas formation, which is early Miocene in age. The sand bodies are believed to be distributary channel sands deposited in a lower coastal plain. Most of the channel sands have lateral and some vertical communication. The field was developed mainly in the period from 1936 to 1950. The ultimate primary recovery is estimated at 18.1 percent of STOIIP, of which 12.5 percent STOIIP had been produced by the end of 1963 from some 540 wells. By then the average pressure had decreased from an initial value of 960 to 550

Jan 1, 1970

By T. D. Parris

Within a 30-mile radius of Sudbury, Ontario, the Ontario division of the International Nickel Co. of Canada, Ltd., operates nine underground mines and two open pits. Prior to 1966, ore removal from the underground workings was mainly by scraping, using 125-hp electric slushers in blasthole and caving mining, and 30-hp or less air-operated slushers for fill- method mining. The development and evaluation of low-profile, diesel-powered front-end loaders, designed specifically for underground operations, was followed with keen interest during the period 1964-1965.

Jan 6, 1969

By G. E. Perry, R. D. Seba

The steam soak process is the most widely applied and most successful thermal supplemental recovery process in use today. This process, which consists of injection of steam in various quantities into a well that is subsequently used as a producer, can be repeatedly applied to the same well. The treatment, if successful, will increase the production rate of the well for a given period. For this reason, the steam soak process should be considered primarily as an acceleration method rather than as a method for producing oil that could not be recovered by primary means. However, application of this process to reservoirs containing viscous crude may produce oil that otherwise would have zero present value or would be produced at such low rates by primary methods that it would not be profitable. A major necessity for the successful application of the steam soak process is a source of natural energy. This source of energy may exist in one or more of the following forms: expansion of fluids, compaction, gravity drainage. This study was undertaken to determine the effect of gravity drainage on steam soak production and recovery and to develop a method for predicting production rates, ultimate recovery, and heat requirements when applying this process to a reservoir in which the primary driving mechanism is gravity drainage. (For additional discussion of some of the factors that affect recovery of viscous oils by the steam soak process, see Ref. 5.) Model The model developed for analyzing the performance of the steam soak process in gravity drainage reservoirs is a modified form of that developed by Matthews and Lefkovits.1,2 The model used in this study consists of two concentric cylindrical volumes of productive formation with cap and base rock of impermeable material. The inner cylinder, extending from the radius of the wellbore to some predetermined radius, is referred to as the hot zone; and the outer cylinder extending from the exterior boundary of the hot zone to the outer drainage boundary will be referred to as the partially heated zone. Fig. 1 is a schematic cross-section of the model, indicating its thermal properties. Horizontally the reservoir is divided into two regions, the upper region having a high gas saturation and residual liquid saturation, and the lower region containing a high liquid saturation. The boundary between these two regions is referred to as the free surface, which drops and changes shape as fluid is produced from the reservoir. The model was used to analyze the steam soak process by applying the following assumptions: 1. The horizontal formation is composed of one or more homogeneous noninterconnected layers containing incompressible fluids. Each layer is treated separately. 2. The hot zone is heated by the initial soak and is maintained at the steam temperature corresponding to atmospheric pressure by resoaking. 3. The height of the free surface at the producing well is at all times zero (producing well is pumped off). 4. The position of the free surface at the start of

Jan 1, 1970

By R. L. Fullman, J. W. Cahn

A method is derived for obtaining the distribution of sphere diameters and plate thicknesses (or pearlite spacings) from size distribution functions obtained along randomly oriented lines. The method described offers several advantages over existing methods. IT is often of interest to measure not only the average particle size but also the distribution in sizes of particles dispersed in an opaque material, such as an alloy. Particles of a single size will exhibit a range of apparent sizes on the plane of polish. If the sample contains a distribution of particle sizes, the distribution of observed sizes on the plane of polish is related to the distribution of actual particle sizes by an integral equation, For example, a distribution of sphere sizes, N(r),

Jan 1, 1957

By L. E. Lynd

Ilmenite concentrates prepared from beach and dune sands of Florida, India, Australia, and other localities deviate to varying degrees from the theoretical composition of pure ilmenite, and generally contain more TiO2 than the formula FeTi03 would indicate, and some Fe2O3. This variation in composition is due to different amounts of altered ilmenite in the concentrates. (') The alteration has been shown to take place by oxidation and leaching of the iron in ilmenite, to yield residual products consisting of amaphous or very finely crystalline titanium dioxide and ferric oxide. Other studies(2) have stressed that ilmenite may alter first to an amorphous titanium-iron oxide, or to an amorphous mixture of titanium and iron oxides, with later alteration to crystalline titanium dioxide. There has been some speculation as to whether alteration of ilmenite is usually caused by weathering or by hydrothermal pr0cesses.(3,4,5 Three main methods of approach were used to study the mechanism and rate of ilmenite weathering. (1) Etch tests were used to learn the relationship between crystallographic orientation and the rate of chemical attack on ilmenite. (2) Altered ilmenites in New Jersey sands were compared with ilmenite in rocks of the area from which the sands were probably derived. (3) Experiments were carried out to determine at what rate, and by what agents, originally fresh ilmenite could be altered in the laboratory under simulated weathering conditions. ETCH FIGURES ON ILMENITE CRYSTAL PLANES A rather large ilmenite crystal from Sannidal, Norway, about 11/2 x 11/2 x 2 inches, bearing a few natural faces with imperfect surfaces, was prepared for etching by grinding and polishing crystal planes whose Bravais-Miller indices (referred to hexagonal axes) had been determined by X-ray methods to be 00.1, 10.0, and 10.2. Metallographic examination showed that the crystal was homogeneous, and uniform in orientation, except for the presence of twin lamellae which were oriented with their broadest surfaces parallel to the 00.1 planes of the large crystal. Since ilmenite is rather inert to standard etch reagents, strong solutions of hot sulfuric acid and cold hydrofluoric acid were used. Preliminary tests established that better etch figures were obtained using HF than with H2So4 solutions. The etch figures produced on 00.1 and 01.0 planes of ilmenite by etching 3 to 5 minutes in cold concentrated HF solution are shown in Figures 1, 2, and 3. The figures on 00.1 are "etch hillocks", rather than pits, and according to Honess(6), this phenomenon

Jan 1, 1961

By L. H. Simons, H. H. Spain

I NT RODUCTION Theory indicates that linear water-floods should exhibit scaling and stabilization properties in both oil-wet and water-wet porous media'. Experimental verification of these proper-ties in oil-wet media was obtained some time ago', and more recently, Perkins', employing a high permeability unconsoli dated sand and a low oil-to-water viscosity ratio, has obtained flooding data in a water-wet medium which follow the pattern of scaling and stabilization. Also, Root and Calhoun3 have found similar trends in the displacement of gas by liquids from unconsolidated sands. In general, however, the experimen-tal evidence appears still somewhat conflicting since scaling and stabilization of water floods in water-wet media has not been observed in most of the studies to date'.".'. The purpose of this paper is to provide a more comprehensive picture of waterflood behavior in water-wet media and to corroborate the validity of the scaling and stabiliza-tion concepts. The need for differentiating between the intrinsic nature of water-oil displacements inside a porous medium and perturbating secondary end effects is discussed, pertinent experimental procedures are outlined and the results of flooding tests conducted at low and at high oil-to-water viscosity ratios on con- solidated, water-wet Alundum cores are analyzed. END EFFECTS IN WATER-WET MEDIA Outlet End Effects The peculiar feature of the outlet end effect in water-wet cores is that it not only results in an excessive wetting phase saturation at the outflow face,'.' but that it also retards the moment of water breakthrough: i.e., causes an undue delay in the appearance of water in the produced effluent. Conditions existing in a water-wet core when water first arrives at the outflow face are schematically shown in Fig. 1. At this flooding stage the saturation distribution, Fig. la, has not yet been distorted by the outlet end effect. A water-invaded pore at the outlet core face, opening into the oil-filled void space between the end of the core and the end plate of the flooding cell, is illustrated in Fig. lb. The curvature of the water-oil interface is concave toward the outlet, and the pressure on the water side of the interface is lower than on the oil side by an amount equal to a certain capillary pressure (Fig. lc). Before the water (of this particular pore) can be produced from the core, it is apparent that the pressure in the water phase must exceed that in the oil outside the core. Accordingly. when water first reaches the outlet face, it will not leave the core but will accumulate inside, and only oil will be produced. Eventually, as the water injection is continued. sufficient pressure is built up in the water to reverse the curvature of the interfaces in the outlet pores, and only then will water be produced from the core. This stage, illustrated in Fig. 2. corresponds to the conventional "water breakthrough" observed in flooding experiments. For purposes of discussion, it is convenient to distinguish between the flooding stage at which water first arrives at the outlet face of a core (Fig. 1) and the stage at which water is first produced from the core ( Fig. 2). The former will be referred to as "water arrival" and the latter water breakthrough. In oil-wet media

By E. H. T. Whitten

Fundamental principles of trend surface analysis are briefly reviewed. Recent analyses of several granite masses illustrate the method of estimating the quantitative regional composition and also variation trends. Mineralogi-cal, bulk chemical, and other quantitative variates can be employed. Such results should improve gravimetric surveys (etc.) relying on correction terms based upon estimated rock composition. For the first time realistic estimates of the mean composition of rock bodies can be made with the aid of trend components. Deviation maps hold promise for delimitation of disseminated mineralized areas. Deviations comprise residual "anomalies" after removal of trend components. In petrogenetic work such deviations have been shown to reflect local influences (e.g., contamination of magma by assimilation of xenoliths or country rocks of the envelope). Similarly, mineralization should be reflected as deviations standing out from background trends. Actual examples are given. Present results largely involved projection of data onto a horizontal plane, but if sufficient data are available use of three spatial coordinates increases accuracy for prediction in "unexposed" areas. Regression analysis shows that statistical correlations between variates for regionally distributed specimens are often strong (particularly with percentage data). Hence, in preliminary surveys more easily determined variates may sometimes be employed to gain first estimates of the distribution of components more economically significant. The sampling plan is of paramount importance. Unless a suitable sampling plan is devised and adhered to rigorously during collection of the samples, all subsequent analysis is liable to strong bias and to yield spurious estimates of the truth. For a wide variety of purposes it is becoming necessary to assess the quantitative composition and variability of rock masses and mineralized zones of various types. In petrology and petrography this is something of a break with tradition, because virtually all descriptions heretofore have been qualitative. Although vast amounts of qualitative data have been compiled, little is known about the quantitative behavior of rock components. However, if variability of rocks could be understood and expressed with sufficient accuracy, there would be considerable advantage to both economic and petrogenetic problems. For example, if, on the basis of a limited series of observations, the three-dimensional geometry of a disseminated orebody (or ore-bearing formation) could be predicted accurately, it would be possible to extrapolate with confidence the structure into unexposed areas. The following example based upon the geometry of an igneous cone sheet illustrates the point. The significant rock-unit gives outcrops at a number of points on the exposed surface ABCD (Fig. 1A). Study of n exposures may enable a geologist to predict that the rock-unit is circular in form and thus indicate where the unit can be expected in unexposed terrain. Such observations are made at sites with u, v coordinates, but, if subsurface information is available at sites with uvw coordinates, it will eventually become possible to predict accurately the three-dimensional form of the rock mass. Geologic rock masses or mineral deposits are not generally as regular and symmetrical in geometric form as the body illustrated in Fig. 1A. However, even where the broad limits of a major rock mass are known, a more acute problem arises in attempting to obtain detailed quantitative estimates of its internal variability; it is difficult to predict the variability of a mass from a limited number of observa-

Jan 1, 1961

By H. C. H. Darley

The principal causes of unstable boreholes* have been known for many years. For example, in a paper published in 1938, Halbouty and Kaldenbachl listed nearly all the classes of troublesome shales that we know today. Much progress has been made since then, but we still have difficulty in identifying the cause of hole problems and consequently in selecting the right technique for combatting them. This subject was discussed at length by Kelly2 in a recent paper and much valuable information was provided. In the Shell Laboratory it was felt that there was a need for a better understanding of the basic mechanisms underlying the various forms of hole instability, and accordingly, a study was undertaken. This entailed reviewing the physics and chemistry of shales, developing tests for characterizing shales, and constructing a model to study borehole stability. Full simulation of underground conditions was not attempted; for example, we used reconstituted shale specimens, and often simple solutions or suspensions rather than practical drilling muds. Such procedures are necessary in order to establish mechanisms, but caution must be exercised in applying the results to the more complex conditions in the field. Part I — Compaction Swelling and Dispersion of Shales Compaction of Sediments Shales are formed by the compaction of sediments. Water is expressed as sediments are buried by subsequent layers, and the degree of compaction is proportional to the depth of burial, provided the water is able to escape easily to permeable formations. The younger sediments, when recovered at the surface, soften and disperse when mixed with water; the older shales have undergone diagenesis (alteration of the clay minerals, secondary cementation, etc.), and remain hard and unaffected by water. In this paper we shall apply the term "shale", as the drilling engineer does, to everything from clays, which are highly reactive to water, to completely lithified materials such as claystones and slates, which are completely inert. However, because these materials behave quite differently when encountered in the drilling well, it is desirable to develop a simple means of characterizing them. The feature that distinguishes the various shales is their dispersibility in water —soft clays disperse readily, harder shales disperse slowly when agitated, the lithified materials will not disperse at all unless milled. In the dispersibility test used in our study, a weighed amount of dry shale particles in the 12-to-100-mesh-size range was placed in a tube filled with distilled water and allowed to soak for 16 hours. The

Jan 1, 1970

By C. G. Cheriton

The Bathurst mining district is located near the town of Bathurst on the north shore of the province of New Brunswick in eastern Canada. This region is the northeasterly extension of the Appalachian mountain system. The area is heavily forested and covered with a fairly uniform layer of glacial gravels and boulder clay. The relief is low, increasing to a few hundred feet in the north west and several hundred feet in the extreme southwest. Thus, rock outcrops are scarce and occur mostly along streams and ridges. These factors made conventional prospecting difficult and resulted in the district becoming an excellent place for the application of regional geological studies, air photo techniques, airborne and gound geophysical methods and soil sampling. Late in 1952, The Brunswick Mining and Smelting Company discovered a very large base metal sulphide body presently known as their No. 6 orebody. Early in 1953, another and even larger sulphide body was discovered and is known as the No. 12 ore body. Mining exploration was very active thoughout 1953 and in 1954, The American Metal Company discovered the orebodies which have been developed into The Heath Steele 'Mines. The location of these deposits is shown on figure one. In the latter part of 1954, The Anaconda Company became interested in the district and decided to make a geological study of this new and promising base rrietal area. The following paper presents the story of the application of modern geological and geophysical techniques, leading to the discovery of a large sul-nhidp liodv. The Bathurst district is generally underlain by the Tetagouche group of Middle Ordovician rocks.

Jan 1, 1961

By H. S. Price, J. C. Cavendish, R. S. Varga

A numerical formulation of bigh-order accuracy, based on variational methods, is proposed for the solution of multidimensional diffusion-convection-type equations. Accurate solutions are obtained without the difficulties that standard finite difference approximations present. In addition, tests show that accurate solutions of a one-dimensional problem can be obtained in the neighborhood of a sharp front without the need for a large number of calculations for the entire region of interest. Results using these variational methods are compared with several standard finite difference approximations and with a technique based on the method of characteristics. The variational methods are shown to yield higher accuracies in less computer time. Finally, it is indicated how one can use these attractive features of the variational methods for solving miscible displacement problems in two dimensions. INTRODUCTION The problem of finding suitable numerical approximations for equations describing the transport of heat or mass by diffusion and convection simultaneously has been of interest for some time. Equations of this type, which will be called diffusionconvection equations, arise in describing many diverse physical processes. Of particular interest here is the equation describing the process by which one miscible liquid displaces another liquid in a one-dimensional porous medium. The behavior of such a system is described by the following parabolic partial differential equation:

Jan 1, 1969

By George A. Roberts, Arthur H. Grobe

H. H. Hausner (Sylvania Electric Products Inc., Bayside, N. Y.)—I tested the 18-8 stainless steel powder described by Grobe and Roberts and the results were excellent. The powder was compacted and sintered in purified hydrogen to a very low density (approximately 20 pct and more porosity). Where the other powders show very little ductility when sintered to such a low density, the described stainless steel was surprisingly high in ductility. I would appreciate it if the authors could give us some data on a similar powder but without silicon and also data on 430 stainless steel powder. H. S. Kalish (Sylvania Electric Products lnc., Bay-side, N. Y.)—I should like to verify the excellent results obtained by Grobe and Roberts with their pre-alloyed stainless steel powder. It seems to me, however, that their data are conservative and that even better densities and higher tensile strength can be obtained with these powders by the selection of proper particle size fractions. The accompanying data were obtained by sintering compacts pressed from the powder made by Vanadium-Alloys Steel Co. to 0.8 in. in diam about % in. high in a cylindrical double action die without the use of lubricant. The compacts were sintered for 1 hr at 2370°F in hydrogen purified by passing it through a Deoxo unit (palladium catalyst) and a Lectrodryer. Even the density obtained using the as-received (—100 mesh) compacted at 40 tsi was quite high as was the hardness. This indicates that by increasing the sintering temperature slightly above 2350°F the coining and subsequent secondary sintering operation can be avoided. With the — 325 mesh fraction, however, very much higher density can be obtained and a high hardness. I have no tensile data to present, but I feel that it is not too optimistic to expect. that unusually high tensile strengths and good ductility can be obtained in this sintered material on the basis of the density and hardness data. The sintered stainless steel made by the above-described methods came out of the furnace very bright and clean. The purified hydrogen atmosphere is important in sintering stainless steel to reduce the oxide present in the powder. The result of such a sintering is good densification and good ductility of the sintered product. G. A. Roberts (authors' reply)—We are greatly indebted to Messrs. Hausner and Kalish for their discussions of our paper. Their encouraging remarks are greatly appreciated and their confirmatory data will do much to add to the value of the paper. At the present time it is impossible to provide complete data on a similar stainless steel powder without silicon, but we do know that, in general, the properties are slightly inferior. Data on such a powder and on other stainless steel types are to be published in the near future.

Jan 1, 1952