Search Documents

By James F. Lewis

The hematite iron ore mines east of the Hudson River are confined to a strip of country ten to fifteen miles wide, commencing on the south, near Fishkill, running northeast through Dutchess County, and striking the southeast corner of Columbia County, extending into Litchfield County, Conn., known there as the Salisbury mines, cropping out again at and near West Stockbridge, Lanesboro, and Cheshire, Mass., also in Bennington County, Vermont. The ores are found in beds or deposits and in veins, with a limestone ledge on one side and slate on the other, from six to thirty feet from the surface. Portions of it are quite solid, called rock ore, but

By G. Roderick Eggert

Material criticality reflects the degree to which a raw material is both necessary and subject to supply-chain risks. This study compares the criticality of selected raw materials from the differing perspectives of the manufacturing sectors of China, the European Union, Japan and the United States.

By L. G. Jones

The objective of this paper is to present a better method for determining reservoir reserves of gas from extended drawdown test data. Means are proposed for conducting the test with or without the requirement of prior build-up to a static reservoir pressure. The test is based upon solutions to unsteady-state gas flow developed in an earlier paper and is believed to be superior to previous methods for determining reserves from flow-test data. The limitatiom of such tests are discussed, with emphasis on the subject of test duration. The test can supply useful information if the mechanics and limitations of it are understood. This type of information can be particularly valuable in deciding whether or not to drill offsets to producing wells. INTRODUCTION In the life of any gas or oil field, the engineer in charge must decide if certain wells need stimulation and whether or not to drill more wells in the field. Well production tests can aid in making these decisions. For instance, the value of stimulation can be estimated by determining formation permeability and the amount of well damage from drawdown or build-up tests. These tests are freely used in the oil industry. A less well-known type of analysis uses data from extended drawdown tests for determining reservoir reserves. Into this category fits the "reservoir limit test" of Jones,'-3 and the "reservoir reserve test" presented in this paper. The chief value of these tests is their ability to "see" reservoir boundaries and thus prevent costly dry holes. Methods for determining reserves of both gas and oil reservoirs were proposed in the work of Jones. His proposals appear sound for single-phase liquid flow. The generality of his work on gas flow is open to question, however, since the analysis is based upon solutions to liquid-flow equations. The test for gas reserves presented here is derived from unsteady-state gas-flow solutions presented in a previous paper.' In deriving this reservoir reserve test, it has been assumed that all of the reservoir energy is provided by expansion of the reservoir fluid and that the reservoir is radially symmetrical. No mention has been made of the effects of deviations from Darcy's law such as quadratic (turbulent) flow and well damage. Since the terms representing these effects are function of rate only, since the rate is constant and since the data are plotted in a differential form, it seems reasonable that these terms will not affect the test results in any way. No attempt has been made to estimate the effects of water drive, faults, or other departures from an idealized and isolated reservoir. For a qualitative discussion of some of these effects, the reader is referred to a paper by Jones3

By T. M. Barton

This paper presents several case studies in which a mechanics-based boundary-element program is used to back-calculate the surface subsidence associated with various panels at several northern Appalachian coal mines. The program used in this case study is called LAMODEL, which incorporates a frictionless, laminated overburden into a general-purpose displacement-discontinuity code primarily designed for calculating the stresses and displacements in coal mines or other thin-seam or vein-type deposits. In this paper, the program is used to calculate both the underground convergence and the resulting surface subsidence at five longwall panels and a room-and-pillar section. The fitted subsidence from the model is compared with the field measurements and analyzed. The results from this work show that the LAMODEL program is not as accurate as available empirical subsidence-predictive methods: the expected correlation between the geology and the optimum input parameters is not evident. However, for a mechanics-based program, LAMODEL does provide moderately accurate subsidence calculations, and it is one of a few programs that can even attempt to practically calculate both underground stress and convergence and the resulting surface subsidence.

By K. V. Raj, C. A. Berry, A. T. Weakley, A. L. Miller, P. R. Griffiths, D. A. Parks

Diesel particulate matter (DPM) has been classified as a carcinogen to humans by the International Agency for Research on Cancer. As a result of its potential carcinogenic nature, DPM exposure is regulated by the Mine Safety and Health Administration. Currently, diesel emissions in the workplace are monitored by collecting the aerosol onto filters, which are then sent to a laboratory for thermal-optical analysis using the NIOSH method 5040. This process can take days or even weeks, and workers can potentially be exposed to excessive levels of DPM before the problem is identified. Moreover, the delay involved in getting the loaded filter to the lab inevitably means the loss of some of the more volatile organic carbon. To remedy this delay, researchers from the National Institute for Occupational Safety and Health are seeking to develop a field-portable, real-time method for measuring elemental and organic carbons in DPM aerosols. In the current study, the use of mid-infrared spectrometry was investigated. It is believed that mid-infrared spectroscopy is more suitable for use in a real-time field-portable device than thermo-optical analysis methods. This article presents a method for measuring organic carbon (OC) and elemental carbon (EC) in DPM for a broad range of OC/EC ratios. The method has been successfully applied to laboratory-generated and mine samples.

Infacon 74 The National Institute for Metallurgy, the South African Institute of Mining and Metallurgy, and the Ferro-Alloy Producers' Association of South Africa extend a cordial invitation to participate in the International Ferro-alloys Congress to be held in Johannesburg, South Africa, in April 1974. Obituary RICHARD AUSTIN COOPER On the 31st March 1972, a once very active member of the Institute, R. A. Cooper, died at the age of 82. He had served on the Council for many years and was Vice-President in 1933/34. ALEXANDER JEREMIAH ORENSTEIN CB., CM.G., C.B.E., LLD., D.Se., M.D., M.R.C.S., F.R.CP., F.R.S.H. Dr Orenstein, who died on the 7th July 1972 at the age of 92, was so well known and so widely acclaimed during his lifetime that it is difficult to add to the tributes which have already been paid to him.

By John Anawati

Bauxite residue is an environmentally unfriendly byproduct of alumina manufacturing, produced worldwide in large quantities. This material contains 50– 100 ppm of scandium, a critical material for the production of stronger, weldable, corrosion resistant, and heat tolerant aluminum products. Aircraft manufacturers are particularly interested in Al–Sc alloys because of its ability to form weldable alloys that could potentially reduce aircraft weight by 15–20%. Because of its abundance and low cost, bauxite residue has the potential to be used as an efficient feedstock for valorization processes to recover its scandium content, while in turn degrading this problematic waste. In this project, we developed an efficient process to recover scandium and rare earth elements from a Canadian bauxite residue with high extraction yields, using industrially scalable and economical techniques. The process employs a modified version of sulfuric acid leaching and subsequent impurity removal by selective precipitation. This article outlines the current progress and design rationale in the development and optimization of this innovative and sustainable recovery process.

As a token of their gratitude, former pupils of Professor Lucien Denoel, a lecturer on coal-mining in the University of Liege, over a period of forty years, have founded the LUCIEN DENOEL prize. It is awarded, regardless of nationality and diplomas, to the author of a thesis dealing with coal-mining, ventilation, mine-draining methods, the use of explosives in the mine or coal-preparation. The thesis must show a truly scientific approach and considerable originality, and must not have been dealt with in a previous publication. The prize is 60.000 Belgian francs. The articles must be written in French, Dutch, German or English, and contain not more than 10,000 words, figures included. They must be forwarded, in triplicate, to the 'Jury du Prix Lucien DENOEL, c/o A.I.Lg., 22, rue Forgeur, Liege, Belgium,' before 1st July, 1970. The examining board reserves the right to publish the winning thesis in a Belgian magazine.

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