Search Documents

By Werner Goldsmith

This discussion is confined to the first section of Professor Clark's paper entitled 'Elastic and Nonelastic Waves' and its application to wave propagation in rocks. Some published results by other investigators concerning measured decay of peak strain and variation of pulse duration with travel distance are cited, but the principal subject of this section is a detailed mathematical analysis of spherical (and, to a lesser extent, of plane) wave propagation in various homogeneous, isotropic media. The materials selected include an elastic substance, a Voigt and a 'modified' Voigt solid-the latter incorporating the assumption of a constant loss factor (defined here as the product of frequency and the ratio of the imaginary to the real part of the shear modulus)-and a solid friction model. While other possible representations of a solid continuum model are mentioned, corresponding solutions are not presented. The pulse inputs considered comprise exponential decay, a unit step, and steady- state sinusoidal excitation: the process of conversion from one solution to another in the operational plane is indicated. The basic problem relative to the analysis of wave propagation in rocks concerns the specification of the appropriate material behavior, since the equations of motion for a given geometry of the source and medium are well established. In principle, the required information can be deduced from two types of measurement: 1) field observations of underground explosions or seismic disturbances, and 2) laboratory tests on small samples of the pertinent rock. However, in practice both types of approaches are beset with difficulties and subject to valid criticism. Field observations must be carried out under relatively uncontrolled conditions, including, in general, a lack of knowledge of the precise arrangement of the medium, stratification, faulting. inhomogeneities, etc. Specimens investigated in the laboratory can not represent the environment found

Jan 1, 1967

The Ringelmann Smoke Chart fulfills an important need in smoke abatement work and in certain problems in the combustion of fuels. A knowledge of its history and method of preparation is, therefore, of interest to many. Since instructions on its use are not shown on the recent edition of the chart, those included in this revision of the previous Bureau of Mines publication now are a necessary complement to the chart, More detail regarding the use of the chart is included than was given in the earlier version.

Jan 1, 1967

By R. R. French

Investigations were conducted by the lndiana Geological Survey for some dry methods of bene-ficiating low-grade gypsum ore. Seventy-two batch and continuous flow tests were performed with a roller mill, rod mill, pebble mill, electronic color sorter, electrostatic separator, and an air separator. Approximately 650 size analyses and 550 chemical analyses were performed during the investigation. Batch samples were tested by the Survey, and most of the continuous flow tests were handled by com-mercial laboratories. INTRODUCTION Although many companies producing gypsum have been interested in dry beneficiation, very little information has been published in the past. One of the most comprehensive investigations was reported by the Canadian Department of Mines.5 Other reports of interest have been published by the Canadian Department of Mines,' the South Australian Department of Mines,2 and the U.S. Bureau of Mines.3 Some beneficiating practices, such as hand picking, selective mining, crushing and screening, milling, and air classifying, were already in use in the Canadian industry prior to MacPherson's investigations.5 MacPherson examined the effect of tabling, air separation, screening, flotation, electrostatic separation, calcining, washing, and various types of milling, in his efforts to eliminate dolomite, limestone, and small amounts of silica and clay from the gypsum ore. RAW MATERIAL Gypsum and anhydrite occur within the lower part of the St. Louis Limestone (Mississippian) in southwestern Indiana. Gypsum in single beds 10 or more ft thick occur at Shoals, Martin County; near Freedom, Owen County; near Bloomfield, Greene County; and in the Devonian strata of LaPorte County. Only the Shoals deposit has been exploited commercially. The western edge of the Shoals deposit is contaminated by a continuous bed of shale, about 1 1/2 to 2 ft thick, near the top of the evaporite and by thin continuous strata or irregular masses of carbonate rock. The low-grade material used in these investigations was obtained from the waste pile of the National Gypsum Co.'s plant at Shoals. The waste ore averaged about 67.4% gypsum and was contaminated by various amounts of shale, dolomite, and limestone. All the material had been previously crushed and screened to minus 1 1/4 to plus 3/8 in. X-ray analyses of powdered, sedimentated, gly-colated, and heat-treated samples of the shale showed that it was composed of slightly structurally disordered illite, Fe-rich chlorite, and very finegrained disseminated silt. The clays and silt were partly cemented with carbonate material. Light-gray argillaceous limestone and gray or brown porous dolomite made up most of the carbonate rock contamination. Crushing Characteristics: The rate unweathered gypsum, carbonate rock, and shale reduce in size was determined by crushing handpicked samples in a rod mill and screening at minus 100 mesh (.0059 in.). Fig. 1 graphically represents the data. These data indicate that carbonate rock should be relatively easy to separate from gypsum, but that shale should be difficult to separate from gypsum by differential crushing. Controlled Samples: Eight rod mill tests of hand-picked samples of mixed gypsum and shale and mixed gypsum and carbonate rock were made. Sample weight and crushing weight were standardized, but the scalping screen size was varied from 3 to 10 mesh (.265 to .0787 in.) in order to obtain the optimum purity and recovery. The data obtained from the tests (see Table I) illustrate that the relatively hard carbonate rock allows appreciable beneficiation (13.6 to 18%), but the low strength gypsum and shale allows only minimal separation (3.1 to 7.5%).* Semi control led Samples: During the investigation considerable slaking of the shale was noted when the ore was subjected to alternate wetting and drying. As demonstrated in laboratory tests, as much as 76% shale can be removed by two cycles of wetting and drying in a 48-hour period with subsequent screening at .265 in. Five ore samples, stockpiled and weathered, were screened, rod-mill crushed and rescreened at .157 or .187 in. with excellent results (Table 11).

Jan 1, 1967

By Richard L. Ash

Usually stresses generated by explosive charges during blasting are considered to propagate through materials in a spherical manner. Although conical stress-waves had been suggested by a few investigators they have never been clearly defined or critically observed. An experimental technique using models and lengths of Mild Detonating Fuse was developed to produce conical stress-waves and study their effects. Several materials with different properties were tested to compare the influence of those properties on wave shape and crater formation. The investigations showed that whenever column charges were end-initiated, conical stress-forms were produced if the explosive's velocity exceeded the energy propagation velocity of the surrounding material. The stress-cone angle significantly affected fragment size and shape as well as the geometry of craters produced.

Jan 1, 1967

By Mark E. Emerson

1967 will be remembered as the year the world realized it was running out of silver. In May of this year a new era in tin the history of the white metal was begun as the power of the U.S. Treasury to regulate world .silver markets was staggered by one of the most intense speculative flurries in recent times. In the past few months the price of silver has risen from $1, 29 to over $1.80 per troy ounce and with it have spiraled the prices of a host of silver-bearing consumer items. Legislation has, been enacted to conserve our coinage and meager silver stocks and today large quantities of non-existent silver are trading hands at a fast pace on the commodity exchanges throughout the world. Now that the price of silver has been set loose, forces are in motion to dull the world's appetite for silver. At the same time there is strong incentive to provide greater quantities of silver to satisfy the world's hunger' for the metal.

Jan 1, 1967

By E. F. Foerster

The Bureau of Mines used a diphenyl ether calorimeter to determine the heat contents of natural and synthetic blast-furnace slags. For solid slags, enthalpy (HT-H298) at 1,448° to 1,680° K ranged from 282.1 to 348.5 cal/g (±3 percent); liquid slag values were 400.1 to 514 cal/g (±3 percent) in the temperature range between 1,626° and 1,890° K. Experimental results for liquid slags were 14.3 to 16.2 percent higher than the enthalpies calculated from thermodynamic tables; values for solid slags ranged from 5.3 percent lower to 3.9 percent higher than calculated values. Experimental results for liquid synthetic slags ranged from 14.9 to 24.3 percent higher than the calculated enthalpies; values for solid synthetic slags were from 1.0 to 3.4 percent higher.

Jan 1, 1967

By T. C. Crosby

"Grinding has a marked effect on concentrator recovery and operating cost, and improved operation through automatic control can enhance plant performance. The correct application of automatic control can only be accomplished after the basic process variables and their effect on mill operation have been defined. This paper develops a control system for a closed circuit grinding system. Typical control systems, such as those for feed, water flow, level and circulating load, are discussed. A method of determining controllability, using the calculated reaction curve technique, is outlined. Various successful automatic grinding control applications are discussed and future possibilities are outlined. Although this paper primarily deals with closed wet grinding• circuits, the principles involved can be readily applied to other processes."

Jan 1, 1967

By Norman F. Schulz

Representative sample s of typical manganiferous iron ores from the Cuyuna Range, Minnesota, were reduced to the metallic iron state at temperatures just short of fusion with an excess of coke in the manner prescribed by the R-N Direct Reduction Process. The magnetic separation characteristics of the reduced ores were determined at several grinds by means of the Davis magnetic tube concentrator. The efficiency of metallic iron concentration by magnetic separation was a direct function of the degree of metallization of the ore sample and of the degree of liberation attained in grinding. Some improvement of liberation characteristics was observed if a small amount of sodium chloride was present in the ore during reduction. Despite the extremely fine -grained structure of the crude ores, iron concentrates containing more than 90% Fe at iron unit recoveries of 90% or more were obtained from well-metallized samples of three of the four ores tested. The rejection of manganese to the tailings corresponding to such separations varied from 80 to 95%.

Jan 1, 1967

By Forrest T. Moyer

General injury experience in the mineral fuel industries was less favorable in 1965. Of the major industry groups, the coal mining and the coke industries had higher injury-frequency and severity rates than in 1964. The oil and gas industries had an increased Tate of occurrence of injuries although the severity rate improved. However, in the relatively small peat and native asphalt industries, the injury experience was more favorable than in 1964. Overall operating activity for the mineral fuel industries, as measured by the average number of men working daily and the total worktime, was higher than in 1964. Moderate declines in employment and worktime in the coal and peat industries were more than offset by gain in the other groups. In addition to the industry classifications included in this chapter, similar injury, employment, and worktime data on the metal, nonmetal, stone and aggregate industries are presented in volume I. Corresponding data for broad classifications of mineral industry groups arc given by States in volume III of the Minerals Yearbook.

Jan 1, 1967

By Bernard G. Carbajal

The dielectric properties and thermal stability of glow-discharge polytnerized films of styrene, chloro-benzene, and other organics and photoresist filrns are presented. Variations in the glow-discharge polymerization procedure are related to thermal stability. Thermal treatwzents to improve the dielectric properties of both positive (A2 340) and negative (KMER) photoresists are given. THIS study of the dielectric properties of thin organic polymer films is part of a program to investigate the electronic applications of organic materials. Thin organic polymer films have been produced in an ac glow discharge,13 with a low-energy electron-beam process,4 and in an rf dicharge. A process has been developed during the course of this program for forming thin organic polymer films in a dc glow discharge. These glow-discharge polymerized (GDP) films exhibit some useful dielectric and insulative properties. Variations in the properties of the GDP films can be induced by changes in the source material (monomer) andor in the conditions of the glow discharge. Thin organic polymer films can also be produced from typical positive and negative photoresists. Some variations in the properties of the photoresist films can be induced by suitable post-application treatments such as thermal crosslinking. The properties of interest in this study are dielectric constant, dielectric strength, resistivity, and thermal stability. POLYMER FILM PREPARATION I) GDP Films. The GDP films were prepared in a dc glow discharge in an atmosphere of organic monomer vapor. The details of the glow-discharge procedure have been presented previously.6 Briefly, under the influence of a dc glow discharge organic vapors can be polymerized in the vapor phase and the resultant polymer (GDP film) deposits on and in the vicinity of the cathode. Substrates placed behind a screen cathode can be uniformly coated over all exposed areas. Pattern definition can be achieved using normal mask stencil techniques. Table I presents the pertinent parameters in the preparation of GDP styrene. Conditions of Forrrzation of GDP Styrene. The rate of deposition for GDP films was monitored using a current-time integrator. The thickness of the deposited film for a given set of experimental conditions was a function of the current-time integral. The fifth row in Table I is a compilation of the deposition rate of GDP styrene films. Other organic monomers behave in a similar fashion. With the exception of the highly halogenated saturated hydrocarbons, all organics subjected to the conditions of the glow discharge form GDP films. 11) Photoresist Films. Photoresists are normally recognized to exist in two forms, a low molecular weight, soluble polymer form (A) and a high molecular weight, crosslinked insoluble form (B). Typical negative photoresists are applied in the low molecular weight form A and converted to the high molecular weight form B by the action of light. (Positive photoresists are applied in form B and converted to form A by the action of light.) Eastman Kodak Metal Etch Resist (KMER) has been chosen as a typical negative photoresist. Polymer film samples were prepared in the same configuration used for GDP films. The KMER was applied in the usual fashion using a 120°C drying step (infrared-lamp heating) and a 180°C post-bake. Care must be taken to completely remove the solvent in the drying step and to properly develop and postbake to exposed photoresist. If the photoresist film is not properly dried and postbaked, severe blistering and peeling will occur during subsequent elevated-temperature exposure. Thickness can be controlled by varying the viscosity of the form A photoresist and the spin rate at application. AZ 340 was chosen as a typical positive photore-

Jan 1, 1967

By C. M. Wayman, S. R. Rauze

It was recently reported' that the initial transformation in (100) iron whiskers grown in the bcc condition can be nucleated on heating in regions of the whisker which are not necessarily the hottest. Although this behavior was observed1 several times, its sporadic nature precluded any documentation other than occasional visual observations during hot-stage optical-microscope heating experiments. The purpose of this communication is to report how, in one instance, the a-y transformation in such a (100) iron whisker was initiated in a region which was clearly at a lower temperature than an adjacent region which did not transform at all.

Jan 1, 1967

By R. W. Parsons, P. R. Chaney

Oil recovery by the imbibition mechanism can be important in fractured carbonate reservoirs with a bottom water drive. Laboratory experiments were performed on water-wet carbonate rocks to model this process. The results are interpreted in light of the applicable scaling laws. The imbibition production behavior of a pillar of rock subjected to a slowly rising water table can be synthesized from data on total immersion experiments on small rock samples. If the rate of water table rise is slow enough and if the physical properties are such that the imbibition zone is of a small vertical extent, then only the final oil recovery on the small samples need be known. Applying laboratory imbibition data to a real reservoir should be done with caution because of possible wettability problems and the unknown behavior of two-phase flow in fractures. INTRODUCTION The spontaneous taking up of a preferential wetting fluid into a porous medium with the simultaneous expulsion of the contained fluid is termed imbibition. One of the many physical imbibition processes is the production of oil by the taking up of water into a water-wet reservoir rock. In certain types of heterogeneous reservoirs, e.g., fractured, it is realized that imbibition alone may be the dominant mechanism. Theoretical studies of imbibition oil recovery have followed several lines of attack. The behavior of certain mathematically tractable and conceptually simple models have been analyzed. Birks1 considered a bundle of vertical capillaries model, while Perotti et al.2 and Kelemen3 used a model of a series of vertically oriented fractures. If the oil production behavior from a block of reservoir rock is known, then the cumulative behavior of a vertical stack of these blocks subjected to a rising water table can be determined by a method outlined by Aronofsky et al. 4 Knowing the capillary pressure and relative permeability functions should uniquely set the imbibition behavior for a given sample. Solutions of the partial differential equations describing the process, however, present some formidable problems. Digital computer techniques have been applied to displacement processes in which the capillary pressure has been included.5"8 More recently, similar techniques have been used by Blair 9 to calculate countercurrent imbibition behavior. The fact that imbibition behavior is described by the usual fluid flow equations (Darcy's Law and continuity equation) implies that the derived scaling laws should apply equally to this situation. Use of model experiments to study imbibition has been employed by Graham and Richardson 10 and Mattax and Kyte,11 who actually tested the scaling law applicability. This study is concerned with scaled model imbibition experiments on water-wet carbonate rocks. Of particular interest is the slowly rising water table experiment. This simulates a highly fractured reservoir with the water table rising uniformly through the vertical fracture system. Experiments were performed on different sized samples and with different rates of rise of the water-oil interface surrounding the rock pillar. Total immersion imbibition experiments on smaller rock samples were performed to test the basic linear scaling laws and to see if these simpler tests could be used to synthesize the rising water table experimental results. THEORY IMBIBITION SCALING LAWS Scaling laws for the flow of two incompressible immiscible fluids in a porous medium have been presented by Rapoport.12 These are derived directly from Darcy's Law for the individual phases and the continuity equation. Results from a laboratory model experiment will duplicate those from some reservoir prototype if a constant proportion is maintained among the three forces acting on the fluids — the capillary pressure gradient, the gravitational gradient and the flowing pressure gradient. Operation of the model (i.e., fluid injection or withdrawal) must be conducted in accordance with specific equations. If the

Jan 1, 1967

By M. C. Flemings, T. Z. Kattamis

Dendrite morphology and grain size were studied in bulk samples of iron and nickel base alloys under-cooled up to 300°C. In the alloys studied, dendrite morphology gradually changes with increasing undercooling from the usual dendritic structure to a structure composed of cylindrical dendrite arms. At a critical undercooling of approximately 170°C in these alloys there is an abrupt transition to a fine-grained structure of spherical morphology. Final dendrite-arnz spacing of melts nucleated at less than the critical undercooling decreases with increasing undercooling and with decreasing distance from a chilled surface. Grain size of melts nucleated at greater than the critical undercooling decreases similarly with these factors. It is concluded that structure coarsening, with reduction of surface area as driving force, is the principal mechanism determining final dendrite-arm spacing in melts nucleated at small undercoolings and grain site in melts nucleated at large undercoolings. IN solidification of highly undercooled melts, nuclea-tion and growth of the solid is accompanied by a rapid release of latent heat of fusion. The heat release is sufficiently rapid in bulk specimens that it takes place essentially adiabatically with respect to the surroundings; hence the specimen heats up ("recalesces") to some maximum recalescence temperature. For bulk liquid metals and alloys, a simple thermal balance shows the undercoolings obtainable in bulk samples are not sufficient to lower the maximum recalescence temperature below the equilibrium solidus temperature. Hence, some liquid always remains at the end of recalescence; this liquid must solidify by rejection of heat to the surroundings. Metals, including iron and nickel base alloys, normally freeze with only a few degrees of undercooling. In general, nucleation begins on impurity particles present in the metal, at a "heterogeneous nucleation temperature" that is quite close to the equilibrium liquidus temperature. Very small droplets of metal can be undercooled considerably before nucleation occurs,1-4 presumably because subdivision of the liquid reduces the probability of solid impurity in a given droplet. Following early work on metal droplets,' a series of studies has shown that it is possible, without great difficulty, to undercool bulk samples of many metals and alloys. Undercoolings obtainable are approximately two tenths of the melting point for nickel, iron, copper, and other metals and alloys.5-9 Specimen sizes have generally been in the range of 50 to 500 g. Studies on undercooling to date have been concerned primarily with techniques of undercooling and with nucleation and growth kinetics. Only limited consideration has been given to dendrite morphology in under- cooled alloys. One such study, in bulk Cu-Ni alloy,8 showed that at small undercoolings a typical dendritic structure was obtained. At undercoolings the order of 175°C the structure of a Cu-40 pct Ni alloy became "ray or starlike", and a dense network of subbound-aries appeared. At still larger degrees of undercooling, of the order of 250°C, the structure became more random and included an ensemble of globular forms. Comparable results were obtained by walker10 on Cu-Ni alloys, and by Kamenetskaya et al.11 in research on Cu-Ni, Sb-Bi, Fe-Ni, and Fe-Cr alloys. In this work, study was made of the processing parameters influencing dendrite morphology and grain size in undercooled alloy melts. Experiments were specifically designed to show the effects of coarsening ("ripening") during solidification on final casting morphology. EXPERIMENTAL PROCEDURE Apparatus used was patterned after that of Walker,10 Fig. l(a). Melting was by induction under reduced pressure of helium. Two different types of crucible materials were used; there was no difficulty in obtaining substantial undercooling in both types of materials. Alumina crucibles were used for the bulk of the runs made, with a glass coating separating the metal from the crucible walls. Other runs were made using fused silica crucibles with or without slag coating. A graphite susceptor was employed around the crucible to aid in melting and to reduce melt stirring. The major portion of the experimental work was on

Jan 1, 1967

By D. R. Uhlmann, T. P. Seward, O. Chalmers

A magnetic field has been used to damp out liquid convection during the solidification of castings. Ingots of Al-2 pct Cu were poured from various sperheats (10" to 65°C) into a room-temperature graphite mold located between the poles of an electromagnet. The effect on ingot structure of varying the strength and time of application of the magnetic field was investigated, and the results are interpreted in light 01- current theories of ingot structure. In the solidification of alloy systems, thermal gradients in the melt and differences in density of the constituents are expected to produce convective motion in the liquid. The thermal convection caused by temperature gradients should be present during the freezing of both alloys and pure materials, while the solute convection caused by density differences should be of concern only for alloys. The possible influence of this natural convection on the structure of castings has been discussed by Ruddle,' Chalmers, and Tiller. More recently, an experimental investigation of such effects was carried out by Cole and Bolling. These workers took stationary molds containing lead- and tin-based alloys and. using a brine solution. quenched them from various superheats above their respective liquidus temperatures. When grids were placed across the mold diameter to suppress convection, the grid region was later observed to contain a larger percentage of columnar grains than other parts of the casting or than similar castings poured without the grids present. From thermocouple measurements, they found the temperature in the grid region to be higher than that in the region outside the grids throughout the freezing process, and concluded that convection should be an important mode of heat transfer in normal castings. In addition to such convection arising from thermal gradients or differences in density of the constituents, a further source of convective motion, which we shall call "momentum convection", may be particularly important in small (laboratory-size) castings. By this term we refer to the convective motion produced in the initial pouring process, a motion which is damped out as the linear momentum and angular momentum of the liquid are dissipated. Besides these sources of convection, which are present to some extent in all castings. external influences may also be effective in producing convective motion. Examples of such "forced" convection which have been discussed in the literature include stirring the melt, shaking the mold, rotating the mold, and so forth The purpose of the present study was to examine the effect of convection on the structure of castings poured into a cold mold. In this work, a magnetic field was employed to damp out the convective motions within the liquid, and hence control the extent of con- The effect of a magnetic field on convective currents in a conducting liquid is to induce eddy currents which oppose and damp out the convective motion. vection. The technique of magnetic damping was previously used by Utech and F1emings" in their study of convective effects on impurity banding and other substructure in alloy crystals grown unidirectionally from the melt. In their investigation, fields as small as 200 gauss were found to be effective in damping the convection currents which influence the banding. I) EXPERIMENTAL TECHNIQUE In the present study, alloys of A1-2 wt pct Cu were poured into a room-temperature graphite mold. The mold, in the shape of a circular cylinder, 2 in. ID, 2-3/4 in. OD, 5 in. inside height, and 5-1/2 in. outside

Jan 1, 1967

By E. A. Breitenbach

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

Jan 1, 1967

By N. J. Olson, G. W. Toop

A regular solution is characterized by a finite molar heat of mixing and entropy due only to ideal or random mixing. The thermodynamic properties of binary regular solutions which exhibit constant atomic coordination number and constant atom-atom intefac-

Jan 1, 1967

By M. A. Jones

Performance of a pilot flood in eastern Kansas indicates improved recovery and accelerated production resulting from mobility ratio control obtained by adding a high moleculur weight polymer to injected water. Ultimate pilot flood recovery of the 75-cp reservoir oil is anticipated to he 350 bbllacre-ft during the six-year flood life. Cumulative recovery to date totals 242 bbl/acre-ft after 30 months of injection. Reservoir parameters and flood performance are disc~rssed in detail to provide a case history relating to the feasibility of polymer flooding. INTRODUCTION The Vernon Polymer flood of the Brazos Oil & Gas Co., located in Sec. 2, 24 S., 16 E., Woodson County, Kans., was begun in late Oct., 1963. Injection of a 500 ppm solution of Pusher* polymer was initiated in a 15-acre pilot to determine feasibility of improving oil recovery by mobility ratio control. Performance to date and the anticipated recovery of 115,000 bbl from the 329 acre-ft pilot indicate efficient displacement of the viscous reservoir oil. Dependence of flood recovery on the displacing fluid-displaced oil mobility ratio has long been recognized. Displacement efficiency of the reservoir oil decreases with increasing mobility ratio. Therefore, economic reduction of this parameter is an attractive approach to improving oil recovery. particularly from viscous oil reservoirs where mobility ratios are often quite high. Waterflood mobility ratios can be controlled by dilute solutions of certain polymers'; which provide a displacing fluid with considerably reduced mobility. Reservoir parameters and flood performance are emphasized in this paper to provide data relating to polymer flood applicability. THEORY The polymer used in the Vernon pilot exhibits the property of unusually high resistance to flow through porous media.'.' A small concentration in injected water effects a mobility reduction considerably beyond that predicted from increases in bulk viscosity as measured in a capillary or rotating cup viscosimeter. Water mobility in the subject flood is reduced by a factor of 7.55 with the addition of 500 ppm Pusher chemical to injected water while the bulk viscosity is increased from 1.03 to only 1.40 cp at 75F. The exact flow mechanism causing resistance to flow has not been established; however, it is observed only in tortuous flow systems and appears to be influenced by brine properties, rock properties and oil saturation. Consequently, mobility reduction is estimated from flow studies for a given sample of reservoir rock and is expressed in terms of a "resistance factor".' This parameter is defined as the ratio of brine mobility to polymer solution mobility measured at residual oil saturation: RF= ?m/?p= (km/µm)/(kp/µp)M.......(1) where RF = resistance factor (dimensionless) k*. and A, = mobility of water and polymer solution, respectively k,,, and k,, = water and polymer permeability, respectively (darcies) µie and µp = viscosity of water and polymer, respectively (poise). Eq. 1 is similar to the expression commonly used to relate the relative mobility of the water bank at residual oil saturation to oil-bank mobility at irreducible water saturation: where M =mobility ratio of water to oil (dimension-less) A, = mobility of oil k, = permeability to oil (darcies) )i* = oil viscosity (poise). From Eqs. 1 and 2 it is apparent that the polymer-solution-to-oil-mobility ratio can be approximated by dividing the brine-to-oil ratio by the resistance factor RF: Mpn=M10/RF...........(3) where M,-, = mobility ratio of polymer solution to oil (dimensionless). This simplified presentation assumes no dilution of the polymer solution by connate water. Also, the presence of a connate water bank separating the polymer solution from the displaced oil is not considered. Displacement of oil by waterflooding is increasingly dominated by viscous forces as the viscosity of the oil displaced increases. When the oil mobility is less than that of the injected water, the oil-water interface is not a piston-like front.3,' Instead, instabilities occur at the flood front in the form of viscous fingers which tend to penetrate the less mobile oil bank. The magnitude of these protrusions, or the degree of channeling, increases

Jan 1, 1967

By H. O. McLeod, J. M. Campbell

Data analysis of pentane and hexane adsorption from natural gas in a fixed bed of silica gel shows that constant length mass transfer zones form, the curvature of the adsorption isotherm controls the growth of the mass transfer zone and surface diffusion of molecules inside the silica gel particle controls the mass transfer rate. Curvature of the hexane isothem is more than the curvature of the pentane isotherm. Because of this curvature the hexane adsorption zones reached a constant length. In contrast, the pentane adsorption zones were always increasing in length during each run. A procedure was developed to obtain correct mass transfer coefficients using effluent curve slopes. These transfer roefficients increase with the amount of hydrocarbon adsorbed on the silica gel particle. The characteristic shape of the hexane effluent curves also show that molecular diffusion inside the silica gel particle controls the adsorption rate of pentane and hexane. INTRODUCTION The purpose of this study was to determine the mechanisms that control the dynamic adsorption of hydrocarbons from a natural gas onto silica gel. Before one can deal effectively with multicomponent adsorption, the transfer mechanisms by which a single hydrocarbon component is adsorbed from the gas stream must be defined. Two principal investigations of this system have been published ' and indicate that diffusion through the gas around the particle controls the adsorption rate. Some of the experimental observations in each study either do not support this transfer mechanism or are inconsistent with the mathematical model used in analysis. In this study surface diffusion of molecules inside the particle controls the mass transfer rate of pentane and hexane. This mechanism is indicated by the effluent curve shape for a constant length transfer zone and by the variation of the mass transfer coefficient with concentration of the adsorbed hydrocarbon. THEORY AND DEFINITIONS — MATHEMATICAL MODELS Mathematical solutions for the isothermal adsorption of a trace component from a carrier gas are derived from three relationships: the mass balance or continuity equation, the equilibrium relationship between the gas and solid phases, and a mass transfer rate equation. The transfer rate is proportional to the adsorbate concentration gradient within either the gas or solid phase. Mathematical solutions of these equations usually give the adsorbate concentration as a function of time and distance from the bed inlet. That part of the bed in which the adsorbate concentration changes from a maximum to a minimum value is called the transfer zone. This transfer zone is directly related to a plot of the effluent concentration vs time which has a characteristic S-shape. This general shape is determined by the continuity equation and occurs in many processes of diffusional transfer.

Jan 1, 1967

By K. P. V. Lei

The Bureau of Mines investigated molten salt electrorefining 'n helium atmosphere cells as a means for the production of ductile vanadium. 'Electrolytes composed of NaCl-VCl2, LiCl-NaCl-VCl2, BaCl2-KCl-NaCl-VCl2 KCl-LiCl-VCl2, and CaCl2-NaCl-VCl2 were tested. Two electrolytes, KCl-LiCl-VCl2 and CaCl2-NaCl-VC12, were selected for laboratory-scale production tests. A molten salt electrorefining process was developed for the production of ductile 99.6 percent vanadium from a commercial 90-percent grade. Using the CaCl2-NaCl-VC12 electrolyte, an 87-percent recovery of ductile vanadium was made; an 85-percent recovery was made with the KCl-LiCl-VCl2 electrolyte.

Jan 1, 1967

By G. W. Nabor, A. S. Odeh

The effect of production history of a well on the results of two-rate flow tests, and conventional build-up analyses was investigated. The effect was examined by means of digital computers and an R-C network model, respectively, for wells with infinite and finite radii of drainage. For systems which behave as infinite, it was found that production history and the duration of production at constant rate prior to the initiation of the test have important effects on the results. During build-up time equal to about one-fourth of the stabilized time, correct permeability-thickness product calculations can be made. For wells with finite radii of drainage, the time was determined during which the straight line can be satisfactorily used for permeability-thickness product calculations in case of drawdowns and build-ups. On build-ups, the dimensionless time (based on the external radius) during which the straight line gives reliable results was detertriined to be about 1/12. This is one-fourth as long as that of the drawdown. The investigation was done theoretically, and subsequently was verified by R-C network model runs. General interpretive rules were formulated which, if not followed, could lead to serious errors. Moreover, a recommended testing procedure is reported. INTRODUCTION The method used by most reservoir engineers for estimating formation characteristics in a producing well is the analysis of pressure build-up data. The method originally devised by Horner' makes use of the point source solution to the diffusion equation. This solution is approximated by a logarithmic function and the superposition principle is employed to arrive at the well known pressure build-up equation:where q, the flow rate, is in reservoir B/D; ft is in cp; kh is in md-ft; At is the shut-in time; and t is the producing time. At and t are in any consistent time units. Ey. 1 is applicable to a well of unlimited drainage radius which produces at a constant rate q from zero to time t and is then shut in. Such a constant production rate seldom obtains in practice. Therefore. a correction term must be applied to Eq. I to account for the varying rate. Two theoretically accurate methods are available for treating the variable rate case. The first, originally derived by Horner,' is based on the application of the superposition theorem. It requires knowlege of production history of the well as a function of time and results in lengthy and laborious calculations. The second t*q* niethod is suited for short production tests and requires that the shut-in time be at least one and one-half times the production time. A third method which is not based on any theoretical justification and which was suggested by Horner as a "good working approximation" is the one used by the majority of reservoir analysts. especially when the well has been producing for a long time and the t*q* method is not practicable. The key to this method is in choosing or determining the t that appears in Eq. 1. Horner suggested using a corrected time t, in place of t. t, is calculated by dividing the total cumulative production by the last established rate. Therefore, a normal procedure of pressure build-up testing is to stabilize the well at a constant rate for at least 24 hours before shut in and to use the stabilized rate to calculate t,. The analysis is then made by plotting either or P. and examining the resulting plot for the expected straight line to calculate kh and the original reservoir pressure. Recently, Russell" proposed a method for determining formation characteristics from two-rate flow tests. His method reduces to pressure build-up if the second flow rate is zero. Russell uses the Horner simplified procedure for calculating a corrected t,. His method also requires the stabilization of the well at a constant rate q which is used to calculate t,. Theoretically, the above procedure is valid for a well with an unlimited radius of drainage or with a limited radius as long as the boundary effect has not been felt by the well. Several authors'' derived formulas which allow the estimation of time during which limited reservoirs behave as infinite ones and, thus, can be treated by unsteady-state mechanics. One equation derived by Swift and Kiel' terminates the application of unsteady-state theory when the drainage radius reaches one-half the reservoir radius. Thereafter, steady-state behavior obtains. Another equation derived by Jones'" initiates steady-state

Jan 1, 1967