Search Documents

By G. W. Nabor, A. S. Odeh

The effect of production history of a well on the results of two-rate flow tests, and conventional build-up analyses was investigated. The effect was examined by means of digital computers and an R-C network model, respectively, for wells with infinite and finite radii of drainage. For systems which behave as infinite, it was found that production history and the duration of production at constant rate prior to the initiation of the test have important effects on the results. During build-up time equal to about one-fourth of the stabilized time, correct permeability-thickness product calculations can be made. For wells with finite radii of drainage, the time was determined during which the straight line can be satisfactorily used for permeability-thickness product calculations in case of drawdowns and build-ups. On build-ups, the dimensionless time (based on the external radius) during which the straight line gives reliable results was detertriined to be about 1/12. This is one-fourth as long as that of the drawdown. The investigation was done theoretically, and subsequently was verified by R-C network model runs. General interpretive rules were formulated which, if not followed, could lead to serious errors. Moreover, a recommended testing procedure is reported. INTRODUCTION The method used by most reservoir engineers for estimating formation characteristics in a producing well is the analysis of pressure build-up data. The method originally devised by Horner' makes use of the point source solution to the diffusion equation. This solution is approximated by a logarithmic function and the superposition principle is employed to arrive at the well known pressure build-up equation:where q, the flow rate, is in reservoir B/D; ft is in cp; kh is in md-ft; At is the shut-in time; and t is the producing time. At and t are in any consistent time units. Ey. 1 is applicable to a well of unlimited drainage radius which produces at a constant rate q from zero to time t and is then shut in. Such a constant production rate seldom obtains in practice. Therefore. a correction term must be applied to Eq. I to account for the varying rate. Two theoretically accurate methods are available for treating the variable rate case. The first, originally derived by Horner,' is based on the application of the superposition theorem. It requires knowlege of production history of the well as a function of time and results in lengthy and laborious calculations. The second t*q* niethod is suited for short production tests and requires that the shut-in time be at least one and one-half times the production time. A third method which is not based on any theoretical justification and which was suggested by Horner as a "good working approximation" is the one used by the majority of reservoir analysts. especially when the well has been producing for a long time and the t*q* method is not practicable. The key to this method is in choosing or determining the t that appears in Eq. 1. Horner suggested using a corrected time t, in place of t. t, is calculated by dividing the total cumulative production by the last established rate. Therefore, a normal procedure of pressure build-up testing is to stabilize the well at a constant rate for at least 24 hours before shut in and to use the stabilized rate to calculate t,. The analysis is then made by plotting either or P. and examining the resulting plot for the expected straight line to calculate kh and the original reservoir pressure. Recently, Russell" proposed a method for determining formation characteristics from two-rate flow tests. His method reduces to pressure build-up if the second flow rate is zero. Russell uses the Horner simplified procedure for calculating a corrected t,. His method also requires the stabilization of the well at a constant rate q which is used to calculate t,. Theoretically, the above procedure is valid for a well with an unlimited radius of drainage or with a limited radius as long as the boundary effect has not been felt by the well. Several authors'' derived formulas which allow the estimation of time during which limited reservoirs behave as infinite ones and, thus, can be treated by unsteady-state mechanics. One equation derived by Swift and Kiel' terminates the application of unsteady-state theory when the drainage radius reaches one-half the reservoir radius. Thereafter, steady-state behavior obtains. Another equation derived by Jones'" initiates steady-state

Jan 1, 1967

By Rex B. McLellan

A simple model for dilute interslitial solid solutions is set up which enables a solubility equation for the equilibrium between the solid solution and the gaseous solute to be deduced. This solubility equation is used to estimate the partial thermody -namic functions of Ag-O and Fe-N dilute interstitial solid solutions from the measured equilibrium between these solutions and the gas phase. It is shown that the partial excess entropy of the solute atom in solution is a vibrational entropy which is due Principally to the oscillation of the interstitial atom and does not arise from electronic or elastic interaction between the solute atoms and solvent matrix. In this work measurements of the equilibrium between dilute solid interstitial solutions and a gaseous phase are used to deduce the important thermody-namic functions of the solid solution with the aid of a theoretical solubility equation formulated from a knowledge of the thermodynamic functions of the gas and a simple model for the solid solution. The quantity of principal interest which is deduced from the equilibrium data is the change in vibrational entropy occurring when a solute atom (u) is placed in an interstial site in the solvent (u) lattice. This paper forms part of a more general study in which the properties of condensed phases are deduced from measurements of the equilibrium between the condensed phase and a gas. Previous investigations have dealt with the equilibrium between krypton and liquid metals,' a pure metal and its vapor,' and dilute substitutional alloys and the vapor,3 and the subject has been popularly reviewed by McLellan.4 In the latter investigation3 the vapor pressures of very dilute Cu-Ag and Cu-Au alloys were measured with a Knudsen cell technique over a large range of temperatures and the analysis of the vapor-pressure data showed that the relative partial excess entropies were large and positive [ds= -Sudk = 30 for Au-Ag alloys and 150 for Cu-Au alloys; these values should be compared to unity, the value assumed for these exponentials in the quasi-chemical theory of solutions]. For many interstitial solutions experimental data is available for the equilibrium between the solution and a gaseous phase, and this can be used to deduce the vibrational entropies of the solute atoms in solution. The general method employed is to calculate the chemical potential of the solute atoms in the gas and equate this to pi, the chemical potential of solute atoms in the solid, which is deduced from a simple model set up for the dilute interstitial solution (in which the solute atoms act as bound oscillators in interstitial sites). This enables a theoretical solubility equation to be set up. In the case of Ag-O and Fe-N solid solutions the solubility c, is found to vary with temperature according to a relation of the type, where PU2 is the pressure of the O2 or Np molecules. Essentially, the theoretical solubility equation enables the parameter A(T), which gives the vibration entropy, and the parameter AH, which gives the energy of solution of a solute atom, to be estimated from the measured variation of solubility with temperature. The solutions treated are dilute enough to enable solute-solute interactions to be neglected so that the partial energy and excess entropy of solution E, and GS are not functions of composition and the configurational entropy is that of a random solution. Large positive values have been found for sxSu this quantity has not been given sufficient consideration in recent work on interstitial solutions and it is felt that the partial thermodynamic quantities of very dilute solutions should be understood before the thermodynamics of concentrated solutions can be fully appreciated. EXPERIMENTAL DATA 1) The Silver-Oxygen System. The measurements of Steacie and Johnson,= using silver foils, indicate that there is a pronounced minimum in the solubility vs temperature curve for O-Ag solutions. Eichenauer and Miiller6 have recently redetermined the solubility of oxygen in silver using a method developed in the investigation of metal-hydrogen systems.7 The solubility determinations were carried out on compact specimens of differing shapes and sizes and the results indicate a much lower solubility than was obtained by Steacie and Johnson with no minimum at 400°C. Subsidiary experiments showed that the erroneous results of Steacie and

Jan 1, 1964

By B. C. Allen

Grain-boundary gvoocing was studied irz chronziu?n. molybdenum, tungsten, and the solid-solution alloys, Cr-35Re. Mo-33Re, and W-25Re at 0.6 to 0.9 of the absolute liquidus temperature under an inevt atmosphere , foveign nzetal zlapor, and liquid metal. The controlling mechanism of matter transport was deduced frolr groove size and shape. kinetics of growth, and cotnparisoz with theoretical relatiotzships derived by Mullins. Grooves in chowzizon formed prittzarily by volume diffusion tkrongk kelium or argon at 1 atm pressure and by e1)apovation-condensatiotz in a static vacuum. Groozles in the other nzetals and alloys forrned priniarily by sirface diffusion, the coefficients of which were calculated. Carbon was found to decrease the acti7:ation energy for surface diffusion on molybdenun. Sitrface diffusizities weve sinzilav for iWo-0.0ZC and Mo-33Re. The presence of about 1 -?nm pressure of silver or tin vapor did not affect the neckanis)n of grooz;e formation. In liquid titz or silz:er, pain bolindary gvooves formed by zolume diffusion in the liquid. BASED on theories of ullins," the mechanisms of grain-boundary grooving have been determined for low-melting metals such as copper in an inert atmospheres or in liquid lead.4 The purpose of this investigation was to study grain-boundary grooving in the refractory Group VI-A metals, chromium, molybdenum, and tungsten, and the bcc solid-solution alloys, r-35e, 0-33e,= and -25e,' at 0.6 to 0.9 of included to see whether Mullins' theories were applicable to them as well. EXPERIMENTAL WORK The materials chromium, molybdenum, tungsten, Cr-35Re, Mo-33Re, and W-25Re were in the form of 0.050-cm-thick sheet assaying 99.98 pct or better. The nominal carbon level was generally 30 ppm. The oxygen + nitrogen + hydrogen content was -50 ppm. Mo-33Re and W-25Re were prepared by powder-metallurgy methods and supplied by Chase Brass and Copper Co. The remaining materials were arc-cast. Chromium and Cr-35Re were prepared at Battelle by arc melting, extruding, swaging, and rolling. Mo-0.003C was supplied by Climax Molybdenum Co., and Mo-0.02C and tungsten by Universal Cyclops Corp. As a source of liquid metals, tin and silver granules assaying 99.98 pct were used. Spectrographic analyses indicated the presence of -100 ppm each of iron and silicon. Vacuum melting reduced the combined gas content to c5 ppm (2.60,, 0.9N2, 0.4H2 in tin and l.602, 0.7N2, and c0.04H2 in silver). Specimens were prepared for grain-boundary grooving. Coupons 0.5 by 1.5 cm with a small hole in one end, or substrates 1.5 by 1.5 cm, were cut from sheet. The specimens were annealed in containers of the same material to provide a vapor-solid equilibrium as well as to protect them from gaseous impurities. Outgassing was done at 2 x 10"5 mm for 15 to 30 min at 1150°C for chromium and Cr-35Re, 1750°C for molybdenum and Mo-33Re, and 2600°C for tungsten and W-25Re. Then recrystallization to a stable grain size of about 0.5 mm was accomplished by annealing at least 1 hr at 1600°C for chromium, 1800°C for Cr-35Re, 2350°C for molybdenum, and 2300°C for Mo-33Re in a rhenium-element resistance furnace under 1 atm of static argon or helium. The atmosphere was gettered by tantalum radiation shielding. Tungsten and

Jan 1, 1967

By A. S. Odeh

This paper analyzes the effect of limited entry to flow at the wellbore on the steady-state productivity of a well. Wells that have been opened to flow along a fraction of their productive interval are termed wells with limited entry. Previous work treated the cases of a partially penetrating well, a well producing from the central portion of the productive interval and a well in which several intervals equally spaced were open to flow. In this paper the open interval can be located anywhere within the productive interval. Thus, in a sense, it generalizes previous work. The finite cosine transform was used to arrive at a solution for steady-state flow of a slightly compressible fluid. The solution was programmed for a CDC 1604 computer. Numerical vaIues for rd = 660 ft, r, = 1/4 ft, and range of sand thickness of 20 to 200 ft are presented in graphical form. The effect of rd and r, values on the result is shown in a table. The correct calculation of skin and damage ratio in the presence of limited entry to flow is explained and illustrated by examples. Moreover, the paper shows how to calculate the net decrease in productivity due to the combined effect of limited entry and perforations. INTRODUCTION In some wells only a fraction of the productive interval is open to flow. Location of this fraction is usually dictated by formation characteristics and reservoir behavior. For instance, if a gas cap exists, the open interval is located away from the gas-oil contact to prevent any possible gas coning. Wells that intentionally have been opened to flow along a fraction of their productive formation are tened wells with limited entry. Obviously, unintentional completions of this type also exist. Limited entry to flow decreases well productivity. Magnitude of the loss depends on the fraction of the formation open to flow, on the thickness of the sand, on the location of the open interval and on the ratio of rd /r, where r , is well radius and rd is the drainage radius of the well. The use of pressure buildup data on producing wells to calculate the condition of the formation around the wellbore is an accepted practice. van Everdingenl and Hurst2 introduced the concept of she skin factor s considered to be due to a thin layer of different permeability immediately around the wellbore. These authors dealt with the case of a well of complete radial geometry, i.e., a well with open-hole completion that completely penetrates the formation. The presence of a low-permeability skin results in a loss of productivity, as does limited entry. Therefore, if pressure buildup data obtained on a well with limited entry are used to establish the presence or absence of skin (i.e., formation damage), and a correction is not made for this loss of productivity, the calculations would result in an erroneous skin value. They might indicate the presence of formation damage when in reality there is none, or they might indicate a value larger than the true value. This could lead to an incorrect basis for planning remedial measures. Muskat3 studied the problem of partially penetrating wells for the case of incompressible flow. He presented equations and figures which allow estimation of loss in productivity. Brons and Marting,4 using equations based on Nisle's work,5 studied the loss of productivity for three cases. The first was for a partially penetrating well; the second was for a well producing from only the central portion of a productive interval; and the third was for a well in which several intervals equally spaced were open to flow. Their work was for steady-state depletion-type reservoirs wherein the well radius of drainage is established and the fluid is considered to be slightly compressible. Considered in this paper is the problem of wells with limited entry in which the open intervals are located anywhere within the productive sand. The finite cosine transform is used to arrive at a

Jan 1, 1969

By D. Turnbull, R. E. Cech

FISHER, Hollomon, and Turnbull have developed a theory for the nucleation of martensite. They first tested the theory on Fe-C alloys and low alloy steels. The major factor influencing nucleation of martensite was considered to be statistical composition fluctuations occurring in small regions at high temperature and frozen-in on quenching. These local regions of varying size and composition serve as nucleation centers. They become supercritical, one by one, as temperature is progressively lowered, resulting in temperature-dependent or athermal transformation. Fisher next applied nucleation theory to substi-tutional solid solution alloys. Detailed predictions were made for Fe-Ni alloys because of the availability of free energy data on this system. It was shown that composition fluctuations that were significant energy-wise did not occur, and nucleation frequencies could be calculated from average properties of the system. Nucleation was predicted to occur as time-dependent and having the functional relationship to give a C curve of nucleation frequency vs temperature. The analysis further predicted that the nucleation frequency was extremely sensitive to composition. Experimentally, it would be found that the transformation in some compositions is so slow that measurable amounts will not form in a reasonable length of time. With other compositions, only slightly different, the nucleation frequency becomes so great that the material becomes transformed while still distant in temperature from the maximum nucleation frequency. On quenching an alloy of such composition, the observed transformation kinetics would be similar to those found in Fe-C alloys. Cech and Hollomon repeated experiments of Kurdjumov n which the kinetics of transformation were similar to those predicted by Fisher for Fe-Ni alloys. The alloy studied in this investigation contained 73.3 pct Fe, 23.0 pct Ni, and 3.7 pct Mn. Fisher,' using an idealized model for the extent of transformation as a function of the number of martensite crystals per grain of parent phase, derived nucleation frequencies from the transformation curves of Cech and Hollomon. Complicating influences such as coupling effects between grains in the polycrystalline specimens were neglected. Nevertheless, excellent agreement was found between the theoretically and experimentally derived nucleation frequencies. These experiments, however, could not provide a critical test of the theory. Experimental nucleation frequencies could vary widely from those calculated, depending upon the extent of deviation from ideal partitioning and the extent of coupling effects. Further, since the compositions of material theoretically analyzed and experimentally determined were different, the free energy changes involved in the experimental work could only be estimated. Also, the effect of heterogeneities on the transformation kinetics was not considered. For these reasons, it was decided that experiments designed to test the validity of the Fisher analysis must be conducted on binary Fe-Ni alloys, which were the ones considered theoretically by Fisher. The martensite transformation in Fe-Ni alloys has been the subject of considerable study. Machlin and Cohen have shown that transformation proceeds in a manner quite unlike that in any other ferrous alloy system. They found that single crystals would transform to a large extent in a single burst. In large grain polycrystalline specimens, frequently more than one grain and sometimes the whole specimen would transform at the same instant in this manner. Results on filings indicated that different particles would undergo the burst transformation at widely different temperatures. These results support the conclusion that the transformation behavior could not be described by a single nucleation frequency as would be the case if the nucleation were homogeneous. It appeared that further work was necessary to define the factors responsible for burst-type transformation, so that the conditions could be altered to favor homogeneous nucleation of martensite if such could be accomplished. This report summarizes the results of some experiments conducted with powdered Fe-Ni alloys for this purpose and the re-

Jan 1, 1957

By H. N. Hall

Reservoir and producing characteristics can govern the decision to use either a one-, two- or three-dimensional model for making predictions for gravity-drainage reservoirs. Examples of conditions requiring one-, two- and three-dimensional calculations are given. In 1961 the author presented a method for predicting one-dimensional gravity-drainage performance. This work has been exrendcd to obtain a three-dimensional model, which is described and for which a sample problem is presented. Introduction Many papers on gravity drainage or gas-cap drive'-3 have been presented for conditions where the reservoir is treated as a single areal segment so that, in effect, only one-dimen-sional flow is considered. Other authors'-Q ave given evidence showing that the one-dimensional concept would be unsatisfactory to predict reservoir performance for certain types of reservoirs. Experience has indicated that specific reservoir and producing conditions govern whether one-, two- or three-dimensional concepts must be used as a basis for reliable reservoir performance predictioils. Previously the author presented a paper outlining a method for predicting the performance of gravity-drainage reservoirs.' That paper treated the reservoir as a single areal segment so that it was applicable for only one-dimensional flow. The work has been extended to obtain a three-dimensional model that can be used to predict performance of gravity-drainage reservoirs. Although capable of being three-dimensional, it can be used equally well for one- and two-dimensional predictions. This paper discusses those reservoir and producing conditions governing the choice of model, and also describes a three-dimensional model suitable for predicting reservoir performance in the more complex situations. Discussion The Choice of Model To predict reservoir performance, it is necessary to represent the reservoir by a physical model. A mathematical model, based on the assumed physical model, then is used to predict the movement of oil and gas within the reservoir. The word "model" is used here to connote the combination of both a physical model and a mathematical model based on the physical model. Certain types of reservoir situations may be predicted suitably with a one-dimensional model, whereas two- and three-dimensional models may be required for other situations. This section illustrates reservoir conditions associated with the use of a specific type of model. Many types of reservoir structures are such that gravity drainage possibly could be an important factor in oil production. Production practices, as well as the nature of the structure. can influence the type of model that should be used. obviously, it is impossible to discuss all combinations of reservoir structure and production practices that might be encountered in gravity-drainage reservoirs. The following are typical examples of conditions where a choice must be made among a one-, two- or three-dimensional model to predict reservoir performance satisfactorily. Conditions Requiring a One-Dimensional Model Fig. la shows contours for an asymmetrical anticline in which gravity drainage might be expected to be important. Fig. Ib depicts a cross-section along Line A-A in Fig. la. Existence of a gas cap is a good clue that gravity drainage will be an important factor in producing oil from this type of reservoir. However, it is not essential since a gas cap will form as oil is produced from the reservoir. If the proper combination of reservoir permeability and withdrawal rates is present, the reservoir will behave like a

Jan 1, 1969

By R. J. Burgess, A. R. Ramey, A. R. Adams

During interpretation of pressure buildup tests on gas wells in a tight dolomite gas reservoir, peculiar behavior was noticed. Two straight lines were apparent. Effective permeability to gas taken from either straight line was about the same, and the Miller-Dyes-Hutchinson dimensionless time check for the straight line was proper for both straight lines. Geological data indicated the likelihood of scattered trending fractures in the reservoirs. Since the first straight Iine yielded permeability values close to the geometric mean permeability from core analyses, it was postulated that the reservoir model was that of an acidized well completed in the tight dolomite, but that widely scattered hairline fractures caused the mean permeability of the reservoir distant from the well to be higher than the matrix permeability. Because all other studies of fractured reservoirs to the authors' knowledge assumed that the fracture matrix was dense enough to communicate directly with the well, no interpretative methods were available. The Hurst line-source solution for a radial change in permeability for interference between oil reservoirs was adapted to pressure buildup testing. The result indicated that the first straight line should yield the proper matrix permeability and wellbore skin effect. The second straight line may be extrapolated to obtain static pressure. The time of bend between the straight lines was used to estimate distance to a fracture. Application to field test data is shown. It is believed that the methods developed and the case history presented will add to present tools available for pressure buildup interpretation. Introduction Since the pioneer studies by Miller, Dyes, and Hutchin-son1 and Horner' in 1950 and 1951, well test analysis has become recognized as one of the most powerful tools available to both production and reservoir engineers. Well test analysis serves as a logical basis for well stimulation and completion analysis, and for long-term reservoir engineering. Since the early 19501s, much effort has been placed on the development of well-test analytical methods. Reservoir and well conditions of increasing complexity have been considered systematically to provide the analyst with a catalog of causes and effects. Matthews and Russella state that some 200 papers dealing with this subject have been published in the last 35 years. Developments in well test analysis appear to have originated in one of two ways. Either a physically realistic field condition was anticipated and analytical solutions for the condition achieved, or anomalous field test behavior was recognized and interpretative methods sought for the anomaly. In recent years, it has appeared that the latter has inspired an increasing number of studies. The analyst today finds an increasing number of known cause and effect studies available for well test analysis, the classic of which is that of finding the specific flow problem that generated the answer — the well behavior. Although it may be impossible to achieve this goal uniquely, the analyst often is able to select a useful interpretation that combines all known performance and geologic data — or to show that various logical alternatives would not significantly affect the interpretation. During a recent reservoir study, we observed gas well test behavior that did not appear to fit behavior described previously. Although it cannot be said that we have found a unique interpretation, we shall present in this paper the peculiar behavior observed, and describe the reservoir and interpretative methods developed. Reservoir Description The subject gas reservoir is a 9-mile-long, narrow dolomite reservoir lying within a limestone of Ordovician age. (See Fig. 1.) The dolomitized rock in the field consists of dark brown to buff, dense to coarsely crystalline, vugular dolomite containing numerous hairline fractures, many of which may have been closed in the reservoir and parted when cores were brought to the surface. Larger fractures are also apparent in core, but usually are filled and sealed with euhedral dolomite crystals. Portions of the north flank of the reservoir are known to be cut by a sealing fault downthrown to the north. Gas wells located near the fault have higher open flow potentials than those more distant from the fault. This is believed to be a result of higher permeability near the fault due to more extensive and open fractures. Detailed coring and core analysis have been performed on several of the wells in this reservoir. Fig. 2 presents permeability variation' plots for both horizontal and vertical

Jan 1, 1969

By P. F. Gnirk, J. B. Cheatham

Single bit-tooth penetration experiments under static load were conducted on six rocks at confining pressures of O to 5,000 psi using sharp wedge-shaped teeth with included angles ranging from 30 to 120° In general, the force-displacement curves for all rocks exhibit an increasingly nonlinear and discontinuous behavior with decreasing confining pressure. The confining pressure at which a rock exhibits a macroscopic transition from predominantly ductile to predominantly brittle behavior during penetration vmies from about 500 to 1,000 psi for the limestones to greater than 5,000 psi for dolomite. The correlation between calculated values of force per unit penetration based on plasticity theory and experimental values is quite encouraging, even at confining pressures as low as 1,000 psi A qualitative correlation between volume of fragmented rock per unit energy input for a single bit-tooth and drilling rate for microbits appears to exist over a confining pressure range of 0 to 5,000 psi. INTRODUCTION Laboratory experiments utilizing a small-scale drilling apparatus have demonstrated that penetration rates are reduced considerably as a result of increasing the confining pressure from atmospheric to a few thousand psi?-3 This undesirable situation can, in general, be attributed to a combination of decreased efficiency of chip removal at the bottom of the borehole, increased rock-failure strength, and a possible change in the mechanism of chip generation and rock fragmentation with increasing confining pressure. To more fully understand the principles underlying the last circumstance, it is the purpose of this investigation to experimentally study the mechanism of single bit-tooth penetration into dry rock at low confining pressures and, in particular, to establish the confining pressure at which the penetration mechanism may undergo a brittle to ductile transition for various rock types commonly encountered in drilling. Confining pressure as used here refers to the differential pressure between the borehole fluid pressure and the formation-pore fluid pressure. EXPERIMENTAL PROCEDURE Using an experimental apparatus previously described: a single, sharp wedge-shaped tool was forced under a "statically" applied load into an effectively semi-in finite dry rock sample subjected to a prescribed confining pressure. To prevent the invasion of the con fining-pressure fluid into the pores of the rock sample during penetration, the exposed surface of the rock was jacketed with a layer of silicon putty.* Electrical instrumentation incorporated into the apparatus yielded a graphical plot of force on the tool as a function of penetration or displacement of the tool into the rock during an experiment. During the course of the experimentation the following conditions were maintained constant: (1) pore pressure — atmospheric (i.e., the rock was dry); (2) temperature — 75F; (3) rate of loading — essentially static (approximately 0.002 in./sec); (4) bit tooth — a sharp wedge-shaped tool loaded normal to the rock surface; (5) rock surface — smooth and flat; (6) drilling fluid — hydraulic oil; and (7) maximum depth of penetration — approximateiy 0.1 in. In addition, each experiment was performed on a different rock sample so the rock surface is free of a layer of cuttings and of any previous indentation craters. The influence of the corners of a borehole was neglected, since each rock sample was cemented inta a section of aLuminum tubing to simulate a semi-infinite body. Indentation experiments were made on a sandstone, two limestones, a marble, a dolomite and a schist with a 60' bit-tooth over a range of confining pressures from atmospheric to 5,000 psi and at a constant confining pressure of 1,000 psi for a variety of bit-teeth with angles ranging from 30 to

Jan 1, 1966

By Hsun Hu, B. B. Ruth

The rates of migration of high-angle boundaries in zone-refined aluminum crystals rolled 20 to 70 pct in the (110)[i12/ orientation were studied. Following a recovery anneal at an appropriate temperature to stabilize the polygonized structure, boundary migration rates of artificially nucleated grains were measwed isothermally at several temperatures. Results indicate that the rate of boundary migration depends strongly on the amount of deformation and on the cell size of the polygonized matrix, and is related to the driving free energy by a power function. The degree of anisotropy in growth 0.f the re crystallized grains nn&apos;th preferred mientation is independent of deformation; the migration rates of the fast-moving and the slow-moping boundary segments of a gowing grain differ by as much as one order of magnitude. The actir\ation energy fm a grain boundary migration, although nearly the same for both the fast-moving and the slow-moving boundaries for a given deformalion, decreases from 45 to 30 kcal per mole with an increase in deformation from 20 to 70 pct reduction. Re crstallization by the growth of the artificially nucleated grains results in preferred orientation. The Percentuge of&apos; grains favorably oriented for growth increases with increasing deformation. None of these grains corresponds to the ideal Kronberg-Wilson orientation relationship. The observed growth aniso-tropy is discussed in terms of boundary structure. The boundary velocity as a function of the cell inter -facial area, or the driving free energy, is discussed in the light of current theories of boundary migration. It is well established that recrystallization with re-orientation occurs by the migration of high-angle boundaries of strain-free grains. The driving force for this process is provided by the free energy stored in the metal during deformation. A quantitative study of the effect of varying driving force on grain boundary migration in deformed metals has not been possible heretofore, primarily because of: 1) concurrent recovery steadily decreasing the available driving free energy for boundary migration, &apos;-3 and 2) in-homogeneity of strain in the deformed metal.4 Aust and Rutter3 studied grain boundary migration in striated single crystals of zone-refined lead. Although the driving free energy in such crystals remains unaltered during annealing, this method does not provide a range of driving free energies over which measurements of grain boundary migration can be made. In the present investigation, the rates of migration of high-angle boundaries in deformed aluminum zone- refined single crystals were studied at various temperatures, after deformation ranging from 20 to 70 pct reduction by rolling at -78°C in the (ll0)[i12] orientation. The boundary migration rates along different crystallographic directions were determined under steady-state conditions, i.e., in the absence of competing recovery processes or impingement of recrystallized grains growing into the deformed single crystal matrix. Simultaneous recovery was eliminated by suitable anneals prior to the boundary migration measurements. The recrystallized grains, which grew a ni so tropically into the homogeneously polygonized matrix, developed flat boundary segments during early stages of growth. These boundary segments subsequently migrated along a direction approximately normal to the boundary plane into the matrix rystal. Increasing deformation over the range employed was estimated to increase the driving free energy for boundary migration by about five times. The kinetics of the boundary migration process, examined under these conditions, indicate that the boundary velocity is greatly affected by a small change of the driving free energy in the matrix crystals. These results were examined in the light of the current theories of grain boundary migration. EXPERIMENTAL PROCEDURES Single crystal strips (9 by 1 by 0.125 in.) of zone-refined aluminum, were seed-grown by the Bridgman method in a high-purity graphite mold (<lo ppm ash) at 1 in. per hr. Precautions were taken to minimize contamination of the metal during crystal preparation and subsequent handling. Spectrographic analysis of the metallic impurities in the grown crystals is Qven in Table I. The crystals were rolled in the (110)[112] orientation at -78°C to various reductions in thickness, ranging from 20 to 70 pct, in 10 pct increments. The desired reduction was achieved by many rolling passes, each being no more than 0.002 in. To minimize surface friction, the crystal was rolled between two thin layers of teflon. For those crystals rolled more than 40 pct, it was necessary to remove the disturbed surface layers by electropolishing at -5" to -10°C at an intermediate stage of rolling. The edges of deformed crystals were removed by a jeweler&apos;s saw while submerged in alcohol at -78° C to obtain samples of about ? by i in. The distorted metal at the cut edges and the surface layers were then removed by electropolishing, with removal of a minimum of 0.004 in. from each surface. The thickness of the crystals prior to rolling was chosen so that the final thickness was 0.025 in. for all samples. These deformed single crystals were each prean-nealed for 1 hr at an appropriate temperature in the range of 130" to 280°C, depending upon the amount of deformation. The purpose of this preannealing was to

Jan 1, 1970

By C. E. Ells, A. Sawatzky

The effect of neutron irradiation on tensile properties of the Zr-2.5 wt pet Nb-0.5 wt pet Cu alloy has been evaluated for integrated neutron fluxes up to 3 x 1020 n per sq cm (E > 1 Mev). Specimen temperature during irradiation was normally 300°C. The material was studied both in annealed and in quenched and aged conditions. When the alloy (fab-ricatedfrom sponge zirconium) is quenched from a temperature =40°C below the a + ß/ß transformation and aged 6 hr at 535oC, neutron irradiation up to the maximum studied has little effect 071 reduction in area, although there is a drop in uniform elongation and the yield strength is increased by = 25 pct. Conversely, irradiation of the alloy in an annealed condition can result in much greater changes in properties, with yield-strength increases of up to 200 pet. The dose dependence of irradiation hardening obeys the saturation equntion proposed by Makin and Minter to explain hardening in copper and nickel: Act = C[1 - exp(-Døt)/1/2 A rnodel is proposed to explain the effect of niobium content and metallurgical condition on irradiatiotz behavior. THE Zr-Nb alloys form a class of material susceptible to marked strengthening by quench and age heat treatment. With niobium concentrations as low as 2.5 wt pet, strengths can be developed which are double those of annealed zircaloy-2.1 Since the neutron-capture cross section of the Zr-2.5 wt pet Nb alloy is nearly identical to that of Zircaloy-2, significant gains in power-reactor neutron economy could be obtained by replacing stressed in-reactor components of Zircaloy-2 with the heat-treated Zr-Nb alloy. When the solution heat-treatment temperature is in the high (a + ß) phase, then the tensile properties of the quenched (and aged at 500°C) material have been shown to be relatively insensitive to neutron irradiation.1-3 Addition of small quantities of copper to the binary Zr-2.5 wt pet Nb alloy gives a useful reduction to the corrosion rate in air and carbon di- oxide,&apos; and a ternary alloy of composition Zr-2.5 wt pet Nb-0.5 wt pet Cu was chosen for one component in the Douglas Point Reactor. This application utilizes both the corrosion resistance to a moist air-carbon dioxide environment and the strength developed by a quench and age heat treatment. By suitable choice of solution heat treatment and aging temperatures, tensile properties of the ternary alloy can be made nearly identical to those considered as optimum for the binary alloy. One effect of the copper addition, however, is to give marked changes in the aging kinetics of the alloy.&apos; For this reason, it appeared that neutron irradiation might promote overaging in the ternary alloy, particularly at the elevated service temperature (=300°C). This paper describes the effect of neutron irradiation on the tensile behavior of the Zr-2.5 wt pet Nb-0.5 wt pet Cu alloy in several metallurgical conditions pertinent to its use as a reactor material. The metallurgical conditions included those for which aging occurs in the unirradiated material, and fully annealed conditions for which a higher degree of thermal stability could be expected. 1) EXPERIMENTAL Two batches of alloy were used, both made from sponge zirconium. The fabricator&apos;s analysis of principle constituents showed no significant difference between the two batches, Table I. The rod material was rolled to 1.5 in. diam, and then swaged to final sizes of 0.5 to 0.375 in. diam at nominal temperatures of 785" and 625°C for the AM and AS batches, respectively. Subsequent hydrogen analyses at CRNL confirmed that the hydrogen concentration of the as-received rod was indeed <20 ppm. However, it was found that the (a + ß)/ß transus for the AM material was approximately 25°C lower than that of the AS material, indicating rather more difference in oxygen concentration than given in Table I. Except for a small difference in yield-point behavior, the tensile prop-

Jan 1, 1965

By A. H. Geisler, F. Keller

Although the subject of precipitation from solid solution appears to be one of the more profitable fields in metallurgy for study with the electron Microscope, few comprehensive studies have yet been made. Occasionally electron micrographs have been published that illustrated alleged precipitation in various alloys, but frequently the apparent size, shape and distribution of the particles do not fully agree with those expected from theory or from observations of the Precipitate when the particles have grown to a size resolvable with the light microscope. The presence of a Widmanstatten pattern for submicroscopic precipitate particles has been demonstrated for only a few aged a1loys.l The initial problem has been the development of suitable techniques for applying the transmission-type instrument to a field that has been inherently associated with reflection-type microscopes. The various techniques have been the subjects of numerous publications, and instead of describing them individually here it will sufficc to say that the oxide film method has proved to be the most satisfactory method yet found for studying the microstructure of aluminum alloys with the electron microscope. This method has the definite characteristic advantage that the actual surface layer of the specimen is examined and not a plastic or silica replica. Doubtless suitable methods for forming and removing the oxide film could be developed for alloys of other metals. The purpose of this report is to point out the characteristics of the oxide film method that have been observed during the studies of numerous aged aluminum alloys and to present the results of these studies. Preparation OF Specimens The oxide film method for preparing specimens of aluminum alloys for examination in the electron microscope has been described in detail previously,¹,² Briefly, the procedure consists of preparing the metallographic specimen, forming the film by anodic oxidation and removing the film from the prepared surface of the oxidized specimen, The specimen is polished accordillg to the usual procedure for microscopic examination.³ he specimen may then be etched to remove flowed metal but etching to attack the alloy constituents or to leave them in relief as in the replica processes is not necessary. Electrolytic polishing is not generally recommended, Deep-etched specimens are used frequently, however, since they provide information that is not revealed by polished specimens and frequently present the same information more clearly than do the polished specimens, The surface preparation of specimens to be deep-etched is not important, since specimens that have been subjected to the first wet PoliShing operation are generally used, but frequently the as-rolled surface is Suitable. These specimens are then etched using a suitable macroetching reagent.

Jan 1, 1947

By F. Keller, A. H. Geisler

Although the subject of precipitation from solid solution appears to be one of the more profitable fields in metallurgy for study with the electron Microscope, few comprehensive studies have yet been made. Occasionally electron micrographs have been published that illustrated alleged precipitation in various alloys, but frequently the apparent size, shape and distribution of the particles do not fully agree with those expected from theory or from observations of the Precipitate when the particles have grown to a size resolvable with the light microscope. The presence of a Widmanstatten pattern for submicroscopic precipitate particles has been demonstrated for only a few aged a1loys.l The initial problem has been the development of suitable techniques for applying the transmission-type instrument to a field that has been inherently associated with reflection-type microscopes. The various techniques have been the subjects of numerous publications, and instead of describing them individually here it will sufficc to say that the oxide film method has proved to be the most satisfactory method yet found for studying the microstructure of aluminum alloys with the electron microscope. This method has the definite characteristic advantage that the actual surface layer of the specimen is examined and not a plastic or silica replica. Doubtless suitable methods for forming and removing the oxide film could be developed for alloys of other metals. The purpose of this report is to point out the characteristics of the oxide film method that have been observed during the studies of numerous aged aluminum alloys and to present the results of these studies. Preparation OF Specimens The oxide film method for preparing specimens of aluminum alloys for examination in the electron microscope has been described in detail previously,¹,² Briefly, the procedure consists of preparing the metallographic specimen, forming the film by anodic oxidation and removing the film from the prepared surface of the oxidized specimen, The specimen is polished accordillg to the usual procedure for microscopic examination.³ he specimen may then be etched to remove flowed metal but etching to attack the alloy constituents or to leave them in relief as in the replica processes is not necessary. Electrolytic polishing is not generally recommended, Deep-etched specimens are used frequently, however, since they provide information that is not revealed by polished specimens and frequently present the same information more clearly than do the polished specimens, The surface preparation of specimens to be deep-etched is not important, since specimens that have been subjected to the first wet PoliShing operation are generally used, but frequently the as-rolled surface is Suitable. These specimens are then etched using a suitable macroetching reagent.

Jan 1, 1947

By K. C. Hong, R. B. Jensen

The problem of determining the optimum set of steam volumes and cycle lengths for a single well undergoing multicycle steam stimulation in order to maximize the cumulative discounted net income has been formulated mathematically and programmed for a digital computer. The mathematical fonnulation of the problem and the method for its solution are discussed in this paper. The oil production performance during each stimulation cycle was simulated by either a constant percentage or a harmonic decline. Using simple analytical expressions for production performance, the cumulative discounted net income and cumulative time of operation were related to the pertinent process and cost parameters and two principal process control variables (cycle length and steam injection volume). The discrete maximum principles was used to transform the equations for cumulative time of operation and cumulative discounted net income into a set of simultaneous equations. The simultaneous equations then were solved by trial and error on a digital computer to determine the set of cycle lengths and stem injection volumes that gives maximum cumulative discounted net income over the project life. INTRODUCTION Steam stimulation is a process for improving the oil recovery rate from wells producing high-viscosity crudes. The process is applied on an individual well basis and is executed in a series of cycles, each consisting of three phases: steam injection, soaking (steam condensation), and production. Significant increases in production rate following the stimulation operation result from heating the reservoir around the wellbore. As heat is removed with produced fluids and dissipated into nonproductive formations, the production rate declines, usually to near the prestimulation value. Typical production responses are given in case histories reported by Owens and Suter,l and are depicted in Fig. 1. Duration of the production phase is equal to the time for the oil production rate to decline to some specified value and is called the cycle length. (Cycle length is defined as the producing portion of the cycle and does not include the downtime required for steam injection and soaking operations.) Termination of the production phase coincides with the start of steaming for the next cycle. The process is continued, cycle by cycle, until it becomes unprofitable. Models 2-4 have been developed to simulate behavior of a single well during one cycle of steam stimulation and have been used to investigate the effects of system and operating conditions on the production responses. However, these investigations have not shown how the models can be used to determine a most profitable set of operating conditions for a multicycle stimulation project. Perhaps these models are too complicated mathematically to be adapted readily to an optimization study. This paper presents a method for optimizing a multicycle steam stimulation project based on simple production performance models. First, the performance of a steam-stimulated well during each cycle was simulated by either a constant percentage or a harmonic decline model. Using these simple analytical expressions for production performance, the cumulative discounted net income (before tax) and cumulative time of operation for each cycle were related to the pertinent process and cost parameters and two process control variables—the cycle length and the steam injection volume. Finally, the profit function to be maximized was formed. The analogy between the above problem and problems of optimizing multistage decision processes that have been reported extensively in the chemical engineering literatures5-10 in recent years has led this investigation to the use of the discrete maximum principle.5-7 This principle was used to transform the equations for cumulative discounted net income and cumulative time of operation into a set of simultaneous equations characterizing the optimum operating conditions in terms of the two control variables. These equations

Jan 1, 1970

By Z. F. Oszter

This paper discusses the materials handling problems encountered in cement plants in the areas of raw and finished product grinding and cement storage. It illustrates by examples taken from actual installations, how optimum feed rate is maintained in grinding circuits with properly applied bin dischargers (vibrating hoppers, Mexican feed boxes and aerated bin bottoms) for materials with high moisture content and/or difficult handling characteristics. Initial operating problems with the equipment will be discussed and maintenance requirements will be indicated for comparison purposes. In a modern cement plant, properly functioning bins and feeders are as essential as process control, for no process control system, however sophisticated, can perform efficiently if raw materials do not flow into the plant processing facilities at the required time and rate. Most materials used in the manufacture of portland cement - such as silica sand, iron ore fines, crushed limestone, gypsum, clinker or cement itself-have difficult handling characteristics under specific conditions. High moisture content or even dampness due to condensation may create bin hang-ups, arching and ratholing, all of which result in improper functioning of the feeders or proportioning equipment. The occurrence of these problems in existing cement plants spurred a search for the corrective devices discussed in this article. The equipment investigated included bin dischargers in the form of vibrating hoppers for moist raw materials, Mexican feed boxes for clinker and aerated bin bottoms for cement. Operating experience with different types of bin dischargers showed that the actual performance of some did not come up to the expected level under prevailing operating conditions. The aim of this paper is to record the results of this experience and there- by aid designers and/or operators in selecting the most suitable equipment for their requirements. BIN DISCHARGERS (GENERAL) A discharger is a device that causes stored solids to flow out of a hopper without metering the rate of flow as a feeder does. The investigation discussed here encompasses the following types of bin dischargers used in a cement plant: (1) vibrating hoppers (horizontal motion), (2) whirlpool hoppers (horizontal and vertical motion), (3) aerated bin bottoms (air fluidization), and (4) Mexican feed boxes. Most of the bins used in the cement industry are of the plug flow type. Raw materials with difficult handling characteristics such as high moisture content may cause flow problems when fed through small openings (12 x 18 in.), even if two walls of the bin are vertical.&apos; The most appropriate solution to these flow problems would be the installation of separate mass flow bins for the different materials such as silica sand, iron ore fines, crushed limestone and gypsum, with each bin correctly designed from test data obtained by measuring the flow properties of representative samples. However, limited storage capacity and plant space usually preclude this approach. Quantitative design of bins for flow based on measured flow properties of representative samples of the bulk solid to be handled is no longer an "engineering art," thanks to the basic research done on the flow of bulk solids in bins during the last decade.2, Where space limitations exist, however, or where improvements are being sought for equipment already installed, bin dischargers offer a satisfactory solution to flow problems. They can be applied effectively under flat-bottom bins or mounted to existing plug-flow bins with very little modification to these facilities. In both the foregoing cases, plug flow will be maintained, but the flow channel will be considerably increased so that it is larger than the discharger&apos;s bin-side diameter, and material flow can be sustained if the discharger is correctly selected. The main purpose of a bin discharger is to draw the bulk solids across the whole area of the bin outlet, to make the opening fully effective. As long as the bin-side opening of the discharger is large enough

Jan 1, 1969

By D. Dever, E. S. Fisher

Measurements of the temperature dependence of the elastic moduli in single crystals of hep ruthenium, rhenium, cobalt, dysprosium, and erbium were carried out for various temperature ranges so as to investigate the correlation, suggested by prior work, between the shear anisotropy ratio, A = C44/C66, and the occurrence of the hcp — bee phase transformation. These measurements confirm the suggestion that A has a positive temperature dependence in that group of hcp metals that undergo thermally induced transfornzations to the bee structure and a negative or zero temperature dependence in the other hcp metals. It appears, from indirect evidence, that the hcp -fee transformation in cobalt is preceded by a highly negative temperature dependence of A. These results indicate that elastic anisotropy is an important factor in determining the temperatures at which martensitic phase transformations occur. Two possible roles for elastic anisotropy are proposed: 1) the shear-modulus ratio within the habit planes must be nearly similar for the two phases and 2) a large anisotropy decreases the stress necessary for the movement of the transfornzation dislocations. FUCHS&apos; calculations of the elastic moduli in the alkali metals have shown that the high shear-anisotropy ratios, c/c&apos;, in certain bee metals and alloys are the result of large negative contributions to c&apos; = (c11 - c12)/2 arising from the nearest-neighbor repulsive forces.&apos; On this basis, and following Zener,2 the association between high c/c&apos; ratios and the martensitic phase transformations that occur upon cooling certain bee metals is explained as a unique characteristic of the bee structure. Recent measurements of the elastic moduli in hcp titanium and zirconium show, however, that a high elastic anisotropy also exists in these hcp structures at the temperatures at which they transform to bee upon heating3. If there exists a direct and consistent correlation between the elastic anisotropy in close-packed structures and the occurrence of transformations to the bee phases then there is evidently a need to re-examine the conclusions of Zener theory. Elastic anisotropy may also be the result of the effects of temperature on the stability of the low-temperature close-packed structures or it may play an important role in the mech- anism of the martensitic transformations. The present measurements were aimed at establishing the phenomenological correlation between elastic shear anisotropy in hcp structures and the thermally induced hcp —bee transformations that occur in the majority of the twenty-three hcp metals, as shown in Table I. Prior to the present work the temperature dependence of the anisotropy ratio, A = c44/c66, was known at various temperatures for nine of the metals listed. The data showed a negative temperature dependence of A for cadmium, zinc, magnesium, and cobalt, which do not transform to bee, and a positive temperature dependence of A for thallium, zirconium, titanium, hafnium, and yttrium. The correlation was, however, ambiguous in that these two groups of metals are also separated with respect to the c/a ratios as noted in Table I. There remained the possibility that the temperature dependence of A is a function of c/a and not consistently associated with the transformation. The measurements of the elastic moduli in rhenium and ruthenium, in which hcp —bee transformations do not occur, were needed to help decide this question. The studies of dysprosium and erbium were undertaken to determine whether the correlation holds for the heavy rare earth metals.

Jan 1, 1968

By L. S. Darken

A quadratic formalism is developed lor the representation of the excess free energy, and of the activity coefficients of each component of a ternary system in the vicinity of a single component selected as solvent. This formalism, in contrast to the e formalism, is thermodynamically consistent at .finite concentrations. Numerous specific systems are used as illustrations, with cross checks where possible; these systems are Fe-C-X systems in the vicinity of 1600°C and include practically all those for which reliable data are available. It is shown that this formalism is usually quite adequate to rebresent available data up to a solute concentration of 20 to 30 at. pet. It is hoped that the formalism here presented may prove useful for the treatment and tabulation of data on ternary and multi-component solutions. SINCE, as discussed in a prior paper,&apos; we do not as yet have an adequate basic understanding of the nature of binary metallic solutions, we cannot hope to be in better position for ternary and multicomponent systems. However, the accumulation of thermodynamic data on ternary and multicomponent metallic solutions has led to an urgent need for a rational formalism for the presentation and use of such data. Notable attempts to represent the isothermal isobaric thermodynamics of ternary solutions were made by Benedict et 01.&apos; and by Wohl.3 These attempts were directed toward developing an analytic expression adequate to encompass the entire ternary triangle, primarily for organic systems; they generally involve the power-series approach, which has been shown&apos; to be inappropriate, in general, even for binary metallic systems. On the other hand, interest in the thermodynamics of ternary and multicomponent metallic solutions has tended to focus either on general thermodynamic relations4-10 over wide compositional ranges, or on cases where one or more components have been at rather low concentration. This latter approach has been followed by Wagner,= by chipman,11 and by Alcock and Richardson,12 and more recently in many publications by the same and other authors. It has already been shown&apos; that the thermodynamic behavior of binary metallic systems is, in general, relatively simple only in the terminal regions. Hence, it would appear that any approach, at this time, to a simple formalism for ternary and multicomponent systems must, a fortiori, be limited to the vicinity of the major component or solvent. The present treatment is aimed to be useful for a given solvent with moderate solute concentrations ranging up to 20 or 25 at. pet and in exceptional cases to considerably higher concentrations. In order to achieve even this modest objective, it is necessary that any appropriate formalism have built into it 1) a thermodynamic consistency (at least within the desired degree of approximation), and 2) an approach at infinite dilution to Raoult&apos;s law for the solvent and to Henry&apos;s law for each solute. Let us first formulate the condition of thermodynamic consistency for a ternary system. Under isothermal isobaric conditions, we may write for the molal value G of any extensive function dG - C1 dN1 + G2 dN2 + G3 dN3 [l] where the G&apos;s are the corresponding partial molal quantities and the N&apos;s are the atom fractions. Taking component 1 as the solvent, this relation may be rewritten in terms of the independent variables N2 and N3

Jan 1, 1968

By William Hurst

Muskat&apos;s depletion performance equation is here derived considering the expansion behavior of the reservoir hydrocarbon system and a simple fractional-flow equation. This nietkod of derivation leads logically to the two extensions that follow. The first of these is concerned with gravity segregation in a depletion-drive reservoir. The second is concerned with including an empirically determined tern? for water influx in the performance equation. The more general equation for gravity segregation when there is a primary gas cap and empirically-determined water influx is stated for completeness. These equations have been found useful in reservoir performance calculations in Eastern Venezuela. A discussion on the methods of solving these equations follows, and considers firstly the effect of taking finite intervals in the numerical integration, and sec-ondly, methods of incorporating the time functions involved in segregation in with the expansion behavior. The paper concludes with a brief general discussion on further extensions to the depletion performance equation. INTRODUCTION The two fundamental sources of energy by which oil is produced from a reservoir result from pressure depletion inside the boundaries of the reservoir and fluid encroachment across the boundaries of the reservoir. The wells in either case form low pressure outlets through which oil and gas may be produced by the expansive force of the reservoir fluids and associated encroaching fluids. When the reservoir pressure is higher than the bubble-point pressure of the oil, so that there is no free gas in the reservoir, these expansive forces are the only ones available for the production of oil. However, when the reservoir pressure is less than the bubble-point pressure of the oil, free gas is vaporized as the pressure falls. With both oil and free-gas phases present, the additional forces of gravity and capillarity may operate on the gas-oil system, as they have previously operated on the oil-interstitial water system. Gravity tends to segregate the free gas from the oil due to their density difference. Capillarity opposes and eventually balances gravity as the more extreme free gas and oil saturations are reached, preventing the independent move- ment of free gas until it is above a certain saturation, and the independent movement of oil when it is below a certain saturation. The type of depletion performance equation chosen for predicting the future performance of a reservoir depends on the amount of past history available. When the reservoir is somewhere past the halfway mark in depletion, some form of decline curve is often used. With less past history, material balance equations which incorporate empirical factors based on the past performance are often used. When, however, the amount of past history is small, the Muskat depletion performance equation will usually be used. The distinguishing feature of this type of equation is that empirical factors based on the over-all or macroscopic reservoir behavior are almost or entirely absent. Each parameter affecting the reservoir performance is ascribed an independent set of values based on measurements made on laboratory samples; that is, incorporating microscopic empirical factors. In establishing Muskat-type depletion performance equations, it is necessary to consider the reservoir as consisting of a number of associated blocks, in each of which the saturations and pressures may be considered uniform, and in each of which all substances have uniform pressure-volume characteristics. Thus, a primary gas cap can usually be considered as one block and an aquifer as another. Gravity segregation may be negligible for practical purposes when the rock and oil properties are adverse and/or the dip or thickness of the reservoir is too small. In this case the whole oil leg may be considered as one block, except in very large reservoirs. In very large reservoirs the fluid and rock properties may vary enough, particularly in the dip direction, for it to be necessary to divide the oil leg into a number of blocks, in each of which the relevant quantities may be considered uniform. When gravity segregation of the oil and free gas is not negligible, it is necessary to consider the space occupied by the initial oil leg as divided into two blocks, a secondary gas cap and an effective oil leg, in each of which saturations may be considered to be uniform. The total volume of these two blocks is thus constant, but the secondary gas cap grows continuously at the expense of the oil leg. Muskat1 derived depletion performance equations for the basic case of an oil reservoir with closed boundaries and without segregation, and for the case of an oil

In the past few years the hydraulic aspects of rotary drilling have received considerable attention. It is generally recognized that accurate prediction of circulating pressures is desirable, particularly in those areas where a delicate pressure balance is necessary to prevent both blowouts and lost circulation. For calculation purposes the standard hydraulic equations for the pipeline flow of Newtonian fluids have been altered to facilitate similar calculations for non-Newtonian or plastic fluids.&apos;.&apos; This latter classification includes most colloidal types of drilling muds. Two properties: plastic viscosity, un, and Bingharn yield value, Y,, are commonly used to define the flow characteristics of a plastic fluid. These are commonly obtained with multispeed viscosimeters such as the Fann V-G meter.3 Normally these measurements are obtained at some surface temperature, and are not corrected to circulating temperatures for calculation purposes. This may result in considerable error, particularly in the laminar flow calculations which comnlonly apply to the drill pipe-borehole annulus. The purpose of this study was to investigate the effect of temperature on the flow properties of some water-base muds. A laboratory model Fann V-G meter was used for measurement of viscosities, yield values and gel strengths. An aluminum water jacket with O-ring seals was machined to fit around the mud cup of the Fann meter. All mud samples were prepared by thorough mixing in a Hobart mixer. Viscosity and gel strength measurements were made with the Fann meter using standard procedures. Flow properties of muds were studied at 80, 120, 160 and 180°F. The compositions of the test samples are given in Table 1. The typical effect of temperature on the flow properties of these drilling muds is illustrated by Fig. 1. It was found that plastic viscosity and apparent viscosity decreased with an increase of temperature. However, the curves were not linear, and did not appear to follow any definite trends or patterns. Yield point data showed much more scattering.&apos; Since a usable method for predicting drilling mud flow behavior with temperature was desired, an attempt was made to find a simple relationship between these variables. The basic approach used by Havenaars for predicting the same effect on turbulent viscosity resulted in Fig. 2, which shows that

By R. D. Taylor, Jim Douglas Jr., H. H. Rachford Jr., P. M. Dyke

A major obstacle to the use of wetting agents in .secondary recovery by water flooding is the adsorption of the agents on the sand. As a result of adsorption, the surfactant always lags behind the floodwater front. Consideration of the chromatographic theory of adsorption indicates that the detergents will not lag as much if used in very high concentrations. An investigation was made of the possibility of using high concentrations economically by flowing slugs of wetting agents followed by normal flood water. The experiments consisted of adsorption studies on Alundum powder and Berea sandstone. Flow rests on a 12-in. Alundum core and 22-in. Berea core were used to determine rate of detergent movement. The results of the flow experiments indicate that the relative rate of surfactant advance is, indeed, sensitive to the concentration of the agent. A 10 per cent slug moved with a rate that war 78 to 95 per cent as fast as the rate of advance of the flood water. By contrast, one with 25 ppm (the number of parts of commercial detergent in a million parts of water on a weight basis) concentration moved less than one-fourth as fast as the flood water, and calculations indicate that in very long porous systems the rate of movement of the lower concentrations will be a small fraction of the rate of advance of the flood front. The results. of the adsorption studies were utilized to calculate the rate of advance of the detergent when only the initial concentration was known. The calculated rates showed substantial agreement with the experimental flow tests in the high concentration ranges. The adsorption results were also used to estimate the cost of the materials for a slug-type surfactant flood in the field. In addition to the faster rates of movement, the concentrated detergent slugs removed much more oil than the dilute solutions. However, the effectiveness of the slug process depends on many variables. The quantity of oil removed is increased markedly by increasing the flooding rate. The efficiency is also influenced by the type of crude, type of reservoir rock and initial water saturation. Therefore, a careful analysis of each reservoir system is required before the economic value of the process can be determined. INTRODUCTION It is well known that the displacement of oil by invading water during water flooding is far from complete. It is generally agreed that the unrecovered oil is retained in the porous medium by the capillary forces which may be relatively large compared to the forces generated by the flowing water. Therefore, it was logical that some early workers should turn to surface-active materials to reduce the capillary forces to facilitate the release of oil. As early as 1927,&apos; a patent was granted for the use of surface-active materials in water flooding. In 1932, when soap solutions were passed through Bradford and Venango sands, it was reported that the results were inconclusive, erratic and that "further investigation is needed to determine exactly the function of the solution and to obtain a clearer insight into the phenomena involved."&apos; Some of the modern scientific reports conclude with a similar statement,&apos; showing that the lack of agreement on the mechanism of oil removal by wetting agents is still very widespread even though several comprehensive studies have been reported.&apos;." Although there is a lack of agreement as to the general effectiveness of the detergents for water flooding, most investigators do agree that all of the common detergents are strongly adsorbed onto the solid surfaces of the reservoir. In the early calculations it appeared that all additives would be lost before reaching much of the formation area which contained the additional oil to be removed. Experiments indicated that if the usual small waterflood concentrations of wetting agents were used, the rate of advance of detergent through the formation would be only a small fraction of the rate of advance of the flood front. Indeed, some investigators4 felt that the use of wetting agents would never be economically feasible because of their adsorption. For example, DunningG estimated that the wetting agent in concentrations of 25 ppm, would advance only 0.05 times as fast as the flood front. Ojeda, et al,&apos; found that a surfactant in a concentration of 10 ppm moved less than 0.01 times as fast as the flood front. It is significant, however, that both investigators found that increased concentrations of wetting agents moved faster, relative to the flood front, than solutions at the lower concentrations. Ojeda showed that an extrapolation of his data indicated a relative rate of 0.5 at 1 per cent concentration, while Dunning6 estimated a relative rate of 0.46 for a 1 per cent concentration. It was obvious that these concentrations could not be used for continuous injection because the cost of the injected detergent would far exceed the value of additional oil produced. Traditionally, detergents are used in very low concentrations for they show good

By Charles Wraith

THE property of the Roan Antelope Copper Mines, Ltd., is situated at Luanshya, north-central part of Northern Rhodesia, and is connected with the main line of the Rhodesia Railways, Ltd., by a 24-mile branch line from N&apos;dola. Most of the material and supplies that are not available in Africa are shipped into the country via Beira, on the East Coast of Africa, a distance of 1473 miles from Luanshya and the nearest ocean port. The output of copper is shipped by rail to Beira and thence by boat to the market. Passenger travel is mostly via Cape Town, a distance of 2162 miles from Luanshya. Excellent passenger trains operate on regular schedules between Cape Town and Northern Rhodesia. Luanshya is 13° south of the equator at an elevation of 4000 ft. above sea level, and has a mild healthful climate. During the rainy season, from October to March, the rainfall averages 53 in. The area is pro-fusely covered with vegetation. An ample supply of water for industrial and domestic purposes is available throughout the year. The more important supplies for smelting purposes, such as wood, coal, limestone, silica, quartzite, and clay for "dolly" mud, are con-veniently available. Excellent grades of coal, coke and firebrick are obtained from the Wankie Colliery Co., Ltd., Southern Rhodesia. The Wankie colliery has just completed the erection of equipment for the manufacture of silica brick. An ample supply of quartzite of good grade is available for this purpose. Silica brick is also obtained from the Belgian Congo and from India. Limestone is quarried from a large deposit a few miles from Bwana Mkubwa, which is 30 miles by rail from the smelter. An abundant supply of clay for "dolly" mud and for lining launders is available within a short distance of the smelter.

Jan 1, 1933