Search Documents

By Harold L. Overton, Leonard B. Lipson

The first paper in this series1 outlined practical methods for applying the theory of steady-state flow of an ideal Bingham plastic liquid through a circular pipe and axially through a stationary concentric annulus to the enginekring analysis of friction loss problems. In the present paper the properties of another useful rheological model, the pseudoplastic generalized Newtonian liquid, are considered. The terminology applied to the rheological models is derived from the following classification of rheological models. The term "generalized Newtonian" is a class designation applied to models which do not have yield points but which exhibit a dependence on shear rate. The term "pseudoplastic" applies to a sub-group of models within this class which exhibit a decrease in "viscosity" with increasing shear rate. The term "power model" refers to a mathematical model proposed for describing the behavior of the pseudoplastic liquid. This model is of interest since certain of the salt saturated, oil emulsion, and "low solids" drilling fluids, inverted emulsion- and oil-base drilling fluids, aqueous gels, gelled crudes and blocking agents employed in hydraulic fracturing operations are of this type. The present treatment considers the equations describing steady-state Poiseuille flow through a circular pipe and a stationary concentric annul us, and Couette flow between concentric rotating cylinders for a particular pseudoplastic model. As in the previous paper' it is demonstrated how these equations can be applied to engineering calculations of friction losses. Graphical procedures for recovering the various flow curves and cxample calculations are included. A method for simplifying turbulent-flow correlations is also proposed. PROPERTIES OF PSEUDOPLASTIC GENERALIZED NEWTONIAN LIQUID In contrast to the Bingham plastic the generalized Newtonian liquid is yield point-free. However, in the flow region the relation between stress and rate of shear is not linear as it is with ordinary Newtonian fluids. If the stress increases with shear rate at a less than linear rate the system is called pseudoplastic. While there is no yield point in this system, the Bingham advocates' coined the term pseudoplastic to distinguish this kind of behavior from Bingham plastic behavior. For the pseudoplastic a plot of flow rate vs friction loss or rpm vs torque, from pipe flow or rotational viscometer measurements, respectively, generally has three qualities: (1) it can be clearly extrapolated toward the origin; (2) there is very little slope at the origin. This accounts for the high apparent viscosities at low mixing and pumping rates; and (3) the flow curve possesses a continuous curvature which gradually diminishes as the shear increases. Rheological measurements on a pseudoplastic liquid can be easily misinterpreted as indicating Bingham plastic behavior. A typical case in point is illustrated in Fig. 1. These data were obtained from measurements on an oil-emulsion mud using a commercially available rotational-type viscometer (Model 35 Fann V-G meter). Applying the "two-point"* method3 gives a plastic viscosity of 17 cp and a yield point of 24 lb/100 ft2; at the same time weighting material settles very easily in this system. This illustrates that it is not possible to distinguish between the pseudoplastic and the Bingham plastic from only two data points. It is possible to determine qualitatively which model is more appropriate by a combination of dynamic and static measurements. Experience indicates that the simplest criterion for the Bingham plastic liquid is that the calculated dynamic yield point should be of the same magnitude as the measured gel strengths. For the oil-emulsion drilling fluid illustrated in Fig. 1, a value of 4 lb/100 ft2 was obtained for gel strengths measured

By Rongdong Deng

A tower-like magnetic system comprising eight permanent magnets was installed at the center of a flotation column to separate chalcopyrite and pyrrhotite minerals. A cyclonic flotation was installed within the underpart of the column to separate the pyrrhotite from the gangues. A three-product magnetic flotation column was fabricated to create a magnetic field and a cyclonic zone. The products were Cu concentrates, pyrrhotite concentrates and tailings. The theoretical magnetic field strengths and the magnetic forces were calculated. The results of the flotation experiments showed that the Cu grade increased from 14.20% without the magnetic field to 26.40% in the presence of the magnetic field and that the S grade in the pyrrhotite concentrates increased in the presence of cyclonic flotation. The separation efficiency of this magnetic flotation column was evident. Therefore, this three-product magnetic flotation column may have a potential application, due to its high separation efficiency.

By John J. Sammarco

The functional safety of programmable electronic (PE) mining systems is an international issue and concern. From 1995 to 2001, 11 PE-related mining incidents in the U.S. were reported by the Mine Safety and Health Administration (MSHA); 71 PE-related mining incidents were reported in Australia. MSHA does not have regulations for formal evaluations of the functional safety of PE mining systems. Hence, the National Institute for Occupational Safety and Health (NIOSH), in partnership with MSHA and the industry, generated the NIOSH safety framework for functional safety of PE mining systems. An overview of the NIOSH framework is given; the key framework elements, the safety life cycle and safety integrity levels are detailed. The safety framework approach has impacted the national and Australian mining industries by enabling the industries to advance from an ad-hoc approach to a formalized and systematic functional safety process. In retrospect, valuable lessons were learned for addressing functional safety and for changing industry perspectives and practices. These lessons continue to benefit mining and are applicable to other industries as well.

By Sean Gallagher, Thomas G. Bobick, William S. Marras

Twelve healthy, experienced underground cool miners performed lifting capacity tests in stooped and kneeling postures using a modified psychophysical procedure. Subjects adjusted weight in a lifting box to the maximum they could handle without undue fatigue in an asymmetric lifting task. Lifting periods were 20 min in duration and the frequency was 10 lifts/min. Tests were performed under a 48-in. roof that restricted the subject&apos;s posture. Psychophysical. physiological, and biomechanical dependent measures included the maximum acceptable weight of lift (MAWL), heart rate (HR), rate of oxygen consumption (VO2), ventilation volume rate (VF), respiratory exchange ratio (R), and integrated electromyography (EMG) of eight trunk muscles. Results indicated that the MAWL was significantly lower when kneeling than when stooped (p < 0.05), Furthermore, metabolic demands were greater in terms of HR (p < 0.005) mid VO2 (p < 0.05) when kneeling. Left and tight erectores spinae muscles exhibited increased EMG activity- in the kneeling posture (p < 0.001). It was concluded that psychophysical lifting capacity is decreased, while the metabolic stress and internal load on the spine are increased, in the kneeling posture. Results of this Bureau of Mines study indicate that it may be advisable to reduce the weight of materials that are handled repetitively in the kneeling posture.

By T. E. Marshall, D. R. Dolinar, T. M. Barczak, T. P. Mucho

Over a portion of the older limestone mine workings at the NIOSH Lake Lynn Laboratory, a mining health and safety research facility, a large sinkhole formed caused in part by the intersection of several tightly spaced joint sets and unusual weather conditions. The anticipated propagation of the associated roof failures threatened to encroach upon a portion of the facility’s secondary escapeway. To protect the escapeway, two cribs walls and over 80 cable bolts were installed. In addition, extensometers were installed to monitor roof movement in the escapeway. However, a large roof fall associated with the propagating roof failures in the sinkhole area partially destroyed one crib wall where the adjacent escapeway was not reinforced by cable bolts. Recent roof movement information indicates that this area is still active. Based on these measurements, steps have been taken to stabilize this portion of the escapeway. This paper discusses roof monitoring, the roof movement, the large roof fall, propagation of the roof failures, and the recent support measures undertaken to stabilize the escapeway.

By Joseph C. Zelanko, Christopher Mark

During August, two eastern Kentucky coal miners died in a roof fall during retreat mining operations. They had just finished extracting one pillar block and were moving to the next one when the intersection they were standing in collapsed. Regrettably, this was not an isolated incident. Since the end of 2000, there have been 28 ground fall fatalities in U.S. coal mines. Of these, 10 occurred during pillar recovery operations. Although pillar recovery has been associated with about one third of recent ground fall fatalities, it accounts for only about 10% of U.S. underground coal production. Statistically, a coal miner engaged in pillar recovery is several times more likely to be killed in a ground fall than a miner on an advancing section of a longwall. Pillar recovery is dangerous because it creates an inherently unstable situation. Once the pillars are extracted, the roof is expected to cave. Safe pillar recovery does not mean preventing roof collapse, it means ensuring that it only occurs after the miners have completed their work and have left the area. Fortunately, there are a number of proven techniques that can greatty reduce the hazards of pillar recovery. These include better mine planning, improved roof support, and safer work procedures. The goal of this article is to describe the most significant "risk factors" associated with pillar recovery, and prescribe control techniques that, taken together, can reduce the overall risk to miners. The conclusions are based on extensive studies conducted by the National Institute for Occupational Safety and Health (NIOSH) and the Mine Safety and Health Administration (MSHA) that have included detailed analysis of every fatal retreat mining incident since 1992 (25 incidents resulting in 30 fatalities). Risk reduction strategies for pillar recovery can be divided into three main groups: Global Stability: Prevention of section- wide pillar failure; Local Stability: Prevention of roof falls in the working area; and Work procedures and worker location: Minimizing exposure to hazardous areas.

By William Adamson

Precise and data driven methods of predicting and modelling blast induced vibration intensity are created and applied under field conditions on a routine basis today. These models are fed with controlled field measurements of vibration pulses generated by controlled charges under known geometric conditions. With such models blasting and geomechanical engineers can estimate vibration intensity at established critical points in mine bench geometry (surface mining) as well as within semi-permanent crown and rib pillars (in underground mining).Damage to bench walls and destabilization of large-scale geomechanical structures can be partially attributed to high intensity blast induced stresses on the rock mass and structures that define important features such as safety catch berms and haul roads. In the underground environment, weakening of critical pillars in block caving or Long Hole Open Stoping(LHOS) operations may lead to ore sterilisation. Blast vibration intensity can be linked with physical alteration in critical rock mass volumes within themine. Often the potential for vibration damage is estimated by the calculation of a critical value of vibration intensity, derived from a consideration of laboratory measured rock matrix properties. Efforts to confirm the accuracy of these estimates of critical vibration level include post blast excavation photography as well as the visual examination of rock mass, pre-and post-blast, using video camera observation of the walls of purpose drilled observation holes. These methods clearly indicate a change in fracture state due to the blast induced vibration, however further modelling of geomechanical stability will benefit if seismic measurement data can be interpreted in terms of geomechanically relevant strength parameters. This paper describes some of the results of previously unpublished research studies conducted in underground operations in Chile and Australia, where cross-hole seismic measurement profiles were interpreted in terms of dynamic modulus. This work dates from a period between 1987 and 1990 and the Author suggests that 30 years on, applying modern day technology, instrumentation and modelling methods, more meaningful, practical and applicable results may be obtained.

By A. J. Duncan, Bath M. D., E. R. Rudolph

Some factors influencing gold recovery by gravity concentration by M. D. BATH,. M.A.Se. (Met.) British Columbia (Visitor) A. J. DUNCAN,. B.Sc. (Hons.) Rhodes, M.S.A.C.I. (Visitor) E. R. RUDOLPH,* M.Se. (Eng.) Rand (Member) Recent studies have shown that the overall extraction of gold from Witwatersrand ores is increased as the proportion of gold reporting to the amalgamation circuit is increased. Testwork has therefore been carried out with a view to improving recovery by gravity concentration. In this paper, the history of gold recovery by gravity concentration is reviewed, and methods and machines at present in use are briefly described. The results of plant and laboratory tests are presented, and the limitation of existing concentrators is discussed. A pilot plant is being built for the testing of improvements to the existing concentrators and of several new devices, particularly the Reichert-cone concentrator.

By L. L. Handy

The most common method of identifying hydrocarbon-bearing strata in a well that penetrates many different formations involves measurement and interpretation of the electrical properties of the formations as determined by electrical logs. Even though this method is used extensively, and even though in a great many instances it is capable of indicating presence of oil or gas, situations arise for which it is extremely difficult, if not impossible, to deduce the presence of hydrocarbons. These situations may involve the following. 1. A thin formation, bounded by highly resistive formations, in which it is impossible to obtain the actual resistivity of the uninvaded zone with existing logging devices. 2. A formation in which invasion has been so extensive that a value for the uninvaded zone resistivity cannot be obtained. 3. A very shaly formation in which the resistivity index, I, is lower than that usually associated with productive formations. 4. Laminated formations comprised of thin productive sands separated by thin shale streaks in which the individual sand and shale streaks are too thin to permit measurement of uninvaded-zone resistivity with existing logging devices. 5. Productive formations in which the water saturation is high. To extend the utility of electric log interpretation to identification of hydrocarbons in all types of formations, there is strong incentive to find a method not subject to these limitations. Some time ago, in connection with research on the wettability of reservoir rock, an investigation was conducted in which the resistivities of cores were measured shortly after they were removed from a core barrel, and again after they had been extracted and restored to their original oil and brine saturation.&apos; The resistivities after extraction were generally lower. Other tests made on the cores indicated that they were more nearly water wet after they were extracted; thus, it was assumed that the observed changes in resistivities were due to a change in wettability of the cores. Other experiments&apos;," have shown that resistivities of rock samples are sharply dependent on wettability. These experiments have shown that oil-wet samples are more resistive than water-wet samples. To obtain an understanding of how the wetting properties of the surfaces of core material affect electrical resistivity, a series of experiments was conducted. Two groups of core samples were prepared for testing. One group contained brine, but no residual oil. The other group was saturated with brine, flooded with oil to a low water saturation, then flooded with brine to a final residual oil saturation. Resistivity measurements were made on each group. Both groups were then flooded with the original brine to which a chemical had been added that renders sand and clay surfaces preferentially oil wet, a so-called reverse-wetting agent. very little change in resistivity was observed in cores containing only water. The group containing residual oil, however, showed resistivity increases of 100 to 200 per cent. These experiments showed that the resistivity of a core containing oil could be altered by changing wet-tability of the core. Moreover, the possibility was introduced that reverse-wetting agents might be employed as the basis for a logging method for identification of oil-bearing strata. Since behavior of a porous rock containing gas and water might be expected to be similar to that of a rock containing oil and water, such a method should also be applicable to identification of gas-bearing zones. In principle the wettability of the invaded zone could be reversed without altering conductivity of the interstitial water or the hydrocarbon saturation therein. Those strata showing significantly increased invaded-zone resistivities would, therefore, contain hydrocarbons; those with no significant change would be filled only with water. Addition of a reverse-wetting agent to a hydro carbon -bearing zone which is, by nature, already preferentially oil wet would not result in an enhancement of its resistivity. It is generally believed, however, that most hydrocarbon-bearing strata are preferentially water wet.

By Jianxin Zhang

A process, in which arsenic in the residue is transformed directly into scorodite and copper dissolves into the solution simultaneously, was studied. The process was carried out in acidic ferrous sulfate solution under oxygen pressure. The optimum conditions of the process were determined as follows: the initial H2SO4 concentration of 10 g/L; liquid/solid ratio of 10; Fe/As molar ratio in the feedstock of 1; the addition of calcium lignosulphonate was 1% of sludge; oxygen pressure, temperature and reaction time were 1.5 MPa, 150 °C and 6 h, respectively. Under these conditions, arsenic transformed to crystalline scorodite in size of 5–10 μm, and the transformation efficiency reached 97.60%; meanwhile, the copper was leached with a high efficiency of 93.34%. The stability of the leaching residue was evaluated by TCLP testing, and the results indicated the concentration of arsenic in the leaching solution was extremely low as 0.04 mg/L and met the standard of safe stockpile.

By F. W. Jessen, N. Lamont

The effect of an existing fracture or joint plane, which may exist in a rock, on the extension of a hydraulically induced fracture through the rock has been investigated in the laboratory. By use of a series of models made from various outcrop rocks, the width and orientation of existing fractures do not alter the extension or direction of the hydraulic fracture. The results, which conform quite well with the Grifith theory of failure, are illustrated by a series of photographs of the various models. INTRODUCTION The process of hydraulic fracturing has been widely used in the oil industry since its introduction in 1948 and has made possible the production of many reservoirs which would have been uneconomical prior to this process. Numerous studies of the mechanics of the process2-h nd of the effects of oriented fractures on recovery&apos;.&apos; have appeared in the literature. The increased recoveries led to attempts to predict the orientation of hydraulic fractures at the wellbore and to the development of methods to control this orientation. Field tests have been devised"." which indicate the validity of theoretical predictions of fracture orientation at the borehole. The direction of extension of hydraulic fractures from the borehole has not received much attention since most studies have predicted, at least implicitly, that a fracture will continue in the same plane unless a change of state of stress in the rock occurs. Since such predictions are based on the assumption of rock homogeneity, the effects of rock heterogeneities are left unknown. Most sedimentary rocks are composed of layers which reflect the changing depositional conditions of geological time. In addition, the more competent rocks frequently are fractured and jointed as the result of structural deformations or tectonic movements. Evidence of such joint systems exists in many surface outcrops, and it can be assumed that similar systems occur in many subsurface rocks, although the individual joints may be rather tightly closed due to overburden forces. This paper investigates the effect that an existing fracture in a rock would have on the direction or orientation of an extending or advancing hydraulic fracture when it intersects the existing fracture. The study was conducted on small rock models under triaxial stress conditions. The results indicate that an existing fracture will have little effect on the hydraulic fracture. THE MODEL STUDY The rock model represents an elemental portion of the earth containing an existing fracture which is located at some distance from a borehole from which hydraulic fracture is being extended. A generalized view of a model is presented in Fig. 1. The models were constructed from cement blocks and natural rocks. The various materials tested and some of the physical properties are listed in Table 1. The rocks were cut into rectangular blocks with dimensions of 11/2 X3Y2 X4 to 8 in. An initial slot was cut into the end face of each block along the longitudinal axis normal to the top face. The slot was filled with plastic aluminum along its outer edge to a depth of 1/16 in. A hole of 3/8-in. diameter was drilled into the sealed slot at the center of the end face to serve as an entry port for the fracturing fluid. Existing fractures of various types were simulated in the models. A "hairline" fracture with essentially zero width was created by cutting the model into two parts with a diamond saw and replacing the two parts together along the cut. A finite-width fracture was created by placing a layer of sand grains between the two faces of the cut just described (see Fig. 12). A large open fracture (similar, perhaps, to a vug) was created by cutting the model in two parts and removing a portion of the downstream block to a depth of % in. See Fig. 13. For a fourth

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

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

By Persifor Frazer

There is always a difficulty in assigning limits to the scope of a science or art, and among the dificulties of that kind besetting practical or stratigraphical geology is that of determining what kinds of natural objects shall be considered as the goals of its search. The term mineral, in its broader sense, includes every solid, liquid, or

By C. H. Mathewson

MICROSCOPIC rnetallography has been exploited quite well enough to bring about a very general understanding that the typical metal or alloy is composed of minute crystalline particles blended into a coherent microstructural mosaic. One does not have to be a specialist in metallography to realize that the properties of such an aggregate are essentially a summation in appropriate form of individual effects derived from the shape, size, placement or orientation and cohesive characteristics of the component particles. It is clearly of great importance to consider the various forms of discontinuity which may occur at the boundaries between these crystalline particles. When a number of particles each possessing the same orderly arrangement of atoms are brought together into a close-fitting system of purely haphazard contacts, not unlike a handful of snowflakes compacted into a snowball, the particles are said to possess random orientation. There are no generally accepted views concerning the arrangement of the atomis or the constitutional reaction between atoms where one crystal meets another, but these contact regions are characterized by strength rather than weakness and it is customary to require the presence of many rather than few boundaries in preparing metal for useful service. Under certain circumstances all of the particles in a metal may be nearly alike in orientation. Other general tendencies of orientation may be associated with particular forms of mechanical and thermal treatment. In contrast with these cases of fortuitous or statistical diversity of orientation, we often find adjacent particles united along a plane which possesses a grouping of atoms belonging equally well to both structures. This, of course, determines a fixed relationship between the two orientations and it is always possible to derive one from the other by some form of rotation or reflection prescribed by the symmetry of the crystal structure under consideration. Particles united in this manner are known as twin crystals although the term refers to the form of association rather than the number of individuals concerned.

By R. W. Snyder, H. J. Ramey

This paper presents the results of applying the Buckley-Leverett&apos; displacement theory to petroleum reservoirs consisting of a finite number of layers. The layers are assumed to communicate only in the wellbores, and the reservoir may be represented as a linear system. Most previous investigations of this nature were limited by assumptions and by inconsistent calculation techniques. This study improves on previous work by applying the Buckley-Lev-erett displacement theory to a noncommunicating layered system where permeability, porosity, initial saturation, residual saturation and relative permeability vary from layer to layer in a logical and consistent manner. Gravity and capillary-pressure effects are neglected. A modification of the Higgins-Leighton calculation method was used in this study. Waterflood predictions were made with all properties varying, and then with only permeability varying using several inability ratios. These results were compared with the Stiles and Dykstra-Parsons predictions. It is shown that the latter methods generally give poor values for the breakthrough recovery and pessimistic predictions for the performance after breakthrough. Similar results were obtained for a gas-displacement case. lNTRODUCTION Field experience with immiscible displacement usually shows constant producing conditions until breakthrough of the displacing fluid. Then oil production continues at increasing displacing-to-displaced fluid ratios until the economic limit is reached. Three different ideal mechanisms are known that will produce this behavior: (1) relative permeability effects as described by Buckley-Leverett frontal advance theory,&apos; (2) vertical stratification as considered by Stiles,2 Dykstra and Parsons5 and others and (3) different path lengths involved in areal (two-dimensional) flow between wells as described by Dyes et al.4 Without question, a combination of these factors modified by formation heterogeneity and other known and unknown factors actually does control the behavior of real systems. This paper presents results of an investigation of certain factors that should affect performance but which have received little attention to date. In 1944, Law5 demonstrated that porosity and perme- ability are often found to have normal and logarithmic-normal distributions, respectively. throughout cored intervals in natural formations. This led to the concept of the noncommunicating, multilayered reservoir model for immiscible displacement. This model assumes that the reservoir is composed of a number of layers that communicate only at the wellbores. Each layer is individually homogeneous, but may be different from every other layer. Stiles&apos; presented one of the earliest applications of this model to waterflood performance. In addition, Stiles assumed that the initial saturations and relative permeabilities were the same for each layer, porosity was the same. displacement was piston-like, fluids were incompressible and injection into each layer was proportional to that layer&apos;s permeability capacity (permeability-thickness product). The last assumption would be true if the mobility ratio for the displacement were unity.21 Dykstra and Parsons" used the same model as Stiles, but rigorously included mobility ratios other than unity for piston-like displacement. Dykstra and Parsons used their general result to produce charts for log-normal permeability distributions between layers. Similarly, Muskat6 Pub1ished analytical solutions for linear and exponential permeability distributions. In 1959, Roberts&apos; described a scheme for calculating water-drive performance for the noncommunicating, layered reservoir model which considered two-phase flow in the displaced region. Roberts used the same model and assumed that the injection rate into a layer was proportional to that layer&apos;s permeability capacity, but that flood front locations could be evaluated from the Dykstra-Parsons results. These assumptions are inconsistent, and a material balance cannot be maintained except for a mobility ratio of unity. At the same time, Kufus and Lynch8 coupled Buckley-Leverett displacement theory with the layered model to provide an improvement of the Dykstra-Parsons method that was consistent. In 1960, Higgins and Leighton9 resented a numerical method for calculating waterflood performance also considering two-phase flow in the displaced region. The result was used to investigate variation in absolute permeability and oil viscosity. An excellent, detailed history of using the noncomrnunicating, layered reservoir model was presented by Nielsen.&apos;" The preceding techniques (and many related ones) were similar in that differences in initial saturations, residual saturations and relative permeabilities from layer to layer were neglected. It is well known that the irreducible water saturation is an important function of absolute permeability. Calhoun11 showed that the irreducible water saturation

By A. M. Skov, L. F. Elkins

Inclusion of anisotropic permeability in mathematical analysis of pressure transients observed during development of the huge Spraberry field indicates a major fracture trend which is in good agreement with that observed by fluid-injection tests spread over a 12- by 17-mile area. Delineation of this trend is important in selecting a pattern of injection for the pending large-scale water flooding in this field. Determination of reservoir parameters yielding best agreement between calculated pressures and observed reservoir pressures in newly completed wells was made using an IBM 650 computer. INTRODUCTION The Spraberry field covering 400,000 acres is a tight sand of less than 1-md permeability cut by an extensive system of vertical fractures. Primary recovery dominated by capillary retention of oil in the fractured sand matrix blocks is less than 10 per cent of oil in place. Strong forces of capillary imbibition of water into the sand, coupled with water flow under dynamic pressure gradient, indicate considerable increase in oil recovery can be achieved through water flooding. Best results will occur if the pattern of water injection is selected to force the water flow across the grain of the major fracture system. Existence of an oriented vertical fracture system in the Spraberry, observed first in cores, was highlighted more recently by the 144-fold contrast in permeability along and at right angles to the major fracture trend required to match relative water breakthrough times in Humble Oil & Refining Co.&apos;s waterflood test there. Spraberry operators since have conducted two gas-injection tracer tests for further areal confirmation of the fracture trend. Re-analysis of early reservoir pressure transients for evidence of anisotropic permeability has permitted many more local determinations of major fracture trend without resort to further field tests. This paper is limited to updating analysis of reservoir pressure transients to include anisotropic permeability as a test for orientation of the major fracture trend in the Spraberry. The reader is referred to Refs. 1 and 2 for information about general Spraberry reservoir performance and to Refs. 3 and 4 for information about significance of fracture orientation in selection of the injection-well pattern for water flooding the Spraberry. RESERVOIR PRESSURE DATA—DRIVER AREA During early development of the Spraberry Driver area, Sohio Petroleum Co. made the extra effort to measure the initial pressure in each of the 71 wells in a 5-mile-long area immediately after completion. Progressively greater reductions in pressure ranging up to 400 psi were observed throughout the six-month development period. Detailed data are presented in Ref. 1. Since the reservoir oil was undersaturated some 300 psi initially, early reservoir performance involving 55 new well pressures is subject to analysis as flow of a single compressible fluid in a porous media. Assumption of uniform permeability in all directions yielded good agreement between calculated pressures and observed pressures of these wells in the earlier study,&apos; but subsequent, additional, mathematical development to include anisotropic permeability in the transient pressure considerations and present availability of electronic computers to perform the much more extensive arithmetical calculations now yield even better agreement. The previous analysis, assuming uniform permeability, consisted essentially of calculating pressure reduction expanding circularly around each producing well and summing these effects at the time and location of each newly completed well for comparison with the measured pressure reduction. Permeability, effective fluid and rock compressibility, and permeability X thickness were varied until the best match with measured pressures was obtained. The present analysis, assuming anisotropic

A Colloquium and General Meeting were held on 17th March 1971, the theme being &apos;Future trends in steel production&apos;. Mr V. C. Robinson (President) was in the Chair. The Colloquium was attended by 121 delegates and was opened by the President at 10 a.m. OBITUARIES THE PRESIDENT: &apos;It is my sad duty to announce the death of the following: F. O. Read, Fellow, who joined the Institute in 1957 and passed away on 13th January, 1971. P. G. M. May, Member, who joined the Institute in 1969 and passed away on 10th August, 1970. R. J. J. Rademeyer, Life Fellow, who joined the Institute in 1942 and passed away on 23rd September, 1970. As a mark of respect to the memory of the deceased and in sympathy with the bereaved I would ask you all to rise and observe a few moments silence. CONFIRMATION OF MINUTES THE PRESIDENT: &apos;May we confirm the minutes of the monthly general meeting held on 21st October and 18th November, 1970, as published in the Journal&apos;. Agreed. MEMBERSHIP THE PRESIDENT: &apos;I have much pleasure in announcing that the names of the undermentioned candidates, having been published in accordance with ByLaw 5.2.2, Council has elected them to membership of the Institute in the following grades: FELLOWS: Romain August Lathioor, Paolo Piga, Renato Ribacchi, Willem Frederick Jacobus Smith, Matthew Churchill, Thomas Anthony James Braithwaite, Hugh Patrick Hart, Eleanor Innocentius Muller, David Alton Smith. MEMBERS: Douglas Adendorff, Petrus Jacob du Plooy, James Valentine Grant Middleton, Karl Anders Ullerstam, Josephus Theodorus Johannes van Wyk, Frank Douglas Abbott. GRADUATES: David Mervyn Gilbart-Smith, John Philip Wilcocks Bennie, Hendrik Theunis Burger, Johann Wilhelm Coetzee, Johannes Stephanus StOols, Leon van Tonder, Brian George Harvey, John Kennedy McLean. ASSOCIATES Richard Cyril Lee Stoyell, Terrence Brian Gouws, Alister Rodney Frederick MacDonald. STUDENTS: Andries Groenewald, Christiaan Johannes Hattingh, William Hendrik Hofmeyr, Stanley Thomas Vincent. I welcome the new members to the Institute and congratulate them on their election. COLLOQUIUM The President welcomed everyone present and mentioned that this was the first of the Colloquia which are to replace the Institute&apos;s monthly meetings. He stressed that discussion at Colloquia would be completely informal and would not be recorded unless the contributor requested publication. In conclusion, he thanked the organising committee, the authors and the contributors, and introduced the author of the first paper. The following papers were presented:- &apos;The future development of the Iron and Steel Industry&apos; by J. P. Coetzee (Published in the Journal March 1971). Among those who contributed to the paper were Dr Bleloch, Dr K. Gebhard, Dr Way and Professor Müller. The main point made by Dr Bleloch related to the reserves of coking coal and of coal for power generation this in country and he deduced that in South Africa we are critically short of coking coal. On the basis that 4 million tons per year of coking coal are consumed for our present steel production, over 16 million tons would be required to meet the anticipated expansion in steel output by A.D. 2000. This must be set against the grim forecast made by the Coal Advisory Board in 1967 that our national reserves of coking coal will be exhausted before A.D. 2000. This gives tremendous significance to the development work presently being carried out by Iscor on alternative methods for coke production. On the other hand the blast furnace may have to be replaced by an alternative process for the reduction of iron ores, e.g., the use of rotary kilns followed by electric arc smelting. Dr Bleloch also expressed grave forebodings about the reserves of bituminous coal, pointing out that with present methods of mining the recoverable coal in any given reserve is in general only 50 per cent of that reserve. He stated that in the past seventy years great quantities of coal have been forever placed beyond the reach of mining by being undermined by extraction of lower seams or lost by being burned. A final point made by Dr Bleloch was the fast rise in the capital investment in a fully integrated steel plant, the figure having risen in the past decade from just over R200 to R350 per ton of finished steel per annum. Dr Gebhard spoke of the direct reduction processes for iron making pointing out the difficulties of securing a high percentage of reduction and the relatively high cost of smelting in electric arc furnaces. As an alternative he suggested briquetting the high carbon sponge iron

By J. P. Proulx, Richard O’Meara, Ruilin Yang

A blast optimization project was carried out at an Open Pit Metal Mine using the Multiple Blasthole Fragmentation (MBF) Model and the Multiple Seed Waveform (MSW) blast vibration model. The aim of the project was to optimize rock fragmentation of the ore blasts to significantly increase metal recoveries while keeping the blast vibration levels at the highwall under control. The project also examined the blasts in the waste rock using the Fragmentation model to expand the blast pattern and select economical explosives for drill and blast cost savings while minimizing any adverse effect on the dig rates. Two sets of near field signature blasthole tests were conducted in both hard ore and relatively soft waste rock. The information from the near-field signature blasthole vibration was used as input to both the MBF and MSW models, as site specific geological information interacting with blasting. The implementation of the full blast design is also input to the MBF and MSW models. A production blast for each rock type was monitored for blast vibration and rock fragmentation. The full blast information and the blast results are used to calibrate and validate the models. The Fragmentation and Vibration model were established for the two rock types at the mine. At each site the Fragmentation modeling was conducted in parallel to the Vibration modeling to ensure blast vibrations were under control while the rock fragmentation was improved. The Fragmentation MBF and MSW are each an integral part of the blast optimization tool. Various scenarios using different explosives, hole delays and blast patterns were modeled for rock fragmentation and blast vibration. The recommendations from the modeling project were to be field trialed at the mine. The implementation of the recommendations was expected to help the mine significantly in terms of the mining economics, and high wall stability and safety.

By W. S. Ayres

The boiler from which these stay-bolts have just been obtained was that of the locomotive Catasauqtia, Lehigh Valley Railroad, built at the company&apos;s shops, South Easton, Pa., in 1864. The iron is Lowmoor, and has been in use nearly ten years. The stay-bolts have become highly crystallized from some cause, and are very brittle, only requiring three blows of an ordinary machinist&apos;s hammer to break them. On the other hand, when annealed they are remarkably tough, as shown by the doubled sample. The majority of the broken bolts were found to be along the curve of the shell from a to c, Fig. 6. All of them were broken just inside the shell, and in no case next the fire-box. The fractures seem to have been of slow procedure, beginning at the top of the bolt, gradually working deeper and deeper until the pressure of steam overcame the tenacity of the remaining parts. Unequal contraction and expansion of the parts of the boiler seems in this case to be the only existing cause for a fracture. We cannot, from the data obtained, prove to an absolute certainty the manner in which the combined forces generated by the unequal contraction and expansion of the parts have arrived at their resultant work ; but we can, from these data and from existing circumstances, arrive at what seems to be a fair supposition, and thus at least anticipate the truth. The contraction and expansion in locomotive boilers has a large range between its extremes, and is very variable, both with reference to time and to their component parts. The shell and fire-box do not, perhaps, at any time, contract and expand at the same rate, nor proportionally; but, on the contrary, owing to the great changes in temperature that take place both in the atmosphere around the shell, and in the fire in the fire-box, their condition is often changed suddenly, and widely from that which they just previously had. This racking tendency produces strains in the different parts, and it is fair to suppose that the continual motion of the two

Today’s world relies upon advanced green technologies that are made of critical elements with unique properties. Examples include rare earth elements (REEs) that are critical in the manufacturing of wind turbines and electric/hybrid vehicle batteries, key components of a greener future that face supply uncertainty and near zero recycling. To tackle their supply challenges, activities have begun for urban mining from waste electrical and electronic equipment (WEEE) that contain considerable amount of REEs, but their current level of recycling is less than 1%. Current recycling practices use either pyrometallurgy, which is energy intensive, or hydrometallurgy that utilizes large volumes of acids and solvents and generates large volumes of hazardous waste. In this study, we developed a novel and sustainable process to recycle REEs from WEEE, NiMH batteries in particular. The developed process relies on supercritical fluid extraction (SCFE) utilizing CO2 as the solvent, which is abundant, safe, cost effective, and inert. The effect of seven operating parameters, namely temperature, pressure, residence time, sample to chelating agent ratio, agitation rate, complex formulation, and methanol addition on REE extraction efficiency was investigated, and optimum operating conditions were determined. This work is the first to utilize a very efficient and safe process that runs at low temperature to extract metals from postconsumer products with minimum hazardous waste generation, while offering about 90% leaching efficiency for REEs. We expect our process to find widespread applicability in urban mining of REEs using green chemistry.