Search Documents

By R. E. Collins, Paul B. Crawford

Potentiometric model studies have been made of ertically fractured reservoirs. All fractures originated at he well and extended into the reservoir for various distances. Studies were made to determine the areal sweep efficiencies for line-drive patterns. The areal sweep efficiency for the unfractured system studied war 70.6 per cent. It was found that if the vertical fracture war paralIel with the direction of flood, the areal sweep efficiency may be considerably reduced, approaching zero in some instances. If the fracture was at right angles to the direction of flood, the areal sweep efficiency may be greater than that for the unfractured system. The sweep efficiency depends on the length and orientation of the fracture and direction of the flood. It is concluded that considerable effort should be made to determine the nature and direction of fractures before initiating a flooding program. Such information is required if the maximum sweep efficiency is to be obtained and channeling is to be reduced to a minimum. INTRODUCTION The recent development and application of fracturing techniques to petroleum reservoirs has served to focus considerable attention on the effect such fractures may have on secondary oil recoveries. The opinion has been expressed that fracturing may not alter appreciably the water flooding possibilities. Others have thought that channeling may occur and that the creation of fractures may serve to substantially reduce the possibilities of conducting a successful water flood in the reservoir The purpose of this study was to obtain quantitative estimates of the effect of vertical fractures on the areal sweep efficiencies of line-drive flooding patterns. In presenting this work on vertical fractures, it is not to be inferred that the present fracturing techniques result in producing a vertical fracture. The types of fracture which are created are not known with certainty in every case. Some are believed to be horizontal, others inclined along bedding planes, and some data indicate vertical fracturing. Each fracture will probably exert a different influence on the sweep efficiency of the secondary recovery program; consequently each fractured well should perhaps be considered as a separate problem, requiring a careful analysis. Preliminary work on the effect of fractures on producing capacity of a well indicated that for practical purposes many fractures may be treated as having a nearly infinite permeability This assumption appears to be valid in some instances where the permeability of the fracture is only 10 or 100 times greater than the permeability of the matrix. It is known to depend on the type, orientation and length of the fracture. In realization that fractures of only 0.01 in. width have calculated permeabilities near 5,000 darcys, it is believed that many fractures may be considered as of near infinite permeability without incurring serious error. Primarily for this reason this study was made with fracture analogs of very great permeability compared to the matrix. FRACTURED SYSTEM ANALOG The pattern sweep efficiencies reported here were obtained with the aid of a potentiometric model. The application of this type model for solving reservoir problems has been described in the literature.1,2,3 The electrolyte reservoir was 20 in. wide by 30 in. long by approximately % in. deep. The liquid depth was constant, which corresponded to a constant permeability-feet product from the injection to the production well. The electrolyte reservoir was built to correspond to a unit in the line-drive pattern shown in Fig. 1. This figure shows part of an infinite array of wells having alternate rows of fractured and unfractured wells. The area within the dashed rectangle shows the area generally studied in the model. In using the model to represent a fractured system, a thin strip of copper was soldered to the well or wells and the length and orientation of the strip was made to conform to the vertical fracture under consideration.

Jan 1, 1955

By H. Kurushima, S. Tsunoda

An investigation of the hydrometallurgy of Cu-Zn concentrates by The Dowa Mining Co. of Japan is described and a detailed flowsheet of the Kosaka hydro-metallurgical plant is given. Sufficient operating experience has shown that the Dorrco FluoSolids system has proven successful for recovery of both zinc and copper by roasting in a fluidized bed. At present, 2400 metric tons of concentrates per month are being treated and recoveries of 93 pct Cu and 65 pct Zn are being effected. A SURVEY of the Japanese chemical industry in 1949 revealed, among other things, that sulphur-bearing blast furnace gas from the Kosaka smelter of The Dowa Mining Co. Ltd., Kosaka, Japan, was being exhausted to the atmosphere as waste. In addition, it was costing that company an average of ten million yen ($28,000) annually for damages to the surrounding property caused by the sulphur dioxide in the gas. It was recommended by this survey that the gas be utilized for the production of sulphuric acid by the contact process. Mr. Kurushima visited the United States in May 1950 to investigate the engineering, construction, and operation of a contact acid plant. Any acid produced at Kosaka would necessarily have to be 98 pct H,SO, to prevent corrosion in railroad tank cars transporting it to the point of use. However, the production of 98 pct acid requires constant volume and SO, content of the gas. As the converter gas at Kosaka fluctuates in both volume and quality, it was decided to evaluate the smelting process with relation to other methods of recovery. This evaluation revealed a relatively new type of roaster, the Dorrco FluoSolids reactor, which was reported to be capable of roasting the copper sulphide concentrate in an excess of air with the production of an SO1 gas strong enough for acid production and a water-soluble copper sulphate calcine. The copper sulphate then could be leached from the roasted ore and electrolytic copper recovered from the leach solution in much the same manner as practiced by the Chile Exploration CO. for treating naturally occuring copper sulphate ores. Use of a reactor of this type would eliminate the need for a blast furnace and converter and, in addition, would effect savings in labor and fuel. The Kosaka ore is a microscopically fine mixture of copper, zinc, and lead sulphides. The lead and a small amount of the barite present can be separated by flotation, but any separation of copper and zinc proved extremely difficult. In blast furnace operation at Kosaka, zinc had been slagged off until a stockpile of about 12,000,000 tons containing 5 pct Zn had been accumulated. Mr. Kurushima felt that if any of the zinc in the concentrate was sulphat-ized and leached with the copper, it also could be recovered electrolytically. The Dorr Co., sponsors of fluidization technique in the metallurgical field, suggested that sulphate roasting be tested out on a laboratory scale and that, if successful, a commercial Dorrco FluoSolids system be constructed. They indicated that the results from such laboratory tests could be transposed to a commercial scale without pilot planting the roasting operation. It was therefore decided to send samples of flotation concentrate to the West-port Mill, laboratories of The Dorr Co., for testing. The objective of Mr. Kurushima's visit to the United States thus had changed from investigation of sulphuric acid plants to a fundamental examination of the hydrometallurgy of copper and zinc. Testing Program In February 1951, Mr. Tsunoda went to the United States to witness the application of the FluoSolids system to sulphate roasting at The Dorr

Jan 1, 1956

By J. L. Lytton, T. E. Tietz

The recovery of flour stress of recrylstallized are-cast molyhdenum after 10 pet prestrain at 150°C was studied for recovery temperatures of 725o 801o and 901oC. The yield phenomenon observed after recovery was such that the yield drop diminished and disappeared as recovery proceeded. A fractional flow-stress recovery parameter fr, was developed and it was shown that, independent of recovery ternperatzwe, equal f, values of about 0.16 produced superposition of stress-strain curves except for the initial yield phenomenon. The activation energy for low degrees oj recovery was found to he 32,600 cal per mole, in close agreement with a reported valzie for difjzdsion of carbon in molybdenum. For values of fr greater than about 0.15 where recovery was thought to predominate over strain aging, the activation energy varied from 43,600 to 82,500 cal per mole and increased with increasing recovery and recovery temperature. A 10 pet pre-strain shifted the ductile-brittle transition temperature, as indicated by tensile elongation, from 10°C. to about -22"C. A 100-hr recovery anneal at 801oC was shown to shift the transition temperature to about —5oC, the major shift occurring during the first hour of recovery. ThE recovery of cold-worked metals is accompanied by changes in the arrangement of the dislocations generated during cold work. In recent years, the use of transmission electron microscopy has allowed direct observation of some of the dislocation configurations formed, such as "tangles", networks, cells, subboundaries, and loops. These observations have recently been reviewed by Keh and Weissmann.' Accordingly, some understanding has been gained concerning the nature of the dislocation rearrangements involved. The relationship between recovery and the decrease of mechanical strength, however, has only recently been explored. Previous efforts of the authors were directed toward investigating the decrease of tensile flow stress of aluminum and Al-1 pet Mg during recovery2 and the effects of concurrent elastic and plastic strain on the flow-stress recovery of aluminum.3 In addition, an investigation of the effect of four different solute additions on the tensile flow-stress recovery of aluminum has been made.4 Many important metals, however, have the bee crystal structure, e.g., iron, molybdenum, tungsten, vanadium, columbium, tantalum, and chromium. These metals might be expected to behave differently during recovery from those with other crystal structures for two reasons. First, the previous and present investigations of aluminum and its alloys have shown that tensile flow-stress recovery is strongly affected by foreign atoms in the lattice. In the bee metals, interstitial impurities are usually present and might be expected to have various effects upon both flow behavior and recovery. For example, heterogeneous yielding might occur during prestraining. The diffusion of these interstitials can also relock dislocations during recovery so that a yield point might be obtained upon restraining (strain aging). Second, some of the bee metals (e.g., iron, molybdenum, and tungsten) undergo the familiar structure-sensitive ductility transition and exhibit brittle fracture at suitably low temperatures. It is of considerable interest to determine the effects of recovery on the transition temperature of such metals. The purpose of the present investigation was to extend our knowledge of flow-stress recovery to the bee metals and to study the effects of recovery on the ductile-brittle transition temperature. Unalloyed recrystallized are-cast molybdenum was chosen as a test material, since it exhibits a well-defined ductility transition near room temperature. As in the previous studies, the degree of recovery was measured in terms of the decrease of tensile flow stress. Prestraining and restraining were carried out at 150°C to avoid the possibility of brittle fracture. EXPERIMENTAL PROCEDURE The are-cast molybdenum sheet used in this investigation was obtained from the Climax Molybdenum Co. of Michigan. It was 0.060 in. thick by 10

Jan 1, 1964

By J. H. Scaff, H. C. Theurer, E. E. Schumacher

The microwave region of the radio spectrum was effectively utilized in radar designs during the recent war and it has become of increasing interest in the field of communications. Work in this field has led to an important use for silicon—that of the point contact rectifiers†1 in the frequency converter of microwave (radar or radio) receivers—and has stimulated considerable interest in the electrical properties and preparation of silicon and its alloys. Silicon is an electronic semiconductor. Its conductivity at room temperature results principally from the presence of certain impurities. While for metals an increase in impurity content increases the resistivity, for semiconductors such as silicon the opposite occurs and, in general, the addition of impurities lowers the resistivity. Silicon materials may be classified into one of two groups depending upon the manner in which the impurities contribute to electrical properties. These have been termed for convenience p-type or n-type. P-type silicon develops a very large positive thermal emf against metals, has a Hall coefficient of positive sign and a low resistance direction in point contact rectification with the silicon positive with respect to the point. N-type silicon, on the other hand, develops a negative thermal emf against metals, has a Hall coefficient of negative sign and a low resistance direction in rectification with the silicon negative with respect to the point. Impurities which produce n-type silicon are called donors inasmuch as these elements contribute to electrical conductivity by donating electrons to an unfilled energy band in the silicon. On the other hand, elements which produce p-type silicon are known as acceptors as these impurities contribute to the electrical conductivity by accepting electrons from a filled energy band permitting what is known as conductivity by "holes" in which the sign of the carriers appears to be positive. A general treatment of the mechanism of conduction in semiconductors from the viewpoint of modern band theory has been given recently by Pearson,2 by Becker, Green and Pearson,3 and by Torrey and Whitmer.4 In this investigation boron and aluminum have been found to be acceptors, and phosphorus, arsenic, and antimony to be donors in silicon. Data on the effect of boron and phosphorus on the electrical properties, when present singly and in combination, have been acquired. These data are discussed in the present paper. Raw Materials Used and Methods for Adding Seeond Constituents Silicon from two sources was used in this work and these will be referred to as A and B respectively. Silicon A is a material of high purity supplied by the Electro Metallurgical Co. It is prepared by chemical purification of "commercial" silicon obtained from the arc furnace reduction of SiO2. It contains 99.8 pct silicon, minimum, with small amounts of calcium, iron, aluminum, boron, and phosphorus as the principal impurities. This material was extensively employed in these studies as well as in the commercial preparation of rectifier materials. Silicon B is a material of high purity from the du-Pont Co. It is prepared by a pyrolytic reduction of Sicl4 and is free of analytically detectable amounts of boron and phosphorus. Its use permitted study of the effect of boron and phosphorus individually on the properties. This material contains, however, spec-troscopic traees of a number of metals, some of which affect electrical and rectification properties. Typical analyses for the two grades of silicon are given in Table 1. To make controlled additions of boron to the charge it was necessary to employ a low boron content master alloy since the quantity of boron to be added was usually only a few thousandths of one percent. The alloy containing nominally one percent boron was prepared by melting chemically pure boron with silicon A using the melting techniques for 320 g silicon

Jan 1, 1950

By J. Orkiszewski

A method i.s presented which can accurately predict, with a precision of about 10 percent, two-phase pressure drops in flowing and gas-lift production wells over a wide range of well conditions. The method is an extension of the work done by Grifith and Wallis and was found to be superior to five other published methods. The precision of the method was verified when its predicted values were compared against 148 measured pressure drops. The unique features of this method over most others are that liquid holdup is derived from observed physical phenomena, the pressure gradient is related to the geometrical distribution of the liquid and gas phase (flow regimes), and the method provides a good analogy of what happens inside the pipe. It takes less than a second to obtain a prediction on the IBM 7044 computer. INTRODUCTION The problem of accurately predicting pressure drops in flowing or gas-lift wells has given rise to many specialized solutions for limited conditions, but not to any generally accepted one for broad conditions. The reason for these many solutions is that the two-phase flow is complex and difficult to analyze even for the limited conditions studied. Under some conditions, the gas moves at a much higher velocity than the liquid. As a result, the down-hole flowing density of the gas-liquid mixture is greater than the corresponding density, corrected for down-hole temperature and pressure, that would be calculated from the produced gas-liquid ratio. Also, the liquid's velocity along the pipe wall can vary appreciably over a short distance and result in a variable friction loss. Under other conditions, the liquid is almost completely entrained in the gas and has very little effect on the wall friction loss. The difference in velocity and the geometry of the two phases strongly influence pressure drop. These factors provide the basis for categorizing two-phase flow. The generally accepted categories (flow regimes) of two-phase flow are bubble, slug, (slug-annular) transition and annular-mist They are ideally depicted in Fig. 1 and briefly described as follows. RUBBLE FLOW (FIG. 1A) The pipe is almost completely filled with the liquid and the free-gas phase is small. The gas is present as small bubbles, randomly distributed, whose diameters also vary randomly. The bubbles move at different velocities depending upon their respective diameters. The liquid moves up the pipe at a fairly uniform velocity and, except for its density, the gas phase has little effect on the pressure gradient. SLUG FLOW (FIG. 1R) In this regime, the gas phase is more pronounced. Although the liquid phase is still continuous, the gas bubbler coalesce and form stable bubbles of approximately the rame size and shape which are nearly the diameter of the pipe. They are separated by slugs of liquid. The bubble velocity is greater than that of the liquid and can be predicted in relation to the velocity of the liquid slug." There is a film of liquid around the gas bubble. The liquid velocity is not constant—whereas the liquid slug always moves upward (in the direction of bulk flow); the liquid in the film may move upward but possibly at a lower velocity, or it may move downward. These varying liquid velocities will result not only in varying wall friction losses, but also in a "liquid holdup" which will influence flowing density. At higher flow velocities, liquid can even be entrained in the gas bubbles. Both the gas and liquid phases have significant effects on the pressure gradient. TRANSITION FLOW (FIG. 1C) The change from a continuous liquid phase to a continuous gas phase occurs in this region. The liquid slug between the bubbles virtually disappears, and a significant amount of liquid becomes entrained in the gas phase. Although the effects of the liquid are significant, the gas phase is more predominant.

By H. Biloni, B. Chalmers

Metallographic studies of the segregation of solute in and near the surfaces of chill-cast alloys have shown that, under sufficiently severe conditions, nucleation on the cold mold surface is followed by the formation of a disc of solid of the same composition as the liquid. This disc is surrounded by two concentric rings in which the concentration first decreases and then increases. The discs and rings are surrounded by a region in which arms of low solute concentration separate regions of high concentration. These arms start at the periphery of the outer ring and are initially radial; they make a gradual curved transition from the radial direction to a crystallographic direction. A normal dendritic substrcture is established beyond the position at which the arms have become crystallographically aligned. When the cooling is less severe, the disc becomes less irregular in shape and has a concentration minimum at the center. The term "Predendritic" is proposed for these transitions. A qualitative explanation is outlined for the phenomena observed at the surface and in the equiaxed region of the ingots. It is well-established1-' that the normal mode of solidification for alloys under most conditions is dendritic. A characteristic of dendritic solidification is the growth of rods or plates in specific crys-tallographic directions. The remaining liquid solidifies later, usually with a higher concentration of solute than the original liquid. Very little is known, however, about the transition from a crystal nucleus until the dendritic pattern is established. It is clear that there must be a transition because the nucleus is much smaller than the characteristic spacings found in the "steady-state" structure. And it also seems clear that the initial conditions must be important since the undercooling existing during the nucleation of the grains has a remarkable influence upon the grain size4 and the segregation pattern. The present paper reports a study of the micro-segregation in regions where crystals have nucleated, and the conclusions that can be drawn regarding the early stages of growth, in connection with the details of the solidification process. Aluminum of 99.993 purity, A1-Cu alloys containing 0.5 to 20 wt pct Cu, and Pb-Sn alloys with 5 wt pct Sn were cast and examined metallographically at the as-cast surface, close to the surface, and in the interior after casting. The alloys were cast from controlled temperatures against surfaces at different temperatures ranging from room temperature to approximately the liquidus temperature of the alloy. The surfaces were, in various experiments, graphite, stainless steel, copper, and aluminum. The methods used for metallographic examination were as follows. 1) The Nomarski phase contrast interferometer6 was used to obtain information on the topography of the as-cast surface. 2) The distribution of any solute present in A1 99.993 pct and copper in the A1-Cu alloys was studied by the anodic film technique originated by Lacombe and Mouflard7 and previously applied in the study of the solidification of aluminum by Biloni.8 The electropolished surface of the specimen is anodized under controlled conditions so that the thickness of the oxide film and therefore the interference colors under white light are sensitive to the local composition. This technique was applied to the as-cast surfaces and to metallographically prepared surfaces from the interior. 3) If in Al-Cu alloys, the surface is anodically oxidized for several minutes, instead of seconds as in the color-film technique, the copper segregation is outlined by a different level with respect to the matrix. This method is less sensitive but complementary to that for the color films. 4) The distribution of copper revealed and mapped by the metallographic techniques was confirmed by electron probe. 5) In samples with high concentration of copper, evidence for the copper distribution could be seen optically after electropolishing.

Jan 1, 1965

By F. W. Jessen, H. Kazemi

A washing technique has been developed to form a spherical cavity in massive salt. The technique is, basically, a process of controlling the fluid motion in the cavity, the concentration distribution of brine, the rate of dissolution at the walls of the cavity, and the particular settings of the wash pipe, or pipes (inlet and outlet pipe) employed. For a particular arrangement of tubing and casing, fresh water enters from the casing annulus at a position above the bottom of the tubing, and the dissolved material leaves the cavity through the tubing near the bottom of the projected hole. The initial configuration of the cavity is cylindrical, representing the drilled hole situated at the axis of the projected cavity. As washing progresses the original hole is enlarged radially (any direction perpendicular to the axis of the hole). An inert fluid (generally a hydrocarbon, natural gas, or air) which serves as a blanket is injected through the annular space between the original hole and the casing. The controlled downward motion of this blanket maintains the upper edge of the brine region that is expanding on the surface of the projected sphere. A large spherical glass container and a radial model were used to investigate the flow patterns and the concentration distribution of brine occurring during the progress of solution of the salt. Aluminum flakes and dye were injected and flow patterns were observed and recorded on motion picture film. Based on the results of these model studies and other work carried out in this laboratory, the controlled washing technique was formulated. INTRODUCTION Underground cavities leached in salt beds and salt domes have been used extensively for storing hydrocarbons. Little attention has been focused on the shape of the cavities as long as large volumes were obtained and the cavities were structurally stable. The general process involves circulating fresh water and removing brine. Two distinctly different circulating systems can be employed, each of which yields characteristically different shapes. By far the most generally accepted method of washing storage caverns in salt domes is the direct circulation system. Water is introduced through the innermost pipe of a concentric pipe system, and brine is withdrawn through the annulus. In the reverse-circulation method, water enters through the annulus and is discharged through the smallesc tubing string. The size of the operation, and hence the size of tubing and casing used, depends on such factors as cavity volume, length of hole, availability of water, disposal of brine and power requirements. Underground caverns are from 100 to 1,000 ft or more in length, and when completed vary from 40 to 100 ft in diameter, providing storage of 100,000 bbl or more. The volumes of water required per volume of salt removed vary from 6 to 7 for the reverse circulation method to 10 to 11 for the direct circulation method. Up to 80,000 BWPD are often circulated. In the direct circulation method, because water enters near the bottom of the projected cavity, solution of salt is greatest at this point, particularly in view of the fact that when the washing is begun turbulent flow undoubtedly exists in the immediate vicinity of the discharge of the inlet water. As the diameter of the hole increases, flow rates also are increased, so perhaps the condition of fluid turbulence is maintained for some time. During the later stages of leaching, when capacity of equipment may be the limiting factor in maintaining flow rates, there may still exist some turbulence adjacent to the point of influx; however, for the most part laminar flow conditions prevail. There appears, therefore, rather good evidence to indicate that the state of fluid motion is responsible for the typically shaped cavity which results from such washing — namely a cylindrical form, generally somewhat larger at the bottom than at the top. Depending on the number of so-called stages which might be employed, the walls or surfaces of the salt diminish rather uniformly from bottom to

Jan 1, 1965

By P. Datta, L. L. Handy

For a wetting phase displacing a nonwetting phase from a porous medium the distribution of the residual fluid may depend on displacement conditions. Although this subject has been debated in the literature, only a few experiments have been cited to support the various conclusions. Experimental results presented in this paper show that fluid distributions are dependent on imbibition pmcedures. Results agree qualitatively with predictions from the pore doublet model. If the rate of water imbibition is restricted, the nonwetting phase is trapped preferentially in .the larger pores as expected. But if the rate of water imbibition is unrestricted, trapping occurs somewhat more in the smaller pores. These conclusions were deduced primarily from relative permeability measurements. INTRODUCTION Relative permeabilities are known to depend on the saturation history of the porous medium. For either continuously increasing or continuously decreasing wetting phase saturations, however, they have normally been assumed to be single-valued functions of saturation. Severat investigators have compared relative permeabilities measured by different methods and have reported acceptable agreement.1-5 Studies of simple models of the displacement process suggest that the distribution of residual fluids can be influenced by the displacement method. If this is so, relative permeabilities also should depend on the method of displacement. These predictions have not been supported by experimental data. The object of this paper is to investigate experimentally some of the predictions of these model studies. The particular question of importance is whether steady-state and unsteady-state methods should be expected to give the same values of relative permeability. Much of the reported data showing agreement between steady-state and unsteady-state methods have been obtained on un consolidated sand. Johnson, Bossler and Naumann, however, compared steady-state and unsteady-state results for Weiler sandstone and found no significant differences.5 Levine made a detailed study of pressure and saturation distributions for a laboratory waterflood in a consolidated system but made no attempt to calculate performance from steady-state relative permeability data Bail and Marsden in a somewhat similar set of experiments did attempt a comparison of observations with predictions but the results were inconclusive.7 Excluding gravity, two forces affect the distribution of wetting and nonwetting phases during an unsteady-state, immiscible displacement: viscous and capillary forces. If relative permeabilities are indeed the same when measured by steady- or unsteady-state methods, the fluid distributions at the same fluid saturations must be similar for both methods. Furthermore, the implication is that the relation between capillary and viscous forces is the same for both methods insofar as the effect on microscopic fluid distributions is concerned. Unsteady-state displacements are normally run at high rates to maximize the ratio of the viscous to capillary forces. The objective is to reduce the effect of capillary forces on macroscopic fluid distributions. For floods in which the phase which wets the porous medium displaces the nonwetting phase, capillary forces compete with viscous forces in determining fluid distributions.8 The residual nonwetting phase is discontinuous. Competitive aspects of viscous and capillary forces and the resulting effect on water-oil distribution in porous media have been illustrated in an elementary way using a pore doublet model. This model and predictions from it have been discussed at considerable length by Rose and Witherspoon, Rose and Cleary and by Moore and Slobod.9-13 Rose and Witherspoon show that so long as water invades the model at a rate equal to or greater than the free imbibition rate, the water will move through the larger capillary of a doublet first and trap oil in the smaller capillary. On the other hand, if a negative pressure is applied to the water

Jan 1, 1967

By Donald A. Brobst

The minerals barite (BaSO4-barium sulfate) and witherite (BaCO3-barium carbonate) are the chief commercial sources of the element barium and its compounds whose many uses are nearly hidden among the technical complexities of modern industrial processes and products. Barite, the major ore mineral, is extremely vital to the petroleum industry which in 1980 consumed at least 90% of the world's production of 8.8 million tons as a major ingredient of the heavy fluid, called mud, that is circulated in rotary drilling of oil and gas wells. The remaining 10% of the barite production was consumed chiefly in the manufacture of barium chemicals and glass and as a pigment, filler, and extender. Barite occurs throughout much of the world and is available from three major geologic types of deposits-vein and cavity filling, residual, and bedded-in sufficient quantity at competitive prices to meet current demands. The world's demand for barite is expected to increase, and geologic circumstances are favorable for the discovery of new deposits of commercial value. Witherite is much less common and abundant than barite, although it is more desirable in many ways as a raw material for the production of barium chemicals. The United States has not produced witherite since about 1950; England is currently the chief producer of witherite. End Uses Most of the world's barite production since 1926 has been used as a weighting agent for the muds circulated in rotary drilling of oil and gas wells. The muds fundamentally are mixtures of water, clay, and barite in different proportions that vary according to local reservoir conditions. Muds with a specific gravity as great as 2.5 are used. The mud is pumped down the hollow drill stem, passes through the bit at the bottom of the hole, and rises to the surface in the space between the drill stem and the wall of the hole. In the course of this circulation the drill bit is cooled, cuttings are removed, the drill stem is lubricated, the walls of the hole are sealed, and the hydrostatic head of the column of weighted fluid helps to confine high oil and gas pressures. The latter feature aids in prevention of gushers, thus reducing both environmental pollution and waste of oil and gas resources as well as conserving the natural reservoir pressure for greater production and rate of recovery of the products contained in the rocks. Barite is particularly well suited for drilling mud because it is clean, easy to handle, soft (nonabrasive), heavy, virtually inert chemically, and relatively inexpensive compared to many other available heavy materials. Barite is used in the manufacture of glass in continuous tanks. The addition of barite homogenizes the melt and gives greater brilliance and clarity to the finished glass. Ground barite, both unbleached and bleached by sulfuric acid, is a common industrial filler, extender, and weighting agent. The rubber industry is a major consumer of barite as a filler. Barite also is added to bristolboard, heavy printing paper, playing cards, rope finishes, brake linings, clutch facings, plastics, and linoleum, to name but a few uses. Bleached barite has long been an extender in white lead paint because of its weight. The low index of refraction of barite, however, makes its ability to cover marks poorer than some other substances, but its low capacity for absorption of oil is a good feature. Off-color or unbleached barite can be used as a filler in colored paints. In the construction industry, some lump barite is used in concrete aggregate to weigh down pipelines buried in marshy areas and to shield nuclear reactors. Because barite absorbs gamma radiation well, its use reduces the amount of expensive lead shielding otherwise

Jan 1, 1983

By G. Sandoz, D. A. Meyn

The fracture surfaces of specimens of titanium a1loys which exhibited susceptibility to subcritical crack growth in a wide variety of environments, including aqueous solutions, alcohols, hydrocarbon gases, carbm tetrachloride, and dry air, were examined. The dominant and characteristic fracture mode was cleavage. which became mixed with an increasing porportion of ductile fracture modes (dimples. and so forth) as the applied KI was increased. The cleavage plane in all the alloys was oriented at 15 deg from (0001) of the a phase. No indices could be assigned to the 15 deg cleavage Plane because of uncertainty as to which zone it belonged to. SLOW crack growth and delayed fracture take place in certain titanium alloys at stress intensity (KI) lower than the critical values for complete fracture (KIc) found in short-time tests for plane strain fracture toughness.' Such crack growth may take place in very dry air on sustained loading if the alloy contains a sufficient amount of hydrogen and if the stress intensity factor is higher than a threshold value designated KIH.2 Similarly, a number of environments such as water, alcohols, and hydrocarbons may cause slow crack growth provided the stress intensity factor exceeds the threshold value KISCC. In this paper it is intended to examine some fracto-graphic and crystallographic details of this subcritical crack propagation process as related to alloy composition, environment, and stress intensity. It has been previously reported that subcritical crack growth takes place in titanium alloys in aqueous environments by cleavage.2"6 The observations to be described here indicate that in the alloys investigated the mode of subcritical crack propagation is more dependent on the stress intensity during cracking than on the environment. The fracture mode is illustrated with photographs of the fracture surfaces at a wide range of magnifications. The crystallographic character of the cleavage mode is also presented and the results discussed. EXPERIMENTAL APPROACH The studies were made on four titanium alloys— Ti-8A1-1Mo-lV, Ti-5A1-2.5Sn, Ti-7Al-1Mo-lV, and Ti-0.350. Of these four alloys the first two were mill annealed and fine grained. The latter two were noncommercial and coarse grained, having been prepared for other experimental purposes. The crystallographic studies required a large-grain size. Consequently the fine-grained alloys were an- nealed in a vacuum at 1950°F and furnace cooled to produce samples for these studies. This treatment did not alter the microstructure significantly except for grain size. All the aluminum-containing alloys showed primary a with a network of ß-phase particles within the a grains; Figs. 1 and 2 show this microstructure for the specific case of fine-grained Ti-8A1-1Mo-1V. The cracks in the figures will be discussed later. Fractographic studies were made on all the alloys, but the electron fractographs displayed herein are from the Ti-8A1-1Mo-1V alloy in the fine-grained condition. Emphasis was placed on this commercial alloy because considerable information was available on the stress-corrosion cracking characteristics in a variety of environments. Sufficient supplementary studies of the alloy after the anneal to produce large grains were conducted to establish that the stress-corrosion cracking susceptibility and fractographic characteristics were not changed significantly. Thus it could be inferred that the crystallographic data obtained exclusively on large-grained specimens would apply also to fine-grained specimens. Fracture surfaces for study were obtained using precracked cantilever beam specimens loaded as described by Brown. The precracked section of a specimen was immersed in the environment of interest. Values of initial stress intensity factor Kli were chosen slightly higher than KISCC in order to obtain a fairly extensive region of slow crack growth. Fracture surfaces representative of a variety of stress intensity levels could then be examined because in the cantilever beam test KI increases as the crack propagates.

Jan 1, 1970

By Donald A. Brobst

The minerals barite (BaSO4 barium sulfate) and witherite (BaCO3 barium carbonate) are the chief commercial sources of the element barium and its compounds whose many uses are nearly hidden among the technical complexities of modern industrial processes and products. Barite, the major ore mineral, is extremely vital to the petroleum industry which in 1973 consumed at least 80% of the world's produc¬tion of 4.4 million tons as a major ingredient of the heavy fluid, called mud, that is circulated in rotary drilling of oil and gas wells. The re¬maining 20% of the barite production was consumed chiefly in the manufacture of barium chemicals and glass and as a pigment, filler, and extender. Barite occurs throughout much of the world and is available from three major geologic types of deposits-vein and cavity fill¬ing, residual, and bedded-in sufficient quantity at competitive prices to meet current demands. The world's demand for barite is expected to increase, and geologic circumstances are favor¬able for the discovery of new deposits of com¬mercial value. Witherite is much less common and abun¬dant than barite, although it is more desirable in many ways as a raw material for the pro¬duction of barium chemicals. The United States has not produced witherite since about 1950; England is currently the chief producer of witherite. End Uses Most of the world's barite production since 1926 has been used as a weighting agent for the muds circulated in rotary drilling of oil and gas wells. The muds fundamentally are mixtures of water, clay, and barite in different proportions that vary according to local reservoir condi¬tions. Muds with a specific gravity as great as 2.5 are used. The mud is pumped down the hollow drill stem, passes through the bit at the bottom of the hole, and rises to the surface in the space between the drill stem and the wall of the hole. In the course of this circulation the drill bit is cooled, cuttings are removed, the drill stem is lubricated, the walls of the hole are sealed, and the hydrostatic head of the column of weighted fluid helps to confine high oil and gas pressures. The latter feature aids in prevention of gushers, thus reducing both environmental pollution and waste of oil and gas resources as well as conserving the natural reservoir pressure for greater production and rate of recovery of the products contained in the rocks. Barite is particularly well suited for drilling mud because it is clean, easy to handle, soft (nonabrasive), heavy, virtually inert chemi¬cally, and relatively inexpensive compared to many other available heavy materials. Barite is used in the manufacture of glass in continuous tanks. The addition of barite ho¬mogenizes the melt and gives greater brilliance and clarity to the finished glass. Ground barite, both unbleached and bleached by sulfuric acid, is a common industrial filler, extender, and weighting agent. The rubber in¬dustry is a major consumer of barite as a filler. Barite also is added to bristolboard, heavy printing paper, playing cards, rope finishes, brake linings, clutch facings, plastics, and li¬noleum, to name but a few uses. Bleached barite has long been an extender in white lead paint because of its weight. The low index of refraction of barite, however, makes its ability to cover marks poorer than some other sub¬stances, but its low capacity for absorption of oil is a good feature. Off-color or unbleached barite can be used as a filler in colored paints. In the construction industry, some lump barite is used in concrete aggregate to weight down pipelines buried in marshy areas and to shield nuclear reactors. Because barite absorbs gamma radiation well, its use reduces the amount of expensive lead shielding otherwise

Jan 1, 1975

By K. H. Ribe

Water often enters an oil reservoir during completion or workover operations on a well and forms a partial "water block" to oil production. A mathematical study of radial two-phase flow, neglecting capillary effects, has been employed to study the formation of such a water block and subsequent re-moval by production from the well. The effects of reduced oil permeability about the well on the well productivity were studied. The fluid saturation distributions about the well during formation and removal of the water block have also been computed. Several relative permeability relations and viscosity ratios were employed. If water has invaded the formation, its influence through relative permeability effects alone can cause the following. 1. Oil productivity will be depressed for extended periods after production is resumed and will build up only gradually as the water is removed. 2. Oil injected for treatment of water blocking will delay rather than promote restoration of full well productivity by enlarging the region invaded by water. Thus, unless the specific action of chemicals contained in the oil is needed, oil injection appears undesirable. INTRODUCTION During oilwell workover operations, water may enter the oil-bearing formation from the wellbore. When production is resumed, oil must flow through the region invaded by this water. The presence of this region can cause both well productivity and oil production rate to be low and oil to be produced with high water-oil ratio for some time after production is initiated. This situation is sometimes described as a water block. The introduction of water into the formation may result in other actions which also lead to reduction in well productivity and which are also usually included in the connotation of the broad term, water block. Often considered, for example, are the possibilities of clay swelling by contact with fresh water and the formation of emulsions with the formation oil. If it is suspected that such specific actions have taken place, remedial treatments are undertaken which usually involve the injection of chemicals in oil. Since the introduction of water, even in the absence of specific interactions with the formation or oil, will cause a temporary water block (which might be misinterpreted as evidence of a more severe situation), it is of importance to evaluate the magnitude and duration of this blocking which results purely from the reduction in relative permeability to oil in the vicinity of the wellbore. It is also of interest to evaluate the effect of oil injection on the productivity of a well blocked by water in this manner. Inasmuch as this unfavorable condition may persist for some time, it may lead to premature condemnation of a workover or premature abandonment of a potentially productive pay zone. A quantitative evaluation of the influence of water entry on the oil productivity through changes in relative permeability was made by solving a radial form of the Buckley-Leverett equation. The distribution of water saturation around the wellbore during the entry of water was calculated and was followed by a similar calculation of the saturation distribution during the period of resuming production. At any stage in the removal of the invading water, knowledge of the distribution of its saturation permitted calculating the attendant loss in oil productivity. The influence of the shape of the relative permeability relationships was also evaluated by carrying out the calculations for two hypothetical cases. Further, the effect of the oil-water viscosity ratio was examined by repeating the calculations, for several ratios of unity and greater, with the same relative permeabilities. Fi-nally, results are presented to show how the length of time a well must be swabbed to resume production depends on the length of time it has been subjected to invasion by water. STATEMENT OF THEORY Differential Equations Water is assumed to enter a producing formation which is initially at the connate-water saturation and contains no gas. The water and oil are treated as in-

By C. B. Manula, H. Gezik

In planning modem-day mining operations, management needs to pass from the area of subjective decision-making to an area of objective decision-makirlg. Planning procedures currently being practiced by many mine managers depend primarily on common sense and trial-and-error methods which oftentimes never reach optimum solutions. For this reason, an attempt is made to provide an understanding of linear programming to solve production planning problems. More specifically, a linear programming model is proposed to solve the production scheduling problem as it was found to exist in the crushed stone industry. Optimal production scheduling is of prime concern here, since the industry has to cope with seasonal demands for its products. The solution, as determined through the application of the model, provides management not only with an optimal production schedule but also with the following information regarding the improvement of existing facilities: (I) the location of bottlenecks which limit production and profit; (2) the directions in which the company i fixed equipment should be expanded or modified; and (3) the evaluation of unused resources. Future trends in stone production will quite likely parallel those in the construction industry. Continued record expenditures for structures such as dams, highways, and buildings are the reasons for this accelerated growth. Production forecasts for the next decade are predicted to be around one billion tons with an approximate 7% annual increase. The problem of increasing production from 570 million tons in 1966 to 1 billion tons by 1976 offers no small challenge. This situation is further compounded by other problems and economic factors. Population expansion, government regulation, labor agreement, and changes in technology wd cause managers and operators to be confronted by many more complexities than characterize present-day operations. Some of these. problems and economic factors are: (1) In some metropolitan and other urban areas, easily accessible sources of stone are almost depleted. Transportation costs will be a prime concern of the industry. (2) Population growth has forced residential expansion to a point where many quarry operators must conform to urban and suburban regulations. (3) Government regulation in regard to air and water pollution, noise, and reclamation as well as zoning developments necessitate that crushed stone producers keep abreast of changing conditions and requirements. (4) The crushed stone industry operates under pressure of a highly competitive situation. Every company's objective is to retain a favorable marketing position by exercising great attention to production and marketing costs. This is rather difficult, since the crushed stone industry is one segment of modern mining which has to cope with a seasonal demand for its products. Among the problems posed here, the one of exercising great attention to production and marketing costs is of current interest to management. Under seasonal demand conditions, this problem creates further difficulties for managers who need to plan and forecast. STATEMENT OF THE PROBLEM Various methods of Operations Research have been developed to handle a problem whose genesis is explained above. It is very seldom, however, that these quantitative means are used in practice by the mining industry. This paper describes the use of Operations Research methods to solve a production planning problem when demand is seasonal. The seasonal production planning problem is an important one because of the substantial costs associated with changing production, employment, and inventory levels to adjust to seasonal sales patterns. Most companies try to operate within reasonable levels of employment changes and hold inventories within storage capacities. In general, little or no attempts were made to develop minimum cosl schedules or to provide objective means to compare various production policies. When demand varies within a fixed time period, 11:-oduction requirements may be met by following one

Jan 1, 1970

By T. D. Mueller

Many investigators have used the response of the "dimensionless aquifer" to a unit pressure drop or a unit fluid-withdrawal volume to calculate the performance of an aquifer in supplying water influx to an oil reservoir. In the past, these response functions have been calculated with the aid of the Laplace transform. With the advent of ultra highspeed digital computers, it becomes practical to solve for the response functions with finite-difference techniques. The computer method also permits extension of the dimensionless-aquifer concept to include the nonhomogeneous aquifer wherein the permeability and other properties vary as a function of the space co-ordinates. This paper gives results of calculating the response functions for a series of nonhomogeneous aquifers. Response functions are presented for both linear and radial aquifers whose thickness, permeability-viscosity ratio and porosity-compressibility vary. These functions are new and should prove useful to the petroleum engineer in analyzing the behavior of nonbomogeneous aquifers. Results are presented in the form of charts that can be easily used by the field engineer. INTRODUCTION Aquifers which surround many oil and gas reservoirs have the ability to supply water influx to such reservoirs as oil and gas are withdrawn. This water influx, called natural water drive, provides one of the most effective driving mechanisms for the production of oil and gas. In producing a reservoir, therefore, it behooves one to make the maximum use of natural water drive. To achieve the maximum use, the reservoir engineer must be able to predict the performance of an aquifer under a variety of production schemes that may be proposed for the reservoir. Unfortunately, the physical properties which dictate aquifer behavior often are known only within limits. Seldom do wells penetrate the porous strata of the aquifer. Even when they do, quantitative information regarding porosity, permeability and water compressibility is seldom available. It is known, however, that the water efflux from most aquifer systems is governed by a single, relatively simple, linear, partial differential equation. Also, the general physical location of the aquifer boundaries often are known. A technique originall proposed by Hurst1 and van Everdingen and Hurst2 has been found useful in analyzing reservoirs in this situation. The idea was later expanded by van Everdingen, Timmerman and McMahan3 to include the mathematical technique of least-squares fitting. This latter method will be referred to as the VTM method. The basic assumptions of the VTM method include the following. 1. The location of the physical boundaries of the aquifer are known. 2. The flow conditions at these physical bound-daries are also known. 3. The aquifer is homogeneous; e.g., thickness, permeability, water compressibility and porosity are constant throughout. In the VTM method, a material balance is made on the fluids entering and leaving the reservoir. In the balance equations, the water-influx term is represented as the product of the water influx from an arbitrarily-selected, dimensionless aquifer system times an unknown conversion number. This balance can be formed as many times as there are data points in the history of the reservoir. Each time, the conversion number can be evaluated. If the reservoir engineer has picked the correct dimensionless aquifer to represent the real aquifer, the conversion number remains constant for all balances that have been made over the history period. If such a situation occurs, the reservoir engineer can then predict the performance of the reservoir for any type of production scheme by using the function associated with that particular dimensionless system and the derived conversion number. These functions will be referred to as the response function of the aquifer. If the aquifer is nonhomogeneous (e.g., if the porosity, permeability, thickness, porosity, or

By William Hurst, A. F. van Everdingen

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By A. F. van Everdingen, William Hurst

For several years the authors have felt the need for a source from which reservoir engineers could obtain fundamental theory and data on the flow of fluids through permeable media in the unsteady state. The data on the unsteady state flow are composed of solutions of the equation Two sets of solutions of this equation are developed, namely. for "the constant terminal pressure case" and "the constant terminal rate case." In the constant terminal pressure case the pressure at the terminal boundary is lowered by unity at zero time, kept constant thereafter, and the cumulative amount of fluid flowing across the boundary is computed, as a function of the time. In the constant terminal rate case a unit rate of production is made to flow across the terminal boundary (from time zero onward) and the ensuing pressure drop is computed as a function of the time. Considerable effort has been made to compile complete tables from which curves can be constructed for the constant terminal pressure and constant terminal rate cases, both for finite and infinite reservoirs. These curves can be employed to reproduce the effect of any pressure or rate history encountered in practice. Most of the information is obtained by the help of the Laplace transformations, which proved to be extremely helpful for analyzing the problems encountered in fluid flow. The application of this method simplifies the more tedious mathematical analyses employed in the past. With the help of Laplace transformations some original developments were obtained (and presented) which could not have been easily foreseen by the earlier methods. INTRODUCTION This paper represents a compilation of the work done over the past few years on the flow of fluid in porous media. It concerns itself primarily with the transient conditions prevailing in oil reservoirs during the time they are produced. The study is limited to conditions where the flow of fluid obeys the diffusivity equation. Multiple-phase fluid flow has not been considered. A previous publication by Hurst' shows that when the pressure history of a reservoir is known, this information can be used to calculate the water influx, an essential term in the material balance equation. An example is offered in the literature by Old² in the study of the Jones Sand, Schuler Field, Arkansas. The present paper contains extensive tabulated data (from which work curves can be constructed), which data are derived by a more rigorous treatment of the subject matter than available in an earlier publication.' The application of this information will enable those concerned with the analysis of the behavior of a reservoir to obtain quantitatively correct expressions for the amount of water that has flowed into the reservoirs, thereby satisfying all the terms that appear in the material balance equation. This work is likewise applicable to the flow of fluid to a well whenever the flow conditions are such that the diffusivity equation is obeyed. DIFFUSITY EQUATION The most commonly encountered flow system is radial flow toward the well bore or field. The volume of fluid which flows per unit of time through each unit area of sand is expressed by Darcy's equation as where K is the permeability, the viscosity and aP/ar the pressure gradient at the radial distance r. A material balance on a concentric element AB, expresses the net fluid traversing the surfaces A and B, which must equal the fluid lost from within the element. Thus, if the density of the fluid is expressed by p, then the weight of fluid per unit time and per unit sand thickness, flowing past Surface A, the surface nearest the well bore, is given as The weight of fluid flowing past Surface B, an infinitesimal distance 6r, removed from Surface A, is expressed as

Jan 1, 1949

By K. H. Coats

This paper describes an approximate technique for handling the problem of percolation of evolved gas upwards through the oil column in computer simulation of natural depletion. This technique has been incorporated into a general model for simulating three-phase flow in one, two or three dimensions. The mathematical model for perfonning these calculations is described in detail in a separate article. 1 While vertical gas percolation occurs during natural depletion in a reservoir of any configuration, it is especially pronounced in the pinnacle reel or bioherm. The reel may have an areal extent less than one section with a thickness of up to 800 It. The upward percolation of evolved gas often limits the time step in calculations as little as one day. Using this small time step, computer expense for a single 30- to 40-year simulation on even a one-dimensional basis has exceeded 92,000. The method described here for handling the percolation allows time steps of 60 days or more (depending on reservoir size and production rate), resulting in a considerable reduction in computing expense. The method also allows calculations in which secondary gas caps build up in tight zones of the oil column below the main gas cap. The validity of the method is indicated in connection with an example pinnacle reef field. Calculated results using a 2-day time step (where the method in question is not needed and not invoked) are compared with results using the method and a 60-day time step. The comparison shows good agreement. Results from one- and two-dimensional simulations of the reef are presented along with the corresponding computing times on the CDC 6600 computer. A three-dimensional simulation was also performed and the required computing time is given. INTRODUCTION A general analysis for simulating three-dimensional, three-phase flow in reservoirs has been developed. 1 In applying this model to reservoirs undergoing natural depletion, a time-step restriction was encountered due to the flow of evolved gas upwards through the oil column and toward the gas cap. A remedy for this problem has been incorporated in the analysis and is described here. The time-step restriction is encountered only in calculations which include flow in the vertical or near vertical direction. The restriction occurs to some extent in reservoirs of any geometrical configuration, but it is especially pronounced in the pinnacle reef or bioherm where the ratio of thickness to areal extent is unusually large. This paper describes the problem, a method of handling it, incorporation of the method in the three-dimensional, three-phase model, and, finally, a test of the method's validity in an application to an example pinnacle reef reservoir. THE PROBLEM During the early stages of reservoir depletion, pressure falls below bubble point in progressively lower regions of the oil column. Continuing production results in evolution of dissoived gas throughout the oil column. This gas then percolates upwards toward the top of the reservoir. If the gas encounters a sufficiently low permeability zone in its travel upward, then it can accumulate and form a secondary gas cap. If no zones are encountered which are tight enough to hold gas against the gravity forces, then the gas travels on until it reaches the top of the sand or the main gas cap. Numerical simulation of multiphase flow in reservoirs is generally performed with time steps such that the flow of a fluid into or out of a block in one time step is a fraction (considerably less than one) of the total amount of that fluid present in the block. Use of a time step where the time step's

Jan 1, 1969

By J. L. Huitt, S. R. Darin

The use of a partial monolayer of propping agent to obtain a high flow capacity for a hydraulically induced fracture is discussed. From the results of laboratory work it was shown that a modified form of the Kozeny-Carman relation could be used to describe the flow in the partial monolayer propped fracture. With equations presented in the paper, the density pattern of the propping agent (number of particles per unit of fracture surface) that results in the maximum flow capacity for the fracture can be determined. The maxitnum flow capacity obtained with a partial monolayer is often an order of magnitude greater than the flow capacity obtained in greater width fractures containing multilayers of the propping agent. INTRODUCTION One of the predominant factors controlling the success of a hydraulic fracturing operation is the propping of the fracture. The trend in hydraulic fracturing recently has been an increase in the ratio of propping agent-to-fluid. The primary purpose of this increased ratio is to sustain a propped fracture of greater width. It was shown in a recent study' that for some formations low concentrations of propping agent would result in "closing" or "healing" of the fracture. In some cases, in an effort to insure sufficient propping agent concentration, the fracturing operation is designed to obtain a pack of the propping agent in the fracture. The placing of a pack of propping agent in the fracture provides a fracture of maximum width; however, the width alone does not control the flow capacity of the fracture. The flow capacity is dependent on the permeability of the fracture, as well as the fracture width. Thus, increasing the permeability to obtain larger flow capacities is as important as obtaining a greater width fracture. By increasing the concentration and particle diameter of the propping agent, fractures of sufficient flow capacities are obtainable for effective well stimulation in a number of formations. However, in some of the recovery techniques, the fracture capacity obtained by commonly used practices in fracturing is inadequate. One such recovery technique is gravity drainage via a horizontal fracture placed near the base of the productive zone. For this recovery technique, the fracture capacity desired is usually one to two orders of magnitude greater than the fracture capacity normally obtained in fractures propped with the size of sand commonly used (20-40 mesh). By the use of particles of larger diameter, the permeability of a pack2 is increased; however, even a multilayer pack results in a fracture flow capacity lower than that desired for some formations. For example, a fracture of 1/4-in. width, packed with 8-12 mesh sand of about 2,000 darcies permeability, would result in a fracture flow capacity (permeability X fracture width) of only 40,000 md-ft. If the fracture flow capacity needed is an order of magnitude greater than that cited in the example, an approach different from the formation of a pack in the fracture is needed. One approach is a partial monolayer of large-diameter propping agent. In working with this type of fracture system, it was found that the flow problems were unique to the system; therefore, the study discussed in this paper was made. THEORY AND RELATED STUDIES Previous studies, although related to that described in this paper, have not pertained directly to the flow in a fracture propped with a partial monolayer of propping agent. However, studies such as flow between parallel plates,' flow in packed columns,' fracture flow capacity' and flow in simulated fractures reduce the scope of this flow problem considerably. A fracture propped with a partial monolayer of dropping agent may be considered to be between two extremes — between an open fracture and a bed of particles packed In a conduit of fracture geometry. Thus, it seems appropriate to review the extreme cases.

By I. R. Kramer

The stress associated with the high-dislocation layer at the surface of deformed aluminum crystals was measured by progressively polishing the specimen and determining the change in the initial flow stress when the specimen was reloaded. This surface stress increased linearly with strain in the Stage II and Stage III region. The depth to which the high-dislocation layer extends was found to be 0.0025 in. The influence of the surface layer on Stage I is discussed and it is shown that Stage I ends when, on the secondary slip system, the difference between the applied stress and the opposing stress due to the surface layer is equal to the critical resolved shear stress. In previous papers'-3 it was shown that the surface exerts a pronounced effect on the plastic-flow characteristics of metals. For fcc metals it was possible to increase the extent and decrease the slope of Stages I and II and decrease the work-hardening coefficient of Stage II by removing the surface of the specimen continuously during the deformation process.2,3 More recently,L it was shown that the activation energy for plastic deformation and the activated volume were also influenced by the surface. These changes were related to the existence of a high concentration of dislocations in the region near the surface of the deformed specimen and it was concluded' that the effective stress, T, acting on a dislocation was T = Ta-Ti-Ts [1] where 7, is the applied shear stress, Tj is the internal stress, and TS is the stress associated with the surface region containing the high concentration of dislocations. From Eq. [I] it is seen that the effective stress, 7, acting on a dislocation is a function of 7,. Therefore 7, becomes important in any description of the work-hardened state. In this paper the value of 7, as a function of strain and its variation with depth are reported for aluminum single crystals. In particular, it will be shown that the end of Stage I occurs when the net stress, 7, on the secondary slip system is equal to the critical resolved shear stress. To measure Ts, single crystals of aluminum were deformed in tension to various strains, whereupon the load was removed and a known amount of metal was removed by electrolytic polishing. The differ- ence between the final flow stress before the specimen was unloaded and the initial flow stress upon reloading is designated as ?Tp and is the decrease in the stress due to the removal of the surface layer in whole or in part. The value for Tp is taken equal to -Ts. EXPERIMENT PROCEDURE The details of the tensile apparatus, method of electrolytic polishing, and crystal preparation were the same as those used previously,1"3 The single crystals were 1/8 by 1/8 in. in cross section with a 3-in. gage length. The initial purity of the aluminum was 99.997 pet. The specimens were held in vacuo at temperatures 50°C below their melting points for 16 hr and furnace-cooled. Just prior to the application of the tensile stress the specimens were electrolytically polished to remove 0.004 in. from the thickness. This latter operation was conducted after the specimen was placed in the tensile apparatus. The tensile tests were conducted at a strain rate, i, of 10"5 sec-' and a temperature of 3°C. The aluminum crystals used for the determination of 7, are designated A1-3-12 and A1-116. The orientation of the various crystals used throughout the investigation is shown in Fig. 1. The ?Tp values for specimens A1-3-12 as a function of Ax were determined mainly at four strain levels, i.e., ? = 0.0217, 0.0381, 0.0593, and 0.075, where Ax is one half of the total reduction of the thickness after polishing and ? is the re-

Jan 1, 1965

By A. J. Ardell, O. D. Sherby

The elevated-temperature steady-state creep behavior of polycrystalline a Zr was studied in vacuo under constant tensile stress. The experiments were conducted from 660° to 845°C over the stress range 870 to 5400 psi. The steady-state creep rate, is, of a Zr was found to be related to the applied stress, a, and temperature, T, by the equation where C is a constant. Two interpretations of the data are considered, each requiring different values of the parameters ß, Q, and . It is shown that this formula can be derived using Gilman's phenomeno logical equation relating dislocation velocity to the shear stress for the case where dislocation mobility is limited by the nucleation of kink motion. From the fundamental dependence of strain rate on dislocation velocity and the density of mobile dislocations, the factor exp(ßo) can be related to the stress dependence of the dislocation density, and Q to the energy of kink formation. The two interpretations considered are as follows: Case I—the parameters are adjusted so that the equation above describes the data over the entire range of stress and temperature, implying that creep is controlled exclusively by dislocation glide; Case IIthe parameters are adjusted to allow for the possibility that creep is diffusion-controlled at high stresses but glide-controlled at low stresses. For both cases the values of the parameters are within the spread permitted by the data. Factors favorable to the observation of the unusual stress -dependence of i, for a Zr are discussed and it is shown that the anomalously rapid rates of self-diffusion and the low modulus on the principal glide system are primarily responsible for the anomalous creep behavior. Metal-lograflhic observations are presented and discussed with reference to the overall creep behavior. PREVIOUS investigations of the creep behavior of a Zr have been conducted from room temperature to about 500C The results of these investigations are primarily of engineering significance because of the variety of conditions under which they were carried out, and the low testing temperatures employed. Little information is available concerning the fundamental nature of the elevated-temperature creep behavior of a Zr. This investigation of the creep behavior of a Zr above half the melting temperature (the melting temperature, T,, of a Zr is estimated to be about 1860"K was initiated because there are two significant anomalies of the physical properties of a Zr: a) the rate of self-diffusion in polycrystalline a Zr is anomalously rapid: b) both Young's modulus for polycrystalline a Zr and the shear modylus as-sociated with the prismatic slip system { 1010)( 1120) (the principal slip system in a Zr 7) have anomalously large temperature variations The importance of these quantities to the understanding of creep behavior has been demonstrated, and it was hoped that this investigation would provide a critical test of the predictions based upon the analysis of creep data for other pure metals. However, as the investigation progressed it became increasingly clear that the unusual creep behavior of a Zr could not be analyzed in a conventional manner. 1) EXPERIMENTAL EQUIPMENT AND PROCEDURES Two grades of zirconium, designated A and B, were used in this investigation. The analyses, as supplied by the vendors, are presented in Table I. Material A, with a nominal purity of 99.84 pct, was supplied by the Wah Chang Corp. in the form of cold-rolled sheet from which specimens were machined without further fabrication. Material B, with a nominal purity of 99.95 pct, was purchased from the Nuclear Materials and Equipment Corp. in the form of iodide crystal bar.

Jan 1, 1968