Search Documents

By H. D. Hodges, J. K. Godbey

In order to better understand the fracturing process, bottom-hole pressures were measured during a number of typical fracturing operations. A recently developed system was used that allows simultaneous surface recording of both the bottom-hole and wellhead pressures on the same chart. The results from six fracruring treatments are summarized on the basis of the pressure data obtained. Al-though no complete analysis is attempted, the value of accurate pressure measurements is emphasized. Important characteristics of the bottom-hole pressure record do not appear at the wellhead because of the damping effect of the fluid-filled column. In four of the six treatments described, the formations apparently fractured during the initial surge of pressure with only crude oil in the well. The properties of the fluids used during the treatments are given and the fluid friction losses are obtained directly from the pressure records. This technique is also shown to be adequate for determining when various fluids, used during the process, enter the formation. INTRODUCTION Hydraulic fracturing for the purpose of increasing well productivity is now accepted in many areas as a regular completion and workover practice. Numerous articles have appeared in the literature discussing the various techniques and theories of hydraulic fracturing'. In general, three basic types of formation fractures are recognized today. These are the horizontal fracture, the vertical fracture, and fractures along natural planes of weakness in the formation'. Any one or all three of these fracture types may be present in a fracturing operation. However, with only the wellhead pressure record as a guide, it is difficult at best to determine if the formation actually fractured, and is almost impossible to determine the type of fracture induced. These difficulties arise in part because the wellhead pressure record, especially when fracturing through tubing, does not accurately reflect the pressure variations occurring at the formation. Several factors contribute to this effect and preclude the possibility of using the wellhead pressure as a basis for accurately calculating the bottom-hole pressure. These factors are: 1. The compressibilities of the fluids which damp the pressure variations. 2. The changes in the densities of the fluids or apparent densities of the sand-laden fluids. 3. The flowing friction of the various fluids and mixtures, which is dependent on the flow rates and the condition of the tubing, casing, or wellbore. 4. The non-Newtonian characteristics of a sand-oil mixture and its dependence upon the fluid properties, the concentration of sand, and the mesh size used. 5. The unknown and variable temperatures throughout the fluid column. Because of these reasons it was determined that in order to obtain a more accurate knowledge of the nature of fracturing, the bottom-hole pressure must be measured along with the pressure at the surface during a fracturing treatment. Even with accurate pressure data, a reliable estimate of the nature of fracturing is still dependent upon knowledge of the tectonic conditions. However, the hydraulic pressure on the formation is basic to any approach to a complete analysis. In order to accomplish this objective a system was developed to record the wellhead and bottom-hole pressures simultaneously at the surface. By recording both pressures on a dual pen strip-chart recorder, it was possible to greatly expand the time scale so that rapid pressure variations would be faithfully recorded. By such simultaneous recording, time discrepancies inherent in separate records are eliminated, thus overcoming one of the most difficult problems associated with bottom-hole recording systems. This paper illustrates the results obtained by using this system during six typical fracturing operations. All of these tests were taken in wells that were treated through tubing. By a direct comparison of the wellhead and bottom-hole pressures, the importance of obtaining complete pressure information during a fracturing treatment is emphasized. THE INSTRUMENTATION AND PROCEDURES The bottom-hole pressure measuring instrument consisted of a pressure-sensing element, a telemetering section, and a lead-filled weight or sinker bar. The pressure-sensing element used was an isoelastic Amerada pressure-gauge element. By using an isoelastic element, no temperature compensation was necessary in the tests described, since the temperature was believed to be well below the maximum temperature limit of 270°F. The rotary output shaft of this helical Bourdon tube element was coupled to a precision miniature potentiometer. The rotation of the pressure-gauge shaft thus changed the resistance presented by the potentiometer

By R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

By R. P. Agarwala, M. S. Anand

Using residual activity technique, the diffusion of palladium, silver, cadmium, indium, and tin in alunzinum has been studied in the temperature range of 400" to 630°C. The diffusivities (in units of square centimeters per second) have been expressed as: IMPURITY diffusion in aluminum,1-9 silverand lead5 for cases of low solid solubility of the impurity in the host metal has yielded frequency factors in the range of l0-6 to l0-9 sq cm per sec whereas the activation energy is practically half the self-diffusion activation energy value. From the observed values of frequency factor, activation energy, and entropy of activation, it has been suggested' that these solutes are not diffusing by vacancy or interstitial mechanisms but by a mechanism which should be consistent with such low values of the diffusion parameters (Do and Q). However in spite of extensive work on these types of systems, the mechanism of diffusion is still not well understood. The present investigation on the diffusion of palladium, silver, cadmium, indium, and tin in aluminum has been carried out to throw further light on the diffusion mechanism in systems, where the solid solubility is very low (except for the case of silver). The results are discussed on the basis of solid solubility and the structural changes involved owing to the presence of the solutes in aluminum solid solution. An attempt has also been made to apply the existing theories of charge5-8 and size8 difference between the solute and the solvent. EXPERIMENTAL PROCEDURE Specimens (1/2 in. diam by 3/8 in. high) were machined out of pure aluminum (99.995 wt pct) rod obtained from Johnson Mattheys. They were sealed under vacuum in quartz tubes and annealed at 620° C for several hours; the grains thus developed were sufficiently large to eliminate the possibility of diffusion along the grain boundaries. The flat ends were prepared carefully after polishing as described previously,10 Radioactive nitrates of cadmium, indium, and tin and chloride of palladium containing, respectively, cd115, 1n114, sn113, and pd103 were dissolved in distilled water and mixed with 30 pct acetone. By means of a micropipet a drop of this solution was placed on a smoothly polished and lightly etched surface of the specimen. Due care was taken to see that the solution spread uniformly on the surface of specimen without trickling down its sides. Radioactive silver was elec-trodeposited using a AgCN-KCN bath. The amount of sample deposited in all the cases was not more than 0.1 µ thick. The samples were then sealed in quartz tubes in vacuum. The cadmium samples were sealed in a purified argon atmosphere to avoid evaporation. The samples were then diffusion-annealed. The temperature of annealing varied between 400° and 630°C and was controlled to ±5°C. On heating to -400°C,the deposits of cadmium, indium, and tin, which were of the order of 0.1 p in thickness, were converted to their respective oxides. The contribution of oxygen present in the lattice of aluminum due to these oxides has been calculated and found to be less than 10 ppm in all cases. Oxide method has already been used by other workers11'12 in diffusion studies without any controversy on the issue. However, in some of these investigations, metallic deposition was also tried. The diffusivities calculated from these measurements were found to agree very well with the diffusivities obtained by using the oxide method. Thus it is assumed that the measured diffusivities represent true diffusion coefficients. Since palladous chloride decomposes at about 500°C, the deposited samples which were to be diffusion-annealed below 500°C were heated in vacuum for a very short time at 500°C to allow the decomposition of palladous chloride to palladium metal. Time taken in decomposition of nitrates to oxides and chloride to metal was negligibly small as compared to the period of the diffusion anneals. The residual activity technique13 was used to study the diffusion profiles where thin layers from the specimen surface were removed by grinding it parallel to a flat surface on a 600-grade carborundum paper. The specimen was washed, dried, and weighed, the differ -ence of the weight being the measure of the thickness of the layer removed. After each such abrasion and weighing, the total residual activity on the surface of the specimens was measured by counting 0.656, 0.94,

Jan 1, 1968

By K. A. Jackson, G. A. Chadwick, A. Klugert

The diffusion field ahead of a lamellar interfnce is analyzed using an electrical analog. A self-consistent solution is obtained for the shape of the interfnce and the diffusion field by an iterative process. The solutions presented here are for a 50-50 eutectoid or eutectic, The shape of the interface is found to he independent of growth velocity and lamellar spacing, and to depend on the relative values of interfacial free energies at the phase houndaries . The mode of growth of lamellar eutectics and eutectoids has been a subject of much interest for many years.1-4 Mehl and Hagel 1 have shown photomicrographs taken by Tardif when he attempted to determine experimentally the shape of an advancing pearlite interface; the results are completely ambiguous. Brandt' and schei13 have made approximate calculations of the composition ahead of a lamellar growth front. The shape of the advancing front and the composition distribution ahead of the front are difficult to calculate because one depends on the other. It is the purpose of the present paper to describe a method by which this calculation has been done. Lamellar-eutectic growth usually occurs under conditions where the growth is fairly rapid, and the interface temperature is close to the eutectic temperature. The growth rate is usually determined by heat flow. Eutectoid growth, on the other hand, can best be studied by quenching to some temperature, and allowing growth to proceed isother-mally. In both cases the growth is believed to be controlled by diffusion* rather than by the atomic kinetics of the transformation. This being the case, a single treatment of the diffusion equation will apply to both cases, provided the region of the interface in a eutectic may be considered to be isothermal. If a part of the interface could appreciably change its thermodynamic driving force by advancing ahead of or lagging behind the mean interface, then the two cases would not be similar. Eutectics normally grow in temperature gradients of the or-der of a few degrees per centimeter. The normal eutectic spacing is the order of a few microns. Part of the interface would have to extend many lamellar spacings ahead of the mean interface before it experienced sensibly different conditions. The interface temperature is usually a few tenths of a degree below the eutectic temperature so that temperature differences of the order of one-thousandths of a degree (a displacement of one lamellar spacing) would be unimportant. Protrusions large compared to the mean spacing do occur when one phase only grows into a eutectic liquid. This is usually a dendritic type of growth, and easily distinguishable from the lamellar mode of growth. A single treatment of lamellar growth will apply equally well to both eutectic and eutectoid decomposition. At the interface, which as shown above is essentially isothermal, the difference between the equilibrium eutectic temperature Teu and the actual interface temperature Ti, can be divided into two parts: 1) the composition varies across the interface, so that the local equilibrium temperature is not Teu; and 2) the interface is curved, so that the local equilibrium temperature is depressed according to the Gibbs-Thompson relationship. This undercooling can be written as Teu-Ti =?T = mAC(x) + a/r(x) [1] where ?C(X) is the departure of the composition at a point x on the interface from the eutectic composition, see Fig. 1, r(x) is the local radius of curvature at a point x on the interface, m is the slope of the liquidus line on the phase diagram, and a is a constant given by where s is the interfacial free energy, TE is the equilibrium temperature, and L is the latent heat of fusion. The calculations in this paper will be made only for the case where the phase diagram is symmetric, that is, the eutectic occurs at 50 pct, the liquidi have the same magnitude slope m at the eutectic temperature, and C,, the amount of B rejected when unit volume of a freezes, see Fig. 2(a), is the same for both phases. As shown in Fig. 2(b), the composition ahead of the a phase will be rich in B, the composition ahead of the ß phase will be rich in A. The composition at the phase boundary is the eutectic composition. The difference between the local liquidus temperature and the actual tempera-

Jan 1, 1964

By L. D. Clark

Pay-out, the ratio of capital expenditure to annual cash return can be very revealing when employed to gage the merits of a proposed mining venture. It is a quantity, numerically but not physically equal to a factor, by which cash return may be discounted to equal required investment. The fewer the pay-out years, the greater the average annual earning power of the project. Considerable material has been written in the past about the most important and difficult problem mining management faces in evaluation of a property: the optimum production rate for the venture. What daily capacity should be adopted for the proposed operation that it may be the most rewarding to those underwriting it? H.C. Hoover has said, "The most important objective is the least cost per ton mined; and minimum working costs can only be gained by the most intensive production." Theoretically and practically, the quicker the ore is removed, the more profitable should be the venture because of the saving in fixed charges, the consequent lower unit costs and the interest unearned by the profit from the unmined ore. Yet, such performance will be tempered by policy and a restraint imposed by the ratio of capital expenditure to annual cash return. * The 'Annual net profit after taxes but including depreciation and depletion allowances. former will be based upon an analysis of the economic climate anticipated, with its many ramifications including the demand for the product, prices, tax laws, the conduct of labor and related wage scales. The variables and intangibles involved in estimating the optimum production rate for a mine do not readily lend themselves to mathematical relationships or equations and the problem can become complex indeed. The more readily applicable and relatively faster appraisal method will be used which, in addition to being reasonably adequate in most cases, will serve as a guide, should a more expanded and detailed estimate be desired. The first question is what is company policy with respect to the rate of investment return? This will depend on that significant period of time referred to as the pay-out, which can be defined as the estimated interval during which capital investment is returned through application of annual cash return to its recovery. Expressed another way, it is the ratio of capital expenditure to the annual cash return. For example, if the capital investment is $C and the annual cash return is $A (after Federal tax), then $C/$A = the pay-out period in years. That is to say, if $C = $1,500,000 and $A = $214,000, then which is the mathematical expression for the period of time it takes to recover the funds originally invested in a mining enterprise. The relationship yp = $C/$A can be useful — an expression which requires only the significant tern yp to be fixed; in other words, the long range investment policy. For what is the pay-out period but an expression of that policy in terms of the speed, within practical limits, with which an individual or company may want to recover capital expenditure? A policy by which pay-out will be fixed is imperative. This ratio must lie within relatively narrow limits governed by what is economically feasible. It cannot be too low, rendering the ratio impractical or absurd, nor too great, delaying the return of capital for reinvestment beyond a commensurate interval. What pay-out will govern this policy? Federal corporate taxes would seem to present the best measure. Consider mining in Canada where new mines are exempt from Federal income tax for 3½ yr from date of first production. This exemption expressly allows The Act states: "A corporation . . . i s not required to include the profits . . . for the period of 36 months commencing with the day on which the mine came into production in reasonable commercial quantites." a company to recover capital expenditure during the early stages of production. With such an incentive, it is logical to strive for as much capital redemption as possible during this period. Thus, 3½ yr would be a major contribution to the pay-out in Canada. Furthermore, it generally takes from 1½ to 3 yr to develop

Jan 1, 1963

By R. F. Krueger

Results are presented from tests of dynamic fluid-loss rates to cores from clay-gel water-base drilling fluids containing different commercial fluid-loss control agents (CMC, polyacrylate or smt,ch), organic viscosity reducers (quebracho and complex metal lignosulfonate) and oil at several different levels of concentration. In the dynamic system the most effective individual additives to the clay-gel drilling fluid, based on cost-equalized concentrutiom, were found to be starch and the viscosity reducers. These results do not conform with the rankings determined by API fluid-loss rests, which indicate CMC, polyacrylate and starch to be the most effective and comparable. Generally, minimum dynamic fluid-losr rates were attained at cost-equalized concentrations of additive (including thinner) of about $1.00/ bbl, or less. For chernically treated clay-gel drilling fluids, both the standard and the high-pressure API filter-loss tests were found to he inaccurate indicators of trends in dynamic fluid-loss rates under the test conditions used, particulurly for drilling muds containing viscosity reducers. From a practical field viewpoint, restrictions on the applicability of the API fluid-loss test are such that it is open to question whether or not results of this test can be used routinely with confidence as an indicator of control of down-hole fluid loss under field treating conditions. INTRODUCTION The petroleum industry spends large sums of money during drilling operations to control the fluid-loss properties of drilling fluids based on the standard API filter-loss test,' which is a static filtration system. Laboratory studies' ' of dynamic filtration have shown that in a given time period filtrate loss from a circulating mud stream is greater than from a static system and that it is a function of linear mud velocity, pressure and the properties of the drilling fluid. Ferguson and Klotz' and Horner, et al," observed that (I) the dynamic fluid-loss rates for the drilling fluids used were not related to the extrapolated API filter loss and (2) the drilling fluids with the lowest API filter losses did not have the lowest dynamic fluid-loss rates. However, there has been no published information on the relative effects on dynamic fluid-loss rate as a given drilling fluid is treated with increasing amounts of chemical additive to reduce the API filter loss. Such information is economically important because drilling-fluid costs rise rapidly as chemical requirements increase. This paper presents the results of a study of dynamic filtratioi rates to cores from a clay-gel water-base drilling fluid treated with various commercial viscosity reducers and chemical fluid-loss control agents. The dynamic fluid-. loss rates to cores are compared with the standard API filter-loss values at several different levels of additive concentration. Dynamic filtration rates were obtained in each experiment under two different simulated wellbore conditions: (1) filtration just above the bit through a new mud cake laid down dynamically on a freshly drilled formation and (2) filtration up-hole through a mud cake formed by deposition of a static filter cake on top of the initial dynamically formed cake. The latter case corresponds to the bottom-hole conditions existing above the bit when mud circulation is restarted after a stand of pipe has been added or a round trip has been made to change the bit. Except for the short-duration, high-rate filtration beneath the bit where no mud cake can form, these two conditions probably represent the two extremes of dynamic filtration. Because thickness of a dynamic mud cake formed on freshly exposed formation is limited by the shearing action of the mud stream, the filtration rate for this condition is high. On the other hand, once circulation is stopped and a static mud cake forms on top of the dynamic cake, re-starting circulation has only a small effect on the cake properties and filtration rate is much lower thereafter. A discussion of the mechanics of mud-cake deposition and dynamic filtration is outside the scope of this paper but may be found in more detail in publications by prior investigators. APPARATUS AND EXPERIMENTAL CONDITIONS The test equipment used to simulate the dynamic flow conditions existing during drilling was a modification of that described previously by Krueger and Vogel: A schematic flow diagram is shown in Fig. 1. In general, a power-driven, high-pressure mud pump capable of delivering up to 60 gallmin was used to circulate drilling fluid parallel to the faces of 1-in. diameter sandstone cores mounted in a 2 3/4-in. ID high-pressure test cell. Pump rates were controlled by means of a magnetic clutch to maintain an average axial fluid velocity of 110 ft/min in the annular space between the cell wall and a 1 1/2-in. rod positioned on the center line of the cell. The core specimens were Berea sandstone plugs sealed with plastic inside 1 1/8-in. OD tubes and were fluid-saturated prior to use. Burettes were used to accumulate fluid discharged from the cores. The mud sump shown was used for treatment and storage of the drilling-fluid samples during a particular test. The valve arrangement permitted either (1) circulating drilling fluid through the by-pass line while treating with

By W. J. Slatosky

As a means of effecting better control of endpoint temperatirres at the Jones & Laughlin basic oxygen furnace plant, a set of mathematical equations has been developed. The eqlutions are the product of a themlochemical anaysis of the process and aye designed to calculate the required scrap, lime, and hot metal additions in terms of a number of independent variables. Results of test heats have warranted adoption of this technique by the Prodrrction Department. BECAUSE of the autogeneous nature of the basic oxygen steel-making process, bath temperature can be controlled without an external fuel supply by charging the furnace with additions that are thermally balanced. The thermal requirements of the charge materials are such that, during the refining process, they throttle the heat generated by the metallurgical reactions in a manner designed to result in a speci-fied temperature at the completion of the heat. In the past, operating personnel at the basic oxygen furnace plant of Jones & Laughlin's Aliquippa Works relied on their experience and technical knowledge of the process to determine the quantities of charge additions needed to result in a finishing temperature in the range 2880"to 2920" F. (The charge consists primarily of 93 tons of scrap and hot metal plus an amount of lime sufficient to maintain a basicity ratio of 2.8 to 3.2). Estimates of these materials are based on a consideration of the effects on finishing temperature of 1) iron silicon content, having a variation of 0.8 to 1.8 pct; 2) iron temperature, ranging from 2250°to 2600°F; and 3)any excessive cooling of the furnace due to a production delay. The end temperature of the preceding heat also serves as a guide in that, if a heat was within the specified temperature range, the succeeding heat could be charged with materials of nearly the same proportions, provided the hot metal used in each of the two charges was of approximately the same temperature and composition. On the other hand, if a heat was outside the specified tapping range, or if the hot metal used in that heat was of different analysis and temperature from that of the iron to be charged, an adjustment in the proportion of additions is in order for the following heat. Due to the complex thermochemical behavior of the process and to the inexact and subjective nature of the described method of determining charge additions, consistently accurate temperature control was not to be expected. Therefore, those heats out- side the specified tapping range necessitated subsequent adjustments by either reblowing the cold heats for a suitable length of time so as to elevate the bath temperature to the desired level, or cooling hot heats with a proper amount of scrap. Because extra time is required to make these adjustments, production is delayed. In an attempt to devise a method for improving temperature control, an analysis of the thermochemistry of the process was undertaken. This, in turn, led to the development of a set of mathematical equations which enable the calculation of the quantities of scrap, lime, and hot metal needed to result in any specified tapping temperature range. The analysis was not intended to be a repetition of work done by others such as McMulkinl or ~hilbrook.' It was meant to be an extension of their work so that charge additions could be calculated not in terms of silicon alone but, rather, as a function of all independent variables. This paper presents the derivation of these relationships, their effectiveness in controlling bath temperatures, and a method of utilizing them on an operational basis. The Heat Balance—The first step undertaken in the analysis of the problem was the enumeration of the pertinent variables. A list is presented in Table I where it is noticed that these quantities have been separated into the following three categories: important variables, variables considered as constants, and variables to be neglected. The breakdown was an arbitrary one designed to facilitate the analysis; otherwise, the mathematical treatment would have been exceedingly cumbersome and complex. Fortunately, experience has shown that these simplifying assumptions do not seriously impair the accuracy of the calculations. These variables along with the limiting assumptions listed in Table n were then used to write a heat balance of the process by applying the equation of continuity, Rate of Rate of Rate of Increase = Income - Outgo PI ] of Heat of Heat of Heat.

Jan 1, 1962

By Francis R. Conley, John A. Sievert, John N. Dew

A brief review is presented of the past performance of a number of large, thin, highly permeable reservoirs with low dips in the Bolivar Coastal fields of Venezuela. The performance of these reservoirs indicates that the fluids are segregated vertically within the sand section by gravity. With this assumption, equations are developed which describe the performance under pressure maintenance operations. Methods of solving these equations and results of simple example calculations are presented. Example calculations indicate that under pressure maintenance conditions injected gas tends to {tow preferentially along the top of the sands and that encroaching water has a tendency to flow preferentially along the bottom. The expected performance of segregated fluids is discussed and compared with that of fluids which are uniformly distributed in any sand section. INTRODUCTION The field performance of a number of reservoirs described herein indicates that the fluids are segregated vertically within the sand section by the force of gravity. Thus, it was felt that any method used to predict the future performance of these reservoirs should consider the effects of segregation in the sand sections. A study of the available information on some of the reservoirs suggests that increased recoveries may be expected if the reservoir pressure is maintained by either crestal gas injection or flank water injection. To predict future performance of reservoirs under pressure maintenance operations, a method of analysis was needed which would account for segregation of fluids in the sand sections. Several methods of analysis have been developed to take into account the segregation of fluids in the reservoir as a whole1,3 To our knowledge no method considers segregation of the fluids within sand sections in the manner indicated by the past performance of several reservoirs in the Mara-caibo Basin. This paper outlines part of the work done in studying some of these reservoirs and contains a description of their performance characteristics. The analysis presented is restricted to pressure maintenance conditions, since space limitations prevent a full discus- sion of the development of the mathematical relations and the various methods of solving them. Only a sketch of the development of the mathematical relations is given. The various methods of solving these relations are pointed out, but the actual determination of solutions to various problems are omitted except for one example. It is felt that these results may aid in the study of reservoirs outside the Maracaibo Basin. Some concepts on which the present analysis is based are outlined by D. N. Dietz. However, the analysis presented herein includes several factors not considered by Dietz: (1) variations in permeability and in the cross section, (2) various shapes of the cross section, and (3) the production of fluids. Also, the mathematical development presented by Dietz differs considerably from our corresponding analysis. RESERVOIR CHARACTERISTICS The reservoirs under consideration, large and thin, with low dip and high permeability, are large to the extent that they contain from 0.5 to over four billion bbl of oil initially in place. The section thicknesses vary from 100 to 400 ft and the lengths from 10,000 to 40,000 ft. Thus, a length-wise cross section of the producing formation appears to be long and thin. The dip angles vary from 0 to 6 degrees and the average permeabilities vary from 0.5 to over three darcies. The reservoirs contain from 20 to 50 per cent shale inter-laminated with the productive sands. Most shale breaks within the major sands do not correlate from one well to the next, and correlation is often difficult, even between wells drilled from the same location. However, it is often possible to correlate the shale breaks between major sands for some distance. Herein we are concerned with the performance of reservoirs containing oil of gravities over 20" API. Past reservoir performance indicates good pressure communication and, in cases where pressure sinks have developed, large amounts of fluid migration have occurred. Few of the reservoirs have had initial gas caps, but also, few have been found to be highly under-saturated. Due to gravity segregation, secondary gas caps usually form before the pressure has been reduced more than 25 per cent of its initial value. The effect of gravity has not only been apparent on a reservoir-wide basis, but has also caused segregation of the reservoir fluids in productive sections. Proof of this is found from the results of selective well tests, workovers, and electric and various other types of logs. Evidence that

By N. A. Gjostein

The formal descnption of the decay of an isolated scratch can be written in terms of an appropriate Fourier integral. With the application of certain approximations, this description leads to the second-moment analysis developed by King and Mullins. It is shown that the second-moment analysis does not adequately describe the behavior of a decaying, isolated scratch, and that it is difficult to obtain internally consistent values of the surface self-diffusion coefficient, D,,from this analysis. This inconsistency arises from contributions of higher-order moments to the decay process. Two new analytical methods, based on the Fourier integval representation of a scratch profile, have been developed for determining values of D, from the decay kinetics of isolated scratches. In the profile-matching method, a given initial profile is resolved into its Fourier components, the components are allowed to decay exponentially over an annealing time interval, At, and finally they are synthesized to form a decayed profile. The predicted profile is adjusted to fit the experimental profile by varying only the parameter, D,At; the contribution from volume diffusion is directly taken into account and the integrals involved in the Fourier analysis are evaluated by numerical integration on an electronic computer. In the kinetic method, the widths, W, of a series of predicted profiles are determined for various volume-diffusion contributions, and the results are plotted in the form W vs DSAt. For a fixed ratio of surface to volume diffusivities, R, these plots are linear; the slope of the lines increases as R decreases. A simi-lar analysis can be performed on the depths of the predicted profiles. In this case, plots of (l/df3 vs D,At are found to be nonlinear. The proper value of Ds is found by superimposing the data for the experimental widths and depths on Plots of this kind and adjusting Ds until a fit is obtained with the line having the appropriate value oj R. An analysis of the slopes of these plots for R —= (surface diffusion only) reveals that they do not give the correct value, 2276, for width tf second-moment standard profile. This finding allows one to modify the second-moment analysis in such a way as to obtain internally consistent values of Ds. ThERE are several processes, involving a flow of matter along a surface, due to gradients in curvature of the surface, that may be used to determine the surface self-diffusion coefficient, Ds. These processes are commonly known as 1) thermal grooving, 2) single and 3) multiple scratch decay, 4) blunting of field-emitter tips, and 5) sintering. In this series of papers, the first three of these processes are examined in detail in an attempt to subject the phenomenological theories governing these processes to a critical test and, furthermore, to make a critical comparison of the values of D, obtained from these techniques. Although the thermal grooving technique has been established'' as a reliable method for determining D, , this has not been the case for the process of scratch flattening. Mullins and coworkers374 have proposed phenomenological laws governing the growth of the width and the decay of the amplitude of an isolated scratch, but only in one study5 has the width-growth law been supported by experimental evidence. In paper I of this series detailed experimental evidence is presented that points up the deficiencies of the second-moment analysis; the origins of these deficiencies are also discussed. Following this development, some new analytical procedures, based on representing the scratch profile by Fourier integrals, are presented as a means of obtaining reliable values of D, from the kinetics of the scratch-flattening process. Further applications of the Fourier integral approach are reported in papers 11' and 111,~ where the method has been useful in interpreting the anisotropy of flattening rates of isolated scratches (paper 11) and in describing the multiple scratch decay process (paper 111). THEORY OF SCRATCH DECAY Following ~ullins" approach we consider an arbitrarily shaped scratch profile, Fig. 6, at some time, t. The profile is designated by a function z(x,t) which can be expressed in terms of a Fourier integral as follows:

Jan 1, 1967

By K. A. Natarajan, I. Iwasaki

Platinum electrodes are not inert as often thought to be. The reactivity of platinum electrodes can explain their erratic behavior in many electrochemical measurements of metallurgical interest, e.g, in flotation systems, streaming potential measurements, contact-angle measurements, and in leaching systems. The anomalous behavior of platinum electrodes in redox potential measurements in aqueous systems was studied through Eh and pH measure ments in water-oxygen, iron-water-oxygen, and manganese-water-oxygen systems. Stability relations between Fe++ and Fe (OH), and between Fe (OH), and Fe (OH), were studied to judge the correspondence between experimental and theoretical equilibrium lines. The practicality of redox potential measurements in estimating ferric-ferrous ratios in aqueous systems was investigated along with their suitability as indicators in leaching operations, e.g., the removal of iron by aeration from manganese leach solutions. Platinum electrodes have often been used in the measurement of dissolved oxygen concentrations and of redox potentials (Eh) in a variety of fields, e.g., analytical chemistry,' corrosion," geology and mineralogy,,'' biology,"' sewage treatment,' * hydrometallurgy,"I" and flotation."la The effectiveness of Eh-pH diagrams, first reported by Pourbaix' in 1949, has contributed much towards the theoretical understanding of numerous problems encountered in the metallurgical industry. Not many references are available in the literature, however, wherein attempts have been made to confirm Eh-pH diagrams from experimental measurements. One reason might be that, in spite of the apparent simplicity of the electrochemical technique, the direct measurement of Eh involves complex practical problems.' Factors such as the purity of the solution, the type of electrodes used, the history of the indicator electrode, and the type of atmosphere (namely, oxidizing, reducing, or inert) do have effects on the measured Eh values. The influence of mixed potentials cannot be underestimated. The poisoning of platinum electrodes by organic and inorganic impurities present in the solution may lead to erratic results. Platinum, commonly thought to be an inert electrode material, is not really so, as attested by a number of previous investigators who advised caution concerning the anomalous behavior of platinum electrodes in various electrochemical measurements.'" In the present article, a few pertinent experiments related to Eh-pH measurements in systems of interest in the metallurgical and water pollution fields are described in an attempt to correlate such information with what is already known, especially in the electrochemical literature. Iron-water and manganese-water systems were selected with a view of studying the correspondence between experimentally observed and theoretically established equilibrium lines. The work included an investigation of the behavior of platinum electrodes with respect to pretreatment and adsorption characteristics, the measurement of dissolved oxygen concentrations and their relation to Eh, the determination of the electrode potential of the ferric-ferrous couple at different pH, and the measurement of oxidation potentials in iron-manganese leaching systems. Experimental Procedure A rotating platinum electrode was used in many of the measurements to study the effect of rotation on measured Eh values. The electrode made by the Pine Instrument Co. consisted of a stainless steel rod with a platinum disk soldered to the end. It was covered with a Teflon insulation along the sides, so that only the circular tip of the electrode was exposed to the solution. Prior to its use, the platinum surface was brightened on a metallurgical polishing wheel with alumina as an abrasive, unless specified otherwise. The electrode was rotated with a Sargent cone-drive motor at 350 rpm. The contact of the electrode with the external circuit was made by filling a notch at the top of the stainless steel shaft with mercury and by dipping a copper wire into the mercury pool. The performance of the rotating platinum electrode was compared with the performances of a Beckman platinum button electrode and a platinum wire electrode. All the potentials were measured with respect to a saturated calomel electrode. A saturated KC1-agar bridge was used to minimize the liquid junction potential. A Beckman Zeromatic pH meter together with a Beck-man electrode switch was used to measure both the Eh and pH. A double-walled, all-Pyrex jar with a capacity of about 1 liter and themostated by circulating water of constant temperature was used for a reaction cell. Four equally spaced ports in the cover provided access for a glass electrode, a salt bridge connecting the saturated calomel electrode, a dispersion tube for bubbling gases into the cell solution, inlet and outlet tubes for passing the desired gases over the solution, and a graduated burette for introducing reagents from outside. The rotating platinum electrode was inserted through an opening in the top-center of the cover, and a positive gas pressure was maintained inside the cell to prevent air from entering into the cell compartment. A magnetic stirrer was used to mix the solution inside. For the determination of dissolved oxygen in the test solutions, the polarographic techniquex was used.

Jan 1, 1971

By R. Ebeling, M. F. Ashby

The paper is in two parts. The first develops the formulae and method needed to calculate the size, nu)nber, spacing, and volume fraction of hard or inert particles in the interior of a specimen from measurernents made on an extraction replica; throughout, the spread of- particle sizes in a single specimen is taken into account, and the extraction efficiency of the replica is considered. A method for testing whether the particles are randotnly distributed in space is described, and the spacing of randomly distributed particles is discussed. The second Part describes the application of the method to copper containing spherical particles of silica; the description includes an estimate of the errors involved, and a comparison of this with alternative methods. All parameters except uolume fraction can be determined satisfactorily, and the method offers advantages over other electron-,rzicroscopic techniques when dealing with hard particles in a softer matrix, or inert particles in a chem-ically less inert matrix. 1 HE experimental measurement of sizes and numbers of discrete particles of a second phase embedded in a matrix material from opaque, plane sections has been considered by many authors, notably by Fullman.' It is assumed that, when the specimen is sectioned, the particles are also sectioned, so that the diameters of the circles of intersection of particles with the plane of the cross section ("surface diameters") are measured. It is difficult to derive the distribution of true particle diameters ("volume diameters") from measurements of the distribution of "surface diameters" on a plane section. The measurement of size and numbers of small second-phase particles by transmission electron microscopy through thin foils of matrix containing the particles has been described by Cahn and ~uttin~' and by ~illiard.~ Because of overlap and the fact that particles which intersect the foil surface may be sectioned whereas those in the bulk of the foil are not, it is again difficult to measure the true distribution of particle sizes. The use of replicas for this sort of measurement has not been analyzed in detail. A replica made from an opaque cross section, thus replicating the circles of intersection of particles with the plane of section, is analyzed by Fullman's method and suffers from the same difficulties. Under certain circumstances, however, extraction replicas can be used, and offer a special advantage in that they make the determination of the true distribution of particle diameters particularly simple. It is the purpose of this paper to describe this extraction replica method. Consider the example of an alloy consisting of SiOz particles, whose mean diameter is about 1000A, dispersed in a copper matrix. When the alloy is sectioned and mechanically polished, the particles are not sectioned, because they are much harder than the matrix, and so will be either left in the surface or pulled out whole, and because the particle size is smaller than the abrasive size of even the finest abrasive. The particles are chemically inert in solutions which chemically etch or electropolish the matrix, so that chemical treatments do not section the particles either. (Experimental measurements of particle shape in this alloy confirm this—truncated particles are never seen.) Small, hard or inert, particles are therefore not sectioned when the specimen is sectioned and polished. An extraction replica taken from this polished surface shows the true "volume" diameter of each extracted particle; the distribution of diameters seen on the replica is the surface distribution of true "volume" diameters, and the volume distribution of true "volume" diameters is easily derived from it. The extraction replica method is not subject to the difficulties associated with sectioning of particles, and avoids the problem of overlap from which the transmission method suffers. The method is well-suited for studying hard particles in a soft matrix. Examples are provided by internally oxidized alloys, oxide-dispersion strengthened metals like T.D. Nickel, and carbides, nitrides, and strong intermetallic compounds in metals. It should be possible to adapt polishing techniques in alloys containing softer but inert particles to avoid sectioning the particles (say by electropolishing or etching). The method is tested and illustrated in Sections 8 and 9 by applying it to a dispersion of spherical SiO2 particles in a copper matrix. I) SYMBOLS Symbols are defined where they first appear in the text. Volume fractions of particles are denoted by the symbol f, particle diameters by x, center-to-center particle spacings by D, and numbers of particles by N. Mean particle diameters and standard deviations about these means are denoted by x and a. Both arithmetic (subscript A) and geometric (subscript G) means are employed. In this paper we measure the distribution of true diameters of particles which touch or pass through a random-plane section of the specimen. We may want to know the distribution of particle diameters in the

Jan 1, 1967

By R. Al-Hussainy, P. B. Crawford, H. J. Ramey

The effect of variations of pressure-dependent viscosity and gas law deviation factor on the flow of real gases through porous media has been considered. A rigorous gas flow equation was developed which is a second order, non-linear partial differential equation with variable coeficients. This equation was reduced by a change of variable to a form similar to the diffusivity equation, but with potential-dependent diffusivity. The change of variable can be used as a new pseudo-pressure for gas flow which replaces pressure or pressure-squared as currently applied to gas flow. Substitution of the real gas pseudo-pressure has a nrtmber of important consequences. First, second degree pressure gradient terms which have commonly been neglected under the assumption that the pressure gradient is small everywhere in the flow system, are rigorously handled. Omission of second degree terms leads to verious errors in estimated pressure distributions for tight formations. Second, flow equations in terms of the real gas pseudo-pressure do not contain viscosity or gas law deviation factors, and thus avoid the need for selection of an average pressure to evaluate physical properties. Third, the real gas pseudo-pressure can be determined numerically in term of pseudo-reduced pressures and temperatures from existing physical property correlations to provide generally useful information. The real gas pseudo-pressure was determined by numerical integration and is presented in both tabular and graphical form in this paper. Finally, production of real gas can be correlated in terms of the real gas pseudo-pressure and shown to be similar to liquid flow as described by diffusivity equation solutions. Applications of the real gas pseudo-pressnre to radial flow systems under transient, steady-state or approximate pseudo-steady-state injection or production have been considered. Superposition of the linearized real gas flow solutions to generate variable rate performance was investigated and found satisfactory. This provides justification for pressure build-up testing. It is believed that the concept of the real gas pseudo-pressure will lead to improved interpretation of results of current gas well testing procedures, both steady and unsteady-state in nature, and improved forecasting of gas production. INTRODUCTION In recent years a considerable effort has been directed to the theory of isothermal flow of gases through porous media. The present state of knowledge is far from being fully developed. The difficulty lies in the non-linearity of partial differential equations which describe both real and ideal gas flow. Solutions which are available are approximate analytical solutions, graphical solutions, analogue solutions and numerical solutions. The earliest attempt to solve this problem involved the method of successions of steady states proposed by Muskat.' Approximate analytical solutions' were obtained by linearizing the flow equation for ideal gas to yield a diffusivity-type equation. Such solutions, though widely used and easy to apply to engineering problems, are of limited value bemuse of idealized assumptions and restrictions imposed upon the flow equation. The validity of linearized equations and the conditions under which their solutions apply have not been fully discussed in the literature. Approximate solutions are those of Heatherington et al.. MacRobertsl and Janicek and Katz.' A graphical solution of the linearized equation was given by Cornell and Katz. Also, by using the mean value of the time derivative in the flow equation, Rowan and Clegg' gave several simple approximate solutions. All the solutions were obtained assuming small pressure gradients and constant gas properties. Variation of gas properties with pressure has been neglected because of analytic difficulties. even in approximate analytic solutions. Green and Wilts8 used an electrical network for simulating one-dimensional flow of an ideal gas. Numerical methods using finite difference equations and digital computing techniques have been used extensively for solving both ideal and real gas equations. Aronofsky and Jenkins"I " and Bruce et al.11 gave numerical solutions for linear and radial gas flow. Douglas et al." gave a solution for a square drainage area. Aronofsky 13 included the effect of slippage on ideal gas flow. The most important contribution to the theory of flow of ideal gases through porous media was the conclusion reached by Aronofsky and Jenkins" that solutions for the liquid flow case'" could be used to generate approximate solutions for constant rate production of ideal gases. An equation describing the flow of real gases has been solved for special cases by a number of investigators using numerical methods. Aronofsky and Ferris 10 onsidered linear flow, while Aronofsky and Porter 17 considered radial gas flow. Gas properties were permitted to vary as linear functions of pressure. Recently, CarteP 18 proposed an empirical correlation by which gas well behavior can be estimated from solutions of the diffusivity equation using instantaneour values of pressure-dependent gas

Jan 1, 1967

Jan 1, 1911

By I. R. M. Chatson

Radial layout's guiding principle is the grouping under one roof of the operating processes which require continuous supervision. Separate treatment sections are isolated by outside stockpiles fed with radiating conveyors. Hence the name "radial" layout. Many plants are designed with "linear" layout where successive operating stages are dotted like beads along a string of conveyors.

Jan 11, 1978

By Gary B. Williamson, Barton K. Cross

Digital simulation is a modern analysis tool which allows the mine operator to evaluate the result of making a decision without actually changing the organization or buying the equipment. The purpose of this paper is to demonstrate the use of digital simulation in evaluating a decision of whether or not to introduce dispatching in an open pit mine. The analysis shows that appreciable savings can be realized by dispatching. In certain cases, an operator can reduce the number of trucks required to move a specified tonnage. A second analysis evaluates the alternatives of using one large shovel or two smaller units. In this case, the operating cost of the haulage equipment necessary to support the shovels turns out to be the critical factor in this decision.

Jan 1, 1969

By H. P. Allen

THE Achotla mill of the Cia. Minera de Peñoles is of special interest in that it is one of the few still employing the chloridization process. The ore milled comes from the Suriana mine and consists of both oxide and sulphide, from neither of which can the gold or silver be extracted by a simple cyanide treatment. The mill was built in 1927, after several years of study in the development of the process1, and operated from December of that year until October 1930. Operations were resumed in February 1935, and have been continuous since, treating approximately 5000 tons per month, while the mine is getting up to full tonnage of 6500 tons per month.

Jan 1, 1937

By T. H. Janes

NATURAL abrasives of some type are found in all countries of the world. In order of their hardness the principal natural abrasives are diamond, corundum, emery, and garnet, which are termed high grade, and the various forms of silica, including pumice, pumicite, ground feldspar, china clay and, most important, sandstone. The properties qualifying materials for use as abrasives are hardness, toughness, grain shape and size, character of fracture, and purity or uniformity. For manufacture of bonded grain abrasives such as grinding wheels, the stability of the abrasive and its bonding characteristics are also important. No single property is paramount for all uses. Extreme hardness and toughness are needed for some applications, as in diamonds for drill bits, while for other purposes the capacity of the abrasive to break down slowly under use and to develop fresh cutting edges is of greatest importance, as with garnet for sandpaper. In dentifrices, soaps, and metal polishes, of course, hardness and toughness are objectionable.

Jan 10, 1954

By S. A. Shoop, J. J. K. Daemen

Site-specific blast vibration prediction can be more accurate than using a universal prediction equation. Vibrations from 37 blasts were monitored using three seismographs. Data analysis showed that square root, cube root, and site-specific scaling predicted ground motion with equal reliability. Imprecise prediction is caused by geologic and blast variability, particularly, inaccurate delay times. Vibrations at some mines were consistently lower or higher than vibrations predicted from US Bureau of Mines (USBM) equations. Vibration restrictions based on USBM studies could then be too conservative or too lenient. Predominant vibration frequencies were generally less than 50 Hz, but several blasts had significant energy up to 125 Hz. In some cases, energy was concentrated at frequencies corresponding to delay intervals.

Jan 1, 1985

By J. V. Emmons

THE processes of hardening and tempering steel tools within the past 15 or 20 years have been so developed that the forward strides of the industry can scarcely be followed by the average observer. No small part of the credit for this phenomenal advancement is due to the improved methods of measuring the high temperatures necessary. Pyrometry has made possible the discovery and use of the many special alloy steels the wonderful properties of which are revolutionizing the industrial world. So sensitive are these steels to changes of temperature that often a variation of only 10° F. at a critical point of the heat treatment will cause the failure of the tool. In the tool-manufacturing industry today, accurate control of the many heating operations is imperative. The principal processes in the manufacture of tools that require accurate high-temperature measurements are annealing, carburizing, hardening, tempering, and laboratory research. Several other operations, such as forging and straightening, use rougher and correspondingly less accurate methods. These operations are usually performed in the following approximate ranges of temperature: DEGREES F. DEGREES F. Annealing 1000 to 1700 Tempering 200 to 1200 Carburizing 1400 to 1000 Forging 1400 to 1800 Hardening 1400 to 2500 Straightening 300 to 1300 Laboratory research covers the entire range. Annealing and carburizing are usually done in muffle or semi-muffle furnaces, heated with coal, coke, oil, gas, or electricity. Many of them are of large size, holding several tons of steel at one charge. Both indicating and recording pyrometers are used. The recording instruments are generally of the thermoelectric type, while for the indicating, both optical and thermoelectric are used.

Jan 9, 1919

By A. C. Emery

The physical and chemical properties for the metalliferrous feed materials to the Imperial Smelting Process blast furnace are defined and the means of achieving these properties by the sintering or briquetting routes discussed. Modern sinter plant performance data is included with a plant flow sheet showing draughting patterns. The necessity for ensuring an adequate SO2 content in sinter machine gases when roasting sulphide materials to enable conversion to sulphuric acid in a contact acid plant is described. The ability of the process to consume limited quantities of oxide materials, residues from other processes and other secondary sources is discussed. The argument for a briquetting route is developed out of this limitation and the merits of cold and hot briquetting examined. Development work on the hot briquetting route for this purpose is briefly described and several flow sheets outlining possible applications for hot briquetting to different raw material situations including the preparation of low lead charges for the furnace are presented.

Jan 1, 1977