Search Documents

By R. M. Rose, K. A. Jones

Resistive measurements ulere made on superconducting niobium single crystals in transverse magnetic fields. Crystals were grouln in both high and ulfrahigh vacua, doped with stnall quantities of oxygen, and annealed in ultrahigh vacua at temperatures from 900' to 1350°C. Maxima in the fluxoid pinr~ir~g force were observedfor the [100] and [210] directions. The increase of the critical current density with dissoved oxygen was seen to saturate at relatively low concentration. This behavior, and the variation of J, and residual resistivity (as a function of oxygen content and annealing) of both clean and dirty gassy) crystals, is seen as a consequence of the formation of clusters with (210)~b habit planes and (200)ffo common boundaries. These suboxide clusters are apparenlly stable at high temperatures and very low oxygen concentrations. It is generally known that dissolved oxygen at concentrations well below the solubility limit sharply increases the critical current density of superconducting niobium.' Such a result would not be expected if Nb-0 solutions were random, not even as the result of segregation on the scale usually encountered in solid solutions. According to the Kim-Anderson theory of flux creep,' the critical transport current density is determined by the pinning force Fp associated with some feature of the microstructure. As the Lorentz force. a or (J X R), of the transport current in the magnetic field approaches Fp, flux creep and then flux flow occur. The pinning force is expected to be quite small if the heterogeneity involved is significantly smaller than tte range of coherency which for niobium is circa 250~.~ Although the strain field of an individual oxygen atom in the niobium lattice is relatively large, it does not seem plausible that individual atoms could be so effective as to account for the data. Also, ~ramer~ found that straining and strain aging specimens containing small amounts of oxygen on occasion decreased the residual resistivity. A possible explanation for the above behavior may be found in the work of van Ooijen and van der Goot,= who found that annealing changed J, for niobium foils containing oxygen, and they attributed the changes to ordering of the oxygen atoms. This phenomenon has been observed by other investigators using X-ray and PROCEDURE Cylindrical 4-in.-diam niobium single crystals with axial orientations [ill], [110], [112], and [oo~] were grown from pure niobium, see Table I for analysis, in an ion-pumped ultrahigh-vacuum system using the electron beam floating zone technique. Total pressure during growth was lo-'9 torr: the oxygen partial pressure was less than lo-" torr, as determined by a mass spectrometer. These crystals will be referred to as the "clean" crystals. Similar crystals were also grown in another system at a pressure of 5 X 10"6 torr. These crystals will be referred to as the "dirty" crystals. Specimens 14 in. long were cut from the as-grown crystals and one end was polished down to no. 600 carbide paper smoothness. A back-reflection Laue photograph was then taken in order to fix the crystallo-graphic directions transverse to the crystal. The critical current density J,, the current density necessary to produce a voltage drop of 1.0 pv between voltage leads placed a in. apart, was determined as a function of the magnitude and direction of the transverse magnetic field and of the oxygen concentration. The technique and experimental arrangement was essentially similar to that of Tedmon et al.' In the remainder of this paper, the azjmuthal angle $ for the field H is referenced to the [ I101 direction in all four crystals, and oxygen was added in small equal increments by anodization and annealing.7 The resistivity ratio, RR, of the resistance at room ternperature to that at 4.2"K (H =7.00 kG) was determined as well.

Jan 1, 1970

By Harold A. Watson

A method for measuring leakage of vapors into a chemical-cartridge half-mask facepiece during simulated use is described. The respirator, fitted with high efficiency cartridges, is worn in a dichlorodifluoromethane (CC12F2)--air test mixture while a series of timed activities (facial movements, varied breathing, exercise, etc.) is performed. Complementing other tightness tests, this method is useful particularly for studying face fit and the effect of facial movements on leakage. The concentration of CC12F2 within the facepiece is monitored continuously to observe the immediate effect of physical activity on leakage. A composite sample is collected simultaneously in a plastic bag as a measure of average leakage during the 18-minute test period. Leakage, reported as percent, is calculated from the concentrations of CC12F2 within the facepiece and in the test atmosphere.

Jan 1, 1970

A new test coke oven of 7 cu. ft. capacity is described. The oven was designed to simulate conditions in a central section of a battery oven and a special feature is the achievement of unidirectional heating from theflues on either side of the charge.The method of charging permits close control of charge conditions such as bulk density and moisture content.These conditions, as well as coking conditions (e.g. flue temperatures and coking time) must be varied to simulate conditions in different batteries. The effects of several variables on coke quality are discussed.Using test oven conditions corresponding to those of Newcastle steelworks batteries, good agreement was obtained between strength indices of test oven and battery cokes. For 17 tests, the average ASTM stability and hardness indices for the test oven coke were 38· 5 per cent and 69·7 per cent, while those for the battery coke were 37· 4 per cent and 67·6 per cent.Observed relationships between coke yield and productivity, and volatile matter and charge moisture are given.The use of a test oven to classify test coals into prime-, medium-, poor-, and non-coking coals is discussed.This may be achieved by observing the effects of the test coal in blends with other coals. The results of a series of blending tests are presented.The dangers of attempting to predict coke quality from a single coal property are pointed out, since this can be done only for very restricted ranges of coals.A diagram is given showing the established trends in the effects of blending various types of Australian coals on coke strength.

Jan 1, 1970

Jan 1, 1970

By W. G. Davenport, M. Paneni

Experiments in sheets of mercury and water have confirmed that the dynamics of bubble rise are simi-lur in Liquid metal and aqueous systems. RECOGNITION that bubble/metal interactions are important in pyrometallurgical refining processes (e.g., steelmaking, vacuum degassing, deoxidation of copper) has led to a number of investigations into the dynamic behavior of bubbles rising in liquid metals. Terminal velocities of bubbles rising in mercury (20°C)' and silver (1000" to 1050"c)1,2 have been measured. Good agreement between aqueous and liquid metal systems has been obtained. Shapes of air bubbles rising in mercury have been determined experimentally by Davenport, Bradshaw, and Richardsonl who used a system of electric probes to measure maximum horizontal and vertical dimensions, Fig. 5. In the present work the motion and shape of nitrogen bubbles rising in a thin sheet of mercury (i.e., a "two-dimensional" environment) have been studied with the view of confirming dynamic similarity between aqueous and metallic systems. The use of thin sheets of liquid permitted direct visual and photographic comparisons of bubble behavior in clear liquids and in mercury. EXPERIMENTAL METHOD The thin sheets of mercury or water were held in a clear acrylic plastic container, Fig. 1. The sheets of liquid were 80 cm high by 30 cm wide and 0.47 cm thick. Single bubbles were formed by injecting the desired volume of nitrogen into an inverted semicircular cup (stainless steel). The cup was then rotated into the upright position using an external handle thus releasing a single bubble of measured volume. Nitrogen was injected into the cup from a nitrogen cylinder through the external handle and a hollow rotating axle. The volume of a bubble (volume range 0.01 to 10 cu cm) at any instant was obtained by connecting a soap film meter to a small gas space above the liquid. By this technique any change in the volume of a gas bubble in the liquid (i.e., during cup filling or during rise) produces an equivalent displacement from gas space A into the buret B. This displacement is determined by measuring the movement of a horizontal soap film in the buret. Bubble shapes were determined using single-frame photography. Rising velocities were determined by cinephotography with a 1/100-sec stopwatch in the field of view. In both cases dimensions were referenced to a rule attached to the container surface, also in the field of view. Linear dimensions were measured on projections of the single-frame negatives. Rising velocities were calculated from measured heights of rise (10 to 20 cm) and times of rise. The effects of disturbances caused by cup rotation and by bubble release were minimized by photographing bubbles at a point not less than 25 cm above the cup. MATERIALS Mercury: Technical-grade purified: minimum 99 pct Hg; Johnson Matthey and Mallory Limited, Toronto. Nitrogen: Technical-grade purified nitrogen: 99.99 pct Nz, 0.01 pct 02, trace argon; Air Liquide Canada Limitee, Montreal.

Jan 1, 1970

Jan 1, 1970

By E. B. Brauer, W. Hurst, J. D. Clark

Because of drilling, completion, and workover practices, the permeability around a wellbore generally is different from the permeability of the formation. The zone with the altered permeability is called "skin" and its effect on the pressure or flow behavior of the well is called "skin effect". Hawkinsl has shown that the radius and the permeability of this zone are related to the skin by: If the permeability in the skin zone is less than that of the formation, the skin is positive; if it is more than that of the formation, the skin is negative. If the two permeabilities are equal, s is zero; that is, there is no skin. Van Everdingen2 and Hurst3 have given mathematical solutions to the case of a zone of reduced permeability around the wellbore. This skin effect is illustrated in Fig. 1. They treat the positive skin as a zone of reduced permeability of infinitesimal thickness around the wellbore. When applied to a well with negative skin, however, their solutions lead to the calculated flowing well pressure which is smaller than the formation pressure, i.e., an injection situation. This is a physical contradiction. In this paper, it is shown that this mathematical difficulty can be overcome by assuming an effective wellbore radius larger than the actual wellbore. Existing solutions2"4 are modified to include the effect of an enlarged wellbore radius, and thus enable the engineer to deal with negative as well as positive skins. Theoretical Development Van Everdingen and Hurst,4 through the use of Laplace transformation, obtained the following expression for transient fluid influx into a wellbore: The values for QtD, reported extensively by van Everdingen and Hurst,4 are obtained from the following integral:

Jan 1, 1970

A rotating disc geometry was used to provide a uniformly accessible surface for a kinetic study of copper cementation with zinc in sulphate solution. For copper concentrations of the order of 10 p.p.m., a porous, adherent deposit is formed on the zinc surface and this has a significant influence on the reaction rate. When the specific mass of deposit is less than about O· 35 mg/cm2 the surface may be considered as smooth and the rate is described by simple convective counter diffusion of zinc and copper ions through a laminar diffusion boundary layer to the deposit-solution interface. In the presence of additional deposit, the diffusional mass transfer rate is increased by increased surface roughness.The rate under these conditions is first order, with the degree of enhancement independent of further increases in the specific deposit mass. The enhancement in rateis, however, directly related to the initial mass flux and the species involved in the reaction.Preliminary work showed the same mechanism to apply to cementation from chloride solution, solutions of high zinc concentration, and cementation of copper with iron..

Jan 1, 1970

A combined reflection seismic and magnetic survey off the south coast of New South Wales has revealed six distinct reflection horizons within the Lower Triassic-Upper Permian sequence, and has enabled the upper part of the lllawarra Coal Measures to be mapped in considerable detail. Large magnetic anomalies in the eastern half of the survey area appear to originate in the Gerringong Volcanic Sequence which forms the magnetic basement in the survey. Small residual anomalies in the western half of the area can be correlated with small scale faulting and folding in the upper part of the llawarra Coal Measures, as revealed by the seismic reflection survey.A synclinal axis striking at approximately 55° has been mapped near the western margin of the survey area, and areas of faulting and folding have been outlined by combined analysis of seismic and magnetic results.

Jan 1, 1970

Plant based on the New Jersey refluxing process was installed at the Cockle Creek Works of Sulphide Corporation Pty. Ltd. to refine portion of the zinc output from the Imperial Smelting furnace.Refining of impure zinc is effected by distillation in two stages in vertical columns comprising a series of superimposed silicon carbide trays designed such that ascending zinc vapours and descending liquid metal are in intimate contact. In the first stage, portion of the zinc and its cadmium content are distilled from impure feed metal, condensed, and while still molten, fed to the second column in which zinc is boiled under controlled refluxing conditions to remove the cadmium content.The resultant zinc has a purity in excess of 99·995 per cent zinc.The availability of refined zinc has made possible the production of a number of zinc alloys for marketing in addition to the basic grades of zinc.

Jan 1, 1970

By R. M. Sellwood

The No. 4 I.S.F. Smelter Complex at the Avonmouth Works of Imperial Smelting Corporation (N.S.C.) Limited commenced operation at the beginning of 1968. The furnace rating is 120,000 tons zinc and 40,000 tons lead per annum. Zinc and lead concentrates are unloaded onto a 1,000 tons per hour capacity conveyor belt and automatically weighed and sampled before receipt into a 48,000 ton capacity raw material store. Distribution is by a reversible boom stacker into 16 x 3,000 ton bays and recovery is by a front end loader. Materials are fed to 16 proportioning bins for feeding via continuous weigh feeders to a 132 square metre effective updraught area sinter machine with a material flow rate of 240 tons per hour over the machine. The required quantity of output lump sinter at up to 800ºC is fed to 2 x 600 ton storage bins on the I.S.F. plant. By-product SO2 rich gas is fed to a conventional wet gas single stream acid plant rated at 600 tons per day utilizing a dry N.T.P. gas volume of 100,000 cubic metres per hour. The hot sinter and preheated coke (800ºC) are automatically neighed and charged to a 27.3m2 shaft area I.S.F. fitted with 2 x 5.5 m. wide lead splash condensers. Grade 4 zinc is tapped and transferred to the zinc refinery. Conventional type lead bullion and slag are tapped from the bottom of the furnace through a separation forehearth. The refinery consists of 6 lead and 3 cadmium columns to produce some 52,000 tons per annum of S.H.G. metal. S.H.G. metal and Grade 4 metal is cast into various shapes as required.

Jan 1, 1970

By W. T. Malone, J. R. Williams, R. L. Tiner, J. M. Tinsley

Hydraulic fracturing methods for production stimulation have become a common procedure in the oil and gas industry. Fracturing treatments are performed on wells of various potentials to help increase the production rate. In addition to the desire for increased production is the need to predict what increase might be expected. This knowledge is useful in economical treatment planning to accomplish the desired production goals for the well. Theoretical production increases due to horizontal radial fractures may be calculated using equations developed from mathematical analyses and electrical models.1,2,6 It also is possible to predict theoretical production increases due to vertical fractures using curves developed from electrical3 and mathematical models.8 However, these studies assumed that the vertical fracture height and the formation height were equal, which may not be true." In addition, the final fracture height through which fluid is produced may not be equal to the created height. The development of vertical fracture proppant transport relationships has made unrealistic the assumption that propped fracture height is equal to created fracture height. Through the use of published proppant transport equations,4'12 it is possible to calculate the height of a bed of propping agent deposited in a fracture of a given height. Assuming that the fluid flowing in the fracture moves through only the propped portion of that fracture (i.e., that the un-propped part of the fracture heals) it becomes essential to know how production increase varies with the ratio of propped fracture height to formation height. It should be recognized that future investigation may determine that the fracture does not heal completely and that a highly permeable channel exists over the deposited bed of proppant. The effect of this channel was not considered for the present study. The following analysis of experimental data allowed development of curves to determine theoretical production increases considering fracture height as a variable. Discussion In designing a fracturing stimulation treatment a number of variables must be considered. The size of the treatment, type of proppant used, proppant concentration, treatment fluid properties and well conditions are all factors affecting the dimensions and permeability of the fracture. Production increase due to a vertical fracture is affected by propped length, height, width and permeability of the fracture for a given set of formation conditions. Production increase may be defined as the ratio of production rate after fracturing to the production rate before fracturing at the same pressure differential, pe — pw. The producing rate in an

Jan 1, 1970

By J. J. Jonas, G. Gagnon

SHG zinc was extruded in the temperature range 110" to 350°C and the strain rate range 0.05 to 5 sec-1 The strain rate/flow stress/temperature results were analyzed using a power sinh stress relationship. When temperature changes due to the heat of deformation and to heat losses are neglected, the exponent of the sinh function of the stress is 5.6, and the apparent activation energy of deformation is 28 kcal per mole. When these changes are taken into account, the exponent decreases to 4.7 and the activation energy to 23 + 2 kcal per mole. The corrected stress exponent and activation energy are in very good agreement with published values obtained from creep experiments, and suggest that the hot extrusion of zinc is controlled by a mechanism involving self-diffusion. When the extrusion and creep data are compared using a Zener-Hollomon parameter and either a sinh or an exponential stress term, an appreciable offset is observed, which may be due to the difference in impurity content. For a given set of extrusion conditions the ram speed, maximum pressure, and initial temperature can also be correlated using a Zener-Hollomon parameter and a sinh pressure term. THE deformation of metals at temperatures over about one-half the absolute melting temperature has been extensively studied at creep strain rates. By contrast, relatively little work has been carried out on the behavior of metals at hot working strain rates. Most of the latter investigations have been performed using simulative methods, such as hot torsion and hot compression, in which the friction conditions and temperature rise during deformation may differ appreciably from those existing under industrial conditions. Recently, however, Wong and Jonas1 used a scaled-down industrial process to determine the stress and temperature dependence of the strain rate during the extrusion of aluminum. In such tests, the effects of friction and adiabatic heating are closer to those produced in industrial operations. Also, with regard to the testing of materials of limited ductility such as zinc, hot compression and hot torsion do not permit the attainment of true strains as large as the deformations usually applied commercially. The present study was undertaken to investigate the extrusion behavior of Special High Grade (SHG) zinc. The detailed objectives were: 1) to establish the stress and temperature dependence of the strain rate with and without a consideration of adiabatic heating, 2) to study the pressure and temperature dependence of the ram speed, and 3) to investigate the microstruc- tural changes occurring during the deformation. The last aspect of the investigation will be covered in a separate paper. The treatment described below differs from that of Wong and Jonas in that the adiabatic temperature rise during deformation is taken into account, and the calculation of the mean equivalent strain rate is based on the redundant as well as on the homogeneous work. EXPERIMENTAL PROCEDURE Rods from two shipments of continuously cast SHG zinc* were used in the investigation. The composition *Supplied by courtesy of Cominco Ltd. range of the impurities present, as given by the supplier, was: Pb: 0.0013 to 0.0023 pct, Fe: 0.001 pct max Cd: 0.0001 to 0.0006 pct, Cu: 0.0002 to 0.0005 pct, Ti: 0.0001 pct max; thus, by balance, zinc valued from 99.9963 to 99.9966 pct. The as-received rods were machined into billets having a nominal diameter of 1.56 in. and a nominal length of 1½ in.; longer billets up to 4 in. in length, were also used to investigate special aspects. The as-machined rods were annealed at 400°C for 24 hr and slowly cooled. This treatment produced a columnar grain structure, with a grain size of about $ by 2½ cm which was appreciably larger than the as-cast one. A 150-ton, direct extrusion, vertical press was used. Ram displacement and force were recorded continuously against time. A constant flow control valve permitted the maintenance of a range of preselected ram speeds up to in. per sec. The selected speed was held constant, irrespective of the required force, as long as the load developed was below the maximum available. Strain gages were used to determine the force; the gages were calibrated before and after each testing period with a 200,000-lb capacity load cell. Further details of the experimental equipment can be found in Ref. 2. The billets were preheated for 90 min in the extrusion container, which was well insulated so as to minimize temperature gradients. This period was sufficient for the billet to reach a uniform temperature at all temperatures between 110" and 350°C. A square-shoulder die having a 0.290-in. diam central hole was used. The extrusion ratio was 30 to 1. This is equivalent to a reduction in area of 96.7 pct, an elongation of 2900 pct, and a true strain of 3.4. The ram speed was varied over two orders of magnitude from 0.0027 to 0.26 in. per sec. The ram was water-cooled during most of the tests, although some experiments were conducted with a preheated, uncooled ram. All extrusions were run without lubricant, which resulted in conditions of sticking friction. EXPERIMENTAL RESULTS Stress Dependence of the Strain Rate Neglecting Adiabatic Heating. The maximum force required to extrude is given in Table 1 for each of the five initial

Jan 1, 1970

By R. G. Aspden, D. A. Colling

THE {120}(001) texture has been reported in near 50 pct Ni-Fe alloys <0.004 in. thick1,2 and in copper 0.014 in. thick.3 This texture is one of two types of secondary recrystallization textures observed to replace a cube texture matrix and is desirable in Fe-Ni alloys because of its unique magnetic characteristics.1,2 However, the conditions for and mechanism of the {120}(001) formation are not well understood.3 During an investigation of the influence of alloying elements on the magnetic characteristics of 49 pct Ni-Fe alloys with a sheet thickness of 0.014 in., we have observed the growth of {120}(001) secondaries when strong deoxidizers such as aluminum and calcium were present. It is believed that this observation in unique and provides additional insight with regard to the mechanism for {120}(001) texture formation. Three alloys were prepared with electrolytic iron and nickel and contained a nominal 49 wt pct Ni. Alloy A was a straight binary Fe-Ni alloy, whereas intentional additions of 0.013 wt pct Ca to alloy B and of 0.59 wt pct A1 to alloy C were made. Ingots were forged at 1100 C, hot-rolled to 0.150 in. gage, pickled, and cold-rolled to 0.014 in. Standard U-shaped magnetic test specimens were punched, deburred, sprayed with MgO, and annealed in ~50°C dew point hydrogen. Textures were identified with the aid of etch pits, the

Jan 1, 1970

By Richard A. Dick

Underground blasting in permanently frozen gravel was studied in a 7-by 12-ft drift by running a nonreplicated 2 X 2 X 2 factorial experiment on the effects of type of cut, delay, and explosive on the fragmentation subsystem. Comparisons were made between a V-cut round and a burn-cut round, millisecond delays and half-second delays, and 40 percent special gelatin and 60 percent high-density ammonia dynamite. The number of large blocks in the muck pile) loading time, shuttle car operating time lost because of cleanup needs, length of advance per round, and length of muck pile were recorded and statistically analyzed. None of the blast parameters had a significant effect on the length of advance or length of muck pile, although the half-second delays tended to give a shorter muck pile than millisecond delays. The burn cut required less secondary breakage, less loading and cleanup time, and gave better overall fragmentation than the V-cut. The ammonia dynamite shots required less secondary breakage and less loading time than those using special gelatin. The half-second delays gave better overall fragmentation and required less cleanup time, but the millisecond delays gave a better loading rate. None of the differences between delays and explosives were notably significant.

Jan 1, 1970

By R. J. Watson, A. H. Juch

About half a million sucker-rod pumps are installed in oil wells in the U. S. alone. In Venezuela, too, the system is widely used; some 5,000, or 90 percent, of Cia. Shell de Venezuela&apos;s wells produce a total of 400,000 BOPD by this method. Most of the crude pumped is viscous, gassy and contains relatively small amounts of free water. API gravities vary from 90 to 16O; dead oil viscosities vary from 10 to 10,000 cp at the pump; and GOR&apos;s range from 100 to 2,000 cu ft/bbl. Pumps are set at depths ranging from a few hundred to 6,000 ft; the average setting depth is about 2,500 ft. From the literature, it appears that although recent years have seen significant progress toward understanding the rod pumping system as a whole, relatively little attention has been paid to those key items, the sucker-rod pump and the subsurface gas separator or gas anchor. Perhaps because of proration in the U. S., pump manufacturers have seemingly concerned themselves more with the durability of their product, achieved mainly through robust design and metallurgical know-how, than with its working efficiency. In 1963, CSV commenced a major rehabilitation and expansion of its heavy oil fields along the Bolivar Coast. The advent of the more efficient gravel-packed completion and of large-scale thermal recovery (mainly steam soak) operations resulted in well potentials *Now with B. I. P. M., The Hague, The Netherlands. greater than pump capacities. Hence, an urgent need arose to improve over-all pump efficiencies and to attack those special problems related to the production of both the extremely viscous and the hot, steam-laden crudes that formed an increasing proportion of the total stream. There followed a detailed investigation and testing program, carried out in close cooperation with a major pump manufacturer. This culminated in the introduction to the oil industry of a new line of sucker-rod pumps. Early Investigations and Experiments Field Tests Extensive field tests were carried out in three primary Bolivar Coast wells with an in-situ dead oil viscosity ranging from 400 to 6,000 cp. In these installations provision was made to (1) monitor bottom-hole pressures below, in and above the pump by means of pressure transducers, and record casing- and tubing-head pressures continuously; (2) measure accurately the annulus gas, as well as measure the tubing gas and liquid production from the well; (3) take special PVT samples of the well fluid; (4) measure power consumption of the prime mover; and (5) measure the rod loads at the polished rod. A typical test installation is illustrated in Fig. 1. It was established from the pressure transducers that with viscous oils, there is considerable flow resistance across the pump and gas anchor assembly. Pressure surges were observed over each pump cycle,

Jan 1, 1970

By J. Douglas, W. T. Ford, G. E. Henderson, I. F. Roebuck

An implicit numerical method is presented for simulating the differential and algebraic relations governing one-dimensional three-phase flow in porous media. The method is based upon compositional representation of the hydrocarbon system. Variable physical properties and water-oil capillary forces are included in the formulation; however, the effect of gravity is ignored. The effects of changing composition and mass transfer are considered through the use of known phase behavior concepts and correlations. The differential and algebraic equations and the numeric approximations to these equations are presented, and the computing algorithm is discussed in detail. This method is used to simulate a single-phase, two hydrocarbon-component gas displacement laboratory experiment. The method is also applied in simulating a solution gas drive and a gas injection problem using a reservoir oil represented as a nine-component mixture. For comparative purposes, the latter two problems also are examined with a one-dimensional volumetric model. The results of the simulation of the laboratory experiment demonstrate the relevance of this method to problems wherein the effects of mass transfer significantly enter the displacement process. The comparison of the results of the volumetric and compositional methods applied to a reservoir system illustrates not only the advantages and benefits derived from examining the hydrocarbon system composition, but also the necessity of this approach in application to a general reservoir problem. INTRODUCTION The advent of large memory, high-speed digital computers enabled reservoir engineers and mathematicians to pool their knowledge in the development of sophisticated techniques for the prediction of oil and gas reservoir performance. These techniques have progressed from volumetric and trend analysis calculations to the development of "mathematical models", or reservoir simulators. These latter models, in general, utilize finite difference approximations to the rather complex partial differential equations that mathematically describe the physics and thermodynamics of fluid flow in porous media. Many forms of this "model" type are available today, and much of the recent literature has been devoted to either improving known methods or developing new techniques for use in these prediction methods. The more advanced techniques in use at this time consist of the solution of differential and algebraic relations that describe a reservoir geometry in one, two or three dimensions and that incorporate reservoir hydrocarbon and water fluid systems that obey specified physical laws, having properties that might easily be measured in the laboratory. This type of model utilizes formation volume factors and other volumetric data, expressed directly as functions of pressure, to represent the reservoir fluid system. Experience gained with this type of simulator, commonly known as a "beta-factor model", has indicated that the approach has many outstanding defects, most of them due to the assumption that those properties measured in the laboratory completely describe the behavior of the fluid system throughout the life of the project under study. The effects of varying hydrocarbon composition and mass transfer between the phases of this hydrocarbon system must, by definition, be ignored. These defects, which may often be ignored in

Jan 1, 1970

By R. G. Bentsen, D. A. T. Donohue

The cyclic steam injection process has become the most widely applied and most successful thermal recovery technique in use today.1-4 Normally, steam stimulation is repeated several times during the life of a project. The problem is to select the optimum number and size of steam treatments so that the profits realized over the entire life of the project are maximum. Obviously, operating costs rise directly with the frequency and amount of steam injection. Operating efficiency, on the other hand, can be adversely affected if the amount and frequency of stimulation are insufficient. Thus arises a multistage decision problem; that is, a decision must be made at each sub-period (day, for example), as to whether steam should be injected, and, if so, how much. It is desirable that these decisions be made in such a way as to secure maximum net profit before taxes. This can be achieved by dynamic programming. To use dynamic programming, we must be able to predict not only the production rate-time curve resulting from a given stimulation, but also how changes in stimulation policy will affect this curve. In this study, the effect of steam injection on well performance is estimated on the basis of a simple model that takes into account both steady-state and transient flow. Oil production rate is calculated for both types of flow, and the larger of the two values is used. This allows the calculation to change to radial steady-state flow at such time as this mechanism becomes more effective than the transient mechanism. We propose to show how dynamic programming can be used to optimize the steam soak process with respect to net profit. For this, a mathematical representation of the steam soak technique is necessary. A further objective, therefore, is the formulation of a mathematical model capable of realistically simulating the physical process. Oil Production Response Model The advantage of using a mathematical model is its generality and ease of manipulation. However, the model must be an accurate representation of the physical process, for a solution derived from a model can be no better than the model itself. The question now arises as to how well the model should be expected to simulate reality. In general, the simplest model should be tried first, for additional accuracy is apt to require additional cost and time. Also, the search for a high degree of accuracy may yield only greater complexity and difficulty without yielding significantly better results. With this in mind, we kept our model as simple as possible, yet generally representative of the actual physical process. The number of state variables is limited to three, and stochastic variations in the parameters are ignored. As a consequence, use of this model should be restricted to cases in which the underlying assumptions are not severely violated. In

Jan 1, 1970

By Robert E. Wood

The temperature dependency of the low-density viscosity coefficients, ?°T, of helium and of nitrogen, 100° to 1,000° K, is correlated with empirical equations, and the temperature and composition dependencies of the low-density viscosity coefficients, ?°T, of helium-nitrogen mixtures are correlated with the Chapman-Enskog expressions by using the Lennard-Jones (6:12) potential function. The residual viscosity, ?T,P - ?°T (the difference between the viscosity of a compressed gas and the dilute gas at a given temperature), was found to be a function of the thermal pressure coefficient, (?P/?T)V, and two parameters, a and 1, which are characteristic of gas composition. The equation ?T, P = ?°T + a [(?P/?T)V]ß represented 1,340 higher pressure experimental viscosity values with a mean absolute deviation of 0.86 percent, and this equation was used to compute viscosities of the gaseous helium-nitrogen system for a range of temperatures and pressures including those encountered in Bureau of Mines helium purification processes and in a crude-helium conservation gas reservoir. Tabular viscosities for helium, nitrogen, and 13 helium-nitrogen mixtures are presented for 110 temperatures, 133° to 740&apos; K, and 49 pressures in the range 1 to 240 atmospheres. The accuracy of the tabulated viscosities varies because correlation parameters were obtained from experimental data of unequal reliability. It is estimated that uncertainties in the computed viscosity values are ±5 percent for the region 325° to 740° K, ±2 percent for the region 183° to 325° K, and ±5 percent for temperatures below 183° K. The latter uncertainty may rise to y10 percent as critical conditions are approached.

Jan 1, 1970

An assessment of methods currently available for the accurate determination of low concentrations of platinum in ores and mineral concentrates indicated that prior treatment of the sample by fire assay procedures is necessary. To enable rapid determinations to be carried out, instrumental methods applied to the metallic bead produced from the fire assay treatment have been evaluated, considering maximum flexibility for concentration variations, minimum chemical operations, and adequate sensitivity.Atomic absorption spectrometry proved suitable, giving solution determination limits of 0&#183;4 p.p.m. Pt, 0&#183;1 p.p.m. Pt, or 0&#183;06 p.p.m. Pt under various extraction conditions using aqueous and organic agents. Thesedetermination limits related to 1 assay ton of sample (32&#183;67 g) produced easily readable and accurate readings within the range of platinum contents (0&#183;03-2&#183;0 p.p.m.) of economic interest in deposit evaluation and prospecting. The technique shows promise of also being suitable for determination of palladium and gold.

Jan 1, 1970