Search Documents

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 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 S. G. Gibbs, A. B. Neely

Jan 1, 1967

Over half of the molten aluminium produced at the 52,000 tons/ year eapacity plant operated by Comalco Aluminium (Bell Bay) Limited in Tasmania is cast into saleable produets by the direct chill process. Considering the importance of direct chill casting to the aluminium industry it is surprising that so little published literature on the subject is available, although Easton (1947), Smart (1954), and Uppert (1940, 1944, 1945) have contributed much toward lifting the veil of secrecy surrounding the process.The direct chill, or D.C. process as it is commonly called when related to aluminium, consists basically of feeding molten aluminium or its alloys into; water-cooled mould from the bottom of which solid aluminium is continuously extracted. Although the casting procedure is continuous in principle, the equipment is normally operated in an intermittent manner. Castings up to 20 ft long are formed during a cycle. At the end of the cycle casting is stopped and the pit into which thc castings are dropped is cleared before a new cycle commences. Depending on the size of the casting in relation to the cross sectional area of the pit and the support mechanism (commonly a hydraulic ram), 'drops' of up to 15 tons are not uncommon. This method of D.C. casting IS termed vertical D.C. casting and is the one most commonly employed in the aluminium industry.Horizontal casting has come more to the fore in recent years with the casting being extracted from the mould in the horizontal plane and is a truly continuous casting method since flying saws can be simply used to cut the continuously cast shape into required lengths as it emerges from the mould...

Jan 1, 1967

By R. Simon

Maurer's review1 summarizes quite thoroughly the various theoretical developments and experimental findings that contribute to the knowledge of rock mechanics in drilling. This discussion supplements Maurer's paper by including generalities that either have been incompletely treated in his review, or formulated from new considerations by the discussant and from hitherto unpublished experimental data. These generalities are associated with the observation that all mechanical methods of rock drilling essentially comprise a sequence of indexed, concentrated loadings on the bottom of the hole being drilled. A concentrated loading is defined as a mechanical loading by a tool that contacts at any one time only a small fraction of the surface of the body being worked and penetrates only a small fraction of the body thickness. The surface in this application is the bottom of the hole. The deformation of ductile material by concentrated loading has received considerable attention because of the importance of indentation tests for measuring hardness of metals. Some of the generalities involved in concentrated loadings on ductile materials apply also to concentrated loadings on more brittle materials such as rock. Among these generalities is the applicability of similitude considerations to scaling laws and to force-displacement and energy-volume characteristics for tools of various shapes, as summarized by Tabor.2 Vork on the flow and fracture behaviors of glass loaded with pyramidal and spherical indenters is reviewed by Douglas.3 Maurer refers to the work of Evans and Murrell on wedge loadings on coal and on mixtures of pitch and sand. It is evident from Maurer's paper that a considerable portion of the studies of rock mechanics in drilling, such as crater formation experiments, has involved concentrated loadings on rock.

Jan 1, 1967

By Thornton L. Neathery

Geologic investigation of a mineralized zone in eastern Alabama indicates that large quantities of talc and soapstone are scattered over a wide area close to the surface. Laboratory experiments show that this material may be upgraded to a high-quality talc. A market survey indicates that the talc deposits may be of current and long-range interest. The anthophyllite asbestos which is associated with the talc does not appear to occur in sufficient quantities to be considered as a primary commercial source.

Jan 1, 1967

By P. T. Allsman

INTRODUCTION The honor of being the recipient of the Mining Society's Daniel C. Jackling Award for 1966 imposes on me the responsibility to search for the subject that I would be best qualified to discuss. It seems that my experience in the mineral industry — nearly 15 years in private and 30 years in public employment — should uniquely qualify me to comment upon the relationships between the mining industry, mineral engineering as a profession, and government, and the changing economic, technologic, political and social changes that dictate their relationships. Such conditions, the increasing knowledge of the world mineral resource base, the greater demands of an expanding population with a rising standard of living, shall be evaluated with respect to our profession — mineral engineering. In all of these forty or so years, I have been in close touch with many phases of the mineral and energy resource industries, especially those of the United States. Particularly during the past 10 years, I have had a ringside seat at the "goings on" in Washington. During all of this period I have observed and analyzed many of the problems of the industry and its relations with the government. During the past few years, I have had increased opportunity for doing research and have found a need to learn something about the mathematical techniques that are undergoing such a revival since electronic computers were developed. One of these techniques is multivariate analysis, which is a system for correlating a large number of complex variables. By applying this technique to the seemingly chaotic present-day situation, there emerges a course of steady improvement in the minerals resource industries that wars, depressions, and political strife have failed to obstruct for long. In short, I advocate a need for new concepts in mineral engineering to cope with the changes, pos- sibly adding many more specialized areas on one end of the spectrum, and broad engineer-generalist approaches at the other. The meeting place shall be where the intricate mineral resource problems at the frontier of knowledge merge into the socio-techno-logical problems of the world. HISTORICAL INVOLVEMENT OF GOVERNMENT IN MINERAL RESOURCE PROBLEMS Involvement of government in mineral resource problems undoubtedly is as old as government itself. But passage of the Mining Act of 1866, superseded by the Act of 1872, was the first major entry of our Federal Government into the problems of mining. Passage of the law was necessary to resolve many conflicting claims and disputes in mining. The General Mining Law of 1872 and the Mineral Leasing Laws of 1920 apply to millions of acres of public lands and to certain privately owned land on which the mineral rights have been reserved to the Federal Government. The extensive involvement of government in the mineral industry business today often comes as a shock to the mineral engineer who has begun to climb the ladder to the top in his company. He is often poorly prepared to cope with the situation and becomes frustrated and even hostile. Better knowledge of the role of government perhaps would ease his situation. In our political system, private industry is responsible for producing minerals. To make certain that the production is adequate for a peace-time economy, and to guarantee the Nation's security in an emergency, the Federal Government finds it expedient to develop policy and engage in regulatory, educational, and research activities. A few examples will illustrate the varied role of government in relation to mineral problems. In varying degrees many agencies of the Executive Branch of Government have responsibilities that influence the functioning of the mineral and energy industries. In addition, of course, Congress, the Judiciary, and State and local governments exercise great influence. The Department of the Interior is primarily responsible for activities to promote the development and conservation of mineral and fuel resources, and the

Jan 1, 1967

By Robert C. Bates

Variations in water percolation rates through saturated cohesionless hydraulic backfill (mill tailings) have been studied in relation to the many physical aspects of the system to determine the most important factors in causing these variations. A multivariable least-squares regression analysis of 135 separate tests shows that a reliable estimate of the water percolation rate through typical hydraulic backfill can be made if the bulk density (unit weight), specific gravity, and grain-size distribution are known or can be estimated. Equations for predicting percolation rate, coefficient of permeability, and seepage velocity are developed both for the range of data analyzed and for extrapolation. The prediction equations for use within the range of data studied are quite accurate, while the ones developed for extrapolation are less accurate but still give reasonable approximations of flow rates. The equations given are valid for cohesionless materials that are free of clays and micas, exemplified by classified mill tailings. The equation developed for predicting percolation rates is as follows: 1nP20 = 11.39147 + 2.853422 1n(VR X D10) + 0.1747436 VR X Cu -178.8039 D10 X D50 + 311.7034 (DSO )3, (1) where P20 = percolation rate, inches per hour at 20° C; correction factors for other temperatures are given in the text, VR = void ratio, dimensionless, D10 = effective size (diameter in millimeters of the 10-percent-finer grain size), Cu = coefficient of uniformity (ratio of the diameter of the 60-percent-finer grain size over D10), D50 = average size (diameter in millimeters of the 50-percent-finer grain size), and In = natural logarithm (log to the base e). The prediction equation for percolation rate is rather complicated, so for practical use a nomograph that simplifies the solution of the prediction equation was developed.

Jan 1, 1967

By K. Leutwyler

The key to realistic casing stress analysis in thermal recovery installations is accurate knowledge of the temperatures involved. Much information leading to prediction of heat losses from tubing string to wellbore has been published. This .same information permits calculation of casing temperature under certain conditions. This paper reviews the basis for most of the published work, the line source solution to the diflusivity equation and examines some of its application criteria in the light of unsteady-state conditions. An alternate finite differences method is presented for correlation, and experimental data pertaining to the heat transfer mechanism in the annulus is reported. In the light of the information presented in the paper, it appears that the assumption of a line source for the tubing-wellbore system is valid even for reasonably short periods of heat flow, .such as "steam soak" operations. The empirically modified equation presented earlier matches computer derived results of the finite differences approach lor periods as short as 50 hours. Previously predicted cas ing temperatures and the associated stress levels are thus verified based on these correlative computations and the reported experimental information. Fall out from the computer output resulting in some studies of the influence of insulating cement properties on casing temperature, and experimentally studied effects of nitrogen placed in the annulus at various pressures., are described. INTRODUCTION Generally speaking, predictions of average casing temperatures have been based on an idealized model involving a centralized tubing string at uniform constant temperature radiating energy towards the casing under steady-state conditions.'. ' Heat is transferred away from the casing to the formation by unsteady-state conduction across thermal barriers consisting of cement and the mud sheath. Heat balances written on the casing permit prediction of instantaneous casing temperatures. During early periods of injection, casing temperature varies appreciably as a function of time; however, for long injection times this variation decays as heat flow from the wellbore to the formation stabilizes, and it becomes virtually negligible as conditions surrounding the wellbore become quasi-steady-state (Fig. 1). Carslaw and Jaeger -escribe a solution to the diffusivity equation for a line source (Fig. 2). A constant supply of heat is introduced into the medium from the line placed into the origin of the x-y coordinate system and parallel to the z-axis. Heat spreads outward in the infinite cylindrical solid and the temperature increase on a line going through point (x, y) parallel to the z-axis at a distance r can be predicted by: The following criterion must be considered: the Carslaw and Jaeger equation is obtained on the basis of a large

Jan 1, 1967

By I. C. Edwards, G. E. Agar

INTRODUCTION Many research efforts over the past years have been devoted to the study of energy consumption in comminution and the characteristics of the broken material. The old "laws" have been superseded by more general methods of analysis1,2,3 but these methods have not received much attention outside the laboratory. The work to be described here was done in order to obtain scale-up information on wet rod milling of sylvinite. EXPERIMENTAL A 30 x 51 in. Denver rod mill with a scoop feeder was used for these tests. Ore was fed from a variable speed screw conveyor into the feed box along with saturated brine. The mill was filled to approximately 45% of its volume with 2400 lbs of a mixed rod charge. The minimum possible feed rate was 2.5 tph and the practical maximum for this mill was 8 tph. Individual tests were run by setting the solid feed rate, then adjusting the brine flow rate to give the desired pulp density. After steady state, as evidenced by constant discharge pulp density, had been reached, samples of the discharge were taken at 15-minute intervals. These samples were filtered, washed with alcohol, dried, then screened on Tyler screens using a Ro-tap shaker. The sylvinite feed was regular plant feed that had been scalped at 20 mesh and was tapped from the coarse ore bins. It had been previously passed through a roll crusher set at 3 inches and after scalping was essentially minus 1 in. plus 20 mesh. RESULTS The screen analysis of samples from each series of tests was plotted on double logarithmic scale paper to determine the reliability of the series. In some cases, a faulty sample or screening procedure dictated that the results of one test be discarded. The criterion is a simple one — if the slope of one plot was noticeably different from the others in the series, it was discarded. The only cases that caused difficulty were those in which the brine was not thoroughly displaced from the solids. When this occurred agglomerates were easily identified on the screens and the whole of the sample was retained on 3 or 4 screens. There were no border line cases. The number N given in Table 1 is the number of size analyses that were pooled to form the data shown in Fig. 1. While examining the reliability of the sizing and sampling by plotting the data, it became apparent that there was no significant difference between grinds made at 50 and 60% solids. For the analysis given here, these tests were regarded

Jan 1, 1967

By C. N. J. Wagner, E. N. Aqua

The bee alloys of the terminal solid solution of rhenium in niobium were investigated by X-ray diffraclion methods. The analysis of the broadening of the powder pattern peaks from the niobium-rich alloys, cold-worked by filing, showed evidence for faulting on (211) planes. The alloys with larger rhenium concentrations were brittle, and fractured when filed. The analysis of the shapes of the diffraction profiles from powders produced by hammering with steel plates) of these brittle alloys indicated that the broadening was due predominantly to platelets bounded by clibe planes. This result was consistent with the Laue patterns from single-cleazlage facets which showed evidence for clearage fracture on {100} planes. The influence of faulting on the mechanical properties of face-centered cubic (fcc) metals and alloys, i.e., how faulting affects the various dislocation configurations in fcc metals, is well-known.' The important question is whether or not faulting occurs in body-centered cubic (bcc) metals and alloys, and, furthermore, if faulting does occur, to what extent it affects the mechanical properties. Theoretical calculations indicate that the stacking-fault energy is rather high in most bcc metals.2 The results of sev-eral studies using transmission electron microscopy have indicated that faults are present in well-annealed and lightly deformed bcc metals and alloys, such as tungsten,3 vanadium,4 and Mo-Re alloys.5,6 These faults, as have been observed, may be stabilized by the presence of interstitial (C,N,O) atom segregation to the extended dislocations. In addition, there is considerable X-ray diffraction evidence for faulting in bee metals,7-8 using the same techniques that were successful for the study of faulting in fee metals and alloys.7 This X-ray diffraction method is currently the only technique applicable for the investigation of faulting in metals with high stacking fault energies, i.e., when severe amounts of deformation are required to produce sufficient amounts of faulting to be detectable by X-rays. The thinned films of heavily deformed metals will become opaque in the electron microscope due to the high dislocation densities present. In a previous investigation of faulting8 in bee metals of group Vb and VIb, the authors observed that niobium was ductile when filed at room temperature, and showed the highest density of faults on the (112) family of planes. A niobium-base alloy system was therefore chosen to study the effect of alloying on the occurrence of faults in a bcc alloy. There is a large range of terminal solid solubility of rhenium in the alloys of the bee refractory metals. All of the Group VIb-rhenium alloys, i.e., Cr-Re, Mo-Re, and W-Re, exhibit solid solution hardening without loss of ductility.9 In fact, the ductility reaches a maximum for alloys with compositions near to the limits of solubility. This enhanced ductility has been attributed to the lowering of the stacking-fault energy in the alloy; i.e., the increased amounts of twinning provide an additional mode of deformation, resulting in the observed high ductility.' Because of this demonstrated effect of rhenium in bee alloys, it was decided to study a series of alloys of Nb-Re. Using the same X-ray diffraction methods as described in the previous investigation,' one is able to identify the possible causes of broadening of the powder pattern peaks, e.q., particle size, microstrains, stacking faults, and/or twin faults. The separation of these various contributions to the measured peak shape is accomplished by the Warren-Averbach analysis.7 I) EXPERIMENTAL PROCEDURE A series of seven bee solid solution alloys containing from 2.5 to 21.5 at. pct Re (remainder niobium) were prepared at the General Electric Research Laboratory. Buttons, weighing 25 g, were alloyed from electron-beam-melted niobium chips and random rhenium sheet by arc melting in a water-cooled copper crucible, under an argon atmosphere. The buttons were inverted and completely remelted to promote homogenization. The nominal alloy compositions are listed in Table I. Samples of the cold-worked alloys were obtained by filing to produce a powder for the X-ray study. Alloys too brittle to be filed were hammered with steel plates. The resultant filings or brittle impactings were screened through 150-mesh classifiers and compacted to the required shape for the diffractometer holder. The X-ray diffraction peaks were automatically recorded by fixed time counting at uniform 28 intervals, or by continuous registration with a ratemeter. In all cases, the (110)-(220) and (200)-(400) pairs of peaks and the (211), (310), (222), and (321) reflections were measured on a GE-XRD5 diffractometer using CuKa radiation with nickel filter or MoKct radiation and a LiF monochromator in the diffracted beam. The first part of the data reduction, completed with the aid of a series of programmed computations,10 was the analytical resolution of the Kol component of the diffraction profile. Subsequent computations were then made with reference only to the resolved Kol peak shape. These included the measurement of the

Jan 1, 1967

Jan 1, 1967

By G. L. Schroeder

Investigations were conducted at two underground locations, foreman's room and 5702 area, of the Kermac Nuclear Fuels Corp. uranium mining installation, Ambrosia Lake, New Mexico. The rock environment is stratified sandstone. Pressures were varied from minus 0.4 an Hg to plus 3.8 cm Hg with respect to barometric pressure, and were maintained for periods ranging from several hours to 21 days. One cm Hg overpressure (-1.5% of normal barometric pressure) produced a 20-fold decrease in Rn flux into the foreman's room as compared to the Rn flux observed during dilution ventilation. The 5702 area experienced a 5-fold decrease in Rn flux for a similar overpressure. During 21 days of continuous 1 cm Hg overpressure in the 5702 area this reduced Rn flux was maintained. The smaller reduction in flux for the 5702 area results in part from the proximity of the foreman's room to a non-pressured area and the presence of an impermeable shale lens overlying the ceiling of the 5702 area. A 0.4 cm Hg underpressure of the 5702 area increased the Rn flux into the area by 3-fold over the non-pressured value. A barometric pressure influence on Rn flux into the mine and on air Rn concentrations in non-pressurd areas of the mine was measured. A falling barometric pressure causes an increase in the flux of interstitial gas into the mine; a rising barometric pressure decreases this flux. These effects are not seen in overpressured mine areas. Radon concentrations in rock interstices show, during overpressure, decreases of all locations monitored except in the immediate vicinity of the shale lens above the 5702 area ceiling. These reductions are the result of air flowing from the mine into the rock interstices under the influence of the pressure gradient induced by the overpressure. The relative humidity of the 5702 area was reduced after 2 days of 1 cm Hg overpressure from a pre-overpressure value of 96% to 81%. The control of harmful or noisesome gases which emanate into underground mining areas from surrounding rocks is often a costly necessity. In uranium (U) mines this problem is presented by the high interstitial concentrations of radon-222 (Rn), a radioactive noble gas of 3.82-day half-period, formed in the decay scheme of U238, and the relative ease with which this noble gas diffuses from its parent rock into the interstices of the rock and then to a working area. Long-term exposure to moderate concentrations of Rn (3000 µµc per 1) and its short-lived alpha-ray emitting daughter products, RaA (Po218) and RaC (po214), has been found to result in an increased incidence of bronchial carcinomaa2 The acceptable limit of Rn daughter products in a mine atmosphere is presently set at a concentration of short-lived radon daughters which would be in equilibrium with 100 ppc of Rn per 1 of air. Ventilation of mines is conventionally accomplished by continuously diluting the mine air with fresh air brought in by fan or natural convection from above ground surface. The diluted air is allowed to vent freely to the atmosphere. This process does nothing to retard the rate at which gases (including Rn) diffuse into the mine from surrounding rock. It merely reduces the concentrations of these gases within the mine. It is suggested that in a sufficiently permeable rock environment a slight pressurization of the mine air may induce air flow from the mine into the interstices of rock surrounding the mine. The convective flow of comparatively fresh air will reduce the concentration and the concentration gradient of undesirable gases in the rock interstices adjacent to the mine. Radon will diffuse counter-current to this convective flow at a rate proportional to its concentration gradient. In time, a steady state of convection and counter diffusion will be established. The Rn flux into the mine will never be reduced to zero, and will only asymptotically approach zero as the overpressure (and hence the convective velocity of purging air) approaches infinity. During 1963 and 1964, an investigation into the effect which such applied pressures have upon the Rn characteristics of a mine atmosphere was conducted at the Kermac Nuclear Fuels Corp. uranium mines, Ambrosia Lake (near Grants), New Mexico.

Jan 1, 1967

By R. W. Parsons

The over-all apparent single-phase permeability of fracture-rock systems was studied using two different two-dimensional models. In a strict sense the results are applicable only to these models, yet may give some insight into real reseruoir behavior. The regular fracture-matrix model can have any number of regular fracture sets. All fractures within each set have the same orientation, width and spacing. The gross behavior is amlytically equivalent to that of an anisotropic permeable medium. In the heterogeneous fracture systems, individual fracture conductivities vary within a two-dinensioml network pattern. A computer program calculates the overall permeability after the conductivities are randomly placed in the network. The calculated permeability depends upon the distribution of these conductivities and the particular fracture pattern into which they are placed. Although detailed information on actual reseruoir fracture systems is not available, this study shows how some of the possible variables may affect obserued gross reservoir characteristics. INTRODUCTION Characteristics of reservoir rocks are determined by direct (core analysis) or indirect (logging and production tests) examination. Both methods, singly and together, are widely used in reservoir evaluations. For fractured reservoirs a problem exists in relating characteristics of cores to the in situ reservoir properties. Knowing fracture characteristics can be of value in determining production mechanisms. Because of the limitations of observation, the intimate details of a fractured rock system will never be known. One approach in analyzing fractures is to study the behavior of some conceptually simple models. The two models analyzed contain features which may be realistic, yet are mathematically tractable. Flow is considered in two dimensions only, with the plane of interest at right angles with the fractures. This concept represents the horizontal flow in a reservoir system with the fractures oriented vertically. Only single-phase, laminar flow is considered. The regular fracture-rock system permeability is calculable from a simple analytical expression. The heterogeneous system requires a computer program for calculating overall permeability.

Jan 1, 1967

By K. H. Coats, M. H. Terhune

Analysis and example applications have been performed to compare the accuracy and computing speed of alternating-direction explicit and implicit procedures (ADEP and ADIP) in numerical solution of reservoir fluid flow problems. ADIP yields significantly greater accuracy and requires about 60 per cent more computing time than ADEP, not 300 or 500 per cent more as reported elsewhere. l,2 INTRODUCTION Several recent papers 1-3 discuss an alternating-direction explicit difference approximation (ADEP) to the diffusion equation. Example applications of ADEP and ADIP4 were reported to support conclusions that ADEP is comparable in accuracy to ADIP and requires one-fifth to one-third the computing time of ADIP. Applications of ADEP in calculation of two-phase flow in reservoirs was also proposed. 3 This study was performed to compare further the relative merits of ADEP and ADIP in simulation of two-dimensional flow of one and two fluid phases in reservoirs. Since two-phase flow equations are often essentially elliptic rather than parabolic, the efficiency of ADEP in solving the elliptic equation was also examined.

Jan 1, 1967

Jan 1, 1967

Table 1.-Salient statistics of carbon black produced from natural gas and liquid hydrocarbons in the United States [ ]

Jan 1, 1967

Major activities of the Bureau of Mines Explosives Research Center during fiscal year 1966 (July 1, 1965 to June 30, 1966) are reviewed briefly. Part 1 summarizes significant accomplishments of the projects that were active during the report period. Part 2 presents short abstracts of the publications that appeared in the Bureau series and in other media during fiscal year 1966. Publications issued after June 30, 1966, are noted where pertinent in foot- notes but are not abstracted. ACKNOWLEDGMENTS The Bureau of Mines gratefully acknowledges the support given through cooperative agreements by the following organizations to investigations conducted by the Explosives Research Center during fiscal year 1966: Manufacturing Chemists' Association. The Sulphur Institute.

Jan 1, 1967

By K. E. Stanfield

Oil-yield results are presented for 65 oil-shale sections in Piceance Creek Basin, Colorado, sampled during the period 1957 through 1963. Including previously published data, oil-yield results for 203 sampled sections have now been released by the Bureau of Mines to aid in appraising the shale oil potential of the basin. Oil yields of core samples (135 sets) were regarded as reliable bases for appraising the Mahogany-zone interval, but only 5 sets of cores sampled the deeper, lower-zone oil shales. Drill cuttings (68 sets) sampled the lower-zone oil shales, but their oil yields are principally indicative and should be confirmed by analyses of cores. More complete core sampling of the lower-zone oil shales is needed.

Jan 1, 1967

By J. W. Givens, P. B. Crawford

A potentiometric model study has been made of the effect of vertical fractures existing in the matrix of the reservoir on the flooding or cycling performance. Fractures can have unusual flow characteristics. Fluid entering one side of a fracture can emerge on either the same or opposite side of the fracture, depending on the particular streamline. Fracture orientation has a great influence on the sweep efficiency. However, sustained fluid injection may still permit large areas to be swept. INTRODUCTION One can find considerable information available in the literature on the effect of fractures originating at the wellbore on sweep efficiencies of water-flooding or gas cycling programs.1-5 However, there arc few, if any, quantitative data reported in the literature on the effect of vertical or horizontal fractures existing out in the reservoir matrix on secondary recovery performance. Several papers have been presented on the effect of vertical and horizontal fractures on the productivity or conductivity of waterflood and gas injection patterns, but in all cases the fractures initiated at the well and simulated commercial fractures.6,7 It is believed that natural fractures may exist out in the reservoir matrix and would effect the displacement performance. The purpose of this report was to study the effect of a few isolated vertical fractures existing out in the reservoir matrix on the performance of waterflood or gas injection patterns. DESCRIPTION AND EQUIPMENT Since the potentiometric model was described by Lee8 it has been widely used to study numerous types of fluid displacement problems. The potentio- metric model can be used for fluid displacement studies when the following assumptions are valid: (a) steady-state conditions exist; (b) the mobilicy ratio is one; and (c) the capillary and gravitational effects can be neglected. Five-spot and direct line drive square patterns were studied using a copper strip of the desired length and orientation to simulate a vertical fracture in the reservoir matrix. The 20 x 20-in. model was considered to represent only one element of an infinite array of similar patterns. DISCUSSION AND RESULTS FIVE-SPOT PATTERN — ONE FRACTURE BETWEEN INJECTION AND PRODUCING WELLS Fig. 1 shows a quadrant of a five-spot pattem with the vertical fracture existing out in the matrix along a line connecting the injection and producing wells. Length of the fracture was equal to 35.4 per cent of the distance between injection and producing wells. In this particular pattern it is noted that fluid breaks through from the injection well to the producing well when the dimensionless time is equivalent to 26.9. Dimensionless time is defined as volume of fluid injected divided by the volume of displaceable fluid in the pattem expressed as a per cent. The cross-hatched areas of this

Jan 1, 1967