Search Documents

By R. E. Cook

This paper presents results of a study to determine to what extent errors in estimated free gas saturation affect the results of static pressure calculations from build-up curves in two-phase systems. Use is made of the method of pressure build-up curve analysis developed by Miller, Dyes and Hutchinson1 , with modifications for two-phase flow advanced by Perrine2. Cornhined fluid compressibility of gas-oil mixtures is included, with special attention given to the discontinuity that occurs in fluid compressibility at the saturation pressure and its effect on static pressure calculations from build-up curves. The sensitivity of static pressure calculations to errors in gas saturation is considered over wide ranges of formation and fluid properties and reservoir pressure conditions. It is concluded from this study that for analyses of two-phase pressure build-up curves (I) the use of under-saturated crude compressibility data, even when producing GOR's are at or near solution values, will result in very significant negative errors in calculated static pressrrre, and (2) except for large values of buildup curve .slope and/or very low crude gravities, rather large errors in calculated gas saturation can he tolerated without .significant effect on calculated static pressures. The conc1usions arc qualified as referring only to the effects of gas saturation and excludc possible errors in the nuumerous other factors that enter into the ancllysis of pressure build-up curves. INTRODUCTION As part of the analysis of basic data for several major Eocene reservoirs in the Bolivar Coastal field (Western Venezuela), it was necessary to analyze a large number of pressure build-up curves for oil permeability and static reservoir pressure. From performance data it became apparent that gravity segregation is very active in these reservoirs. As a result of this performance it was recognized that the free gas saturation within the oil leg of the reservoirs would be small but would be subject to certain errors in calculation due to the inaccuracy of relative permeability data in the low gas saturation region. Since all of the various methods of build-up curve analysis were originally derived for sin- gle-phase, constant compressibility systems, and since the presence of a small but rather indeterminate gas saturation would greatly affect the reservoir fluid compressibility, a study was undertaken to investigate the effect of gas saturation on the results of two-phase pressure build-up analyses. TWO-PHASE BUILD-UP CURVE ANALYSIS The method of analysis used in this investigation is basically that originally developed by Miller, Dyes and Hutchinson1, with modifications for two-phase flow advanced by Perrine2. Miller, et al, showed that the original equation for single-phase permeability away from the well bore, held for two-phase flow except that the calculated permeability is the effective permeability to oil away from the wellbore. This finding was acknowledged by Per-rine, who also presented the equation for effective gas permeability* as Since Eqs. 1 and 2 contain no terms influenced by fluid compressibility, no distinct problem should develop in calculating permeability from two-phase build-up curves. In regard to the calculation of static reservoir pressure from build-up curves in single-phase systems, the modified Miller, et al, method as presented by Perrine is briefly as follows. 1. From build-up curve data, PVT data for the reservoir fluid and the well spacing and sand thickness, calculate the climensionless shut-in time corresponding to the highest value of actual shut-in time (t, hrs) on the straight-line portion of a semi-log build-up curve. Dimenrsionless shut-in time is related to actual shut-in time by the relation, t- 0.000264 kt 2. From a graph of dimension less shut-in time vs dimensionless pressure difference (Fig. 1 1 ) find the value of dimensionless pressure difference corresponding to the value of 7 calculated under (1) above. Dimen-

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 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 Pol Duwez

The effect of the rate of cooling on titanium, zirconium, and thallium has been measured. For titanium-molybdenum alloys, it has been shown that for molybdenum concentration up to 8 pct the rate of cooling has no effect on the transformation temperature, but this temperature decreases progressively with increasing molybdenum content. HIGH strength to weight ratio and remarkable corrosion resistance have often been advanced as the main reasons for the present widespread interest in titanium. For the physical metallurgist, the existence of an allotropic transformation in titanium is another important factor. As in the case of iron, an allotropic transformation may be the fundamental property upon which the response of Ti-base alloys to heat treatment will depend. It is therefore important to study the mechanism of the allotropic transformation in pure titanium and the effect of alloying elements on this transformation. The existence of an allotropic transformation in titanium was reported in 1931 by Schulze.' Later, De Boer et al.' determined the transformation temperature from a discontinuity in the electrical resistivity at 882"42O°C. According to recent work," the transformation temperature is 882.5"&1°C. The crystal structure of the room temperature a phase is hexagonal close-packed and the high temperature fi phase is body-centered cubic.' It is well known that the high temperature /3 phase of pure titanium cannot be retained by quenching to room temperature. Several investigators'-' have shown that in alloys of titanium with metals which are soluble in the @ phase (for example, nickel, molybdenum, chromium, or manganese), the /3 structure can be retained, providing the concentration of the alloying element is greater than a certain critical value. Below this critical value, the alloys transform during quenching, at least partially, into a supersaturated solid solution having the a. structure. This behavior is very similar to that found in some iron binary alloys (Fe-Ni and Fe-Mh, for example) and is an indication that the transformation might occur by means of a shear mechanism analogous to that observed in the austenite-martensite reaction. The present investigation is mostly concerned with the effect of the rate of cooling on the /3 to a transformation. In the first portion of this paper, the results obtained on pure titanium are described. Two other pure metals, zirconium and thallium, were also investigated, because they both have an allotropic transformation like titanium, from an hexagonal close-packed to a body-centered cubic structure. For purposes of comparison, the effect of rate of cooling of the y to a transformation of iron was included in the experimental program. The second part of this paper deals with Ti-Mo alloys. The equilibrium phase boundaries involving the solid state reactions on the Ti-rich side of the system were determined and a study was made of the effect of the rate of cooling on the transformation of the high temperature /3 solid solution. Experimental Procedure Since the main object of the present investigation was the determination of the effect of the rate of cooling on a transformation temperature, the best technique available was that used by Greninger" in his study of the martensite arrest in carbon steels. This method consists of heating a very small sample of metal to which fine thermocouple wires are .attached, cooling the sample rapidly by means of a blast of helium, and recording the temperature vs. time curve on an oscillograph. The apparatus used in the present experiments was very similar to that described by Greninger, except that the heating coil was made of molybdenum wire (0.040 in. in diam) and the assembly was placed in a bell jar under vacuum of 10" mm Hg or less. The specimen con-

Jan 1, 1952

By S. T. Yuster

An analysis has been made of the factors responsible for tilted fluid contacts in petroleum reservoirs. The factors are static and dynamic with the former being controlled by those variables responsible for the capillary rise of fluids, and the latter include rate of formation tilt and hydraulic circulation. The explanation that rate of formation tilt may be responsible for tilted contacts appears untenable. The role of some of these factors in the migration, fractionation. accumulation and possibly the prospecting for oil is indicated. INTRODUCTION While it is generally considered that the segregation of gas, oil, and water in a petroleum reservoir results in approximately horizontal planes separating these phases. there are factor; which modify this simplified picture to a very great extent. Some of these factors are static. and the departure of the fluid interfaces from horizontal may represent a static equlibriun. Some of these factors may be dynamic and the interface positions may be non-equilibrium and represent a steady .state condition. Of course combinations of these two situations may also exist. The reasons most commonly given for this behavior are hydraulic circulation of fluids, change in character of the formation, and a rate of tilt of the formations more rapid than the rate at which the fluids may re-establish an equilibrium. Russell's Principles of Petroleum Geology' contains an extensive discussion of tilted fluid interfaces and their importance in the migration and accumulation of petroleum. Russell cites examples of tilted water tables at Lance Creek, Wyo.; Rock River, Wyo.: Cut Bank. Mont.; Cushing, Okla.: and Graham; Okla. Also given was the Wasson Field of West Texas which has saucer shaped oil-water and gas-oil interfaces with the saucer being deeper in the case of the oil-water contact. The Bradford Field in Pennsylvania has a tilted water table. San Ardo. Coalinga None and some other California field.; have tilled contacts. A recent paper by King Hubbert discuses the effect of hydraulic circulation of the water underlying the oil in a given stratun and it- implications in the accumulation of and drilling for the oil. The foregoing would indicate that tilted fluid interfaces are not uncommon and that it would be of importance to underbtand the forces which are responsible. Further. it would be of interest to examine certain idealized models of situations which may cause tilting with the thought in mind of deriving relatinnships between the important variables. These relationships may only hold exactly for the model* assumed but the results may apply semi-quantitatively to actual formations and at least indicate trends. The discussion will be divided into the effect of static factors and the effect of dynamic factors. STATIC FACTORS It has already been mentioned that one reason given for tilted contacts is a variation in formation properties. Specifically this is a change in the pore size of the formation capilaries. If the internal surface of the reservoir is preferetially wet by water, capillarity will cause the water 10 rise to a greater extent above the free water table in a small capillary than in a large one. Consequently the water-oil contact will be higher in the lower permeability sands than in the higher permeability sand.;. The free water table may be defined as the position the oil-water interface would take if capillarity were absent. Since the phenomenon is the capillary rise of water in an oil environment, a derivation of this shold give a relationship between the variables. Fig. 1 is an idealized model of the situation with a uniform bore circular capillary ubstitttting for the tortuous non.circular formation pore. The upward force.; on the column of water will be equated to the forces downward to give the equilibrium situation. forces upward = force downward where r is the radius of the capillary in ceritimetera. a is the interfacial ten4on oil-water in dynes per centimeter, 0 is the

Jan 1, 1953

By Robert K. Waits

The resistivity, temperature coefficient of resistance, stability at 200°C, and structure of annealed Si-Cr films have been studied as a function of film cowposition. Colorimetric analyses of the films indicated 20 to 60 wt pct Cr over the resistivity range of 1 to 10-4 ohm-cm. Resistors were delinealed by photolithographic etching, contacted with aluminum, and mounted in TO-5 packages. The average temperature coefjicient of resistance, a, between 50" and 150°C is related empirically to the film resistivity, p, by a (ppm per °C) = -2900 - 1140 log p (ohm-cm). Near a resistivity of 3.5 x 10-3 ohm-ctn (30 at. pct Cr) , the temperature coefficient changes from positive to negative. The temperature at which the resistance maximum occurs increases with decreasing resistivity. The conduction mechanism between —196" and +360 C appears to have two components—an activated conduction dominating a1 lord chromium content and higher temperatures and a linearly temperature-dependent conduction dominating at high chromium content and low temperatures. Electron diffraction patterns showed the films to consist of Cr3Si and amorphous silicon. Resistors thicker than 150A typically drifted less than 0.6 pct after 5000 hr stor-age at 200°C. Resistors as thin as approxinzately 80°A that had a silicon overlayer showed similar stability. Such films had sheet resistivities of 10 kilohm per square and temperature coefficients near —600 ppm per "c. THIN-FILM resistors having sheet resistivities greater than 1 kilohm per square are required in many integrated circuit applications. A process for producing such resistors should be compatible with present circuit fabrication and packaging, including exposure to maximum temperatures of 300° to 570°C during contact alloying, die-attach, and ceramic package sealing. In addition, the resistors should have a low temperature coefficient of resistance and change less than 1 pct during storage at temperatures up to 200°C for 1000 hr or longer. Reports of work beginning in 1949 at Battelle Memorial Institute indicated that, among other materials, the silicides of chromium, Cr-Co, molybdenum, tantalum, and tungsten were worthy of further study as high-stability resistive films. Later work on Si-Cr films vacuum-deposited by flash evaporation from a tungsten boat1 led to the conclusion that Si-Cr had "excellent properties for development into high- stability ... resistive elements". More recently, Si-Cr films have been vacuum-deposited from a tantalum boat,2 from a BeO-coated tantalum boat,3 by flash evaporation from an electron-beam-heated tungsten Source,4 and by triode sputtering.5 This paper describes a study of Si-Cr films vacuum-deposited by electron-beam heating of Si-Cr alloys. Resistivity, temperature coefficient of resistance, stability at 200°C, and structure were studied as a function of film composition. EXPERIMENTAL PROCEDURE Film Deposition. Si-Cr sources were prepared by vacuum-melting pressed pellets of mixtures of powdered chromium (99.9 pct, -100 mesh, Lunex Co., Pleasant Valley, Iowa) and silicon (99.99 pct, -325 mesh, United Mineral and Chemical Co., New York). A permanent-magnet focused electron-beam gun having a water-cooled copper crucible was used to heat the pellet until complete melting occurred. The mixture was then kept molten for an additional 5 min. The alloy thus formed was used as the source material for Si-Cr film depositions using the same electron-beam gun. Pieces of semiconductor-grade silicon ingots were used as sources for silicon depositions. The films were deposited within a standard vacuum evaporator having an 18-in.-diam Pyrex bell jar, a liquid-nitrogen-cooled trap, and a Consolidated Vacuum Corp. 6-in. diffusion pump containing Dow Corning DC-705 pump fluid. The substrates were thermally oxidized* silicon wafers 2 to 2.5 cm in diam. The

Jan 1, 1969

By P. B. Mee

The secondary re crystallization behavior of two vacuum-melted, nominally 3 pct Si-Fe alloys has been studied, one alloy being produced from pure electrolytic iron and semiconductor-grade silicon and the other from commercial-purity hot band. The vacuum-melted and vacuum-remelted alloys were found to develop {110}(hhl) and {lll)(hkl), respectively, on annealing in high vacuum and {100}(0kl) and {11l}(hkl), respectively on annealing in pure hydrogen of dew point, < — 60°C. H2S doping of the hydrogen or annealing in commercial argon was found to effect a reversal for the vacuum-remelted alloy from (111) to (100) secondary recrys-tallization, while prolonged vacuum annealing at 2 x 10-6 Tow resulted in a tertiary re crystallization from (111} to (110). The growth selection of the {lll}(hkl) component was found to be governed predominantly by the {111}/planar misorientation, the maximum deviation of the resulting secondaries being 5 deg. It is considered that the driving force for {111}(hkl) secondary recrystallization arises from the preferential lowering of ?111, possibly by phosphorus or nitrogen diffusion, just as oxygen and sulfur under certain conditions have been found to preferentially lower ?100. When heavily cold-rolled metals are annealed at high temperatures, primary recrystallization usually occurs rapidly and may be followed by considerable grain growth. However, in the presence of a fine dispersion of inclusions1"3 or with a strong matrix preferred orientation4&apos;&apos; the matrix may be stabilized and a secondary recrystallization may then be effected. Also, in vacuum-melted Si-Fe it has been shown possible to develop {100}(001) and {110}(001)6 surface energy driven secondary recrystallization textures. Further, it has been found that, on annealing very pure metals, grain growth normally ceases when a grain size approximately twice the strip thickness is acquired, due to the formation of V-shaped thermal grooves which result in an equilibrium being set up between surface and grain boundary stresses. If, however, a surface energy difference ??s, exists between adjacent grains (of the order of a few percent of the grain boundary energy ?b), then according to Mullins8 a secondary recrystallization may be effected, since the grains having lowest ?s would grow at the expense of their higher-energy neighbors. For clean bee metal surfaces ?110 should be lowest; however, for 3 pct Fe-Si (100) secondary structures have frequently been developed and have been attributed by Walter and Dunn9 to an alteration of the relative surface energies due to the preferential chemisorption of oxygen on (100) planes. Kohler,10 however, indicated that this reversal may be achieved by employing minute quantities of highly polar compounds in the annealing atmosphere, i.e., sulfur compounds, in par- ticular H2S being beneficial in (100) growth. This work has been verified by Kramer" using {100)/{110} bicrystals, while Mager12 demonstrated the importance of sulfur diffusion within the strip. Detert13 has intimated that carbon may be responsible for the preferential lowering of ?110 although no direct evidence has appeared yet to support this supposition. It must , be appreciated that the above-mentioned examples of surface energy driven secondary recrystallization were obtained with very pure vacuum-melted materials with oxygen contents (which is now known to be a critical factor") less than 35 ppm. During the course of a textural study on Si-Fe it has been found that under certain conditions (111) secondary recrystallization can be obtained by what is believed to be a surface energy growth selection mechanism. The purpose of the present paper is to briefly describe this result and to try and postulate the possible causes. EXPERIMENTAL Two vacuum-induction-melted, nominally 3 pct Si-Fe alloys were produced, one from pure electrolytic iron and semiconductor-grade silicon designated A and the other from commercial-purity hot band, designated B.

Jan 1, 1969

By Alton B. Cook, G. B. Spencer, F. P. Bobrowski

In estimating production gas/oil ratios and oil recoveries from reservoirs containing highly volatile oils it is highly important to include condensate that may be recovered from the gas produced from the reservoir. The volume of the hydrocarbon liquid condensed from the production of the solution gas that has been liberated. in some reservoirs may equal or even exceed the volume of recoverable stock-tank oil as estimated by present methods. The usual procedure1, 2, 3 in estimating future reservoir performance includes use of relative permeability data of the reservoir rock to gas and oil and includes material-balance calculations utilizing data on the volume of the reservoir fluids, oil and gas production and reservoir pressure decline, and data on the physical properties of the initial reservoir fluids. The inaccuracy in the procedure results from the false assumption that all of the free or liberated gas that enters the well bore of reservoirs containing highly volatile oils remains in the gaseous phase as it is produced. This paper presents a method for estimating future reservoir performance and oil recoveries based on special labora. tory analyses of reservoir-oil samples and recognition of the additional volume of hydrocarbon liquid that will be recovered at various stages of depletion from a solution.gas.drive reservoir. The method also includes calculations of the volume and composition of the hydrocarbon liquids that can be recovered by processing the produced gas in a natural-gasoline plant. The development of the method resulted from a study of laboratory analyses of reservoir samples and of field data obtained from two reservoirs containing oil of the highly velatile type. Relatively high pressures and temperatures prevail in oil reservoirs discovered at the greater depths of present-day drilling. The contained fluids may have high saturation -Dressures and large volumes of gas in solution. The oils and gases under these conditions may be almost identical in compositions and densities.4,5 Also, engineers have noted that, for a number of oil reservoirs and especially the deeper reservoirs, the API gravity of the stock-tank oil increaees as the production gas/oil ratios increase. They attributed this increase in API gravity of the stock-tank oil to the formation of condensate as the liberated gas in the reservoir is being produced. SOURCE OF DATA FOR CALCULATIONS Subsurface oil samples and well-test data were obtained from the E. Constantin-Lambert Well No. 1 in the Elk City Field, Okla., to demonstrate the need for and a method of calculating future reservoir performance when consideration is given to the recoverable hydrocarbon liquid from the production of solution gas that has been or will be liberated in a reservoir. The Lambert well is located in C, SW 1/4, SW 1/4, Set- 18, T-l0-N, R-20-W, Washita County, Okla. At the time of completion it was on the east edge of the development and Was producing from conglomerate through 30 perforations in the 5 1/2-in. casing at a depth of 9,859 to 9,864 ft. The well Was completed on July 9, 1949, and subsurface oil samples were obtained July 29 and August 1, 1949. Production gas/oil ratio data were obtained during the operation in series of two Separators at high and low-pressures of 314 and 44 psia, respectively, and at an average temperature of approximately 60°F. The flowing bottom-hole pressure at the time of sampling was 3,934 psis. The static bottom-hole pressure was not obtained at the time of sampling, as the well could not be shut in long enough for testing because gas produced from the well was used for fuel at an offset drilling well. However, the initial reservoir pressure was calculated to be 4,364 psia from knowledge of pressure gradients in other wells in the field and the depth of the casing perforations in the Lambert NO. 1 well.

Jan 1, 1951

By H. N. Mead

In a reservoir when gas comes out of solution and rises, additional pressure is created because of the change in position of this gar in the bounded volume. If this pressure effect is not taken into account as a pseudo-influx ill mate rial balance calculations on reservoirs in which there is evolving gas, an error is introduced that is directly proportional to the length of oil column and amount of gas evolved and inversely proportional to the reservoir pressure. In reservoirs with small oil column this pseudo-influx is of little importance. The larger the oil column the more important the effect becomes. The most important application, however, will be for reservoirs with declining pressure in which increasing water entry has been indicated by material balance calculations. It could be that this indication is the result of the cited pseudo-influx and not water entry at all. This paper attempts to explain this phenomenon in relation to material balance calculations and prevents an approximate method of determining tile pseudo-influx. INTRODUCTION Only recently has the pressure-volume relationship of gas rising in a liquid appeared in engineering literature, although its effect has been discussed for several years—especially among engineers engaged in drilling operations. The writer published a paper in 1957 in which this phenomenon was examined in relation to blowouts in drilling wells.&apos; Also, in 1957 Stegemeier and Matthews&apos; authored a paper in which the relationship was presented with respect to pressure build-up tests when shutting in producing wells. In this paper the effect will be called "pressure inversion" which is descriptive of the phenomenon. Fundamentally, the effect just mentioned is the result of changing the position of a quantity of energy contained within the gas itself (as a result of its compression) within the bounded volume in which the gas exists. If the gas were allowed to expand "beyond all bounds" there could be no pressure inversion effect. This is because all of the energy contained within the gas as a result of compression would have been re- PRESSURE INVERSION A phase that has not been discussed in the literature is the effect of pressure inversion upon material balance calculations in reservoirs in which the reservoir pressure drops with reservoir withdrawal and is below the bubble point. This action releases solution gas which, after the critical gas saturation has been reached, will tend to rise and collect at the top of the reservoir. The pseudo-influx caused by this effect could be great enough in reservoirs of large oil column and high permeability to render invalid any analyzer studies made on them — if pressure inversion has not been taken into account. Most suspect would be those reservoirs with large oil column in which increasing water entry has been predicted but the bottom wells have not yet produced water. Example 1 Simply stated, this phenomenon exists in oil reservoirs because of the difference in density and expansibility between gas and oil. For example, if a bubble of gas exists at the bottom of a 1,000-ft column of oil having a gradient of 0.30 psi/ft with the pressure at the top of the column zero, the bubble of gas will be compressed to a pressure of 300 psi. If temperature and the combined volume of the oil column and the gas bubble remain constant, as the gas bubble rises in the oil column it will have the same volume and consequently the same pressure at any point. Under these conditions when thc bubble reaches the top of the oil column, the pressure will be approximately 300 psi at the surface and approximately 600 psi at the bottom of the oil column. If it is assumed that the oil is incompressible the pressures will be exactly 300 and 600 psi, respectively.

By F. Paredes, W. R. Holman, R. W. Crawford

The diffusion coefficient for hydrogen in the ß titanium alloy containing 13 pct V, 11 pct CY, and 3 pct A1 was measured over the temperature range 20° to 500°C. Results fit the expression: D= 1.58 x lgs exp-------==-----) cm2,sec THE diffusion coefficient for hydrogen in the metastable ß-phase titanium alloy containing 13 pct V, 11 pct Cr, and 3 pct A1 was measured as part of an investigation of the effects of hydrogen on the mechanical properties of titanium alloys. Several theories on the mechanisms of hydrogen embrittle-ment include diffusion of hydrogen at ambient temperature as an important step in the embrittling process.&apos;-5 In hydrogen-contaminated a-ß titanium alloys, it has been postulated that the rate-controlling step may be the rate at which hydrogen diffuses from the ß phase to stress- or strain-nucleated brittle hydride platelets at the a-ß interface.&apos; However, embrittlement of the all-ß alloy investigated in this work may proceed by a different mechanism, since this alloy can be embrittled without the formation of detectable hydride platelets.6 Troiano has suggested that in this case the mechanism of delayed failure or embrittlement is similar to that proposed for high-strength steel, i.e., repeated crack initiation caused by a decrease in bonding energy as hydrogen diffuses to points of high tri-axial stress at the root of a notch. In both types of alloy the rate of embrittlement is believed to be controlled by diffusion of hydrogen in the ß phase. Accurate hydrogen-diffusion data in the vicinity of room temperature should therefore be of some assistance in understanding the mechanisms of embrittlement in these two classes of titanium alloys. Diffusion of hydrogen in pure titanium has been studied by Wasilewski and Kehl7 at temperatures above 500°C. Their measurements on diffusion in the ß phase were extended below the ß-a transformation temperature (882°C) to approximately 600°C. This was possible because of the 8 stabilizing effect of hydrogen. No results at lower temperatures have been published and there are no data in the literature on hydrogen diffusion in the ß -phase alloy used in this program. Hydrogen diffusion in metals is generally investigated by measuring the rate of transport of hydrogen past a solid-gas interface, i.e., by absorption, permeation, or degassing-rate measurements. These methods produce accurate diffusion data when bulk diffusion is the rate-determining step in the over-all process. However, at the lower temperatures of interest in hydrogen-embrittlement studies or under conditions where surface films are present, the measured rates may be surface-controlled and the calculated diffusion coefficients will be in error. Surface effects were eliminated in the present work by using a modification of the common welded-couple technique m which two bars with uniform but different compositions are joined before the diffusion anneal." The same type of step change in the initial concentration vs distance curve was produced by electrolytically hydrogenating thin-strip specimens which had been electrically insulated over half of their length. Concentration vs distance curves were determined, after appropriate diffusion-annealing treatments, by vacuum-extraction analyses of transverse samples sheared from the strips. Diffusion coefficients were calculated by the methods of Grube.M This technique was feasible because of the large concentrations of hydrogen which can be charged into p -phase alloys without the formation of hydrides15&apos;16 and because of the extremely low degassing rates observed when hydrogenated titanium alloys are heated in air. EXPERIMENTAL WORK The diffusion couples were cut from a 0.022-in.-thick mi ll-annealed sheet.* Strips 8 in. long by 1/2

Jan 1, 1965

By R. Lagneborg, R. Gullberg

An X-ray method for determinig volume fractions of phases in textured materials has been developed. The method involves measurements of the integrated intensities of reflections of each phase for one single orientation of the sample. The method was checked by evaluating the phase contents in samples from two stainless Cr-Ni steels containing austenite and ferrite in known amounts. In the present case measurements of about five reflections from each phase were required. The agreement between calculated and known volume fractions was found to be very close. AN X-ray diffraction method is often employed for determining the amounts of two or more phases in a material. The details of the method have been described at several instances.&apos;" In short, the method is based upon the fact that the integrated intensity of a reflection is directly proportional to the volume fraction of the phase considered. In principle, the method only requires the measurement of the intensity of one reflection for each phase. The most severe shortcoming with this method is that it disregards influences from preferred orientation in the material. In a diffractometer, where this type of analysis is most conveniently carried out, the geometry is such that the diffracting planes for a specific Bragg peak are all parallel to the surface of the specimen. This implies that it must be required that the occurrence of the diffracting planes parallel to the specimen surface must be as frequent as for any arbitrarily oriented physical plane of the sample. This requirement will generally only be fulfilled for a material free of texture. Since all materials in reality always exhibit a preferred orientation to some extent. the deficiency of the method is quite serious. In cold-deformed material the texture is appreciable, but also annealed material will in many cases display a considerable texture. The problem of quantitative phase analysis in textured material is one of obtaining the integrated intensity averaged over all orientations of the specimen with respect to the X-ray beam. Lopata and ~ula~ have shown that this can be accomplished by integration of intensities over the whole pole figure, or by random rotation of the specimen in the X-ray beam, or by a combination of both. Hence, in this method the average integrated intensity for a specific reflection is obtained by directly determining the intensity for a sufficient number of physical orientations of the sample. However, since an excess of diffracting planes in a certain surface plane necessarily must be accompanied by a deficiency of another type of diffracting plane, it appears to be possible to obtain the correct average intensity for a specific diffracting plane by performing an appropriate averaging of intensities belonging to a large number of diffraction planes of different hkl when all intensities have been produced from one single orientation of the specimen surface. In the present work such a method has been attempted and tested experimentally for two stainless Cr-Ni steels containing austenite and ferrite. 1) THEORY The expression for the integrated intensity of a Bragg reflection diffracted from a phase in a polycrys-talline sample in a diffractometer is: where C is a constant only dependent on the incident X-ray power and the diffractometer geometry, v the volume of the unit cell, F the structure factor, m multiplicity of the diffracting plane, L Lorentz factor, P polarization factor, e-2 temperature factor, p average linear absorption coefficient of the irradiated vO1umej and V volume fraction of the phase &apos;Onsidered. Hence, the factor R in the abbreviated form of Eq. [I] can be calculated and depends only on 0 and the diffracting plane {hkl}. For an alloy containing two phases, a and y, we can write

Jan 1, 1967

By F. P. Bobrowski, G. B. Spencer, Alton B. Cook

In estimating production gas/oil ratios and oil recoveries from reservoirs containing highly volatile oils it is highly important to include condensate that may be recovered from the gas produced from the reservoir. The volume of the hydrocarbon liquid condensed from the production of the solution gas that has been liberated. in some reservoirs may equal or even exceed the volume of recoverable stock-tank oil as estimated by present methods. The usual procedure1, 2, 3 in estimating future reservoir performance includes use of relative permeability data of the reservoir rock to gas and oil and includes material-balance calculations utilizing data on the volume of the reservoir fluids, oil and gas production and reservoir pressure decline, and data on the physical properties of the initial reservoir fluids. The inaccuracy in the procedure results from the false assumption that all of the free or liberated gas that enters the well bore of reservoirs containing highly volatile oils remains in the gaseous phase as it is produced. This paper presents a method for estimating future reservoir performance and oil recoveries based on special labora. tory analyses of reservoir-oil samples and recognition of the additional volume of hydrocarbon liquid that will be recovered at various stages of depletion from a solution.gas.drive reservoir. The method also includes calculations of the volume and composition of the hydrocarbon liquids that can be recovered by processing the produced gas in a natural-gasoline plant. The development of the method resulted from a study of laboratory analyses of reservoir samples and of field data obtained from two reservoirs containing oil of the highly velatile type. Relatively high pressures and temperatures prevail in oil reservoirs discovered at the greater depths of present-day drilling. The contained fluids may have high saturation -Dressures and large volumes of gas in solution. The oils and gases under these conditions may be almost identical in compositions and densities.4,5 Also, engineers have noted that, for a number of oil reservoirs and especially the deeper reservoirs, the API gravity of the stock-tank oil increaees as the production gas/oil ratios increase. They attributed this increase in API gravity of the stock-tank oil to the formation of condensate as the liberated gas in the reservoir is being produced. SOURCE OF DATA FOR CALCULATIONS Subsurface oil samples and well-test data were obtained from the E. Constantin-Lambert Well No. 1 in the Elk City Field, Okla., to demonstrate the need for and a method of calculating future reservoir performance when consideration is given to the recoverable hydrocarbon liquid from the production of solution gas that has been or will be liberated in a reservoir. The Lambert well is located in C, SW 1/4, SW 1/4, Set- 18, T-l0-N, R-20-W, Washita County, Okla. At the time of completion it was on the east edge of the development and Was producing from conglomerate through 30 perforations in the 5 1/2-in. casing at a depth of 9,859 to 9,864 ft. The well Was completed on July 9, 1949, and subsurface oil samples were obtained July 29 and August 1, 1949. Production gas/oil ratio data were obtained during the operation in series of two Separators at high and low-pressures of 314 and 44 psia, respectively, and at an average temperature of approximately 60°F. The flowing bottom-hole pressure at the time of sampling was 3,934 psis. The static bottom-hole pressure was not obtained at the time of sampling, as the well could not be shut in long enough for testing because gas produced from the well was used for fuel at an offset drilling well. However, the initial reservoir pressure was calculated to be 4,364 psia from knowledge of pressure gradients in other wells in the field and the depth of the casing perforations in the Lambert NO. 1 well.

Jan 1, 1951

By G. R. Purdy, J. A. Malcolm

The isothermal dendritic precipitation of y in p CuZnSn is compared with the isothermal formation of Widmanstatten a plates in p CuZn. Although both precipitates are evidently bounded by semicoherent interfaces, the dendritic shape of y evolves through smooth "sinuoidal" shape perturbations, while a-p interfaces tend to be strictly planar, or rnicroscopically faceted. An estimate of the average inter-facial free energy of cube-oriented y in p through the Mullins-Sekerka analysis indicates that it is low (of the order of 50 ergs per sq cm). Interface mication is apparently controlled by volunle dij&apos;usion, and the observed anisotropy in growth probably results from anisotropy of inter facial free energy. Studies of the grain boundary precipitation of y in p indicate that the precipitate can grow in both $ grains, but becomes unstable with respect to the grain with which it shares a mobile semicoherent interface. Facets observed at the tips of Kurdjumov-Sacks-oriented Wid-manstatten a plates correspond to variants of an irrational hahit, near (231)0. Comparison with Hu and Smith&apos;s determination of relative interfacial free energies of a-(3 boundaries suggests that this system possesses an extremely deep o cup, which results in an ej&apos;ective low-o zone. The possible effects of such a cup on the stability of a-p interfaces are discussed, a needles, grown during rapid continuous cooling, possess microscopically smooth interfaces, and are characterized by a high degree of&apos; damage in the region of the boundaries. These boundaries, induced to move on tempering at 200"C, become faceted as equilibrium is approached, which seems to indicate that boundary structure and mobility in this system are dependent on the local degvee of departure from equilibrium MULLINS and sekerkal have elucidated the phase boundary stabilizing effect of surface tension. a factor not taken into account in previous treatments of morphological development. Their analysis applies to the isothermal growth of a pure spherical precipitate (planar in the limit of infinite radius) in a slightly supersaturated solution, under conditions of inter-facial local equilibrium and isotropy of interfacial free energy. More recently, several other factors which might tend to stabilize planar interfaces have received attention. These are: i) interface diffusion;4"8 and ii) slow interfacial reaction. cahn7 included a slight anisotropy in interfacial free energy, and showed that it has no stabilizing effect, but forces shape perturbations to reflect the anisotropy. shewmon4 noted that dendritic morphologies are rarely observed in solid-solid systems, and, in particular, postulated that slow interfacial reaction occurs at all boundaries between proeutectoid ferrite and austenite, except those which propagate parallel to {11 l)y planes (for Kurdjumov-Sachs [K-S] related phases). Previously, Aaronson~ had proposed interfacial reaction at the broad faces of Widmanstatten plates, and expected that the lateral motion of high-mobility incoherent ledge edges would produce planar precipitates bounded by law-mobility semi-coherent interfaces. The problem of interfacial reaction during the dif-fusional formation of precipitates has often been ignored, although it is apparent that a small component of chemical reaction must exist at any moving interface." Further, all solutions of the diffusion equations for boundary conditions corresponding to isothermal growth predict infinite growth rate at time zero, so that elements of interfacial reaction would be most easily detected at early times. However, reported deviations from diffusion control are comparatively rare," and the local-equilibrium diffusion-control model appears adequate for the description of many isothermal growth systems.&apos;2"9 Cahn14 has treated the problem of migration of an interface which is able to reach positional energy minima, and concluded that such an interface, if subjected to a high local ther-modynamic "driving force", may migrate normal to itself without the benefit of ledges. As local equi-

Jan 1, 1968

By G. Krauss, J. M. Thompson, H. Y. Kumagai

Boron-doped p-type and arsenic-doped n-type source materials were used to deposit thin silicon films on amorphous fused quartz substrates by cathodic sputtering in argon atmospheres. All as-sputtered films were found to have high electrical resistivities of about 104 ohm-cm, despite controlled variation of the substrate temperature from 50O to 400°C. The high resistivity of the n-type films persisted even after postdeposition heat treatments as high as 1000°C, while p-type films showed a sharp decrease in re-sistivity to values on the order of 1 ohm-cm or lower after heat treatments above 650°C. X-ray and electron diffraction, together with transmission electron microscopy of the p-type films, revealed that the improved electrical properties resulted from two distinct processes. A primary crystallization or recovery within the pains, which were approximately 200Å diam in the as-sputtered films, accompanied the initial sharp drop in resistivity at 650°C. Following this process, a secondary recrystallization or discontinuous grain growth was nucleated in films heated at 1000°C. The resulting significantly increased size of certain grains and the corresponding decrease in grain boundary area accompanied an additional decrease in resistivity after heat treatment at 1000°C. The properties of thin silicon films are of interest because of their possible application in device technology. Although evaporated silicon films on amorphous substrates1-3 and both sputtered and evaporated germanium films4-&apos; have been studied in some detail, sputtered silicon films have not been widely investigated. It has been observed that evaporated silicon films deposited at substrate temperatures below approximately 600°C are amorphous.1,3,9,10 At substrate temperatures higher than this, the films are polycrys-talline, sometimes exhibiting a fiber texture with the fiber axis normal to the substrate surface.&apos; Similar results have been obtained with germanium films, although in this case the substrate temperature for the amorphous to crystalline transformation is reported to be as low as 150OC.5 The mechanism of film deposition by cathodic sputtering differs markedly from that of thermal evaporation. It might be expected, therefore, that the properties of sputtered and evaporated silicon films would also differ. The purpose of our study was to determine the effects of deposition conditions and heat treatments on the electrical resistivity, structure, and composition of thin silicon films sputtered on fused quartz substrates. EXPERIMENTAL DETAILS Apparatus. A schematic diagram of the sputtering apparatus is shown in Fig. 1. The vacuum chamber was enclosed by an 18-in.-diam by 12-in.-high Pyrex bell jar mounted on a stainless-steel base plate. Access holes in the base plate were used to supply argon gas to the chamber, to allow for various electrical feedthroughs, and to provide mechanical motion to the shutter between the source material and substrate. The chamber was evacuated by a pumping system consisting of a 720 liter per sec oil diffusion pump backed by a 13 cu ft per min mechanical pump. Ultimate pressures of less than 1 x 10-6 torr could be achieved. The substrates were securely mounted on the stainless-steel substrate holder by quartz tabs held in place by stainless-steel screws. Temperatures of the substrate holder were measured by a thermocouple inserted within the holder. A quartz lamp was used to heat the holder and temperatures as high as 420°C could be maintained. Deposition Conditions and Procedure. The substrates used were polished fused quartz slides 116 by 1-12 by 3 in. and were outgassed at 200°C for 1 hr at a pressure less than 1 x 10-5 torr before film deposition. This was followed by a presputtering period of 2 hr during which the movable stainless-steel shutter was so positioned that deposition on the substrate was prevented. At a sputtering voltage of 5000 v dc and an argon pressure of 0.03 torr, the cathode cur-

Jan 1, 1967

By J. Szekely, G. P. Martins

A formulation is given and a general analytical solution is presented for the problem of spherical phase growth in multicomponent systems. The results are valid for describing the rate of bubble growth in the absence of translational motion and, therefore, should be relevant to heterogeneous nucleation processes. Typical worked examples are given, which include bubble growth due to the simultaneous diffusion of hydrogen and nitrogen in molten steel, and the growth of CO bubbles in steel, due to the reaction: C + 0 = CO. The calculations predict that the presence of additional components would enhance the growth rate and also increase the absolute rate at which the individual components are transferred into the bubble. These findings should be relevant to degassing processes. THE problem of spherical bubble growth in the absence of translational motion is appropriate to the description of the heterogeneous nucleation and subsequent growth of gas bubbles while attached to a solid surface. Such problems are of considerable practical relevance to various vacuum degassing and steelmaking processes where bubble nucleation and subsequent growth play an important part in the overall reaction mechanism. A definitive formulation of spherical phase growth has been given in the now classical paper by Scriven1 and the analogous problem for solid phase growth has been discussed by Horvay,2 furthermore a thorough review of the relevant literature has also become available.3 All this previous work was concerned with one component or binary systems. The purpose of this paper is to show that the technique developed by Scriven may be generalized to alow the consideration of multicomponent systems. The results of this generalization would be particularly relevant, as gas bubble nucleation in the metallurgical systems previously mentioned frequently involves the simultaneous diffusion of several components. In the following a general formulation is given and typical worked examples are presented. FORMULATION Consider the growth of a spherical gas bubble due to diffusion in a quiescent liquid of infinite extent, as sketched in Fig. 1. The formulation of such problems requires the statement of the following conservation equations relating to the liquid phase: Equation of continuity Equation of motion Conservation of the diffusing species—Diffusion Equation. The following assumptions are made: a) The fluid is incompressible, quiescent, and of infinite extent b) The bubble is stationary, i.e., it undergoes no translational motion c) Spherical symmetry is observed throughout d) The mass diffusivity is constant e) Chemical equilibrium exists at the surface of the bubble, and any reactions occurring there are fast compared to the diffusion process f) Inertia and surface tension effects are neglected. Equations derived on the basis of these assumptions by previous investigators were found to provide a reasonable representation of bubble growth while attached to a surface in organic and aqueous systems4,5 For gas bubble growth in molten metal systems assumptions a) to d) chould hold reasonably well. The exclusion of chemical kinetics as a possible rate limiting step may not be universally valid, but should hold for many systems operating at high temperatures. The effect of hydrodynamic inertia may be neglected for moderate growth rates; finally, the assumption neglecting the effect of interfacial tension on the bubble growth was made for the sake of mathematical convenience; its inclusion would have required the use of quite cumbersome numerical techniques. While surface tension effects could be appreciable on the initial values of the overall growth rate, especially at low absolute pressures and in case of small bubbles, they are not thought to modify the relative effects present in multicomponent systems, the study of which is the objective of this paper. Within the framework of the above assumptions

Jan 1, 1970

By H. E. Wilborg, G. V. Henderson

This chapter deals with those sands employed by foundries for the manufacture of cores and molds used in the casting of such common metals as steel, gray iron, ductile iron, aluminum-based alloys, and copper-based alloys. Most metal is cast in "green sand" molds in which the sand mixture, consisting primarily of silica, clay, and water, forms a plastic material. This plastic material is rammed around a pattern in a split box called a flask, then the box is separated, the pattern is withdrawn, and the metal is poured. These green sand molds are broken apart when the metal has solidified. The sand is then remixed and a new mold made for the next casting. To a lesser extent silica sand is mixed with resins or oils and, after forming, the mixture is hardened by either heat or chemical reaction. In any case, the foundry mold or core is used just once before being broken up and remixed for manufacturing another core or mold. By far the greatest consumption in the foundry is of those types of sand high in silica, low in clay (under 0.5%), and occurring in numerous locations thoughout the United States. Consumption End Uses Foundry sands fall into four distinct categories: silica sand, lake sand, bank sand, and natural molding sand. Silica sand is a general term used to describe washed, graded, dried, and cooled clay-free sands which are found in many parts of the United States. They generally contain 90 to 100% silica but vary greatly in chemical analysis and grain structure. Lake sand is the sand dredged from Lake Michigan or dug from the dunes bordering that lake. Bank sand is composed primarily of silica grain but containing up to 5% clay. It is shipped in both the damp and dry state but the normal processing of bank sand consists only of drying and scalping. Occasionally two different grades of bank sand are blended prior to or during loading. Natural molding sand is a term used to describe a sand containing sufficient clay as mined to be suitable for molding purposes. Clay content may vary from 10 to 24% with the balance of the material consisting primarily of silica sand. These sands range in fineness from one peaked on the 20 mesh to one which contains more than 50% of the particles passing the 270 mesh. Lake sands are relatively coarse and are used in great quantities for the production of automotive castings. Bank sands and natural molding sands are gradually phasing out as the production of castings becomes more a science and less an art. The vast majority of foundry sands used today are those known as the silica sands and lake sands. Foundry sands are used in two basic ways. The majority are used in the green or uncured state with the balance used in the cured or set state. Green sand is grain bonded with clay and contains varying amounts of water to plasticize the clay in order to bond the grains together. Green sand is either naturally bonded or "synthetic" sand. The naturally bonded molding sands are generally found in the Hudson River Valley of New York State, Cumberland County in New Jersey, and in Tennessee and Mississippi. Synthetic sand is molding sand compounded in the individual foundry. It is made up primarily of graded silica sand, western or southern bentonite, water, and carbonaceous and/or cellulosic materials. Green sand is used primarily in making molds. The use of naturally bonded green sand has been declining for many years, primarily because naturally bonded sands are more difficult for the foundryman to find and control. Synthetic green sand, on the other hand, can be

Jan 1, 1983

By C. Bezemer, I. Havenaar

An investigation was carried out on the dynamic liltration behavior of drilling fluids. In a set-up consisting of a porous pipe through which the drilling mud was circulated, dynamic liltration rates and cake thicknesses were determined as a function of flow rate and filtration time. Eventually, a state of equilibrium was reached; i.e., both cake thickness and filtration rate attained constant values. Concurrently, dynamic filtration tests were carried out in a simple dynamic filter apparatus where mud flow was similar to that in a rotational vismmeter. In both arrangements the relation between equilibrium liltration data and rate of shear at the cake surlace was studied. Filter pipe and dynamic filter apparatus were lound to give identical relations. It was therelore concluded that the latter apparatus, which is suitable lor routine measurements on small amounts of mud, could be used to determine the dynamic liltration behavior of mud under borehole flow conditions. Some data are presented on the effect of filtration pressure and temperature on equilibrium filtration rate and cake thickness. INTRODUCTION In oil well drilling, the drilling fluid used generally forms a mud cake on the permeable sections of the borehole wall while filtrate is lost into the formation. As to the conditions under which this filtration process takes place, a distinction should be made between static mud filtration (mud is not circulating) and dynamic mud filtration (mud is circulating). In recent years several investigatorsl-5 have studied the process of dynamic mud filtration about which little had previously been known. This type of filtration is the subject of this paper which confirms several findings of previous investigators. It differs from already published works, however, in the following respects. 1. It contains a description of a simple laboratory apparatus developed for dynamic filtration measurements on small amounts of mud. 2. The effect of rate of shear at the cake surface on dynamic filtration has been studied. As a result of this study, dynamic filtration data pertaining to the borehole conditions of mud flow can be found from laboratory tests. 3. Some data are presented on the effect of filtration pressure and temperature on dynamic filtration. The work was limited to a small number of muds. No attempt was made to compare data on API filtrare loss with dynamic filtration rates. Previous workers have clearly shown that no correlation exists between static and dynamic filtration rates. EXPERIMENTAL ARRANGEMENT To establish whether the effect of mud flow conditions on dynamic filtrations is governed by the rate of shear at the cake surface, two different types of mud flow were used. In this connection the experiments were partially conducted by means of a porous pipe through which mud was circulated, partly in a dynamic filter apparatus which somewhat resembles a rotational viscometer. In both cases the flow conditions of the mud were well defined, thus making it possible to-calculate rates of shear at the cake surface. The filter pipe was 200 cm long with an internal diameter of 3 cm. Its perforated wall was covered on the inside with filter cloth on which the mud cake could form. Mud was pumped through the pipe at linear rates varying from 20 to 186 cm/sec. The required filtration pressure was obtained by means of a lack-pressure valve downstream of the filter pipe. In this arrangement the type of flow existing in the borehole is reasonably well simulated. The dynamic filter apparatus (Figs. I. and 2) consists mainly of two concentric cylinders with a clearance of about 0.5 cm. The outer cylinder (rotor) can be rotated at rates up to 1,000 rpm. The outer surface of the inner cylinder (radius about 2.5 cm) is covered with filter cloth or paper supported by copper gauze. The two cylinders are mounted inside a pressure

Jan 1, 1967

By D. de Fontaine

A convolution procedure is described for obtaining an analytical expression for the composition profiles of solute concentration:; in binary solid solutions. The free energy of solutions containing such compositional inhomogeneities can then be evaluated and the stability of spherically symmetric clusters (G.P. zones, for short) us a function of their radii and central compositions can be studied accordingly. Two characteristic zone radii are defined: the minimum radius (RM) corresponding to the smallest zone for which growth is energetically favored over dissolution, and the critical radius RK) corresponding to zones associated with the minimum work of formation. The values of RM and RK coincide at the two classical limits of nucleation and growth and spinodal decomposition. As the spinodal is approached, however, RM and RK diverge increasingly and it is proposed that the initial] decomposition morphology is governed more closely by the minimum, rather than by the critical radius since the former possesses no singularity within the m eta stable phase boundary. Reversion is briefly discussed in the light of the foregoing and it is shown that a set of reversion temperatures located between the coherent miscibility gap and the coherent spinodal must correspond to a given distribution of low temperature zones. SINCE the early works of Guinier1 and preston,2 the kinetics and morphology of GP zones have been studied extensively, see Ref. 3, for example, for a recent review. Nevertheless, a completely satisfying theoretical treatment of the decomposition of meta-stable or unstable binary solid solutions is still lacking due to the almost unsurmountable mathematical difficulties that one encounters in attempting to solve the pertinent diffusion equation which, in principle, describes the whole process. This equation, first derived by Cahn for the initial stages of spinodal decomposition,4 is a nonlinear partial differential equation of the fourth order in anisotropic space, a general analytical solution of which is completely out of the question. From this equation, some prediction can be made, however, concerning the later stages of spinodal decomposition5 and a recent numerical solution6 of the equation in the case of one- and two-dimensional composition modulations of a solid solution inside the spinodal has led to the following considerations: the sinusoidal morphology7 is to be expected only for solutions of average composition c, located near the center of the coherent miscibility gap.8 For solutions with co closer to the spinodal decomposition, depleted shell zones9 should constitute the first recognizable decomposition product, later evolving to iso- lated, enriched preprecipitates in a uniformly depleted matrix as this must, in all cases, constitute the final metastable morphology, barring, of course, loss of coherency and change in crystal structure. In agreement with the qualitative description of Bonfiglioli and Guinier,l0 it thus appears that the various models discussed by previous authors are, in fact, idealized successive stages of the general decomposition process. Unfortunately, corresponding calculations for a solid solution outside the spinodal are not yet available because of two additional difficulties: it is essential that the diffusion equation be solved in three-dimensional space, and furthermore, one must include in the equation an appropriate "source term" describing the statistics of the creation of composition fluctuations required to trigger the nucleation events. At present, a general treatment of this problem does not appear to be tractable. Nevertheless, by using a modelistic approach, i.e., by selecting a priori certain composition variation profiles and by evaluating their associated free energies, one can draw certain conclusions relative to the stability of these isolated composition modulations or "zones", both outside and inside the spinodal. This is the object of the present investigation. THEORY We shall now describe a standard four-parameter spherically symmetric GP zone, evaluate its free energy according to the Cahn and Hilliard11 formula (Eq. [5], below), and plot the corresponding energy surfaces as a function of zone amplitude and radius as

Jan 1, 1970

By R. E. Collins, Paul B. Crawford

Potentiometric model studies have been made of ertically fractured reservoirs. All fractures originated at he well and extended into the reservoir for various distances. Studies were made to determine the areal sweep efficiencies for line-drive patterns. The areal sweep efficiency for the unfractured system studied war 70.6 per cent. It was found that if the vertical fracture war paralIel with the direction of flood, the areal sweep efficiency may be considerably reduced, approaching zero in some instances. If the fracture was at right angles to the direction of flood, the areal sweep efficiency may be greater than that for the unfractured system. The sweep efficiency depends on the length and orientation of the fracture and direction of the flood. It is concluded that considerable effort should be made to determine the nature and direction of fractures before initiating a flooding program. Such information is required if the maximum sweep efficiency is to be obtained and channeling is to be reduced to a minimum. INTRODUCTION The recent development and application of fracturing techniques to petroleum reservoirs has served to focus considerable attention on the effect such fractures may have on secondary oil recoveries. The opinion has been expressed that fracturing may not alter appreciably the water flooding possibilities. Others have thought that channeling may occur and that the creation of fractures may serve to substantially reduce the possibilities of conducting a successful water flood in the reservoir The purpose of this study was to obtain quantitative estimates of the effect of vertical fractures on the areal sweep efficiencies of line-drive flooding patterns. In presenting this work on vertical fractures, it is not to be inferred that the present fracturing techniques result in producing a vertical fracture. The types of fracture which are created are not known with certainty in every case. Some are believed to be horizontal, others inclined along bedding planes, and some data indicate vertical fracturing. Each fracture will probably exert a different influence on the sweep efficiency of the secondary recovery program; consequently each fractured well should perhaps be considered as a separate problem, requiring a careful analysis. Preliminary work on the effect of fractures on producing capacity of a well indicated that for practical purposes many fractures may be treated as having a nearly infinite permeability This assumption appears to be valid in some instances where the permeability of the fracture is only 10 or 100 times greater than the permeability of the matrix. It is known to depend on the type, orientation and length of the fracture. In realization that fractures of only 0.01 in. width have calculated permeabilities near 5,000 darcys, it is believed that many fractures may be considered as of near infinite permeability without incurring serious error. Primarily for this reason this study was made with fracture analogs of very great permeability compared to the matrix. FRACTURED SYSTEM ANALOG The pattern sweep efficiencies reported here were obtained with the aid of a potentiometric model. The application of this type model for solving reservoir problems has been described in the literature.1,2,3 The electrolyte reservoir was 20 in. wide by 30 in. long by approximately % in. deep. The liquid depth was constant, which corresponded to a constant permeability-feet product from the injection to the production well. The electrolyte reservoir was built to correspond to a unit in the line-drive pattern shown in Fig. 1. This figure shows part of an infinite array of wells having alternate rows of fractured and unfractured wells. The area within the dashed rectangle shows the area generally studied in the model. In using the model to represent a fractured system, a thin strip of copper was soldered to the well or wells and the length and orientation of the strip was made to conform to the vertical fracture under consideration.

Jan 1, 1955

By Robert E. Green, Kenneth Reifsnider

A microscopic X-ray diffraction technique was employed for the simultaneous study of the behavior of several families of lattice planes, whose local orientatiotz changes are manufestations of internal distortions resulting from specimen deformation. Photographs taken using aluminum were single crystals 1/64 in. In diam are presented which illustrate the appearance of various types of deformation- or growth- induced structural variations. Slip and lattice bending were detected in tensile tests for strains as small as 0.1 pct. Evidence of deformation bands, lineage structure, and void formation wns also shown. AmONG the most important methods of studying internal stresses in metals is the utilization of various X-ray techniques. The general subject of X-ray metallography has been treated in a most excellent manner by Taylor,1 and the use of X-rays for the direct observation of imperfections in crystals has recently been reviewed in an article by Newkirk2 and in the book edited by Newkirk and Wernick.3 Since the inception of metallographic X-ray analysis, there have been numerous technological developments and refinements, each with its particular advantages and disadvantages and, acrording to these, its particular applications. X-ray diffraction techniques can be divided into two general categories, that of X-ray diffraction microscopy and that of X-ray diffraction macroscopy. In the former are included such methods as those of Berg-Barrett, Lang, and Schultz2,3 which are concerned with microscopic imperfections such as dislocations. The second category includes all work concerned with the study of the condition of a crystal lattice structure as indicated by the marroscopic detail in its diffraction patterns. An example of this is the study of the diffraction manifestations of inhomogeneous plastic deformation as manifested by formation of deformation bands and inhomogeneous lattice rotations. The method used in the present investigation was an X-ray diffraction macroscopy technique. The geometry and analytical aspects of the method were similar to those of Julien and Cullity.4,5 The source was collimated with a single slit in a manner simi- lar to the "parallel beam&apos;&apos; method of Barth and Hosemann6 except that a continuous spectrum of radiation was employed. The apparatus produced large-beam transmission Laue patterns of thin cylindrical aluminum-wire crystals, which were deformed either in tension or torsion. Although the conventional X-ray target and interposed slit were not capable of the high resolution obtainable with a fine-focus X-ray source as employed in a microfocus apparatus, resolution was a secondary consideration in our technique, and was sacrificed to obtain sufficient intensity along a sizable length of the specimen. By using a grid indexing technique similar to that of Barth,7 we were able to study several diffraction lines from various planes simultaneously. The results of these tests give additional information with regard to the inhomogeneous deformation of single crystals, and the indexing procedure used provides a correlation between portions of the test specimen and portions of a diffracted beam. EXPERIMENTAL TECHNIQUE The experimental arrangement used is shown schematically in Fig. 1. A source of continuous X radiation from a vertically mounted Machlett X-ray tube possessing a copper target and operated at 20 kv and 10 ma was passed down a horizontal lead tube 80 cm long which terminated in an adjustable vertical brass slit. The slit opening was adjusted to correspond with the diameter of the test specimen. In front of the brass slit it was possible to place a screen wire mesh, one horizontal strand of which was a lead wire which served as a fiduciary mark, the screen itself serving as an indexing system as will be described later. Behind the screen was the test specimen which consisted of an aluminum single crystal wire 1/64 in. in diameter and 3 in. long. These 3-in.-long crystals were cut from

Jan 1, 1965