Search Documents

By A. B. Cook, C. J. Walker, G. B. Spencer

Although water will displace oil from a petroleum reservoir to a greater extent than gas will, there are some reservoirs in which gas rather than water should be used for pressure maintenance. This is indicated by the high-percentage oil recovery in the case history reports on the Pickton field1 in East Texas and the Raleigh field2 in Mississippi. In gas cycling, total oil recovery includes vaporized products from the immobile oil in addition to oil produced by displacement. If a large part of the immobile oil is vaporized, total oil recovery may be higher than that obtainable from pressure maintenance by water injection. However, the amount of immobile oil vaporized may range from almost 0 to 100 percent, and no simple and reliable method has been presented in the literature for calculating oil vaporization. Development of such a method was the purpose of this study. Several papers concerning the calculation of oil vaporization have been published,3 6 and each is based on the concept dealing with K values — equilibrium constants for the various components in reservoir oils. (K for a component in a system of vapor and liquid in equilibrium is the mole fraction in the vapor phase divided by the mole fraction in the liquid phase.) This appears to be the most logical approach. Yet, the problem is too great for a perfect solution, even with modem laboratories to analyze oils and high-speed digital computers to perform calculations; the number of components in reservoir oils is too large. Because of their complexity no two reservoir oils are exactly alike. Also, K values change not only with variations in pressure and temperature but also with composition of the reservoir oil. In addition, during gas cycling, the lighter hydrocarbons tend to vaporize first. Thus, the reservoir oil becomes more dense and less volatile as gas cycling continues. Furthermore, the greatest amount of oil vaporization occurs near the injection well. Therefore, a simplified method is required for calculating oil vaporization because a rigorous method is not practical. The present methods of calculating vaporization generally require knowledge of the reservoir oil composition in which the mole fractions of the lighter components through hexanes are given; the remaining oil is described as C, + (heptanes and others of greater molecular weights). Assigning one K value for the C,+ system provides the simplest solution, but this incorrectly assumes that K for C, + does not change as gas cycling continues. This assumption can cause large errors because the K value may become less than a thousandth of the original after a large amount of oil vaporization. Another way to solve the problem is to determine changing K values for C, + according to how much of the C, + has been vaporized. To determine these values a sample of the reservoir oil is injected into a pressure-volume cell, and dry gas is batched in and out of the cell. Appropriate measurements and tests are made to determine K

Jan 1, 1970

By Michael Tenenbaum

Selected characteristics of nonmetallic phases were determined in low-carbon sheet steels which had been produced by various refining; deoxidizing, and casting procedures. The results showed little influence of the refining process on the quantity, distribution, size, and character of the nonmetallic phases. The quantity of inclusions was governed by the deoxidation procedures. While the mode of casting influenced the internal structure and the pattern of inclusion distribution, it had little effect on the quantity of nonmetallics. The vast majority of inclusions were under 30 µ in diam, and there was little change in size distribution during hot and cold rolling. The composition and structure of the inclusions were expained an the basis of known equilibrium phase diagrams for the appropriate systems. Additional work is needed to explain some of the observations of this paper. HENRY Marion Howe's contributions to the metallurgical profession include his many writings, the wisdom that he imparted to his students, and the inspiration that he gave to the many scholars, teachers, and industrialists who knew him personally. It was Professor Howe's intimate colleague, Albert Sauveur,1 who, in the first Howe Memorial Lecture, could have been speaking for the metallurgist of our generation as well as his own when he wrote, "...how many of us on the point of going astray were called back to the right path by Howe's guiding torch?" More than a quarter of a century later, Francis Foley,2 one of the last in the Howe Lecture series who worked directly with the great man, delivered the ultimate tribute, "...truth, in its time is so hard to identify." While Professor Howe strove to be right, "he was content if he provoked the reader to take the next step along the path to greater knowledge." Like so many of my predecessors, I sought guidance for this lecture by referring to Professor Howe's two classic texts, Metallurgy of steel,3 and Metallography of Steel and Cast Iron: the former book interpreting the facts as he saw them near the start of his career, and the latter—published 26 years later—pointing to many fertile fields for metallurgical investigation. Those of us who harvest the fruits of Henry Marion Howe's contributions to the metallurgical profession now look with wonder at the tremendous legacy he has provided. It is rewarding to study his writings and to recognize the anticipation of such aspects of modern ferrous metallurgy as: 1) Direct reduction 2) Steelmaking reactions 3) The approaches and benefits of degassing 4) Solidification mechanisms 5) Habits and behavior of nonmetallic inclusions 6) Recrystallization phenomena 7) Slip and slip propagation 8) Dislocations 9) Discontinuous deformation 10) Microscopic and macroscopic segregation

Jan 1, 1970

By Rita D. Sylvester

"IntroductionThe Bureau of Mines was established in the public interest to conduct inquiries and scientific and technologic investigations concerning mining and the preparation, treatment, and utilization of mineral substances; to promote health and safety in the mineral industries; to conserve material resources and prevent their waste; to further economic development; to increase efficiency in the mining, metallurgical, quarrying, and other mineral industries; and to inquire into the economic conditions affecting these industries. The organic act of the Bureau, as amended by Congress and approved February 25, 1913, made it the province and duty of the Bureau to ""disseminate information concerning these subjects in such manner as will best carry out the purposes of this Act.""In accordance with this directive, the Bureau reports the findings of its research and investigations in its own series of publications and also in articles that appear in scientific, technical, and trade journals; in proceedings of conventions and seminars; in reference books; and in other non-Bureau publications. The number of these reports, the wide range of subjects they cover, and the variety of mediums in which they appear make the kind of list and index presented in this special publication both necessary and valuable. This issue describes reports and articles published during the period January I, 1965, to December 31, 1969. It supplements the 50-year list of Bureau publications issued from July 1, 1910, to January I, 1960;' the 50-year list of articles by Bureau authors published outside the Bureau from July 1, 1910, to January I, 1960;3 the 5-year list of Bureau publications and articles from January I, 1960, to December 31, 1964.' It includes all the material in the four annual lists of Bureau publications and articles from January 1, 1965, to December 31, 1968, as well as the Bureau publications and articles for 1969."

Jan 1, 1970

By D. Bienstock

A rig was constructed to burn 100 to 125 lb of coal per hour at 4,000° F in a cyclone burner with oxygen-enriched air preheated to 1,500° F to ascertain the problems in the coal-fired MHD generation of power. Potassium carbonate was added to the coal at several concentrations up to 1.23 g-moles/kg of coal (17 lb K2C03/100 lb coal). Boiler-tube materials maintained at wall temperatures of 800° to 1,500° F were exposed to the combustion effluent of 1,800° to 2,500° F. Practically complete removal of sulfur dioxide from the stack gas was achieved with the seed addition. Nitric oxide formation followed theoretical expectations at combustion temperatures with some decomposition of the NO occurring as the gas cooled flowing through the system. Nitrates in the fly ash were less than 1.7 percent of the total fixed nitrogen. With aqueous extraction of the fly ash a maximum recovery of 75 percent of the potassium was obtained. Haynes 25 was unattacked at a wall temperature of 1,500° F in combustion gas at 2,500° F; the stainless steels 310, 316, and 446 were resist-ant at a metal temperature of 1,100° F in gas at 2,100° F.

Jan 1, 1970

By Y. Oka, H. Majima

When a rock particle is subjected to an external force, internal stresses and deformations are experienced by the particle. Assuming that the breakage of a rock is mainly due to tensile stresses generated by the action of a pair of concentrated loads, F, directed in compression towards the loading points at the shortest distance, as illustrated in Fig. 1, the stress and deformation experienced by a particle of irregular shape can be determined with the aid of the theory of elasticity. According to the photoelastic experiments made by one of the present authors,' the stress states in a rock particle of irregular shape are quite similar to that of an inscribed sphere of diameter x, as indicated in Fig. 1, but accompanied by the faint stresses, in the vicinity of the particle outside of the inscribed sphere. Therefore, the stress components in an irregularly shaped particle can be estimated by referring to stress components in a sphere. Applying polar spherical coordinates to an elastic spherical body, as shown in Fig. 1, the stress components are all expressed as multiples of F /x2 and more or less influenced by the loading angle and Poisson's number m. In general, the energy w, acting on a rock particle is defined as the product of loading compression F

Jan 1, 1970

By Peter N. Calder

The orientation of planes of inherent weakness (joints) with reference to an excavation in rock is the critical factor determining stability. The purely frictional 'residual' strength remaining following the breakdown of initial cohesion is the only strength which can be counted on to resist sliding for the long term. A simple method of resolving driving and resisting forces in three dimensions along joint intersections 'undercut' by the excavation is presented. Contrary to the case of slope failure in soils, it is shown how significant changes in the slope angle sometimes have no effect on the degree of stability pre-sent in rock slopes. Wall control by reorientation of the slope, the use of cable bolts and water diversion are briefly discussed.

Jan 1, 1970

By R. C. Earlougher

Jan 1, 1970

By J. R. Hoskins

This is the final report for Research Contract No. 00101300 (Min. 34) granted by the U.S. Bureau of Mines to the University of Idaho Mining Research Bureau. The research project is entitled "Development of Improved Ventilation Technology for Non-coal Mines and Other Underground Excavations.' The duration of the project was six calendar months, from 16 April to 16 October 1970. The budgeted cost was $34,867.00 of which $31,953.00 was to be furnished by the U.S. Bureau of Mines. The sequence of reporting that follows is to state the objectives as contracted, outline the program of investigation, highlight briefly the important technological factors peculiar to vein-type metal mine ventilation, and describe the normal methods of control. This is followed by the results achieved by the survey, and an itemization of the observed deficiencies in current technology. There follows a section dealing with the adopted Federal Regulations.

Jan 1, 1970

By L. A. Price

With increased stripping ratios and with dwindling ore content, the mining fraternity is increasingly turning to larger and larger equipment for profitable open pit operations. The availability of higher horsepower prime movers and increased-capacity electric transmissions using either wheel-mounted or axle-mounted traction motors makes a reality of single-axle, rear dump trucks having 150 to 250 ton capacity. The size of electric loading shovels is inherently geared to the size of haulage equipment. To work effectively with this new generation of trucks, the Harnischfeger Corporation has introduced its P&H Model 2800 with a capacity of 20 to 25 cu yd and a working weight of nearly 1,500,000 lb. The Model 2800, depicted in Figure 1 loading a 200 ton truck, has many new electric and mechanical features designed to give excellent production and low cost operation. Chief among the new features is the ELECTROTORQUE control system. Electrotorque is a completely static, solid state electronic system supplying controlled power both to the armatures and to the fields of the d-c shovel drive motors. It eliminates many major electrical rotating machines and their attendant maintenance. The system was developed jointly by ASEA of Sweden and Harnischfeger Corp.

Jan 1, 1970

By James D. Cooper

Technologic advances in materials utilization during the past two decades have brought is a long way toward the point where we can make the products we need from the materials we have readily available. Notable examples are the synthesis of diamond, zeolites, quartz crystal, mullite, and mica, and the development of phosphate bonded tale. Two of these--zeolites and mullite--involve the use of clays. There is little reason to doubt that equally important advances will be realized in the next decade, and that clays will play an important role in the developing materials technology. Indeed, the stage is already set for many of them. This discussion is concerned primarily with the high quality, relatively high priced clays, kaolin, bentonite, fuller's earth and ball clay. The first three have enjoyed demand growth rates over the past two decades which on the average exceeded the growth rate for the economy as a whole. There are indications that demand for ball clay will be growing at a faster rate in the future.

Jan 1, 1970

By M. H. Yoo

Possible interactions of the perfect dislocations of six slip systems or the c dislocation with the (10i2f (ioii), {ioIi}(ioiZ), {1122}(1123), and {1121}(ii26) type twins in hcp metals have been analyzed from the crystallographic and the energetic points of view. Twenty-six distinct types of possible interactions were identified, and those selected based on crystallographic constraints were examined for their energetic feasibilities by use of the anisotropic energy factors. No long-range elastic interaction exists for a dislocation when its Burgers vector is parallel to the twin interface. Under a suitable applied stress, a screw dislocation can cross slip at the twin interface. For basal mixed dislocations in cadmium and zinc, the interaction with {1012} twins is found to be attractive, indicating that incorporation of these dislocations into the twins is energetically feasible and that twin growth will result. On the other hand, the interaction between both basal and Prism mixed dislocations and the {1012} and (1121) twins is found to be repulsive in Mg, Co, Re, Zr, Ti, Hf, and Be. This indicates that under an applied stress a local stress concentration will develop due to a dislocation pileup at the interface, which may result in a site for either the nucleation of other twins or the formation of a crack, depending on the cleavage strength. WHEN a metal undergoes plastic deformation, a certain configuration of slip dislocations will result in a state of dislocation pileup against an obstacle. The stress concentration thus developed may enhance the process of twin nucleation and also twin growth. Furthermore, once formed and dispersed in the crystal, twins can act as effective barriers against slip dislocations. The degree of such mutual influence or interrelation between slip and twinning is generally known to be pronounced in the case of hcp, metals. It is also known that deformation by twinning occurs more commonly in hexagonal metals than in cubic metals. In fact, under suitable stress states, all hexagonal metals exhibit {1012) <1011> type twinning.&apos; In addition to this common type, deformation by (1151) <1126> type twinning occurs in zirconium, titanium, and rhenium, which show remarkable ductility.&apos; The importance of twinning during general deformation to the ductility of hcp polycrystals has been briefly discussed in recent review works.2&apos;3 The purpose of this paper is to analyze the interaction between slip dislocations and twins in the hcp structure and to discuss the nucleation and growth processes of twinning and the role of twinning in the <"°" noil) o, 1/3[112O] (OOO2) 1/3[1123] Fig. l—-Slip systems in hcp structure. ductility of hexagonal metals. The problem will be discussed from the geometric and the energetic points of view in a manner similar to that of the previous work on zinc.4 Since hcp crystals deform by several slip and twin systems, numerous interactions result as possibilities. The Burgers vectors of six slip systems and the c dislocation shown in Fig. 1 and the four twin systems listed in Table I are considered here. A complete tabulation of the possible interactions is followed by discussion of those that are more likely to occur on the basis of crystallographic constraints and energetic considerations. 1) CRYSTALLOGRAPHY OF TWINNING The crystallographic elements, K1, K2, n1, and n2, for the four compound twin systems are now well established.= A unit cell with the base vectors n1, and n2 is shown in Fig. 2 for each twin system. The unit cell before twinning is shown in solid line, and the corresponding unit cell after twinning is shown in dashed line. Also shown in Fig. 2 are the following crystallographic parameters: S is the plane of shear, d the interspacing of the twin habit planes K1,Ø Iis the acute angle between n1, and n 2, e is a numerical factor, and q is the number of K, lattice planes intersected by 17&apos;. These parameters can be expressed in terms of the axial ratio, y = c/a, as listed in Table 11. The macroscopic shear strain of twinning, s, and the magnitude of a "unit twin dis-l~cation,"4 bt, are also expressed in terms of y and given in Table 11. In Table 11, K1 and q1 are given in both Miller-Bravais and Miller indices. In double lattice structures, shuffling of atoms in addition to a homogeneous shear of the lattice is generally required if the original crystal structure is to be restored after twinning. The extent of current understanding on this problem of atom shuffling is per- Table I. Four Twin Systems in Hcp Structure

Jan 1, 1970

By W. D. Bachman

A study has been made of wet, autogenous grinding of disseminated copper ores, including testing of a large number of samples from Kennecott Copper Corporation&apos;s Chino mine. The efficiency with which the various samples could be ground was found to vary over an extremely wide range, even among samples taken from different locations within the same ore body. These variations in grindability can be qualitatively correlated with the composition of the rock and the fracturing and alteration that had occurred. Because of the large-scale non-homogeneity of disseminated copper ore bodies, extensive and careful testing of autogenous grinding is needed. An accurate correlation of the grinding data from various samples with the geology of the ore body and with mining plans will indicate if variations in grinding mill throughput can be avoided in practice by mixing or blending of mill feed and if, therefore, autogenous grinding can be applied successfully.

Jan 1, 1970

By Horst Huettenhain

The Rowland Coal Preparation Plant is located on Workman&apos;s Creek in Raleigh County, West Virginia, approximately 20 miles Northwest of Beckley and 10 miles Southeast of Ameagle. It is owned and operated by the Rowland Division of Consolidation Coal Company. This paper describes the plant flowsheet, process considerations that were evaluated in its development, and observations made during and since the plant was constructed. The final flowsheet was developed specifically to wash the coal from the Rowland Mine and contains several equipment combinations and applications heretofore never placed in a single, commercial operation. The major process equipment is presented in a functional description of the flowsheet. The coal reserves are leased from the Rowland Land Company and consist primarily of #2 Gas and Middle Eagle coals and various major splits of the Middle Eagle seam. The Rowland Coal Preparation Plant was designed to process 560 TPH of raw coal consisting of 50 percent #2 Gas coal intimately mixed with 50 percent Middle Eagle coal and containing lumps no larger than a 12-inch cube.

Jan 1, 1970

A mechanised sub-level overhead benching method is used at Cleveland, to produce 22,000 tons per month of ore assaying 0.92% Sn and 0.40% Cu. The ore occurs in a series of vertical lenses which can be mined by open stoping. Levels are developed at 60ft., vertical intervals by jumbos. The ore is silled out and the stope block formed between levels drilled out by long hole drilling. Each row of holes is fired and broken ore within safe reach from the stope brow is mucked out by Wagner ST-5A Scoop Trams. The remaining broken ore is recovered on extraction of the next lower level. Ore is hauled by road to the crusher station from the mine. As mining proceeds below the lowest surface portal a spiral decline will be used for ore haulage and mine servicing.

Jan 1, 1970

By G. H. K. Schenck

Weather changes of a daily, weekly, or seasonal period affect mineral production and cause significant losses that are considered unavoidable. Private weather forecasting services are available on a contract basis that can significantly reduce the uncertainty associated with weather-induced changes in the production of and demand for minerals. This paper provides examples of current practice and reviews the opportunities that await those producers who incorporate weather services as an integral part of their operating, scheduling, and planning procedures. Private weather forecasts can be highly accurate because atmospheric trends such as conditions and movement tend to persist and are related to previous meteorological events. Commercial meteorological services are prepared to deal specifically with a mineral producer&apos;s special problems (Table 1). Tactical use of private forecasts is a major advance over the severely restricted value of public forecasts which, of necessity, are general in nature and cover a dispersed geographic territory. Private weather services provided on a consulting fee basis offer the mineral producer the opportunity for tactical and strategic response to weather conditions at an acceptable cost. This contrasts with the too-frequent current tendency to view weather as a random variable the costs of which are unknown and about which nothing can be done.

Jan 1, 1970

By Martin E. Chenevert

The early work of Terzaghi1 showed that the true stress state of a rock can be expressed by the principle of effective stress. This principle states that effective stress is equal to total applied stress minus the pore pressure. Since the introduction of this principle, a great effort has been made to measure the pore pressure of rocks so that the true stress state can be determined. For clean, permeable sandstones, it is usually found that the effective stress principle can easily be applied. Direct measurement of the pore pressure is achieved with a hydraulic pressure device and the effective stress state thereby computed. However, for fine-grained rocks which contain large amounts of clay (referred to herein as shales), pore pressure measurements are usually quite difficult because of the adsorptive forces of the clay platelets, which produce negative internal pressures. No suitable device is known for measuring internal pressures of dry shales; the effective stress approach is therefore not feasible and data are usually presented in another form. As an example, data from compressive strength tests are presented as a function of rock wetness, rather than effective confining stress. The work of Belopoloskii2 and Colback and Wiid3 are examples of such data. It was pointed out by Lambe and Whitman4 that a better way must be found to express the effective stress of argillaceous rocks. The information presented herein is an attempt to fill that need. Various authors have been able to measure moderate negative pore pressures (0-240 psi) within rather wet, argillaceous rocks by means of various pore pressure devices. In 1961 Lambe5 resented information which shows that pore pressures in argillaceous rocks are indeed negative ; he presented values of suction pressure as low as - 14 psi. Likewise,

Jan 1, 1970

By A. Timur

Fluid flow properties of porous media have long been of interest in such varied disciplines as geology, geophysics, soil mechanics, and chemical, civil, and mechanical engineering. This interest has resulted in numerous models of porous media that have been proposed, tested, and found to be useful under, at best, only special circumstances. A critical review of most of these models is given by Scheidegger. 1 The application of each new tool to the study of porous media has helped to test some of the existing theories and to form new ideas based on the new parameters being measured. In the application of nuclear magnetic resonance (NMR) techniques to the study of the properties of fluids in porous media, the theoretical studies of Korringa et al.2 resulted in a model for the relaxation of spin polarization of protons in a hydrogenous fluid in a pore of a solid. Seevers3,4 verified this model, established a method for measuring the surface-to-volume-ratio distribution in a porous medium, and proposed a technique for determining specific permeability in sandstones. In the present investigation, the uses of NMR methods are extended and a simple model of porous media is developed from an analysis of pulsed NMR measurements. In this model, the pore spaces of a porous medium are divided into three groups on the basis of their surface-to-volume ratios. To evaluate predictions of the NMR model, laboratory measurements of spin-lattice relaxation time, porosity, permeability, and residual (irreducible) water saturation were conducted on more than 150 sandstone samples obtained from four oil fields in North America. Applications are discussed for estimating the volume of movable fluid, ?p, and the specific permeability, k, of sandstone samples. The former parameter, which can be considered as producible porosity, is expressed by ?p = [1 -- (Swr/100)]? , .... (1) where + is porosity in percent of bulk volume and Swr is the residual (irreducible) water saturation in percent of pore volume at a capillary pressure of 50 Psi. Fast, nondestructive laboratory methods are described for determining porosity, movable fluid, and permeability of sandstone samples from measurements with a pulsed NMR (spin-echo) apparatus. Procedures are given for adapting these methods to samples of irregular shape and small size. Theory General In the Korringa, Seevers, and Torrey2 (KST) model the observed decrease in the spin-lattice relaxation time, T1, of protons of a hydrogenous liquid contained in the pore spaces of a porous solid is attributed to an increase in the correlation time for the random motion of the water molecules and to the presence of paramagnetic centers at the liquid-solid interface. This KST model predicts the observed relaxation time of a hydrogenous liquid in a pore of a solid through

Jan 1, 1970

By D. S. Pye, J. P. Gallus

Successful hydraulic pressure treatment of oil and gas wells that have long producing intervals, multiple pay zones, fluid thief zones, or natural or induced fractures has long been a serious problem plaguing the petroleum industry. High fluid-injection rates, straddle packers, ball sealers, and either liquid or solid additives to the treating fluids are being used to cope with this difficult and costly problem.l-6 The use of ball sealers and straddle packers is limited, in practice, to wells that have been completed with perforated cemented casings or liners. Well stimulation with ball sealers and straddle packers is usually successful as long as the cement remains competent. The other techniques of diverting are used with all types of completions. However, there are two significant disadvantages of the high rate method: excessive horsepower requirements, and the need of large treating-fluid volume. Consequently, most well stimulation techniques that involve diverting depend upon ball sealers or granular materials. Liquid diverting agents, and some of the solid ones in use today, frequently will cause formation damage. On the other hand, the most popular particulate diverting agents often are so soluble in the treating fluid that they become largely ineffective; their particle size is severely reduced as they dissolve while being transported through thousands of feet of pipe to the formation face down hole. In our attempts to deal with the shortcomings of these diverting agents, we experimented with a large variety of potentially suitable deformable materials the solubility of which could be controlled. Following are the criteria for the diverting agent sought: 1. The material must be effective in reducing fluid loss to thief zones, including existing fractures and very high permeability layers in the formation to be pressure treated. 2. The material must have enough strength to bridge under the treating pressure. 3. The material must be insoluble or slowly soluble in the treating fluids so that the particles retain sufficient size to be effective after they reach the formation face. 4. The material must be completely soluble in the produced oil or condensate liquid so that it can be removed readily; it must not cause formation damage. A number of products consisting of various polymers blended with waxes and other solid hydrocarbons were examined. All these materials satisfied the criteria for a diverting agent and also exhibited some degree of plastic deformability. On the basis of initial tests of strength and hardness, several blends of wax and polymer were selected for extensive testing in the laboratory and the field. These compounds have congealing temperatures ranging from 140 to 190F; they are insoluble in water and soluble in oil. They conformed, therefore, to the criteria specified for a theoretically superior diverting agent. Laboratory and field tests were conducted using these wax-polymer materials to demonstrate that the use of a deformable diverting agent can significantly improve fluid diverting performance. Data in this

Jan 1, 1970

By C. F. Weinaug, P. M. Blair

Many difference equations used to approximate reservoir flow problems treat the phase pressures implicitly but not the mobility-density coefficients. Such difference equations are neither wholly explicit nor implicit, but might be described as mixed. Mixed equations are relatively easy to apply. But the associated time truncation error is relatively large, and when used to solve problems characterized by high flow rates, these equations may be unstable for practical size time steps. This paper outlines the development of a completely implicit difference analogue for reservoir simulation, along with a Newtonian iterative method for solving the resulting nonlinear set of algebraic equations that arise at each time step. While this implicit equation requires two to three times more work than does the mixed equation, it is shown to markedly decrease the time truncation error and to yield a stable solution for much larger time steps than does the mixed equation. INTRODUCTION Calculation of multiphase, multidimensional flow in porous media is generally accomplished by numerical methods that involve approximating systems of partial differential equations by systems of partial difference equations. The early development of difference equation techniques was directed toward the solution of linear differential systems. However, the equations of multiphase fluid flow through porous media are highly nonlinear in that mobility and density often are strong functions of pressure and saturation. Thus, solution of the flow problems by difference equations involves solving sets of algebraic equations whose coefficients change from step to step of the calculations. Current practice is to evaluate most coefficients at the beginning of a time step and then to apply difference methods well suited for solving linear problems. At least two major difficulties arise from the use of this practice. First, evaluation of coefficients at the old time level and evaluation of pressures at the new time level results in larger time truncation errors near displacement fronts than if all quantities in the distance derivative were evaluated at the same time level. Second, evaluation of mobility coefficients at old time levels results in an unstable difference equation in regions of high flow rate. The first of these difficulties often results in optimistic recovery behind the displacement front, and the second requires the use of small time steps or special techniques in solving coning problems1 or problems of gas percolation.2 Evaluation of all quantities in the distance differences at the new time level results in a completely implicit difference system. Such a difference system has a lower time truncation error than equations in which mobilities are evaluated at the old level and pressures evaluated at the new level. The algebraic equations resulting from the fully implicit difference equation are, in this case, nonlinear and require some iterative method for their solution. This paper develops fully implicit difference equations for two-phase flow in porous media, describes their solution by Newtonian iteration, and gives examples of problems more easily solved by the new method than by previous methods. EVOLUTION OF DIFFERENCE EQUATIONS FOR NONLINEAR PROCESSES The completely implicit difference equation used in these calculations results from evolution beginning with explicit difference equations, proceeding through what we will call "mixed" equations, and arriving now at a completely implicit equation. Let us illustrate the evolution by considering different analogs of the nonlinear equation

Jan 1, 1970

By L. B. Davidson

The petroleum literature contains many reports on the relative permeability properties of porous media. However, only recently have studies of relative permeability at elevated temperatures been made.1 -5 Each of these studies presents relative permeability data on water-oil systems in linear sandpacks or sandstone cores under isothermal conditions. The results of these investigations are dissimilar. Edmondsonl found that in Berea sandstone a temperature-dependent effect was present, while Shilobod,2 using the same material, found no temperature dependence after normalizing for initial water saturation. Workers at Texas A&M U.3-5 found that the ratio kw/ko was temperature dependent in an unconsolidated system. The differences in results of these studies are probably a consequence of differences in the material employed, in the experimental procedures followed, and in the properties of the porous systems used. This paper describes an experimental investigation of permeability ratio temperature dependence, conducted to clarify some of the results obtained by previous workers. Each experiment described in this paper involved the isothermal displacement of Chevron No. 15 white oil by either nitrogen, steam, or distilled water and displacement of water by nitrogen. All permeability ratios were calculated using the Welge procedure. Values of the water-oil permeability ratio were obtained over the temperature range 75OF to 540°F. Results indicate the ratio is temperature dependent at low water saturations, probably due to interfacial effects. However, as the saturation is increased, temperature dependence decreases until, above some minimum saturation, the ratio becomes insensitive to temperature. Previously published work indicates that the dependence reappears at very high water saturations, a region of the ratio curve not completely examined in the present study. It appears that water-oil permeability ratios are temperature dependent at the extremes of the ratio curve where changes in interfacial tension with temperature induce variations in residual saturations. Between the extremes, there is a region in which the permeability ratio is insensitive to temperature. The size and location of this region is probably a function of the characteristics of the porous system. Gas-oil permeability ratios calculated from nitrogen-white oil displacement experiments indicate a definite temperature dependence over the range 75OF to 500°F and over the entire gas saturation range studied. Evidence discussed in this report indicates that this dependence is a consequence of molecular slippage in the gas phase. This is an extension of the Klinkenberg effect6 to two-phase systems at high temperatures (the existence of slippage in two-phase systems at room temperature has been demonstrated7). Evidence for a similar source of temperature dependence in superheated steam-white oil experiments

Jan 1, 1970