Search Documents

By Fred H. Poettmann

An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of the compressibility factor of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. Comparison of calculated results with field measured results were made on 20 dry gas wells from a given field and 11 distillate wells from different fields. The agreement between calculated and observed results is good. In addition to their use in the calculation of sandface pressures of flowing wells, the factors can be used in the direct calculation of the static bottom hole pressure of gas wells, the capacity of gas transmission lines, and in the calculation of the theoretical isothermal horsepower necessary to compress a natural gas. Examples demonstrating the use of various equations are given. INTRODUCTION Many of the equations used to calculate the sandface pressures of flowing gas wells from well head data do not take into consideration rigorously the deviation of the natural gas from ideal gas behavior. For low pressure wells this error is not serious. For high pressure wells flowing at high rates this error call be serious. Sandface pressures are used to determine the producing capacity of gas wells as shown in the U. S. Bureau of Mines' Monograph 7.' An equation has been derived for use in calculating the sandface pressure of flowing gas wells in which the variation of compressibility of the gas with pressure is taken into consideration. This variation due to compressibility has been put into both graphical and tabular form. In addition to their use in the calculation of sandface pressures, these "factors" can also be used to give a direct solution in the calculation of the static bottom hole pressure of gas wells, the calculation of the capacity of gas transmission lines, and the calculation of the isothermal horsepower necessary to compress a gas. CALCULATION OF SANDFACE PRESSURE OF FLOWING GAS WELLS The starting point in the derivation of any specific flow equation is an energy balance on the fluid flowing between any two points in the system. This energy balance based on one pound of flowing fluid is expressed by the well known general flow equation: jp~- vdP+ Ah + -— + wt + r. = o . . (i) P, 2 g. where V = specific volume of flowing fluid in cu ft/lb P = pressure in Ib per sq ft ah = difference in height above datum plane in, ft ------ = kinetic energy change of one lb of flowing fluid, v equals velocity in ft/sec and g. = 32.174 Ws = Work done by fluid while in flow (similar to shaft work in driving a turbine) W, = Znergy losses due to the irreversibilities of the flowing fluid This equation contains no limiting assumptions and can be made the basis of any fluid flow relationship by limitations or substitutions. The energy units are usually expressed in term* of foot-pounds of energy per pound of mass. Considering the vertical flow of gas through tubing or annu-lus, the external work W, done by the gas is zero. Let us also assume that the value of the kinetic energy function is small and can be neglected. Equation (1) is reduced to: /P: VdP + Ah+Wt = 0 .......(2) Pi The energy loszes PI can be expressed in terms of the well known Fanning equatinn: 4/Lf: rf = —— ............(3) 2g,D where L = the integrated average velocity in ft/sec of vertical flow in a uniform flow string -~h = L v = the integrated average velocity in ft/sec / = dimensionless correlating function

Jan 1, 1951

By V. Kondic, W. M. Rumball

HUNT and Jackson' have recently suggested that binary eutectics might be classified according to the entropy of fusion of the component phases. Such a classification was in fact proposed2 some time ago by the authors and the purpose of this communication is to present some supporting evidence. As pointed out by Hunt and Jackson' typical metals of low entropy of fusion, ?S, tend to form normal eutectics (using Scheil's terminology3). Where one phase is metallic (low ?S) and the other nonmetallic (high AS) then anomalous structures are formed. These anomalous structures may either be regular (e.g., Bi-Sn)or irregular (e.g., A1-Si). Finally Hunt and Jackson suggest that where both phases are of high ?S regular faceted structures are formed. In the proposed classification2 the ratio of the entropies of fusion is taken as the quantitative parameter, i.e., ?Sa/?SB = E, where e > 1.0. An extensive literature survey has been carried out2,4 for reliable photomicrographs of eutectic structures. The structures, numbering 151, have been classified into two classes, viz., normal (i.e., lamellar and rodlike) and anomalous. The latter group has been arbitrarily subdivided into regular and irregular anomalous eutectics. Details and sources of the structures are given elsewhere.2,4 In the absence of better thermodynamic data the value of E has been determined for each of the systems making the following assumptions and approximations: l) The entropy of fusion of solid solutions may be calculated additively;5 ii) Inter metallic compounds are treated similarly but with the addition of the following term for fully ordered structures.5 -4.573(Ni lOg10 Nl + N2 log10 N2) where Nl and N2 are the atomic fractions. It might be argued that the entropies of fusion are meaningless as the eutectic phases form at some temperature (T,°K) below their true melting point (Ta or Tb). As AS, = Lo/To if La the latent heat of fusion of a, is assumed independent of temperature, the actual entropy change AS& is given by La/Te. The ratio of entropies of fusion then becomes Fig. 1 shows E' plotted against cumulative totals for normal and anomalous systems from which it is seen that normal eutectics predominate at low values of E'; e.g., some 90 pct of normal systems lie below E' = 1.5 with 65 pct of anomalous systems above this value. Fig. 1 also shows that the regular anomalous structures occur mostly at low values of c' possibly suggesting growth more akin to normal than anomalous growth. The choice of c as a classification parameter was based on nucleation concepts, it being observed that metals widely separated in Sundquist and Mondolfo's8 nucleation series were anomalous and those close together normal. The basis of these classification views is really outside the intended scope of this communication but it is hoped to explain them fully in the near future.7 Having arrived at this entropy of fusion classification by nucleation considerations it is interesting to note that Hunt and Jackson's ideas are based on growth considerations. Elsewhere one of the authors4 has shown that an even better classification of the metallographic structures reviewed is obtained from purely conventional heterogeneous nucleation concepts. It has been found that systems having a value of the parameter +b less than 0.28 are almost certainly normal and above this anomalous: V, is the gram atomic volume. This parameter is nevertheless closely related to c and c', these being an approximation of p. The parameter is also related to the other classification3,8,9 which have met with any success in the past. All these classifications are to be reviewed shortly by the authorsU7 A significant thing about the classifications based on

Jan 1, 1968

By John G. Hall

The pumping operations at the Chief mine have been unique in the respect that for many years the entire flow of water into the mine has been disposed of by pumping into natural underground "caverns" or "sinks." Volumes up to 4000 gpm have been pumped into these caverns with no apparent tendency for the caverns to fill with water, and it has never been proven that any of the water returns to the mine workings. Major pumping operations have been carried on during two periods, the first from 1918 to 1927, and the second from 1942 to the present time. Early Pumping Operations CHIEF MINE The Chief Consolidated Mining Co. was organized in 1909, and for the following nine years ore was developed and mined from the 600 level down to the 1800 level, all above the permanent water table, the 1820 level. The ore occurs as extremely-irregular pipes of replacement ore in limestone, and to date a total of approximately 2 1/2 million tons of ore has been produced. The ore was highly oxidized down to the permanent water level, and some highly-oxidized ore runs occur beneath the permanent water table. In 1916, with extensive mining being done on the 1800 level, 20 ft above the water table, there were indications that a zone of secondary enrichment might exist below the water, and in order to follow the ore runs to greater depth, it was realized that arrangements would have to be Made to pump water to the surface. The original pumping installation was designed to pump 350 gpm in five stages against a 1900-ft head to the surface. By January 1919, the water had been lowered to the 1900 level and an excellent grade of silver-lead ore was being produced from the zone of secondary enrichment. At this point. the flow of water equalled the capacity of the pumps, and pumps of an additional 1000 gpln were installed to supplement the original system. By January 1920, the water had been lowered approximately 100 ft, and the total flow of water was 1000 gpm. The water was being pumped from many different places in the mine where ore was being followed downward, with 27 different pumps operating below the 1800 level, mainly Cameron air-driven sinker pumps. These first sinking pumps were cumbersome and inefficient, and each of the larger pumps weighed approximately three tons and required tremeridous volumes of compressed air. During this period, an exploration drift being extended north and west from the Chief shaft on the 1600 level encountered a large natural cavern some 2000 ft from the shaft. This cavern in places measured 60 ft in diameter. Sometime later the disposing of water into this cavern was suggested, and as an experiment to determine the results of such an operation, an attempt was Made to fill the cavern by pumping in a small volume of water, but the water disappeared as fast as it was discharged into the cavern. It was decided then to try the cavern as a drainage outlet for the entire volume of water being pumped to the surface, and the water apparently did not return to the mine workings. Subsequently water up to 3000 gpm was pumped into this cavern. An attempt was made to determine whether any water was being returned from the cavern to the lower workings, and 100 lb of sodium fluorescein salt was added to the water discharged into the cavern, but no trace of it was found in the drainage water in the mine. Results of another such test will be described later.

Jan 1, 1950

By J. J. Taber, W. K. Meyer

Experiments on consolidated sandstones have shown that a variety of conditions are responsible for the wide range of behavior observed when oil and water are displaced simultaneously by a solvent which is miscible with both. The type of displacement is influenced most strongly by the relation between the plait point of the ternary phase diagram and the composition path followed by the system as it achieves mis-cibility. Neglecting viscosity, the displacement becomes less efficient as the distance between the composition path and plait point is increased in either direction so that a very wide spectrum of displacements is possible with different alcohol, oil and brine systems. It is possible by choosing various composition paths as miscibility is achieved to exercise some control over the actual displacement mechanism. Control over the composition path is exercised by adding oil or water in specific amounts to the alcohol prior to its injection into the porous medium. This technique apparently makes it possible to obtain piston-like displacements with systems which are usually characterized by the efficient displacement of only one of the two phases i.e., either oil or water, while the other phase lags behind. However, the simultaneous displacement of oil and water is not as efficient as other miscible displacements at a similar mobility ratio. The evidence indicates that both oil and water "leak" into the alcohol "piston" in compensating amounts. With oils of low viscosities the displacement is inefficient because the low oil saturation in the stabilized bank causes oil to be trapped by water before miscibility is attained. INTRODUCTION The simultaneous displacement of oil and water by solvents such as alcohol has been investigated in the laboratory and in the field. The laboratorv workers investigating all-liquid miscible systems have recognized the great importance that the nature of the traingular phase diagram exerts on the efficiency of the miscible displacement. Gatlin and Slobod pointed out that displacement efficiency will be improved if miscibility can be retained for long periods and that this behavior tends to be followed if the applicable phase diagrams have very low binodal curves. In their unconsolidated sand tests the alcohol was considered an ideal miscible piston displacing both oil and water phases. Later workers have pointed out that an important feature of the phase diagram as far as oil recovery is concerned is the location of the plait point. This is particularly true in the case of consolidated Sandstone.3 Since the original discussion of this effect, the importance of plait-point location has been confirmed and discussed by other author The early workers394 recognized that the position of the plait point could be altered by changing the oil or the alcohol used in an experiment. In addition, recent work has demonstrated that the plait point can be shifted by altering the nature and concentration of dissolved salts in the aqueous phase. However, the entire scope of the relationship between displacement and plait point has not been discussed previously and a technique for changing the composition path as miscibility is achieved in a particular ternary system, i.e., oil, water and solvent, has not been demonstrated. * The importance of oil viscosity on the efficiency of oil displacement by alcohol has not previously been observed, probably because the range of viscosities used in earlier experiments was too restricted. It is the object of this paper to describe the relationships between the ternary phase diagram, fluid properties and the displacement achieved. A technique for controlling the nature of the displacement is also presented. THEORY AND DEFINITIONS PISTON-LIKE THEORY The theory of alcohol flooding has been discussed in some detail.2-4 The concept was advanced that

Jan 1, 1965

By Klaus Detert, Joseph Ziebs

The progress of primary recrystallization has been studied in alloys of nickel with 0.01, 0.1, and 1.0 at. pct Ta. It can be shown that in one particular alloy containing 0.1 Ta and 0.02 02 a large grain structure (0.5 mm diameter) could he developed by primary recrystallization in thin sheetswhich were cold-worked to 97.5 pct reduction in thickness. Electron transmission microscopy demmstrates that this phenonmena was linked to the occurrence of fine precipitates at subboundaries. A discussion of how fine precipitates can exert under certain circumstances a rather distinct retarding effect on nucleation rate is included. IN a former investigation, the recrystallization and texture formation in a series of binary nickel alloys including pure nickel were studied.' Nickel alloys with some tantalum added displayed peculiarities during primary recrystallization which have been studied in more detail by the present investigation. The alloys were prepared by powder-metallurgy techniques similar to the previous investigations.' Pure metal powders were thoroughly mixed to obtain alloys with 0.01, 0.1, and 1.0 at. pct Ta, referred to as alloy A, B, or C, respectively. The metallic impurities were kept below 0.005 at. pct except iron and cobalt with slightly higher values, by using highly purified nickel and tantalum powders.' The mixtures were pressed into bars and sintered at a temperature of 1300°C. The sintering atmosphere was either high vacuum, dry hydrogen, or wet hydrogen. The sintered rectangular bars were rolled 60 pct to a thickness of 3 mm, then annealed under dry hydrogen at 1000°C and again cold-rolled to about 97.5 pct reduction of thickness. After this heavy deformation, the sheet thickness was 0.07 5 mm. In order to study primary recrystallization which occurred in measurable times between 400" and 600°C, samples of the rolled sheets were annealed in a salt bath at various temperatures for various times in order to study the progress of recrystallization by hardness measurements and microscopic examination. For optical observation, the samples were polished and etched in a solution of 65 cu cm HNO3, and 18 cu cm acetic acid in 20 to 100 cu cm H2O, depending of the particular sample and annealing treatment. For transmission electron microscopic observations, thin foils were prepared by the combined window and Bollmann techniques.273 The electrolyte consisted of 23 pct perchloric acid and 77 pct acetic acid held at 15°C. Polishing was performed at a current density of 0.4 amp per sq cm. The examination of the thin films was accomplished by a Siemens Elmiscope-type electron microscope operating with an accelerating voltage of 100 kv. RESULTS The progress of primary recrystallization was very accurately followed in these alloys in order to establish the time when recrystallization started and ended during an isothermal anneal. Table I gives some of the data obtained by these studies. The time for beginning is a mean value, where somewhat less than 5 pct of the matrix has recrys-tallized. The end of recrystallization is the time when at least 95 pct of the matrix has recrystal-lized. Listed are three different B alloys, with 0.1 Ta distinguished by different sintering atmospheres which influence the rate of primary recrystallization. By the above method, the state of the grain structure when primary recrystallization was completed, but grain coarsening had not started, could be clearly defined. The average grain diameter of the recrystallized matrix is also given in Table I. The grain size did not depend on the temperature of the isothermal anneal. With one exception the grain size of the primary recrystallized matrix was found to be rather fine

Jan 1, 1965

By J. L. Moriarty, Irving Johnson, H. M. Feder

The crystallization of ternary alloys has been studied as a coprecipitation phenomenon by measuring the substitutional carrying of trace metals with CeCd11 precipitating from liquid cadmium. The Doerner-Hoskins logarithmic distribution law has been found to apply in all cases; the coprecipitation coefficients obtained were: 1.48 (La), 1.06 (Th), 0.63 (Pr), 0.24 (Gd), 0.17 (Sm), 0.13 (U), 0.10 (Er), 0.09 (Sr). Li, Na, K, Ba, Y, Sc, and Zr were not measurably carried. It has been observed that coprecipitation was most extensive for those tracers which form with cadmium an intermetallic compmdnd isostructural with CeCd11 or closely related thereto. For these tracers an inverse relation between Solubility in cadmium and coprecipitation coefficient was noted. THE partition of impurities between the liquid and solid phases of a metal or a congruently melting intermetallic compound has been extensively studied by the methods of zone melting.' However, when the composition of a liquid alloy and of the solid phase crystallizing from it differ grossly, partition of an impurity cannot be investigated in this manner. Partition of the impurity under equilibrium conditions may be studied by investigation of the appropriate phase diagram and partition under nonequi-librium conditions may be studied by observation of coring in the solidified alloys.2 The nonequilibrium case is of special interest in metallurgical practice because many complex alloys are cooled too rapidly to permit the uniform diffusion of impurities throughout the crystallizing phase(s). With this in mind, the present study was undertaken to obtain quantitative data on the phenomenon of coprecipitation via substitutional solid solution in certain metallic systems. It will be noticed that rather than adopting metallurgical ter- minology we have chosen to refer to the partition of impurities during alloy crystallization as coprecipitation.3 In the terminology of coprecipitation the impurity is referred to as the tracer and the precipitating phase as the carrier. From aqueous solution, coprecipitation may occur as a result of adsorption on the surface of the carrier, by occlusion of the liquid within the crystals, or by solid solution of the tracer in the carrier. By slow cooling of liquid metallic solutions it is relatively easy to form large, regular crystals of the carrier; the effects of adsorption and occlusion may be minimized and the influence of substitutional solid-solution formation may be made to predominate. Previous studies of the crystallization of some intermetallic compounds4 and of copper, silver, and gold5 have shown that the coprecipitation or partition behavior of trace metals is entirely analogous to that observed in aqueous systems when ionic substitution occurs, e.g., coprecipitation of Ra++ with BaSo4. The precipitation of CeCd11 from liquid cadmium was chosen for study. This choice was influenced by the following consideration. If the relative atomic sizes of the tracer atom and the atom for which it substitutes determine the tendency to coprecipi-tate, then a lattice having a substitutional site of well-defined size would be desirable. In CeCd116 (structure type cubic BaHg117) the cerium atom is completely surrounded by a cage of twelve nearly equidistant cadmium atoms. Thus, atoms much larger or much smaller than cerium should be poorly accommodated in this cage. The choice of a lanthanide central atom allowed the testing of the size hypothesis with a number of lanthanide tracers which differed in size but not in valence or electro-negativity. EXPERIMENTAL The general technique of this investigation was to prepare a slightly undersaturated solution of carrier in liquid cadmium to which a tracer was added. The concentration of tracer was always so low that saturation over the entire temperature range investigated was never exceeded. The solution was cooled in a stepwise fashion. After each step a por-

Jan 1, 1964

By Pablo Hadzeriga

Great Salt Lake brine was subjected, in laboratory scale, to conditions simulating solar evaporation. Solid phases and the variation in composition of the liquid phase throughout the potassium salts crystallization zone are presented. These experimental results are discussed and compare with stable equilibrium data of the quinary system Cl-, SO4, Na+, K+, Mg + + and H2O. The Great Salt Lake brine in Utah is a vast reserve of several potential chemicals which has attracted much attention in the recent past. From the standpoint of an ore reserve, it may be said that the lake contains a low grade ore which has to be concentrated by evaporation of water. Because of the fairly favorable weather conditions in the Great Salt Lake area, it seems reasonable to consider solar evaporation as an economical means of concentration. From the view of phase chemistry, the main constituents of Great Salt Lake brine make it belong to the five component system: C1-, SO4=, Na+, K+, Mg++ and H2O. Because of their low concentration, the presence of other elements may be regarded as having little effect on the equilibrium of the five main components. Equilibrium data for the quinary system of our interest has been available for many years. Using these data, it is possible to calculate how the brine of the lake will concentrate when subjected to solar evaporation. In practicing solar evaporation, the weather factor and supersaturation conditions cannot be ignored, but their introduction into calculations using equilibrium diagrams is almost an impossible task. With the experimental work presented in this paper, an attempt is made to compare the stable equilibrium data of the quinary system as applied to the case of the Great Salt Lake brine and the dynamic equilibrium when that brine is subjected to continuous solar evaporation. Sea water brine belongs to this quinary system when only the major constituents are considered. The study of this system has been carried on since the last decades of the nineteenth century. The reason for these intensive studies was to give a geological explanation to the underground potash deposits. AS a result, many theories were presented trying to relate the equilibrium data with the geological formation, paragenesis and metamorphism of these deposits. Two main approaches were followed for the equilibrium studies of this quinary system. One of them, and the best known to us, was the school started in Germany in the last century by J. H. Van't Hoff. The other one, which started about the same time, was the school of N. S. Kurnakov in Russia. Both approaches were different and as a consequence generated opposite theories about the geological formation of potash deposits. The Van't Hoff school studied the quinary system of sea water under careful laboratory control, resulting in stable phase diagrams. During these studies, Van't Hoff recognized the difficulties in crystallizing some double salts.' ". . . . we unexpectedly found that some compounds whose formation in the solutions under investigation was possible at 25°C, nevertheless totally failed to put in an appearance. These, not to mention the calcium salts, were leonite, MgSO4 . K2SO4 4 H2O; kainite, KC1-MgSO4 . 3 H2O; and kieserite, MgSO4 .H2O. Even exceedingly slow crystallization with the addition of the compounds themselves did not remove the condition of supersaturation in the case of these bodies." However, those salts (leonite, kainite, and kieserite) were included in the stable phase diagram at 25°C as shown in Fig. 1 using D'Ans data.3 The figure also shows the sequence of crystallization of the Great Salt Lake brine according to the stable diagram. In later studies by the Van't Hoff school, especially by D'Ans,3 Karsten,4 Autenrieth,5,6 attention

Jan 1, 1968

By A. Gittens

In two recent papers stevens13,14 has compared the methods used by various investigators to calculate the contribution of grain boundary sliding to the total creep strain. In previous work Gittins15 has shown that the strain due to grain boundary sliding Egb may be estimated from the lateral, d, and vertical, v, displacements of a longitudinal marker line as it crosses a grain boundary. He obtained the formula, where nL is the number of grains per unit length measured on the specimens before creep, and \p and ?2 are, respectively, the angle between the grain boundary and the stress axis in the plane of the surface and normal to the surface (see Ref. 13, Fig. 1). The subscript L indicates that observations are made in a longitudinal direction parallel to the stress axis. It was generally found more convenient to measure the grain size of a specimen before creep; this was necessary to characterize the specimen. Stevens13,14 however prefers to define Egb in terms of NL, the number of grains per unit length after creep. Since nL - NL(l + €1) where cr is the total strain, Eq. 1 11 becomes The quantities d and ?, can be determined directly from the marker lines, while v and ?2 are found by sectioning normal to the surface. In the interior when symmetry considerations suggest that v = d and ?1 - ?2 the strain due to grain boundary sliding may be estimated from a suitable internal longitudinal marker line by the equation To determine cgb at the surface using either Eqs. |1| or [2] requires a large number of observations, particularly if many estimates of cgb are required, and also it is impossible to determine the variation of cgh with creep strain on the same specimen because the specimen must be sectioned to determine ?2. Gittins15 found that for a number of materials the sliding vectors d and v were related to each other so that provided one was measured the other could be estimated. Also during creep, grain boundaries tend to migrate to meet the surface at angles close to 90 deg,16 so that the term nL(v/tan?2)L was comparatively small when 4,- 90 deg. By determining the frequency distributions of ?2 on control specimens this term can be estimated without sectioning the creep specimen. In this way an empirical equation based on Eqs. |1| and |2| enables cgb to be determined with reasonable accuracy with fewer observations of the angles and sliding components. For B brass, the relationship between the sliding vectors indicated that ceb = 0.85nLVL . In Table I of Stevens paper this equation is stated to be incorrectly derived. However it should be pointed out that this is an empirical equation derived from Eq. |1|, based on a knowledge of the relationship between d and v and on the frequency distributions of ?1 and ?2 for a particular material under given experimental conditions. It was used only when these relationships were known, to reduce the number of observations required to estimate cb. Therefore the comment "derivation incorrect does not apply. Gittins did in fact make use of Eqs. |1| and 131. Author's Reply I must apologize to Dr. Gittins for the inadvertent misrepresentation of his work. Other slight errors in Table I have been pointed out to me. Langdon (No. 12 Table I), Gifkins (No. 7, Table I), and possibly other workers did not measure the number of boundaries per unit length after creep (A') but the corresponding quantity bejove creep (n). As Dr. Gittins points out, a completely analogous set of equations to those given in the paper can be derived by the same method: using the number of boundaries per unit length before instead of after creep. The basic equations corresponding to Eqs. [2] and 171 of the paper, from which all the others can be derived, are

Jan 1, 1968

By L. Burris, C. H. Stockman, I. G. Dillion

Zone melting is a purification process in which separation of impurities is effected by slowly moving a narrow melted zone through a bar of solid material. Equations are presented which 1—predict the concentration profiles after successive passes of the zone and 2—give the limiting distribution of solute along the bar after an infinite number of passes of the molten zone. These equations take into account the effect of normal freezing in the last zone length. These equations were solved numerically using an electronic computer, and the concentration profiles are presented as a function of the solute distribution coefficient and the ratio of molten zone to total bar length. ZONE MELTING is a purification process in which separation of impurities is effected by slowly moving a narrow melted zone through a bar of solid material. Separation occurs if the solid freezing from a solution at any moment has a composition different from that of the liquid from which it is freezing. Solutes which lower the freezing point will tend to remain in solution and will be moved in the direction of zone travel. Solutes which raise the freezing point will preferentially freeze out, and they will be moved in a direction opposite to that of zone travel. After repeated unidirectional passage of a molten zone through a bar of material, impurities will be concentrated in the ends of the bar. A limiting concentration gradient will be reached in the bar after a large number of passes of the molten zone. Both the magnitude of this limiting concentration gradient and the rate of movement of solute in any pass depend on the zone length and on the equilibrium concentrations of solute in the liquid and solid phases. The basic theory of zone melting was thoroughly covered by Pfann.1 The purpose of this paper is to extend the theory of zone melting by presenting 1—equations which relate the concentration profiles in successive passes of the zone and 2—equations which give the limiting distribution of solute along the bar after an infinite number of passes of the molten zone. These equations take into account the effect of normal freezing* in the last zone length concentration of a given solute in the solid to the concentration of the same solute in the liquid. Analytical methods for determining the concentration profiles after any pass have been developed by Lord2 and Reiss. The method developed by Lord does not take into account normal freezing in the last zone and also becomes very complicated when the number of passes becomes large. The method developed by Reiss enables the direct calculation of the complete concentration profiles (including that in the last zone) after any pass for values of k close to unity (0.9 < k < 1.1). Both methods are useful and provide valuable insight into the zone melting process. The equations given in this article are based on the same assumptions used by Pfann. These are: 1—uniform composition in the liquid, 2—negligible diffusion in the solid, 3—constant distribution coefficient, 4—constant length of the liquid zone (except at the ends), and 5—solubility of the solute in the melt is not exceeded at any point. Equations for Concentrations in Successive Passes For a melted zone of length 1 at the left-hand end of a bar of uniform unit cross-section, the amount of solute in the liquid is J0Cn-1(x) dx [1] where the subscript is the number of the pass. Cn-1(x) represents the concentration (expressed as weight per unit volume) obtained in the previous

Jan 1, 1956

By M. R. J. Wyllie, H. W. Patnode

In the quantitative interpretation of electric log data it is essential to know the formation factor, defined as the ratio of the resistivity of the formation 100 per cent saturated with brine to the resistivity of the brine.&apos; This concept implies that the formation factor for any given rock is constant irrespective of the resistivity of the brine or solution with which it is saturated. It is apparent that this can be true only when the rock matrix is wholly non-conductive. As reservoir rocks may contain clay or shale having an appreciable conductance one may expect these constituents to influence the determination of true formation factor. The formation factor may be theoretically derived from log data by dividing the resistivity of the mud filtrate invaded formation by the resistivity of the mud filtrate. Such a method of calculntion has implicit in it the assumption that the fluids in the invaded formation are wholly replaced by mud filtrate. that the solid matrix is non-conducting, and that the formation factor of the invaded zone is identical with the formation factor of the uninvaded zone. Since the validity of the first two basic assumptions regarding the insulating character of porous matrices and the displacement of the interstitial fluids in porous media do not appear to be fully established, the formation factors of a number of cores and clay slurries have been measured under conditions suitable for verifying the assumptions and the extent of any errors inherent in them. CORE STUDIES Apparatus and Procedure Resistivity measurements were made on short cylindrical permeable cores which had been previously mounted in lucite.2 The cores were clamped by means of insulated bolts between brass plates about 5 cm square and 1 cm in thickness. The plates were provided with a ring seal to prevent leakage between them and the core mount, and had a centrally located orifice so that fluids could be flowed through the mounted core. The brass mounting plates were used as electrodes and any current flow between them necessarily passed only through the constituents of the intervening core. Fluid saturation of the core, and fluid flow through the core, were effected by a centrifugal pump connected to the outflow end of the core, thus providing a differential pressure across the core of approximately one atmosphere. The procedure used was first to evacuate each core by means of the pump and then to saturate the core with an aqueous sodium chloride solution of known resistivity. The resistance of the core was measured and flow of the solution through the core continued for sufficient time to establish a constant resistance reading. After equilibrium had been reached with one solution, another sodium chloride solution of a different resistivity was flowed through the core and the flow continued until a constant resistance reading was again obtained. This procedure was continued until resistance measurements of the core had been made with a series of solutions within a salinity range from 1.0N NaCl to distilled water. The data obtained enabled the apparent formation factor of given porous rocks when saturated with fluids of different resistivity to be computed and also provided information regarding the ability of an invading fluid to displace interstitial water from a core. Results Berea Sandstone — The first core examined was a specimen of Berea sandstone." This core had been cut essentially parallel to the bedding for permeability studies. The core was clamped in the core holder and saturated with 1.0N NaCl solution and the resistance determined. It was subsequently flushed and the resistance measured when saturated in turn with 0.1N, 0.025N, 0.01N, .005N sodium chloride solutions and finally with distilled water. In flushing the core with each

Jan 1, 1950

By J. O. Brittain, E. Lautenschlager

Yield point measurements were made on polycrys-talline specimens of a iron deformed various amounts at 770, 1950, or 273°K, aged, and redeformed at 77: 195°, or 273°K. When the prestrain temperature is lower than that subsequent to aging, the yield point return is suppressed relative to specimens deformed at the same temperature both before and after aging. The degree of suppression is a function of both the magnitude of prestrain at the lower temperature and the difference between the deformation temperatures before and after aging. It is proposed that the suppression is caused by the formation of dislocation barriers during the initial low-temperature deformation. The aging kinetics of polycrystalline a-iron wires deformed at 77o or 303oK were studied at 44 " and 64 "C using internal friction measurements. No change was noted as a result of the difference in prestrain temperature. IN an earlier investigation1 it was observed that the return of the yield point was suppressed when a poly-crystalline specimen of a iron was plastically deformed at 77oK, aged at room temperature, and subsequently tested at room temperature. In order to investigate this observation, a research program was undertaken to measure the effect of prestrain temperature on the strain aging of a iron.&apos; According to Thomas and Leak,3 return of the yield point and internal friction are similar measures of strain aging since both are dependent on solute atoms leaving the solid solution. The decrease in internal friction as a function of aging time is indicative of the removal of solute atoms from the solution, while the return of the yield point as a function of aging time gives the rate of arrival of solute atoms at the dislocations. According to Cottrell and Bilby,4 the solute atoms create an "atmosphere" which locks the dislocations and causes a return of the yield point upon restraining. In this experiment, both methods of measuring strain aging were employed in order to have two independent measurements of the strain aging phenomenon. In a previous study,5 the strain aging kinetics of a similar lot of a iron were measured for specimens which were prestrained and tested at the same temperature over the range of 77" to 273°K. Although some differences were found in the correlation of the two methods employed for measuring strain aging, the kinetics of strain aging were unaffected by the change of prestrain and test temperature. In another study,6 dislocation line tension and the density of the solute atmosphere were found to contribute to the return of the yield point when ingot iron specimens were prestrained at 303°K and tested at a lower temperature, namely 77°K. This paper is concerned with changes in strain aging kinetics produced by pre-straining at 77" to 303°K and testing at either the prestrain temperature or at higher temperatures up to 273°K. EXPERIMENTAL PROCEDURE The polycrystalline a iron specimens used in both the tensile and internal friction tests were 0.033-in. diam ingot iron wires, wet hydrogen treated and re-carburized. A more detailed procedure for processing a 3/16 in. ingot iron rod (nominal wt. pet composition: Fe - 0.013; C - 0.0125; Mn - 0.023; Si -0.005; P - 0.0023; S - 0.11; Cu - 0.004 N) to 0.033 in. diam specimens containing approximately 0.020 wt pet C and no detectable nitrogen with a grain diameter of 0.033 to 0.040 mm is given elsewhere.5 Immediately prior to testing, specimens were solution treated for 15 min at 700o ± 2°C in a static vacuum of less than 50 ?. of Hg and rapidly quenched into an ice-water bath while still enclosed in a 3 mm quartz vial. The internal friction specimens had 8-in. gage lengths and the tensile specimens had 2-in. gage lengths. The wires used as tensile specimens were first run as internal

Jan 1, 1962

By L. S. Darken, H. A. Wriedt

An investigation has been made of nitrogen absorption by the lattice defects in a low-carbon steel afte~ Plastic deformation. Specimens in which defects were distributed by various combinations of cold rolling and vacuum heat treatment were equilibrated with controlled ammonia + hydrogen mixtures at 300° to 450°C to study the relation between the solubility and the thermodynamic activity of nitrogen. The results are interpreted in terms of nitrogen being absorbed upon two different types of defect sites, called A sites and B sites. The concentrations of these sites are calculated for specimens with various histories of straining and annealing; and the thermodynamic relations descrihing the equilibrium distribution of nitrogen among the A, B, and ordinary interstitial sites in ferrite are evaluated as follows: THE generalization that the solute atoms in a solid solution interact with imperfections in the solvent lattice is now universally recognized and the subject of an extensive literature. All types of imperfections—point, line, or surface—exhibit these interactions, but the relative importance of each type depends upon its concentration, the strength of its interaction with the particular solute atom, and whether the interaction is attractive or repulsive. The strength and sign of the interaction vary with the vectorial displacement of the solute atom from the imperfection in the crystal. The concentration of lattice defects increases during plastic deformation and the effect of the imperfections upon the solute atoms is usually perceived only after cold work or irradiation. Two experimental methods have previously been used to observe the effects of interaction between nitrogen atoms and lattice imperfections in a iron. Internal friction, which was first used for this purpose by Dijkstra,1 has been applied more often than the electrical-resistivity method of Cottrell and Churchman.2 The method of the present investigation has not been used before.* Uniform specimens of steel with a metastable defect structure are produced by cold rolling and subsequent heat treatment below the recrystallization temperature. Variations in the defect structure are obtained by varying the amount of strain and the heat-treatment conditions. The specimens are then saturated with nitrogen at moderately elevated temperatures in ammonia + hydrogen mixtures having various thermodynamic activities of nitrogen below saturation with Fe4N. From the nitrogen contents of analyzed samples, the metastable equilibria of interactions between the defects and nitrogen can be established and the concentration of sites for nitrogen atoms in the lattice defects calculated. A steel was used in the present experiments instead of purified iron or iron alloys because the amount of nitrogen absorbed at defects is larger in the steel over a much wider range of experimental conditions and could therefore be measured more accurately. The effect of alloy composition is discussed in Part II of this paper; it is sufficient to say here that impurities and phase heterogeneity in ferritic alloys favor the generation and/or retention of larger defect concentrations during plastic deformation and annealing at 300° to 500°C. MATERIALS AND EXPERIMENTAL PROCEDURE The experimental specimens were prepared from stock of a low-carbon steel, which will be called "X steel" here, with the composition: C, 0.12 pct; Mn, 0.43 pct; P, 0.010 pct; S, 0.034 pct; Si, 0.007 pct; As, 0.086 pct; Cu, 0.022 pct; Ni, 0.009 pct; Cr, 0.003 pct. The stock was all from a single ingot made by the basic open-hearth process, which was hot-rolled to strip 0.080 in. thick.. This hot-rolled strip, given the code designation "XH", was the starting material for preparing a series of cold-rolled sheets. The latter are designated "XC27",

Jan 1, 1965

By W. V. Green

The creep-rupture behavior of commercial powder-metallurgy tungsten rod is reported for temperatures of 2250°, 2500°, 2700°, and 2800°C, stresses up to 7000 psi, and times up to 4 hr. The temperature dependence and the stress dependence of creep rates are evaluated. Rupture time is fitted to the Larson-Miller parameter and compared to lower-temperature and higher stress results reported by pugh.4 THE tensile strength of single-crystal wire1 and of polycrystalline wire2 have been measured by Goucher. More recently Bechtold3 and pugh4 have studied the mechanical properties of tungsten rod between room temperature and 1200°C. Since tungsten has the highest melting point of all metals, it should retain useful strength to very high temperatures, and should find many applications in missiles, rockets, and so forth. Therefore, the mechanical usefulness of polycrystalline tungsten rod was investigated in the 2250" to 2800°C temperature range by short-time creep-rupture tests. MATERIALS TESTED Creep specimens were machined from 1/2-in. round commercial tungsten rod which was manufactured from powder by cold pressing, sintering, and hot swaging. Purity reported by the vendor was 99.95 pct W. Spectrochemical analyses showed "traces" (less than 0.01 pct) of silicon, iron, and titanium, and "faint traces" (less than 0.001 pct) of calcium and manganese. Carbon content was 30 ppm and nitrogen content was 70 to 110 ppm. No oxygen analysis was obtained. Room-temperature stress-strain tests on the as-received tungsten gave badly scattered results. Threaded end specimens broke at the threads. When these specimens were retested using shoulder grips, the specimens broke at the shoulders. No fractures occurred in the gage lengths. The fracture stresses were: 127,000; 137,000; and 154,000 psi for three specimens which fractured at the shoulders. The microstructure of the as-received tungsten rod is shown in Fig. 1. SPECIMENS USED Fig. 2 shows the two types of creep specimens tested. The threaded-end type specimen, which was used in most tests, had a 3/4-in. gage length and 1/4-in. gage diam. The cone-end specimen, which was used primarily to investigate the effects of a smaller gage diameter and the practicability of the modified grips and specimen, had a 3/4-in. gage length and a 1/8-in. gage diam. In both cases, the gage length terminated at 0.005-in. high shoulders which served as fiducial marks for optical strain measurements. Surface finish on the gage length was 16 µ-in. rms. Creep rates and rupture times were not significantly different for the two types of specimens and were not affected by electropolishing the specimens before testing. TESTING PROCEDURE The constant-load creep-testing machine used has been described by Smith, Olson, and Brown.5 In it the specimen is held vertically by grips of similar material and on the axis of a cylindrical tungsten or tantalum heater-tube. The specimen is loaded by means of the specimen grips, a pull rod (which is sealed to the chamber lid by a vacuum-tight bellows), a mechanical lever system, and hanging dead weights. Both specimen and grips are heated by radiation from the heater-tube, which is heated by its own electrical resistance. An optical system views the incandescent specimen through slits in the radiation shielding and heater-tube. The optical system projects an enlarged image of the specimen gage length on a ground-glass screen. Periodic gage-length measurements are made on this image with two cathetometers. Thorium oxide coatings on the specimen threads or conical end surfaces completely prevented diffusion-welding of specimens to grips, and were used

Jan 1, 1960

By M. Ahlers, Horace Pops

Slip features and formation of stress-induced mm-tensite have been studied in single crystals of fl phase Cu-Zn alloys in the composition range between 43 and 48 at. pct Zn. A phenomenological analysis is presented in which the glide planes and directions that me associated with the inhomogeneous shear in the martensite phase are predicted from the experimentally determined habit planes and macroscopic shear directions. It is shozun that the resolved shew component of the applied stress on the secondary shew system is a very important factor governing the formation of stress-induced martensite. A comparison is made between martensitic phases that form spontaneously upon cooling and those phases that arise from deformation. THE ß phase in Cu-Zn alloys is unstable and may transform spontaneously by a martensitic transformation1-3 during cooling. Thermoelastic and burst-type spontaneous martensitic phases have been observed3 in the composition range between 38 and 42 at. pct Zn. Martensite forms in alloys containing more than 42 at. pct Zn only by cold work below a temperature denoted as Md.4 Thus, alloys near the equiatomic composition can undergo a stress-induced transformation upon cold work at room temperature. Although data have been obtained concerning the effect of temperature5 and composition4 on the deformation of polycrystalline Cu-Zn alloys, it is rather surprising that few detailed studies6-8 on 6 brass have been conducted for single crystals where the stresses acting on crystallographic planes can be determined. A thorough understanding of transformation mechanisms requires knowledge of the role of stress and plastic strain. With this objective in mind we have made an investigation of the stress-induced rnartensite transformation in single crystals of 0 brass in the composition range between 43 and 48 at. pct Zn. Particular emphasis is given to the morphology and crystallographic features and a comparison is made between stress-induced and spontaneous martensite. I) EXPERIMENTAL PROCEDURE Master alloys having nominal zinc contents of 43, 45, and 48 at. pct were prepared from spectroscop-ically pure metals (>99.999 pct) in sealed quartz tubes under a partial pressure of helium. These alloys were cut into several pieces and converted into small rods by melting in quartz tubes. Subsequently, single crystals were grown from the rods by the Bridgman technique and homogenized for 3 hr at 820°C. Lattice parameter measurements were made from stress- relieved filings. Compositions, as determined from the reported variation of lattice parameter with atomic percent zinc,&apos; indicated a zinc loss of less than approximately 0.25 at. pct. The single crystals were electroturned using a phosphoric acid electrolyte into tensile specimens each with a center section of 1 mm diam and 1 cm gage length. Cylindrical specimens approximately 3 mm in diam and 10 mm long and spark-machined rectangular blocks approximately 3 by 3 by 10 mm were prepared for compression studies. These blocks were stress-relieved and electropolished to eliminate any stresses introduced by spark machining. Crystal orientations were determined to an accuracy of ±2 deg primarily from back-reflection Laue photographs, and sometimes by the rotating single-crystal method. Figs. 1 to 3 show the location of the tension and compression axes within unit stereographic triangles. All specimens were deformed plastically at room temperature using an Instron machine, the crosshead speed being 0.02 or 0.05 cm per min. A hard steel hemisphere was placed at the upper end of the compression samples to achieve uniaxial loading. Trace normals of slip markings and habit planes of strain-induced martensite were determined to an accuracy of ±4 deg from the rectangular blocks by a two-surface analysis and from cylindrical specimens by measuring the angle between the tensile (or compression) axis and the tangent to the surface trace for different rotations about the deformation axis. Slip

Jan 1, 1969

By Arthur H. Wilson

INTRODUCTION The stratified deposits of the coal measures are strongly laminated and exhibit a high degree of anisotropy. The stronger rocks in the sequence possess joints and shear planes, the coal itself has a closely spaced cleat pattern, and the softer rocks are capable of flowage. To calculate the general distribution of stress in such a medium by analytical means based on the theories of elasticity is difficult, and has not yet been satisfactorily achieved. Nevertheless, some estimate is required to aid the design of supporting pillars and permit the correct siting of main roadways. This is important when laying out districts for longwall extraction, particularly if the results of previous experience are not available in the area to be exploited. Although elastic stress analyses, including finite element analysis, may not be very appropriate, the basic equilibrium conditions will apply. These, combined with logical deductions based on observed fact, can be used to give appropriate solutions which experience has shown to be capable of giving acceptable results in the design of pillars associated with longwall workings. Constancy of cover load is one of the fundamental equilibrium conditions. When mining below ground, the cover is not removed, therefore the average vertical load after mining must equal the initial cover load. If mining produces a local decrease in load in any particular area, it must be balanced by an equivalent increase in load nearby. Bending moments in the strata above are capable of compensating for the turning moments which this redistribution of load may produce, and the existence of a "rear abutment", much discussed between 1950 and 1960, is not an essential equilibrium requirement. If the distribution of vertical stress is plotted three-dimensionally, then the load on any small element of area 6S will be o.6S, where o is the average stress. This also represents the volume of the "stress envelope" which lies above 6s. By making the individual areas very small and summing together, it is possible to show that the total load on any particular area can be found by calculating the volume of the stress distribution envelope which lies above it. A similar consideration applies to the stress distribution in a two-dimensional situation. If the stress distribution curve is drawn, then the area below any portion will represent the load per unit depth of section. For example, Figure 1 represents the possible stress distribution across a caved waste and associated ribside. If AW is the load deficiency [ ] These loads must balance, therefore area AW = AS, and if one can be deduced, the other will then be known. Attempts have been made to measure stress distributions in the ribside and in the waste. Stress meters in soft rocks such as the coal measures have not proved very effective. Those based on the measurement of strain have experienced difficulty in strata capable of flowage, high modulus stress meters suffer from cross sensitivity, and overcoming to obtain absolute stress is almost impossible. Load measuring instruments left in the waste have also proved largely ineffective, mainly as a result of the difficulties of protecting the leads in a caving environment. Certain logical deductions, however, can be made about the conditions in the waste. CONDITIONS IN THE WASTE Over the area of the waste the surface trough which develops is remarkably consistent, despite the geological variations which occur from coalfield to coalfield. In Britain, an accuracy of 210% is claimed for subsidence prediction based on width of extraction and depth of cover (National Coal Board, 1975). The subsidence in the USA is more variable (Adamek and Jeran, 1981), particularly with respect to the point of inflexion, but the same general pattern occurs. At the shallower depths the somewhat harder roof strata involved in the caving in the USA may suffer less compaction, giving reduced values of subsidence. Subsidence at intermediate horizons below the surface

Jan 1, 1982

By R. M. Goldhoff, J. W. Spretnak

IN connection with establishing the mechanism by which boron enhances the hardenability of heat treatable steels, this research work has been undertaken. Spretnak and Speiser1,2 indicated the need for studying high temperature effects in boron steels and, chief among these, the possible nonhomogeneous distribution of boron in iron grains. The occurrence of this effect may be of significance in explaining the mechanism of enhancement of hardenability by boron as well as in explaining some of the observed characteristics of boron steel behavior. Primarily, then, this work was undertaken a) to establish whether or not boron is adsorbed in ? solid solutions, b) to determine the nature of any such adsorption effects, i.e., whether it is positive or negative as well as its temperature coefficient, and c) to establish the magnitude of any such effect. While past literature on this specific subject is rather sparse, a number of authors have postulated that the observed temporary loss of hardenability as well as the occurrence of the characteristic boron constituent in these steels is due to the nonuniform distribution of boron in austenite.3-0 Experimental Techniques Attempts to obtain direct evidence of grain boundary effects in boron steels have been unsuccessful, both by using the technique of Chalmers&apos; in which small concentrations of solute made radioactive by irradiation are detected by autoradiographic means, and by means of studying cleavage fractures with the electron microscope.8 The present study attempts the design of experiments to measure the effect of boron on properties which might be subject to variation because of grain boundary effects. A study of the theoretical aspects of grain growth indicates that such studies may be pertinent. Another valid type of experiment involves the measurement in polycrystalline materials of properties whose values depend on the bulk of the grain, being negligibly influenced by grain boundaries. One such property is the lattice parameter measured by X-ray diffraction. The value of the lattice parameter reflects the average conditions which obtain over a relatively large number of atoms in a fixed crystalline lattice, and is therefore mainly dependent on bulk grain material. Measurement of lattice parameter variation over the austenite range of temperature in iron with and without boron should give evidence of concentration variations in the bulk grain material if adsorption effects are significant. Finally, use is made of the Grange and Garvey3 metallographic test for boron in which a characteristic boron constituent can be developed in austenite grain boundaries. The very occurrence of this phenomenon suggests grain boundary enrichment of boron and a study of the occurrence of this constituent and its varation with austenitizing temperature should help to detect grain boundary effects. Materials One lb ingots of pure iron, an Fe-B alloy, an Fe-C alloy, and an Fe-C-B alloy were prepared. The melting technique included an initial melting of electrolytic iron under vacuum followed by remelting under hydrogen and, finally, by melting under argon with the addition of alloying elements. Each ingot was surface machined between meltings. A quantitative spectrographic analysis was obtained on the pure ingot and, in addition, oxygen, nitrogen, and carbon were determined by the vacuum fusion method. The results of these analyses are given in Table I. Semiquantitative spectrographic analyses were performed on the alloy ingots, while quantitative determinations of the alloy additions were obtained. These analyses are shown in Table 11. The cast ingots were subjected to the following sequence of fabrication treatments: forging to 5/8 in. bars, swaging to 0.300 in. rounds, drawing to 0.025 and 0.010 in. wire. Portions of each ingot were retained in the swaged condition for experimental use. For certain portions of the work typical commercial steels were used. These materials were supplied by the Republic Steel Corp., Canton, Ohio, and correspond to AISI 8640 and AISI 86B40 compositions.

Jan 1, 1958

By T. F. Mueller, P. A. Witherspoon

For the case of an infinite radial system operating at constant terminal rate, the reservoir engineer often uses the "point source" solution of the diffiusivity equation to study pressure interference effects. At early times and at short distances from the inner boundary these solutions are invalid. The amount of error is often not precisely defined because of mathematical difficulties. Work presented here shows that, if dimensionless time is defined appropriately, previous solutions of the pressure equation can be displayed .as a family of curves on one chart. These solutions include the point source solution (referred to in the field of hydrology as the Theis solution) and other solutions obtained with digital computer methods. With these curves, an exact evaluation of the pressure drop within a reservoir or an aquifer can be made by the engineer. Examples of field problem solutions are presented. In most reservoirs the error involved when the Theis solution is employed is often negligible; whereas, in the calculation of interference effects in an aquifer, a substantial error can occur through such an approach. INTRODUCTION Flow equations are used in petroleum engineering to study the behavior of individual wells and reservoirs. In the case of wells, the pressure response at the wellbore face is the major point of interest; whereas, in the case of reservoirs, the pressure response at the interface of the aquifer boundary is sought. To aid in such studies, the flow equations have been solved in terms of the behavior at these two inner boundaries. Only limited work has been published in regard to the pressure conditions away from these points, i.e., within the reservoir or aquifer. Theis and Mortada are among the few who have reported on this problem. The Theis approach employs the exponential integral and is valid for pressure conditions that occur some distance away from the flow disturbance. It is derived from the concept of a point source, as opposed to a flow across some finite area. The Mortada results, on the other hand, are valid at all points within the reservoir or aquifer. They are presented in terms of dimension less ratios of the radius where the pressure is desired to the radius where the flow rate is measured. Their main use, in the past, has been in aquifer studies. The published results are presented in the form of graphs that are limited to a maximum radius ratio of 64. These graphical results are cumbersome to interpolate at non-integral radius ratios, so that one may be forced to utilize a rather involved analytical expression presented by Mortada. BASIC EQUATIONS The solutions of Mortada and Theis are both based on the diffusivity equation as applied to the case of an infinite radial system subject to a constant terminal rate. The equation is obtained by combining the material balance equation with Darcy&apos;s flow equation. The assumptions implicit in the use of this equation are as follows: (1) a single fluid is present that occupies the entire pore volume; (2) the reservoir is horizontal, homogeneous, uniform in thickness, and of infinite radial extent; (3) compressibility and viscosity of the fluid remain constant at all pressures; and (4) fluid density obeys the equation Using the diffusivity equation in situations where the above conditions do not hold will result in errors. These errors (not discussed here) but only the errors which arise in the solution of the equation itself. The diffusivity equation for the homogeneous reservoir conditions cited above can be written in cylindrical coordinates, as To obtain a dimensionless equation, so a single solution can be used for applications of different porosity, permeability and fluid properties, the following transformations are usually made: After these transformations are made, Eq. 3 .can be written in dimensionless form, as MORTADA SOLUTION One solution of Eq. 6 has been given by Mortada.

Jan 1, 1966

By M. Prats

The effect of a sand-filled vertical fracture of limited radial extent and finite capacity (fracture capacity is the product of the permeability and width of the fracture) on the flow behavior of a cylindrical reservoir producing an incompressible fluid through a centrally located well has been investigated mathematically. The shape of the lines of equal pressure near the fracture is essentially independent of the size of the reservoir, provided that the field radius is of the order of the fracture length or larger. For reasonable values of production rates and of fluid, reservoir and fracture properties, the total pressure drop between the end of the fracture and the well is gelrerally negligible compared with the pressure drop in the reservoir. With regard to production response, the effect of vertical fractures can be represented by the production response of an equivalent or effective well radius. For a high-capacity fracture, the effective well radius is a quarter of the total fracture length, decreasing with the fracture capacity. When invasion effects are simulated by decreasing the width of the damaged zone with distance from the well, the effect of formation damage around a fracture on the production response is not so seriolls as indicated by the literature. This suggests that frac fluids with a conventional filter-loss response are better than high-spurt-loss frac fluids, provided the effective permeability of the damaged zone is the same. INTRODUCTION This paper considers the effect of the fracture capacity, as well as the formation damage which can result from fracture treatments, on the productive capacity of vertically fractured wells. Other publications, notably those of van Poollen,l,2 consider these same effects. In addition to providing more general results for vertical fractures than are available from the liter- ature, the present paper gives the equivalent well radius of fractures having different lengths and capacities and, also, includes pressure distributions in and around the fractures. The effect of a damaged zone around a fracture on the production response was not found to be so great as that reported by van Poollen. 2 This difference probably stems from the fact that we consider a damaged zone which is widest (but is still small) near the well and thins out toward the extremities of the fracture, whereas van Poollen considers a damaged zone having a uniform width for the entire fracture length. Simple, but adequate, equations which describe the effect of these variables on production response are presented (in Appendixes A and B). Thus, results can easily be extended to values of the variables not specifically considered here. IDEALIZATION AND DESCRIPTION OF THE FRACTURED SYSTEM It is assumed that a horizontal oil-producing layer of constant thickness and of uniform porosity and permeability is bounded above and below by impermeable strata. The liquid is incompressible. It is assumed that, after the well is fractured, the cylindrical outer "boundary" is at a uniform potential, provided that it is not too near the fracture. The fracture system is represented by a single, plane, vertical fracture of limited radial extent, bounded by the impermeable matrix above and below the producing layer (reservoir).3 Fig. 1 indicates the general three-dimensional geometry of the fractured reservoir described above. If gravity effects are neglected, the flaw behavior in the reservoir will be independent of the vertical position ill the oil sand. This means that the flow behavior in the fractured reservoir can be described by the two-dimensional flow behavior shown in the horizontal cross section of the reservoir, Fig. 2. PARAMETERS The effect of a fracture on the pressure distribution in the fracture itself and in an otherwise undamaged (no skin) reservoir indicates that the pressure distribution is a function of three parameters.

By Charles E. Mann, Douglas J. Lootens, Rufino B. Bomasang

INTRODUCTION The economy of the Philippines, like that of many other developing nations, has been severely strained by the rapid escalation of imported crude oil prices over the past decade. In an effort to reduce its dependence on imported crude, significantly increased emphasis, both in the public and private sectors, is being placed on the rapid development of Philippine coal resources. This paper is designed to present an overview of the Philippine coal industry, its historical development, present status, and future potential. COAL RESOURCES AND RESERVES In 1977, a British consulting firm collated all available data on Philippine coal resources and estimated the total coals potential at about 0.9 to 1.4 billion metric tons (mt). In 1981, using more recent data the Philippine Bureau of Energy Development (BED) estimated the potential at 1.7 billion mt. The BED estimated the proven coal reserves at 257 million mt. These coal resources are scattered throughout the major islands comprising the Philippine archipelago. However, about 52 percent of the currently proven reserves are located in the tiny island of Semirara between the islands of Mindoro and Panay. About 25 percent is located in the Cagayan Valley in northeastern Luzon while 13 percent is located in Mindanao. The remaining 10 percent is scattered in the islands of Cebu, Samar, Mindoro, Negros, Polillo, Batan, and Catanduanes. The bulk of Philippine coal is sub-bituminous in rank. Some bituminous coal with coking properties has been reported in Malangas, Zamboanga del Sur in western Mindanao, and in Catanduanes, but the exact tonnages have not been accurately established. It is believed, however, to be at most 5 million tons. A little semi-anthracite coal exists in Malangas with the remaining unmined reserves estimated at a little more than one million tons. Table 1 shows the quality ranges of Philippine coal deposits.. Figure 1 shows the distribution of known Philippine coal resources PHILIPPINE COAL PRODUCTION In 1981, total coal production was 330,720 mt. Though modest by world standards, this was the highest production figure in Philippine coal mining history. This came from 25 mines, mainly small scale underground mines of capacities varying from 10 to 100 tons per day (tpd). The actual production capability was about 800,000 mt in 1981, but the producers had to scale down production to the level of local demand. Even then, some 50,000 tons of coal remained unused by year-end. The 1981 production came from the following areas: Cebu 224,751 mt (68%) Zamboanga del Sur 57,095 mt (17%) Semirara 13,222 mt (4%) Polillo 17,132 mt (5%) Cebu has always been the main source of Philippine coal production. Coal mining has been carried out in Cebu for the last 50 years, although it has been characterized by a very crude mining technology. Due to the limited size of coal deposits and complicated geology, however, Cebu&apos;s share of the total national production is expected to go down to only 10 percent by 1985. Except for Semirara&apos;s open-pit production, the balance of Philippine coal production comes from underground mines. These underground mines generally use room-and-pillar methods, except for those at Zamboanga del Sur and Polillo which use modified longwall methods patterned after those used in Taiwan. Due to the much higher mining recovery, most underground mines are now planning on switching to longwall. At Semirara, an open pit mine is now being developed for a one million ton per year capacity by 1984, and will then be the largest Philippine coal mine. While present earthmoving work is done with bulldozers, loaders, and trucks, the main coal extraction equipment will be bucket wheel excavators imported from Austria. PRESENT COAL UTILIZATION Existing coal production is sold mainly to coal- fired cement plants and power plants. The balance is sold to a sugar central, a coconut oil mill, a food processing plant, and two plants producing carbide and ferro-silicon alloys.

Jan 1, 1982

By P. Hazebroek, C. S. Mathews, E. E. Allen, H. C. Lefkovits

A rigorous study was made of the behavior of reservoirs composed of horizontal layers, unconnected except at the well and filled with a compressible fluid. The report is presented in two parts. Part I deals with the practical implications of the results obtained; Part 11 contains the mathematical derivation of the solution. The main portion of Part I is concerned with application of the results to reservoirs composed of two layers. Relative rates of depletion of the layers are studied, and it is shown that differential depletion between layers exists in the transient stage of the reservoir, which generally is of an order of magnitude longer than the transient stage in a single-layer reservoir. In the transient stage, the more permeable layer is depleted faster than the less permeable. However, the reservoir approaches a steady-state situation where both layers contribute equally to production. Theoretical build-up curves from such reservoirs are presented, and the influence of a skin effect and flow into casing and tubing are also evaluated. The build-up curves are predicted to have, first, the familiar logarithmic straight-line section and, subsequently, a rising and flattened section. The results found are compared with an extension of an approximate theory, and it is found that this simplified theory is applicable over most of the life of the reservoir. Extensions to reservoirs containing more than two layers are indicated, and examples of the performance and build-up curves for three-layered reservoirs are presented. The theoretical results on build-up curves were applied to field examples of such reservoirs, and satisfactory results were obtained. From the build-up curve, the pemeability-thickness product, the wellbore damage and the static pressure are obtained. However, it apparently is not possible to determine the properties of the individual layers from the combined build-up curve. A large number of reservoirs have been found where the producing formation is composed of two or more layers of differing physical characteristics, such as permeability and porosity, and where also the thicknesses of the layers differ. Since in general a differential depletion should exist between the layers of such reservoirs, it would be advantageous to know the amount of the differential depletion under a given production history. It also would be valuable to know the pressure behavior, because it could be undesirable to produce a reservoir in such a manner that the pressure in one of the layers would have fallen below the bubble point of the oil at a time when another layer might still be virtually undepleted. The determination of static pressures in shut-in wells with slow build-up is another feature of importance in layered reservoirs. If the behavior of the reservoir differs appreciably from the predicted behavior of a single-layer reservoir, the type of analysis for build-up curves suggested for single-layer reservoirs is not applicable. Several authors have treated the problem of finding the pressure and production characteristics of multilayer reservoirs. Horner treated the problem of n layers in an infinite field, the layers being connected only at the well. Tempelaar-Lietz2 has presented an approximate treatment of the characteristics of a bounded reservoir composed of two layers of different permeabilities and equal thicknesses. Homer gave no information on differential depletion between layers, and the treatment of Tempelaar-Lietz was not applicable for prediction of pressure build-up behavior. In addition, it was thought that a more rigorous treatment than the one by Tempelaar-Lietz would be valuable in determining the range of validity of this solution. The present study, therefore, was made of a bounded reservoir composed of two or more layers,