Search Documents

By Raynor Linzey, Karl M. Busen

The formation of thin films of silicon oxide by anodic oxidation of silicon and the subsequent removal of these films by an etch is a process often used for the evaluation of concentration distributions Profiles) in silicon layers by the differential sheet conductance method. The accuracy of the resulting profile is very strongly influenced by the uniformity of the thickness of the reacted silicon. Normally, it would be expected that for a constant number of coulombs passed the thickness would be the same from oxidation to oxidation. Investigations show that in certain electrolytes, for a given number of coulombs passed through an n-type silicon sample, the thickness of the reacted silicon increased with increasing resistivity. Even for the same resistivity the thickness varied sometimes by a factor of 1.5. When an electrolyte was used which consisted of 10 pct water by volume in ethylene glycol with 4.0 g KNO per 1000 ml of solution, anodiza-tion at 5 ma per sq cm led to satisfying results. Short-time anodizations gave oxide layers of a higher apparent density than those experienced from thermally gown silicon oxides. THE functioning and the electrical characteristics of semiconductor devices are based upon the incorporation of "impurities" into a single-crystalline body of suitable material and on the concentration distribution (profile) of this impurity. The incorporation can be achieved by well-known processes as, for example, by diffusion, epitaxial growth, alloying, or ion implantation. Often the profile resulting from these processes is not known. A powerful tool to learn about a concentration distribution is given by the method of differential sheet conductance which employs successive four-probe measurements on a layer subjected to stepwise removal of thin sublayers. Differential sheet conductance vs position (or penetration depth) of a sublayer is plotted and a smooth curve is drawn through the data. From this curve the profile is then calculated. When the semiconductor material is silicon, the sublayers are removed suitably by anodic oxidation of the silicon and subsequent dissolving of the formed silicon oxide by an etch. The accuracy of the resulting profile is very strongly influenced by the uniformity of the sublayer thickness. Normally, it would be expected that, for a constant number of coulombs passed, the thickness would be constant from oxidation to oxidation. Investigations showed that for sublayers several hundred angstroms thick the reproducibility can be rather poor. Therefore, efforts were made to obtain a reliable technique for uniform removal. The present paper describes such a technique and certain phenomena which were encountered during the investigations. APPARATUS AND TECHNIQUE The first report on the investigation of concentration distributions, where for differential sheet conductance measurements thin silicon sublayers were removed from a diffused layer by anodic oxidation, was given by Tannenbaum. The author reports the removal of sublayers which for the most part were 400 thick. More advanced device designs now require much narrower layers. When it was tried in these laboratories to determine profiles within such layers, difficulties were encountered with respect to sublayers which by necessity had to be thinner than the ones reported by Tannenbaum. The apparatus used for the investigations is sketched in Fig. 1. Two cylindrical containers connected by a wide tube are filled with an electrolyte. The left container receives the cathode whereas the right container is closed at the bottom end by a sample support consisting of a Teflon base and a copper pedestal. The silicon sample (1 by 1 cm) is mounted to the pedestal embedded flush in the Teflon base using silver paint (Degussa) for electrical contact. Pyseal is applied to the edges of the sample to protect the copper pedestal from the electrolyte. The base is mounted to the container using a water-tight silicon rubber gasket. Fig. 2 gives a view of the sample support. The copper pedestal is connected to the positive side of a power supply operating at a constant current output (Kepco Model ABC 425M). The electrolyte which has been reported by Duffek et al. consisted of either 2 or 10 pct water by volume in ethylene glycol and 4.0 g KNO3 in 1000-ml solution. The electrolyte is best prepared

Jan 1, 1967

By O. D. Gonzalez, H. A. Wriedt

Alloys ranging from pure iron to pure nickel were saturated with nitrogen gas at 918°, 999°, and 1217°C and analyzed. The solubility of nitrogen at 1-atm pressure was obtained as a function of nickel content for the whole range of iron-nickel alloys and the temperature dependence of solubility is shown for the experimental alloy compositions. The heat of solution and the entropy of solution of nitrogen in the alloys as func-tiom of composition up to 41 pct Ni are given. 1 HE acquisition of more exact thermodynamic information for systems composed of iron with various alloying elements continues to interest ferrous metallurgists. The interstitial alloying elements, especially carbon, hydrogen, and nitrogen, comprise a field of particularly intense study. A recent investigation of carbon solubility in Fe-Ni alloys held at 1000°C in gases of fixed carbon activity showed an unpredicted minimum in the solubility near 72 pct Ni.1 The present investigation was undertaken to see if a similar anomaly in solubility existed for nitrogen in Fe-Ni alloys. The solubility of nitrogen in ? iron at known nitrogen activities has been measured several times.&apos;-&apos; Corresponding information for the solubility of nitrogen in solid Ni and Fe-Ni alloys is not available. Juza and sachsze9 reported that metallic nickel dissolved 0.07 pct N in equilibrium with Ni3N at 445°C. This result was regarded as dubious by Turkdogan and Ignatowicz10 who found less than 0.0004 pct N in nickel specimens equilibrated at 600°C with ammonia-hydrogen mixtures having fugacities of N2 gas up to 35,000 atm. Such fragmentary reports indicate only that nitrogen gas has a very much lower solubility in nickel than in iron. MATERIALS A purified iron from Battelle Memorial Institute was used for solubility measurements. This grade contained more than 99.95 pct Fe, the highest content of any measured individual impurity being 50 ppm for each of As, C, Cu, and Pb. The nickel and iron-nickel alloys, supplied by the International Nickel Co., contained 0.02 to 0.04 pct C, 0.17 to 0.19 pct Mn, 0.04 to 0.08 pct Si, < 0.01 pct P, <0.03 pct S, and <0.20 pct Co as principal impurities. The Fe-0.59 pct Mn alloy was of zone-melted quality, the highest concentration of any individual impurity being 15 ppm Al. The "prepurified" grade of nitrogen gas from Matheson Co. had a nominal nitrogen content of 99.996 pct. An "electrolytic" grade of hydrogen gas from National Cylinder Gas. CO. in which the principal impurity is water vapor was admixed with the nitrogen to prevent oxidation of these specimens. This grade was preferred to a purer grade because the presence of some oxygen helps in suppressing contamination of the specimens, for instance, by silicon from the porcelain furnace tube. EXPERIMENTAL PROCEDURE The specimens were equilibrated in vertical tube furnaces through which a 99 pct N2-1 pct H2 gas was passed at ambient pressure. For the 915", 918", and 999°C runs, a wound-resistance furnace was used, and for the 1217°C runs, a furnace with a tubular silicon carbide element was used. The temperature control system on both furnaces was similar: a Weston "Celectray" controller with a modulator, responding to the electromotive force of a thermocouple with its hot junction in the furnace, switched the furnace current between high and low settings. The temperature distributions in the working zones of the furnaces were determined with calibrated Pt/Pt-10 pct Rh thermocouples. The overall uncertainty in temperature was estimated to be ±3° at 915° to 999°C and ±5° at 1217°C . Consideration of the small temperature coefficients of the equilibria studied, together with the error limits in the method of analyzing for nitrogen, showed that further refinement in temperature uniformity would give no improvement in resultant accuracy. The reactive gas was generated by mixing together "prepurified" nitrogen and "electrolytic" hydrogen directly from cylinders commercially supplied. The flow rates of nitrogen and of hydrogen were fixed by setting with bleeders the pressure drops across the individual calibrated flow meters of constant resistance. For any given alloy, the solubility of nitrogen at 1-atm pressure (% N) was obtained from the analyzed nitrogen content of equilibrated specimens (% Na)according to Sieverts&apos; law, thus

Jan 1, 1962

By W. J. Leas, L. A. Rapoport

As a preliminary, consicleration is given to the conventional definition of relative permeability and to the conditions governing the simultaneous flow of oil and gas through porous media. For the conditions of flow prevailing throughout most of a gas drive reservoir, the oil and gas can reasonably be supposed to be in capillary equilibrium with each other. Under these conditions, and these conditions only, the relative permeability to liquid can be expressed as a function of saturation. The relative permeability to liquid in that case is dependent upon the distribution of fluids which itself is shown to be related to the capillary pressure, and, in turn, to the saturation. As a consequence, relative permeability to liquid can be expressed in terms of the volume and surface area of a network of liquid channels bounded by the rock and the gas phase. While the volume of this network can be evaluated accurately, the surface area cannot. However, for any such volume, maximum and minimum values of the corresponding surface area can be calculated from capillary pressure data. It is then possible to establish for any saturation the limits within which the value of the relative permeability to liquid must lie. As a consequence of the theoretical development, the validity of an experimental method for measuring relative permeability to liquid which utilizes a stationary gas phase is demonstrated. In this method capillary barriers are cemented to the ends of the core sample to permit the maintenance of capillary equilibrium between the two phases. At the same time, this procedure eliminates undesirable secondary phenomena such as end effects, fissure effects, etc., the presence of which adversely affect the results of other laboratory methods. The results obtained by theoretical calculations, and experimentally, are discussed. In view of the overall precision that can presently be obtained in reservoir calculations, the agreement between the calculated and measured relative permeability to liquid data can be considered satisfactory. In conclusion, for reasons of economy and simplicity, the procedure of calculating limiting relative permeability to liquid curves from capillary pressure data is indicated for general engineering purposes. It is shown that the above procedure can easily be extended to the cases where connate water is present. Its use for reservoir studies is particularly recommended in conjunction with the method for measuring relative permeability to gas&apos; which simultaneously yields the capillary pressure data necessary for the calculations. THEORETICAL Definition of Relative Permeabilities — Basic Equations for Heterogeneous Flow The equations by which the relative permeability concept is defined and upon which the formulation of all of the gas-oil flow problems rests at the present time are expressed as: V, = — Grad PL = — Grad PL ....(la) PL µL, k k KG VG = —K - Grad PG k- Grad Pc .... (lb) Mo where ,. and G refer to liquid and gas; V is the volumetric rate of flow per unit gross area. µ the viscosity, Grad P the potential gradient. and k the specific permeability of the porous medium.* (For horizontal flow, Grad P becomes the pressure gradient; otherwise, gravity must be included.) According to these expressions, each of the constituent phases is considered similar to a homogeneous system where the volumetric rate of flow is proportional to the pressure gradient, and for each of which the constants of proportionality, k, and kG, are termed effective permeabilities. by analogy to the specific permeability as defined by Darcy&apos;s law in its original form. In order to obtain a convenient basis of comparison, the effective permeabilities are referred to the specific permeability, k, of the considered porous medium, with the help of the relations: k1. = KI. k kQ = Kn k..........(2) K,. and Kr. are defined as the relative permeabilities to the liquid and to the gas phase respectively, and frequently expressed in per cent of specific permeability. It may be seen that Equations (1); which appear to be a direct generalization of Darcy&apos;s law, correspond to the assignment at any given time of a set of "local" permeabilities to each point of the porous medium, and represent in a differential form the two fluid flow system as a simple superposition of the individual single flow systems. The above interpretation implies that the effective or relative permeabilities are independent of pressure or rate of flow,

Jan 1, 1951

By L. A. Rapoport, W. J. Leas

As a preliminary, consicleration is given to the conventional definition of relative permeability and to the conditions governing the simultaneous flow of oil and gas through porous media. For the conditions of flow prevailing throughout most of a gas drive reservoir, the oil and gas can reasonably be supposed to be in capillary equilibrium with each other. Under these conditions, and these conditions only, the relative permeability to liquid can be expressed as a function of saturation. The relative permeability to liquid in that case is dependent upon the distribution of fluids which itself is shown to be related to the capillary pressure, and, in turn, to the saturation. As a consequence, relative permeability to liquid can be expressed in terms of the volume and surface area of a network of liquid channels bounded by the rock and the gas phase. While the volume of this network can be evaluated accurately, the surface area cannot. However, for any such volume, maximum and minimum values of the corresponding surface area can be calculated from capillary pressure data. It is then possible to establish for any saturation the limits within which the value of the relative permeability to liquid must lie. As a consequence of the theoretical development, the validity of an experimental method for measuring relative permeability to liquid which utilizes a stationary gas phase is demonstrated. In this method capillary barriers are cemented to the ends of the core sample to permit the maintenance of capillary equilibrium between the two phases. At the same time, this procedure eliminates undesirable secondary phenomena such as end effects, fissure effects, etc., the presence of which adversely affect the results of other laboratory methods. The results obtained by theoretical calculations, and experimentally, are discussed. In view of the overall precision that can presently be obtained in reservoir calculations, the agreement between the calculated and measured relative permeability to liquid data can be considered satisfactory. In conclusion, for reasons of economy and simplicity, the procedure of calculating limiting relative permeability to liquid curves from capillary pressure data is indicated for general engineering purposes. It is shown that the above procedure can easily be extended to the cases where connate water is present. Its use for reservoir studies is particularly recommended in conjunction with the method for measuring relative permeability to gas&apos; which simultaneously yields the capillary pressure data necessary for the calculations. THEORETICAL Definition of Relative Permeabilities — Basic Equations for Heterogeneous Flow The equations by which the relative permeability concept is defined and upon which the formulation of all of the gas-oil flow problems rests at the present time are expressed as: V, = — Grad PL = — Grad PL ....(la) PL µL, k k KG VG = —K - Grad PG k- Grad Pc .... (lb) Mo where ,. and G refer to liquid and gas; V is the volumetric rate of flow per unit gross area. µ the viscosity, Grad P the potential gradient. and k the specific permeability of the porous medium.* (For horizontal flow, Grad P becomes the pressure gradient; otherwise, gravity must be included.) According to these expressions, each of the constituent phases is considered similar to a homogeneous system where the volumetric rate of flow is proportional to the pressure gradient, and for each of which the constants of proportionality, k, and kG, are termed effective permeabilities. by analogy to the specific permeability as defined by Darcy&apos;s law in its original form. In order to obtain a convenient basis of comparison, the effective permeabilities are referred to the specific permeability, k, of the considered porous medium, with the help of the relations: k1. = KI. k kQ = Kn k..........(2) K,. and Kr. are defined as the relative permeabilities to the liquid and to the gas phase respectively, and frequently expressed in per cent of specific permeability. It may be seen that Equations (1); which appear to be a direct generalization of Darcy&apos;s law, correspond to the assignment at any given time of a set of "local" permeabilities to each point of the porous medium, and represent in a differential form the two fluid flow system as a simple superposition of the individual single flow systems. The above interpretation implies that the effective or relative permeabilities are independent of pressure or rate of flow,

Jan 1, 1951

By J. K. Elliott, P. H. Kelly

Many engineering functions such as surface metering work and laboratory compressibility check points involve the use of liquid densities of light hydrocarbon mixtures at various pressures and temperatures. However, at the present time, no simple reliable method exists for determining density variation, particularly if the composition of the liquid is unknown. Consequently, a study was undertaken to develop and present a simple and accurate method of predicting density variation of a light hydrocarbon liquid with pressure and temperature, knowing only the density of the liquid at some condition. The experimental liquid compressibility data from API Project 37 by Sage and Lacey&apos; have been considered to be accurate within 0.5 per cent and cover a wide range of pressure (14.7 to 10,000 psia), temperature (100" to 400°F) and molecular weight (up to 150). From these data, a set of liquid density curves, which relate density to pressure, temperature and molecular weight, was developed. These curves make it possible to predict density variation with pressure and temperature. Compared to extensive laboratory compressibility data on a complex, light hydrocarbon liquid, the use of the liquid density curves resulted in an average error of less than 0.5 per cent. Based on the results of this analysis, it is concluded that the set of liquid density curves developed from the data of Sage and Lacey provides an accurate and simple method for predicting the density variation of light hydrocarbon liquids when the density at some condition is known. These curves should be very helpful in many engineering calculations, particularly in the surface metering of light hydrocarbon liquids. INTRODUCTION Many situations arise in field and engineering laboratory work, such as reservoir engineering studies, check of experimentally determined laboratory data and orifice flow-meter formulas, where liquid density factors at various pressure-temperature conditions are required. Also, the need for accurate light hydrocarbon liquid information has become more important with the advent of miscible-type displacements for secondary recovery purposes in oilfield operations. Several reliable methods are available1 - "or determining the density of liquid hydrocarbons if the composition of the liquid is known. However, there is a definite lack of methods for accurately determining the variation of density when the composition of the liquid is unknown. The purpose of this study is to review the various methods for determining hydrocarbon liquid densities and to develop a simple and reliable method of determining variation in density of light hydrocarbon liquids with pressure and temperature when the compositio~n of the liquid is unknown. METHODS FOR DETERMINING DENSITY OF LIQUIDS OF KNOWN COMPOSITION Sage, Lacey and Hicks&apos; have proposed a method to predict the density of light liquid hydrocarbons by using partial molal volumes. Data are available on experimentally developed partial liquid volumes of hydrocarbons over a rather limited range of temperature, pressure and composition. The partial mold volume method has proved satisfactory for determining the density of some hydrocarbon liquids when the composition is known. Within the range covered in the Sage, Lacey and Hicks1 data, the results agree within about 3 per cent of the experimental values. Hanson mentions the limitation of this method to a composition range of approximately 10 per cent by weight of methane, which will not allow this correction to cover most low molecular weight-light hydrocarbon liquids. Standing and Katz2 studied data on light hydrocarbon-liquid systems containing methane and ethane at high temperature and pressure and have presented a method for determining liquid densities, assuming additive volumes for all components less volatile than ethane and using apparent densities for methane and ethane. The compressibility and thermal-expansion curves used by Standing are based on assumptions that compressibility of a hydrocarbon liquid at temperatures below 300°F is a function of the liquid density at 60°F and that thermal expansion of the liquid is affected little by pressure. The information required to use this technique with an example problem is furnished by Standing.&apos; Hanson eports an average error of - 0.5 per cent using the method of apparent densities in calculating liquid densities of several volatile hydrocarbon mixtures. However, as implied, the apparent density method is not applicable for liquid density calculations when the composition of the liquid is unknown. Watson- as presented a method

By Nicholas J. Grant, Michio Yamazaki

Cast copper alloys containing 10, 20, and 30 pct Ni and 0.75 to 0.80 pct Al were machine-milled into chips, then comminuted in a rod mill to fine flake powder utilizing a number of processing variables. The powders here internally oxidized, mostly at 800°C, in a low-pressure oxygen atmosphere. The consolidated powders were hot-extruded into bar stock. Room-tenmperature tension tests, stress-rupture tests mostly at 650°C, but also at 450° and 850°C, and hardness measurements after various annealing temperature treatments to study alloy stability were perfomted. Excellent room-temperature strength, high rupture strength at 650°C, and resistance to recrystallization at 1050°C were obtained. Problems in optimizing conditions for internal oxidation of Cu-Ni base alloys are discussed. THE interesting high-temperature properties of SAP&apos; have stimulated considerable effort in the study of more refractory alloy systems where the potential for high-strength alloys at high temperature is great.2-13 A number of methods have been utilized to produce the desired fine, hard particle dispersions, of which internal oxidation2,9,7 of dilute solid-solution systems offers considerable promise by virtue of the potential for producing ultrafine, well-dispersed oxides. While most of the published works are concerned with pure metal matrices, a number of investigators have studied the effects of solid-solution strengthening.10,19 Use of more complex alloy matrices (for example, aging systems) has been unsuccessful because overaging still occurs at high temperatures in the oxide-containing alloys.14.15 Solid-solution strengthening is, however, effective at very high temperatures9,10 and might be expected to contribute importantly to the strength of oxide-dispersion strengthened alloys. For this study, internal oxidation of solid-solution alloys of copper and nickel, containing small amounts of aluminum, was chosen as the method of alloy preparation. PREPARATION OF ALLOYS Three copper alloys containing about 10, 20, and 30 pct Ni and each containing 0.75 to 0.80 pct A1 (enough to yield about 3.5 vol pct alumina) were prepared as air-cast ingots measuring 2.5 in. diameter by 6 in. high (see Table I for the analyses). Processing steps for all the alloys were as follows (also see Table 11): 1) Homogenization of the ingot at 982°C (1800°F) for 45 hr in an argon atmosphere. 2) Machine milling of ingots into fine chips. Average thickness was about 0.1 to 0.2 mm. 3) Hydrogen reduction of chips at 593°C (1100° F) for 1 hr to reduce copper and nickel oxides. 4) Rod milling of chips to finer powders. 5) Hydrogen treatment of powders as in step 3. 6) Internal oxidation of the powders. 7) Hydrogen treatment of oxidized powders as in step 3. 8) Hydrostatic compression of evacuated powders. 9) Sintering of compacts in hydrogen. 10) Hot extrusion. Variations in processing among the alloys were made in steps 4, 5, and 10 (see Table 11). In the past, two methods were utilized to internally oxidize alloy powders. Preston and Grant3 surface-oxidized dilute Cu-Al powders to obtain the necessary amount of oxygen to oxidize the solute metal (aluminum and silicon), and then permitted the formed copper oxide to diffuse and react with the solute in an argon atmosphere. Bonis and Grant4 exposed Ni-A1 and other nickel alloys to an oxygen pressure derived from the decomposition of nickel oxide at a preselected temperature, in an argon atmosphere. Both methods are applicable and can be modified to generate a range of oxygen pressures for oxidation of the solute but not the solvent metals. Procedure I: Surface Oxidation of Alloy A3, Cu-10Ni-0.76A1. Powders of -20 to +28 mesh were surface-oxidized at 500°C (932°F) to obtain the desired amount of oxygen for oxidation of the aluminum to alumina; the powder was then sealed in Vycor and heated at 900°C (1652°F) for various times up to

Jan 1, 1965

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&apos;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 J. C. Melrose, W. R. Foster, J. G. Savins, E. R. Parish

Many differences can be imagined between gas-oil flow in which the gas is supplied at the face of the core and gas-oil flow in which the flowing gas was originally dissolved in the oil. If capillary pressure characteristics and flow requirements control gas saturation distribution, the gas would be expected to be located at preferred sites within the porous medium as determined by pore sizes. On the other hand, during solution gas drive the gas first appears as bubbles through a nuclea-tion process. Nothing in self-nucleation theory specifies at which sites the first bubbles should be formed. In all probability they will be randomly distributed throughout the porous medium. Furthermore, it is not at all certain that even at low rates of production the gas will redistribute itself after nucleation to the channels normally occupied by gas in simple gas flow. Stewart, et at, have shown that at least for some limestone samples, oil recoveries could not be predicted for all rates of production using any one set of relative gas and oil permeabilities. An important factor in controlling recoveries during solution gas drive was the rate of bubble formation, higher rates giving higher recoveries. Stewart, et al, attributed the increase in recovery to a better distribution of the gas phase in heterogeneous limestone samples than is obtained by simple external gas drive. Differences in recovery from these causes were not reported for sandstone cores. In the experiments to be reported here, oil recovery, pressure and producing GOR history were measured during solution gas drive for a 5-ft sandstone core. The results were compared with predictions from the Muskat method for computing solution gas-drive behavior using external gas-drive relative permeability. The effects of changing the rate of production and oil viscosity were studied. At high laboratory rates of average pressure decline, two observations were made which would not have been predicted by Muskat&apos;s depletion theory: (1) oil recovery increased with increasing rate of production for a given viscosity oil, and (2) oil recovery increased with increasing oil viscosity for a given high rate of production. Both of these observations are explained as consequences of diffusion control of gas saturations superimposed on the normal gas-oil flow requirements, Fur- thermore, discontinuous gas phase flow appears to be significant during solution gas drive. The laboratory tests were performed at rates of average pressure decline many times greater than the maximum possible rate of average pressure decline in an actual oil field. It is, therefore, not possible to draw any direct conclusions regarding the effect of rate on recovery for the solution gas-drive mechanism under actual field conditions. However, at the lower laboratory rates, recoveries were nearly independent of rate and could be predicted by the Muskat method, using external gas-drive relative permeability data. These results suggest that at normal oilfield rates the effect of rate on recovery for the solution gas-drive mechanism is negligible. EXPERIMENTAL PROCEDURES The core material for the pressure depletion studies was Bandera sandstone from an outcrop in the Mid-Continent. This sandstone was selected because of its low permeability (about 10 md), which would permit the development of substantial pressure gradients in the corn at moderate flow rates. The core was 5-ft long and 2-in. in diameter. Its properties are listed in Table 1. Relative gas-to-oil permeability ratios were measured by an external gas-drive method2. The results are shown for a short 2-in. core and for the 5-ft core in Fig. 8. Oils used in the pressure depletion experiments were kerosene and a highly refined white oil (standard white oil No. 3) with gas-free viscosities of 1.8 and 25 cp, respectively. The gas was a naturally occurring methane from Gough field, Inglewood, Calif. The oil viscosities, gas solubilities and formation volume factors are plotted as functions of pressure at 75°F in Figs. 1 and 2. Methane viscosities and compressibilities were obtained from the literature"&apos;. A core mounting was required which could withstand up to 2,500 psi internal pressure. This was obtained by first encasing the core completely in a plastic resin (Scotch Cast, manufactured by Minnesota Mining & Manufacturing Co.) The plastic covered core was then inserted into a steel pipe equipped with screw caps so that the plastic coating could be pressured from

By Craig S. Hartley

Expressions are obtained for the energy changes associated with the reaction of (a& (111) slip dislocations on intersecting (110)planes in anisotropic bcc metals. An energy criterion for assessing the likelihood of dissociation of the products of such reactions is also presented. It is found that the "burrier reactions" which form a(100) dislocations at the intersection of two active {110) slip planes are more energetically favorable in metals which exhibit a high value of Zener&apos;s anisotropy factor, A, than those which have a low value. The results are presented in a form which permits the stacking fault energy to be obtained from a measurement of the separation between par-tials in a dissociated configuration. However, until accurate calculations or measurements of the stacking fault energies involved are available, it is not possible to assess the physical importance of dissociated dislocations. In a recent paper,&apos; the energy changes associated with several types of reactions between two slip dislocations, (a/2)(111){110), in bcc structures were calculated.* Isotropic elasticity and the approxima- tion v = -3- were employed. The purpose of this work is to present calculations of the energy changes for many of the same reactions using anisotropic elasticity. The problem of dissociation of a(100) and a(110) dislocations is also considered, and maximum fault energies for which dissociation will be energetically favorable are calculated for several bcc metals. Two general types of reactions are considered; those for which the reactant (a/2)(111) dislocations have long-range attractive forces and those for which the reverse is true. An example of the former is: (a/2)[lll] + (a/2)[lll]-a[l00] while the latter are typified by: (a/2)[lll] + (a/2)[111] -a[011] Only reactants lying in different slip planes are considered; therefore, the products must lie along (111) or (100) directions, which are the intersection of two {llO} planes. It will be assumed that the reactants and products are infinitely long parallel dislocations, since in this case the energy change associated with the reactions is a maximum.&apos; THEORY The self-energy per unit length of a straight mixed dislocation in an anisotropic medium can be written? where b is the Burgers vector, K is an appropriate combination of the single-crystal elastic constants, and R and ro are, respectively, outer and inner cut-off radii of the elastic solution. The energy given by Eq. [I] does not account for any variation of the core energy with orientation. This could be manifested by an orientation dependence of the core radius or, equivalently, the Peierls width, of the dislocation. However, the energy contribution due to this source is expected to be small, and current models of the dislocation core are not sufficiently accurate to justify such a refinement. It has already been shown that for the isotropic case the energy contributions due to nonzero tractions across the cores of the reactants and products exactly cancel one another in the reaction.&apos; Accordingly, it will be assumed that this contribution to the total energy change in the anisotropic case is small. In the subsequent discussion it is also assumed that the core radii of the reactant and product dislocation are the same and that, where stacking faults are formed, the faulted region is bounded by the centers of the partials. Consequently only changes in elastic energy due to the reactions will be considered. When the dislocation is parallel to either the (111) or the (100) directions, K may be written:375 K = (Ke sin2 a + Ks cos2 a) [2] where K, and Ks are the combination of elastic constants corresponding to an edge and screw dislocation lying along the same direction as the mixed dislocation, and a is the angle between the direction tangent to the dislocation line and the Burgers vector. Eq. [2] should not be confused with the isotropic approximation to the variation in energy with line Orientation.6 It should be noted that the essentially isotropic expression for K is a result of the characteristic symmetry of the (111) and (100) directions and is not, in general, valid for other dislocation directions in anisotropic cubic metals. The energy* change for a reaction in which the re- actant and product dislocations are parallel perfect dislocations can be written: where Ep and E, refer to the self-energies of the products and reactants, respectively. For dislocations parallel to (100) and (111) directions, Eq. [3] becomes:

Jan 1, 1969

By B. G. Matson, M. A. Whitfield, G. R. Dysart

This paper describes the development of u computer program to calculate changes that occur in the length of tubular goods due to temperature and pressure changes during stimulation operations. Due to the numerous variables involved and the uncertainty of all static and dynamic conditions that could exist, it becomes a staggering task for individuals charged with completions to perform the necessary mathematical calculations. The computer program permits advance calculations for several sets of conditions. INTRODUCTION In the Delaware basin of West Texas alone, 50 wells were contracted or drilled to 15,000 ft or deeper in 1965. Deep well activity is continuing in this and other areas on an expanding scale. Many of these deep wells require extensive stimulation for successful commercial production, and during these operations, pressures and temperatures are encountered that have a pronounced effect on the length of tubular goods. This length change during a large-volume, high-pressure stimulation treatment utilizing fluids considerably cooler than bottom-hole temperature can be of such a magnitude that permanent damage to casing and tubing will result unless mechanical design, pressures and fluid temperatures are evaluated and controlled. These pressure and temperature effects can be calculated. However, the process lends itself well to computer solutions because of the mathematical nature of the problem and the calculating hours involved in arriving at an answer. The engineering-hour demand becomes more severe as tapered strings are involved. On initial treatments on a given well, surface pressure and injection rate conditions are unknown, and offset well conditions have not proven to be a reliable method for making predictions. For these reasons, it has become rather standard procedure for operators to compensate for these uncertainties by placing unnecessary pressure and fluid temperature restrictions on stimulation design. On a number of occasions treating fluids have been preheated to as much as 160F as a means of compensating for thermal contmction resulting from pumping cool fluids. The maintenance of packer seals has been treated by Lubinski, Althouse and Logan&apos;,&apos; and the problem of therma1 effects on pipe has been explored by Ramey." These works were expanded and the results made applicable to everyday oilfield terminology before submitting them to computer programming. The pressure and temperature effects on tubing movement previously mentioned occur simultaneously as fluid moves through the pipe. The pressure changes, for purposes of explanation, are categorized here as to the various effects these pressures have on a tubing string. These divisions are (1) the piston-like results of forces acting on horizontal surfaces exposed to pressure, (2) swelling or ballooning of the tubing along its entire length due to the forces of pressure acting against the tubing walls, (3) the elongation of tubing due to frictional drag and (4) corkscrewing of the pipe due to internal pressure. Thermal changes are also of great importance, as their results may be more significant than any of the pressure effects. Steel is an excellent conductor of heat and the earth is a relatively poor conductor. It has been calculated that pipe temperatures at depths of more than 20,000 ft approach within as little as 25" the temperfature of the surface fluid after pumping for 2 hours, or a drop in temperature in some treatments of more than 220F. The equations presented in this paper were developed for computer programming and simplicity of input information; therefore, numerical constants such as Young&apos;s modulus for steel (28 X 10\ si), the coefficient of thermal expansion of steel (6.9 X 10."IF) and Poisson&apos;s ratio for steel (0.3) are included with unit conversion factors. The moment of inertia of tubing cross-sectional area with respect to its diameter was changed to a constant times (D&apos; — d&apos;) where D is outer diameter and d is inner diameter. Units in the equations are length in feet, diameter in inches, density in pounds per gallon, pressure in psi, rate in barrels per minute and time in hours. PISTON-LIKE REACTIONS A change in tubing internal dimensions and the exposure of other horizontal surfaces to different pressures on the inside and outside of the tubing result in a reaction much like a piston under pressure. Such is the case when the internal diameter changes in a combination string of pipe, when seals of a slick joint assembly are subject to pressure and in the end effects of a tubing string. The change in tubing length due to the piston effects of a slick joint packer is affected by the various diameters involved, the tubing pressure Ap,, the casing pressure ,Ap,, length of pipe L, densities of fluid in the tubing before and during pump-

By K. A. Slagle

Hydraulic analysis of the wellbore has become increasingly inzportant for designing cementing operations and selecting equipment, materials and techniques to complenzent modern well-c-ompletion practices. Non-Newtonian fluid technology has advanced beyond the point where former empirical methods of analysis adequately define the hydraulic system and fluid properties. In view of these factors, this paper describes a series of rheological calculations which have been found practical, through field usage, for assistance in selecting a cementing program. A relatively simple laboratory method using standard viscometric equipment is suggested for determination of the rheological properties of slurries, and clrrta are presented on some of the more common cementitrg conzposition.A. A criterion for divergence from laminar-flow characteristics has been proposed. Usefulness of the calculations is indicated by examples of cementing operations where they have been used. INTRODUCTION With the changing aspects of well-completion practices during the past few years, it has been increasingly important to have a relatively simple method of analyzing the flow conditions existing in the well during cementing operations. This is particularly true in view of the improved economics toward which most of the changes have been directed. Rheological characteristics of slurries used for cementing should be a major consideration in the trend toward smaller casing sizes, either single or multiple strings. Receiving increased attention is the practice advocated in 1948 by Howard and Clark&apos; of attaining turbulent flow with the fluids circulated during a primary cementing operation. While there may still be a difference of opinion concerning this technique, most available information indicates that superior primary-cementing results are generally obtained when high displacement rates are employed. Fluid properties of the slurry to be used must be available, as well as calculation methods, to determine what flow rates should be attained and the probable consequences in terms of frictional pressure and horsepower utilization. It would certainly be inappropriate to attempt high displacement velocities if sufficient pressure might be developed to create lost circulation. Since cementing slurries are non-Newtonian fluids, it is not possible to define their rheological or fluid properties by the single factor of viscosity and then make calculations for the quantities just described. Because the shear stress-shear rate ratio is not constant: it becomes necessary to establish at least two parameters for adequate fluid-flow calculations. It is not the purpose of this paper to delve into the mathematical development of non-Newtonian technology, nor to discuss the arbitrary classification system under which a single fluid may resemble two or three different classes depending upon experimental conditions. Rather, the intention is to present a useful series of calculations based on a concept applicable to both Newtonian fluids and to the preponderance of non-Newtonian fluids encountered in the oil-producing industry. Development of this approach was begun some 32 years ago,&apos; and has most recently been brought to fruition by Metzner and his co-workers at the U. of Deleware. Some non-Newtonian fluids encountered in the petroleum industry, other than cementing slurries, have also had the benefit of this method of analysis."&apos; The two parameters required to define the fluid are usually denoted by the symbols n&apos; and K&apos; and, for the purposes of this discussion, are called "flow behavior index" and "consistency index", respectively. These two slurry properties permit calculation of the Reynolds&apos; number and the "critical" velocity, or the velocity at which departure from laminar flow begins. EXPERIMENTAL DETERMINATIONS The two principal instruments used for rheological studies are the pipeline (capillary-tube) viscometer and the rotational viscometer. When conveniently possible, a capillary-tube viscometer (where the pressure drop and flow rate of the material can be measured) is the better method for rigorous determination of the flow behavior index and consistency index for non-Newtonian fluids. With pressure-drop data at various flow rates, it is then possible to prepare a logarithmic plot of shear rate as the abscissa-shear stress as the ordinate. For fluids which do not exhibit time-dependency, these data will usually produce a straight line. The flow behavior index n&apos; represents the slope of this line, while the consistency index K&apos; becomes the intercept of this line at unity shear rate in accordance with the mathematical derivations associated with this concept of rheology. Due to the difficulties anticipated in maintaining a uniform, pumpable cement slurry for the time interval required to obtain measurements from the pipe viscometer, the n&apos; and K&apos; data reported herein were obtained using a direct-indicating rotational viscometer (Fig. 2). The

By Harry C. Swanson

TRANSPORTING concrete from mixer to forms has always been a problem. Twenty-five years ago this task was generally accomplished by means of wheelbarrow or concrete buggy. On large dam jobs, as the number of these projects increased, the gantry crane or highline came into use. Today several methods of handling concrete are employed on smaller surface construction jobs, for example, transit-mix trucks or dumpcrete trucks, which have crawler cranes with buckets for placing concrete into forms. In 1944, during early stages of developing Mather mine A shaft, several large underground concrete jobs were necessary. At this time the Cleveland-Cliffs Iron Co, purchased the first pump-crete machine, introduced by the Chain Belt Co. of Milwaukee. The machine was used to pour approximately 200 cu yd of concrete for a dam, or bulkhead, located 400 ft from the shaft. Concrete was mixed on surface, lowered down the shaft 1000 ft in a 2-cu yd bucket hung under one skip, spouted into the bowl of the pumpcrete machine from the bucket, and pumped directly into the forms. Since the day of the first pipeline concrete in 1944 to the present time, other equipment and other methods have been developed to permit transportation of concrete by pipeline through vertical and horizontal distances totaling 1 mile from mixer to forms. Much of the efficiency in present handling of underground concrete can be credited to the Bethlehem Cornwall mines, where concrete was transported through 6-in. pipe for great distances down an inclined shaft and along levels into forms.&apos; During initial development of Mather mine B shaft, with concrete work under way on two or more levels at one time, the pneumatic concrete placer, Fig. 1, was selected as best adapted for underground concrete transportation. The 3/4-cu yd pneumatic placer is a small machine readily moved from one location in the mine to another. It can be equipped with two sets of mine car wheels, which will permit moving on regular mine tracks. It is therefore possible to send concrete through the pipe at great velocity; the pipeline is clean after each shot except for the film of cement adhering to the inside. With the proper slump in the concrete, it is possible to shoot concrete 2000 ft with this machine, using the mine supply of compressed air at 95 psi. This equipment was first used at Mather mine B shaft to concrete slusher drifts, Figs. 2 and 3, and finger raises located about 2000 ft from the shaft. In several instances there were bends into crosscuts and up vertical distances into the forms. For the first pours two placers were used. The first was located near the shaft where the concrete could be spouted into it from a 2-cu yd concrete bucket on the cage. The second was set on the side of the drift at a point approximately 1500 ft from the shaft. The concrete was shot directly into the second placer from the first unit and from the second machine directly into the forms. After completion of several pours with the two machines, a trial pour with only one placer located at the shaft proved that the second placer could be eliminated. Since then all pours have been successfully completed with only one placer underground. As production of iron ore from the mine increased and the development program expanded, use of the cage for handling mine supplies and concrete became a major problem. This brought about the first attempt at shooting concrete vertically down the shaft for 2600 ft. A 6-in. pipeline with victaulic couplings installed during shaft sinking was used for the trial. One placer was set on surface 250 ft from the collar of the shaft so concrete could be spouted directly into it from the mixer. This machine shot the concrete 250 ft horizontally on surface to the shaft, 2600 ft vertically down the shaft, and 100 ft horizontally into the second placer located near the rib of the shaft station or plat. The second machine shot the batch into the forms, about 2000 ft. Total distance horizontally and vertically was 4800 ft. The entire time cycle for a ¾-cu yd batch of concrete from the mixer on surface to the forms underground totaled about 5 min. During the past two years the two-placer method from the mixer on surface to the forms underground has proved a very efficient means of transporting underground concrete. Advantages of using pipeline concrete are as follows: 1—Interference with normal mining operation is eliminated. When the cage, skips, mine cars, or mine openings are used for transporting concrete and materials used for making concrete, mine operation suffers in one way or another.

Jan 1, 1955

By W. G. Pfann

In zone-melting, a small molten zone or zones traverse a long charge of alloy or impure metal. Consequences of this manner of freezing are examined with impurerespect to solute distribution in the ingot, with particular reference to purification and to prevention of segregation. Results are expressed in terms of the number, size, and direction of travel of the zones, the initial intermsofsolute distribution, and the distribution coefficient. IF a charge of binary solid-solution alloy is melted and then frozen slowly from one end, as for example in the Bridgman method of making single crystals,&apos; coring usually occurs, with a resulting end-to-end variation in concentration. Such coring, or normal segregation, is undesirable where uniformity is an object. On the other hand, for certain systems, it can be utilized to refine a material by concentrating impurities at one end of the ingot.&apos;. &apos; In the present paper a different manner of freezing will be examined with respect to the distribution of solute in the ingot. A number of procedures will be indicated which have in common the traversal of a relatively long charge of solid alloy by a small molten zone. Such methods will be denoted by the general term zone-,melting, while the process described in the preceding paragraph will be called normal freezing. It will be shown that, in contrast to normal freezing, zone-melting affords wide latitude in possible distributions of solute. Segregation can either be almost entirely eliminated or it can be enhanced so as to provide a high degree of sttparation of solute and solvent. A number of simplifying assumptions will be invoked which, while not entirely realizable in practice, nevertheless provide a suitable point of departure for more refined treatments. Moreover, our own experience with zone-melting has shown that, for certain systems at least, the analysis holds quite well. The present paper will be confined to a discussion of principles and a general description of procedures. Comparison with experiment is planned for later publication. Normal Freezing Before considering zone-melting, segregation during normal freezing will be reviewed briefly. If a cylinder of molten binary alloy is made to freeze from one end as in Fig. 1, there usually will be a segregating action which will concentrate the solute in one or the other end of the ingot. If the constitutional diagram for the system is like that of Fig. 2, then the distribution coefficient k, defined as the ratio of the concentration in the solid to that in the liquid at equilibrium, will be less than one and the solute will be concentrated in the last regions to freeze. If the solute raises the freezing point, then k will be greater than one and the solute will be concentrated in the first regions to freeze. The concentration in the solid as a function of g, the fraction which has solidified, can be expressed by the relation: C = kC0 (1-g)k-1 [I] where C, is the initial solute concentration in the melt. Eq 1 is based on the following assumptions: 1—Diffusion in the solid is negligible. 2—Diffusion in the liquid is complete (i.e., concentration in the liquid is uniform). 3—k is constant. Concentration curves representing eq 1 for k&apos;s from 0.01 to 5.0 are plotted in Fig. 3. This equation, in one form or another, has been treated by Gulliver,³ Scheuer,4 Hayes and Chipman5 for alloys and by McFee2 for NaCl crystals. It is derived in Appendix I. It should be pointed out that the k which is calculated from the phase diagram will be valid only in the ideal case for which the stated assumptions are correct. In all actual cases, the effective k will be larger than this value for solutes which lower the melting point, smaller for solutes which raise the melting point, and will probably vary during the beginning of the freezing process. For simplification it will be assumed that the ideal k is valid. Zone-Leveling Processes The processes of this part are designed to produce a uniform, or level, distribution of solute in the ingot. Single Pass: Consider a rod or charge of alloy whose cross-section is constant and whose composition, C2, is constant, although permissibly varying on a microscopic scale." Such a charge might be a rapidly frozen casting or a mixture of crushed or powdered constituents. Cause a molten zone of

Jan 1, 1953

By J. H. Swisher, E. O. Fuchs

Rate equations were derived to describe the kinetics of internal oxidation of cylinders and spheres. The derived equations for cylinders were checked experimentally by means of sub scale thickness and electrical conductivity measurements on Cu-Cr alloy wires. The properties of the internally oxidized samples were examined with conductivity applications in mind. It was possible to produce uniform dispersions of Cr2O3 in copper with an initial chromium content as high as 3 wt pct. While electrical conductivities only a few pct less than that of OFHC copper were obtained, the Cr2Os particle size and spacing were too large for effective dispersion hardening. T.HE process of internal oxidation has been used widely in basic studies of the permeability of gases in metals. In a review article, Rapp1 has discussed the principles of internal oxidation in considerable detail. From a technological standpoint, internal oxidation is often considered undesirable, since it is a means by which inclusions can be introduced into an otherwise clean material. Another important aspect of internal oxidation is its use as a means of dispersion hardening a material. Broutman and Krock2 discuss this and other methods for making dispersion hardened alloys. The only internally oxidized material known to the authors which is commercially available is a Cu-BeO alloy.3&apos;4 This alloy is made from Cu-Be alloy powder, using a so-called Rhines pack. It has a tensile strength of 80,000 psi and retains its strength at relatively high temperatures. The objectives of the present study were to derive rate equations for the internal oxidation of cylinders and spheres, to check the derived equations for cylinders experimentally, and to examine the structure and properties of internally oxidized Cu-Cr alloys. The Cu-Cr system was chosen for this study because uniform dispersions are obtainable at high alloy contents, which is a desirable characteristic in dispersion hardened materials. RATE EQUATIONS FOR VARIOUS GEOMETRIES A number of authors5--9 have derived equations to describe internal oxidation kinetics. These derivations differ somewhat in mathematical assumptions and approximations, and all except one of the derivations deal exclusively with the internal oxidation of plates. The exception is a brief treatment of cylindrical and spherical geometries given by Meijering and Druy-vesteyn9 as a part of a comprehensive paper on the general subject of internal oxidation. These authors did not obtain rate data to check their derivations, although they did show that the hardness profile across an internally oxidized sample is directly related to the rate of interface movement. For cylindrical and spherical geometries, a quasi-steady-state approximation is needed to circumvent mathematical complications in obtaining a solution to the basic differential equations. In using this approximation, we consider the concentration gradient of dissolved oxygen in the internally oxidized zone or sub-scale to be the same as the gradient which would be present if there were no movement of the subscale interface. The steady-state approximation introduces an error of about 1 pct in computing the rate of internal oxidation of an Fe-1.0 pct Mn alloy plate, if the present method is compared to the more exact method of Wagner.7&apos;10 The details of the derivations of the rate equations for cylinders and spheres are given in the Appendix, and only the results of these derivations are given below. The final equations obtained by Meijering and Druyvesteyn9 can be shown to be equivalent to our Eqs. [1] and [2], although the two approaches are somewhat different. Cylindrical Geometry. [2] where r1 is the outer radius of the cylinder or sphere, cm, r2 is the radius of the unreacted core, cm, see Fig. l(a), D is the diffusion coefficient of oxygen in copper, cm2 per sec, %O is the concentration of dissolved oxygen at the surface of the specimen, wt pct, %Cr is the initial chromium concentration in the alloy, wt pct, and t is the reaction time, sec. Plate Geometry. The analogous rate equation for a plate has been derived previously for internal oxidation of Fe-Al alloys.8&apos;11 For Cu-Cr alloys, we may write the same equation as follows: [3] where r1 is the half-thickness of the plate, cm, and r2 is the distance from the mid-plane to the subscale intherate is An analysis of Eqs. [1], [2], and [3] shows that for a plate the rate is completely parabolic. The initial

Jan 1, 1970

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett&apos;, Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,&apos;. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By W. R. Purcell

An apparatus is described whereby capillary pressure curves for porous media may be determined by a technique that involves forcing mercury under pressure into the evacuated pores of solids. The data so obtained are compared with capillary pressure curves determined by the porous diaphragm method, and the advantages of the mercury injection method are stated. Based upon a simplified working hypothesis, an equation is derived to show the relationship of the permeability of a porous medium to its porosity and capillary pressure curve, and experimental data are presented to support its validity. A procedure is outlined whereby an estimate of the permeability of drill cuttings may be made with sufficient acuracy to meet most engineering requirements. INTRODUCTION The nature of capillary pressures and the role they play in reservoir behavior have been lucidly discussed by Lev-rett&apos;, Hassler, Brunner, and Deah12, and others. As a result of these publications the value of determining capillary pressure curves for cores has come to be generally recognized within the oil industry. While considerable attention has been directed toward the subject in an effort to provide a reliable method of estimating percentages of connate water, it has been recognized that capillary pressure data may prove of value in other equally important applications. This paper describes a method and procedure for determining capillary pressure curves for porous media wherein mercury is forced under pressure into the evacuated pores of the solids. The pressure-volume relationships ob- tained are reasonably similar to capillary pressure curves determined by the generally accepted porous diaphragm method. The advantages of the method lie in the rapidity with which the experimental data can be obtained and in the fact that small, irregularly shaped samples, e.g., drill cuttings, can be handled in the same manner as larger pieces of regular shape such as cores or permeability plugs. Based upon a simplified working hypothesis, a theoretical equation will be derived which relates the capillary pressure curve to the porosity and permeability of a porous solid, and experimental data will be presented to support its validity. This relationship aplied to capillary pressure data obtained for drill cuttings by the procedure described provides a means for predicting the permeability of drill cuttings. METHODS FOR DETERMINING CAPILLARY PRESSURES Several techniques have so far been employed in determining capillary pressure curves and these fall into two principal categories: (1) Liquid is removed from, or imbibed by, the core through the medium of a high displacement pressure porous diaphragm (2) Liquid is removed from the core which is subjected to high centrifugal forces in a centrifuge4,&apos;. There are? however, certain limitations inherent in both methods. The greatest capillary pressure which can be observed by method (I), above, is determined by the maximum displacement pressure procurable in a permeable diaphragm which at the present time appears to be less than 100 psi. An even more serious limitation of the diaphragm method is imposed hy the fact that several days may be required to reach saturation equilibrium at a given pressure; hence, the time re- quired to obtain a well-defined curve may be measured in terms of weeks. Furthermore, to date, no suitable technique for handling relatively small, irregularly shaped pieces of rock, such as drill cuttings, has been reported and, therefore, measurements must be made, in general, on cores, or portions thereof. The centrifuge method offers the distinct advantage over the porous diaphragm method of arriving at saturation equilibrium in a relatively short time by virtue of the elimination of the transfer medium for the liquid. The calculation of capillary pressures from centrifuge speeds is somewhat tediousa, however, and the equipment required is fairly elaborate. While there exists the possibility that this method might be adaptable to the determination of the capillary pressures of cuttings, this particular ramification has not been investigated, as far as is known. In view of the limitations of the two principal methods for determining capillary pressures, the apparatus described in the following sections has been devised in order that difficulties previously encountered might be circumvented. MERCURY INJECTION METHOD FOR DETERMINING CAPILLARY PRESSURES Theory The methods described above for determining capillary pressures are characterized by the fact that one of the fluids present within the pore spaces of the solid is a liquid which "wets" the solid, i.e., the contact angle which the liquid forms against the solid is less than 90" as measured through that phase. For these "wetting" liquids the action of surface forces is such that the fluid spontaneously fills the voids within the solid. These forces likewise oppose the withdrawal of the fluid from the pores of the solid.

Jan 1, 1949

By K. J. Gill, P. Chiotti, J. T. Mason

Thermal, metallographic, and vapor pressure data were obtained to establish the pkase boundaries and the standard free energy, enthalpy, and entropy of formation for the compounds in the Y-Zn system. Three coinpounds with stoichiometric formulas of YZn, YZn2, and Y2Zn17 melt congruently at 1105", 1080°, and 890°C, respectively. Four compounds with stoiclziometric formulas of YZn3, YZn4, YZn5, and YZn,, undergo perztectic reactions at 905", 895", 870º, and 685ºC, respectively. Three eutec-tics exisl in this system with the .following eutectic temperatures and zinc contents in wtpct: 875ºC, 23.2 Zn; 1015ºC, 51 Zn; 865ºC, 82 Zn. The YZn, pkase undergoes an allotropic transformation. In the two phase YZn2 -YZn alloys the trans.formation gives a weak thermal arrest at 750°C, whereas in the two phase YZn2-YZn3 alloys no thermal arrest is observed and the transformation occurs over a temperature range below 750°C. At 500°C the free mzergies of formation per lnole vavy from —18,090 for YZn to —53,430 fov YZr11 and corresponding enthalpies vary from -24,050 to -92,080. The free energies and enthalpies per g atom as a function of composition show a maximum for the YZn2 phase; the 500°C values are -9580 and -13,180, vespectively. 1 HE only information found in the literature on Y-Zn alloys was the observation reported by Carlson, Schmidt. and speddingl that Y-20 wt pct Zn forms a low melting alloy. The alloy was produced by the bomb-reduction of YF3 and ZnF2 with calcium in an investigation of methods for producing yttrium metal. The solubility of yttrium in zinc has been determined by P. F. woerner2 and reported by Chiotti, Woerner, and Parry.3 In the temperature range 495" to 685°C the solubility may be represented by the relation In these equations N represents atom fraction of yttrium and T is the temperature in degrees Kelvin. The purpose of the present investigation was to establish the phase diagram for the Y-Zn system and to determine the standard free energy, enthalpy, and entropy of formation for the compounds formed. MATERIALS AND EXPERIMENTAL PROCEDURES The metals used in the preparation of alloys were Bunker Hill slab zinc, 99.99 pct pure, and Ames Laboratory yttrium sponge. Arc-melted yttrium buttons contained the following impurities in parts per million: C-129, N-12, 0-307, Fe-209, Ni-126, Mg-13, Ca < 10, F-105, and Ti < 50. Some of the alloys containing 70 wt pct or more of Zn were prepared from yttrium containing 5000 ppm Ti as a major impurity. Tantalum containers were found to be suitable for all alloys studied and were used throughout the investigation. The pure metals, total weight about 30 g, were sealed in 1 in. diam tantalum crucibles by welding on preformed tantalum covers. A 1/8 in. diam tantalum tube was welded in the base of each crucible for use as a thermocouple well. Welding was done with a heli-arc in a glove box which was initially evacuated and filled with argon. The sealed crucibles were enclosed in stainless steel jackets and heated in an oscillating furnace at temperatures up to 1150°C. Homogeneous liquid alloys were obtained within a half hr at these temperatures except for alloys containing less than 20 pct zinc. The latter alloys were held at 1000º to 1100°C for 2 to 3 days in order to obtain equilibrium. After the initial equilibrations the tantalum crucibles containing the alloys were removed from the steel containers and used directly for differential thermal analyses. Further annealing heat treatments for alloys in which peritectic reactions were involved were carried out in the thermal analyses furnace. After thermal analyses the tantalum crucibles were opened and the alloys sectioned and polished for metallographic examination. In the following discussion alloys referred to as "quenched" were obtained by quenching the sealed stainless steel jacket containing the tantalum crucible and alloy in water. The differential thermal analyses apparatus used was a modified version of the one described in an earlier paper., The graphite crucible was replaced by an inconel crucible, the nickel standard and sampie container were separated by a 1/8 in. MgO plate, no getter was used, and provisions were made to

Jan 1, 1963

By B. L. Landrum, J. Simmons, J. M. Pinson, P. B. Crawford

Patentiometric model data have been obtained to estimate the effect of vertical fractures on the areas swept after breakthrough in water flooding and miscible displacement programs such as gas cycling where the mobility is near one, The data are presented for the case of the fire-spot pattern in which the cemer well is fractured various lengths and orientations, the data indicate that for 10-acre spacing, fractures extetidirrg over 1300 ft in either directior1 from the fractured well may re.srrlt in reductions in sweep efficiencics from 72 to approximately 34 per cent. However. the area swept after break through may be quite largr and only 10 or 12 per cent 1ess than would be obtained if the reservoir were trot fractured. For the specific case when the volume of fluid injected is equivalent to 100 per cent of the pattern vol-unie, the swent area may vary from 80 to 88 per cent, depending on the lenght of the fracture. The former value is that which occurs when the break through or sweep efficiency was orrly 34 per cent and the latter figrrre of 88 per cent is that which is obtained if the reservoir were unfrac-ttm&apos;d. It is pointed out that although the sweep efficiency may he very low in vertically fractured five-spot patterrz.s, the area swept at low water-oil ratios may be only 5 to 10 per cent less than those achieved if the reservoir were unfractured. INTRODUCTION Since the initiation of commercial reservoir fracturing techniques it has been desirable to determine the effect of fractures on the areas swept after breakthrough. Most water flooding or gas cycling projects are continued for substantial periods after the brcakthrough of the injected fluid. Although the sweep efficiency serves as one criterion for rating various flooding patterns. the area swept after breakthrough for various water-oil ratios or percentage wet gas, if cycling. is of perhaps more importance than the sweep efficiency alone. Sweep efficiency data on the vertically fractured five-spot have been presented3. Previous work on the line-drive pattern has shown the effect of vertical fractures on the area swept after breakthrough for the case in which the distance between injection and producing wells divided by the distance between adjacent input wells was equivalent to 1.5 (see lief. 2). The data indicated that for the line-drive pattern it may be desirable to flood or cycle substantially perpendicular to the fractures in order to achieve the greatest recovery for the smallest volume of fluid injected. For this study the center well of a five-spot is assumed as the fractured well. All fractures were assumed to originate at this well and extend into the reservoir for various distances and orientations. All the fractures are straight and are of large permeability compared to the matrix proper. These data are presented to aid the engineer in estimating fractured five-spot pattern performance. ANALOGY The potentiometric model was used in making this study. The model used was 20 20 in. by approximately 1-in. deep. For certain portions of the study one corner of this model was considered to be an injection well and the opposite corner a production well. To simulate vertical fractures a copper sheet was soldered to the wire well and made to conform to the desired length and orientation. In other studies the same model was used except that the four corners of the model might be considered as the corner wells of a five-spot pattern and a fifth well was placed in the center of the model. The well placed in the center of the model was fractured. The total fracture length is L and the well spacing. d. The complimentary fracture angles will be obvious from Figs. 3 and 4. The data obtained on the potentio-metric model assumes the pay to be uniform and homogeneous, the mobility ratio is one, steady-state conditions exist and gravity effects arc neglected. The permeability of the fractures is very great compared to that of the matrix proper. The po-tentiometric model has been used widely both in water flooding and gas cycling projects, and may be used for miscible displacement; how-ever. it is believed that the poten-tiometric model data are more properly applicable to gas cycling than water flooding because the model as-

By D. B. Robinson, C. A. Macrygeorgos, G. W. Govier

Experimental data have been obtained on the volurrletric behavior of ternary mixtures of methane, hydrogen sulfide and carbon dioxide at temperalures of 40°, 100" and 160°F up to pressures of 3,000 psia. The results indicate that the compressibility factors for this system do not agree with compressibility factors for sweet natural gases at the same pseudo-reduced conditions. The deviation increases as the temperature and methane content decrease. Discrepancies of up to 35 per cent were observed. A careful analysis has been made of the existing pUrblished data on compressibility factors for binary systems containing light hydrocnrbons and hydrogen sulfide or carbon dioxide. It has been found that the deviation of actual from predicted compressibility factors for methane-acid gas mixtures is a function of the methane content and the pseudo-critical properties,.v of the mixture. The ratio between actual compressibility factors for methane-acid gas mixtures and compressibility factors for sweet natrlral gases at the same pseudo-reduced conditions has been currelated over a range of pP,, from 0 to at least 7 arid a range of pT, from about 1.15 to at 1east 2 0 with an error not exceeding 3 per cent and over most of the range within I per cent. The validity of the correlation for mixtures containing appreciable hearvier hydrocorbons has not been fully established, but it is shown to be preferable than the use of a corretation based only on hydrocarbons. INTRODUCTION Although a relatively accurate method for predicting compressibility factors of pure materials is provided by charts based on reduced properties and the assumption that the compressibility factor is a unique function of T P and z the determination of the correct values of compressibility factors for gas mixtures is somewhat difficult. Two general methods of dealing with gaseous mixtures have been proposed. The first assumes a direct or modified additivity of certain properties of the mixture in terms of the properties of the individual components. Examples of this method are based on the familiar laws of Dalton and Amagat. The second method averages the constants of an equation of state applicable to the pure components. Both of these methods are of limited value in engineering calculations because the first usually provides reliable answers only over narrow ranges of pressure and temperature and the second is cumbersome to handle. In petroleum engineering practice accurate estimations of the volumetric behavior of natural gases arc frequently required. To fulfill this need, several generalized compressibility charts have been developed.&apos; &apos; Of these, the one prepared by Standing, el al is widely used at present. In the construction of charts of this type a third method for dealing with mixtures has been followed. It is based on correlation of pseudo-critical properties as outlined by Kay and calculated from the critical properties of the individual components in a mixture. Although these charts provide relatively accurate information on the compressibility of dry or wet sweet natural gases, they are less reliable when used for gases containing high concentrations of hydrogen sulfide or carbon dioxide or both. Thus, an experimental program, although time consuming, is the best means now available for the determination of the volumetric behavior of sour or acid gas mixtures. An increased interest in the behavior of these gas mixtures, particularly in connection with some of the fields in Western Canada where the acid gas concentration of the reservoirs may be as high as 55 per cent and where hydrogen sulfide alone may be as high as 36 per cent, provided the incentive for this study. It was the purpose of the investigation to determine the volumetric behavior of selected mixtures of methane, hydrogen sulfide and carbon dioxide over a range of temperature from 40" to 160°F and at pressures up to 3,000 psi. EXPERIMENTAL METHOD The apparatus used in this investigation was basically the same as that described by Lorenzo.&apos;" The amount of each pure component used in preparing the gas mixtures was measured over mercury in a glass-windowed pressure vessel. The pure components were then transferred individually in the desired amounts to a second glass-windowed pressure vessel where the volumetric behavior of the mixture was determined. Volume was varied by mercury injection or withdrawal. The capacity of the cell was about 125 cc. Temperatures in the cells were measured with copper-constantan thermocouples and a Leeds Northrup semi-

By J. B. Conway

Stress-vuptuve data for several of- the refractory metals are frequently found to yield a linear relationship between the Larson-Miller parameter and the logarithm of the applied stress. In such cases linear stress-rupture isotherms result with slopes bearing a definite relationship to the temperature. It also follows that the stress to produce rupture in a certain period of time will be linear in temperature. Data for several refractory metals have been reviewed and excellent linearity is shown in this type of isochronal plot. In addition, the af ore - mentioned lineavity leads to a linear relation between the log of the stress to produce rupture in a certain time and the homologous temperature. This has been illustrated for the Group VI-A metals, tungsten and molybdenum. EXTENSIVE use has been made of the Larson-Miller&apos; parameter in the interpolation and extrapolation of stress-rupture and creep data. In those cases where this particular parametric approach is applicable a convenient and fairly straightforward procedure is made available for the correlation of experimental stress-rupture data. It is quite common to employ this parameter in the form of a master rupture plot in which the parameter, T(C + log tr), is expressed as a function of log stress. In many cases this functional relationship in log stress is linear within acceptable accuracy and hence the following relation results: where P is the parameter, C is the Larson-Miller constant, T is the absolute temperature, t~ is the rupture time, a is the stress, and a and b are constants. Examples of such a relationship are contained in the work of Green, Smith, and 01son2 dealing with high-temperature rupture behavior of molybdenum and in the work of Green&apos; dealing with the high-temperature behavior of tungsten. In addition, pugh4 has shown a similar linearity for some fairly low-temperature data for molybdenum. It can be shown that when the relationship in Eq. [I] is exhibited certain generalizations can be made concerning the form of the stress-rupture isotherms. For example, rearranging yields: For a given material (constant C) at a given temperature the first term on the right-hand side of Eq. [2] is a constant and hence this equation defines a straight line when log stress is plotted as a function of log-rupture time. This is recognized as the standard form usually employed in this type of data presentation. Such linearity then suggests the linear form of the Larson-Miller parameter. Or, in other words, the linear parametric relationship in Eq. [2] results only when the stress-rupture data are linear on a log-log plot of stress vs rupture time. Another interesting observation can be made in regard to Eq. [2]. It can be noted that the slope of the stress-rupture isotherms is given by - T/b and hence a direct calculation of the constant b is available. It also follows that since the value of b is the same for all temperatures the slopes of the various isotherms on the log-log stress-rupture plot cannot be the same. Indeed, the existence of the relationship in Eq. [2] precludes a system of parallel lines on this common stress-rupture plot. As a matter of fact it further specifies that in addition to being nonparallel the slope of these isotherms must decrease (i.e., become more negative) with increasing temperature. Such a condition is indeed found to exist in the case of the stress-rupture data reported for molybdenum.&apos; As a corollary to the above, it may be stated that stress-rupture data which do not lead to a linear log-log stress-rupture plot or whose isotherms do not exhibit a decrease in slope as the temperature increases will not yield the linear relationship of Eq. [I]. Applying Eq. [2] to two different temperatures and solving for C yields: Eq. [3] affords a simple and rapid method for calculating the Larson-Miller constant from the log-log stress-rupture plot. The slope of a given linear isotherm is measured and the value of "b" calculated based on Eq. [2] as: slope = - -Tb Then at an abscissa value of 1.0 hr (making log tr in Eq. [3] equal to zero) read the stress corresponding to rupture for two different temperatures. Substitution in [3] yields: A value of the Larson-Miller constant can thus be calculated from a few simple mathematical procedures employing data read directly from the log-log plot of the stress-rupture data. Of course, it is not to be overlooked that the above reasoning has been based on the linear relationship of Eq. [I] being applicable. However, if as mentioned above the log-log plot is

Jan 1, 1967