Search Documents

By E. L. Dougherty, J. W. Sheldon

Using the numerical techniques shown in this paper it is possible to compute the simultaneous dynamical behavior of multiple fluid-fluid interfaces in two dimensions. Hence, fluid-fluid inter face models of several oil recovery processes can be constructed which allow prediction of the effect of fluid flow in two dimensions on the behavior of these processes while taking into account various physical effects, such as saturation gradients, phase changes and thermal stimulation. Techniques for constructing fluid-fluid interface approximations for several oil recovery processes are reuiewed. The validity of these techniques is established by comparison to experimental results. The availability of a computer program for computing the behavior of multiple fluid-fluid interface problems makes the use of potentiometric models uneconomic. The results indicate the areal sweep efficiencies obtained from potentiometric model studies at high mobility ratios may be somewhat optimistic. The results also indicate that the areal sweep efficiency curves obtained at high viscosity ratios 2—4 using miscible fluids in porous media are misleading and that considerable care must be taken in defining the mobility ratio for a displacement process, especially for miscible fluids. INTRODUCTION Many of the methods used to estimate areal sweep efficiency assume that the oil recovery process can be represented by one or more fluid-fluid interfaces. When the mobilities of the displaced and displacing phases are assumed equal, this quantity can be computed using analytical mathematical methods.l.2 For non-unit mobility ratio, experimental methods based upon the fluid-fluid interface concept include potentiometric and electrolytic blotter models. Scaled models of porous media have also been used to study this problem. Fay and Prats3 and Aronofsky4 used numerical techniques to compute breakthrough sweep efficiency, but neither of these studies were extended beyond breakthrough. Numerical techniques which have been developed for solving a general fluid-fluid interface problem are presented in a companion paper.5 These techniques have been incorporated into a program for a large-scale digital computer which can treat simultaneously the dynamical behavior of as many as six fluid-fluid interfaces. Thus it is possible to simulate the behavior of several oil recovery processes taking into account such effects as saturation gradients. Techniques for representing several oil recovery processes by fluid-fluid interfaces are presented. Processes considered are conventional water flooding with and without a mobile gas phase, hot water flooding, miscible flooding and enriched gas drive. Other processes, whose effects could be studied by this technique, but are not considered in detail here, are underground combustion and steam injection processes and various slug processes such as gas-water injection and alcohol-solvent processes. Theoretical results obtained with single interfaces for mobility ratios of 0.25, 1, 2, 4, 8, 16, 32 and 64 in a repeated five-spot are compared to experimental results for this geometrical configuration obtained using miscible fluids and conductive analogs. Predicted recovery curves for water flooding in a repeated five-spot are compared to experimental results reported in the literature. Example calculations are shown for a conventional water flood in the presence of a free gas phase and for a hot-water flood. MATHEMATICAL PRELIMINARIES We assume that the rock matrix is homogeneous with constant thickness, absolute permeability and porosity. The fluids are incompressible. We also neglect the effects of gravity.

Jan 1, 1965

By L. G. Austin

My first contact with industrial milling was during the time I worked in the electricity generating industry in the United Kingdom. In visits to power stations to investigate either deposits in the boiler furnaces or polluting deposits settling around the stacks. I had to check the performance of the vertical coal pulverizers, since poor pulverization aggravated both problems. Naturally, then, when I came to the USA in 1957 to take a PhD in fuel technology at Penn State, I was put to work to review the science of coal pulverization. After this reviewing, I was completely confused. On one hand, there was a well-developed understanding of stress-strain equations, and a rapidly developing knowledge of how stressed, brittle solids fractured, based on the Grifith crack theory. On the other hand, reading in the grinding literature gave me: -Kick's Law, which was clearly not correct in the light of modem fracture theory; -Rittinger's Law, which was also clearly not correct; -Bond's Third Law of Comminution, which was claimed to have something to do with the Grifith crack theory, but where the connection between the two was made by intuitive pseudo-scientific reasoning I could not accept; -the choice of mill motor power for the most common type of coal mill, the Raymond pulverizer, was calculated from the fan power required to move air through the mill. Although I could accept the empirical connection between the two, it made no sense from the point of view of fracture energy. Even today, most books or review chapters on size reduction start from these laws. Incorrect statements abound in the literature, such as "the Hardgrove Index is based on Rittinger's Law," which it is not, "The Bond theory states that work input is proportional to new crack tip length produced in particle breakage," which is not true, etc. My own test work showed that these "laws" did not fit the data for grinding of coal. At about this time, Epstein (1948) and Broadbent and Callcott (1956), following the original work by R.L. Brown (1941) at the British Coal Utilization Research Association, proposed describing breakage as a series of fracture stages. I took their concepts and developed the basic differential equation for a batch grinding process continuous in time, analogous to a batch chemical reactor. Robin Gardner then joined the project and did his PhD on treating batch grinding in the same way as a batch chemical reactor. He found that the basic equation had already been partially derived by Sedlatschek and Bass (1953) in Ger- many. We confirmed experimentally the validity of the equations for describing batch grinding (1962) and formulated the equation describing steady-state continuous grinding in a fully-mixed mill. At about the time this work was published, Gaudin and Meloy (1962) and Filippov (1961) independently published essentially the same equations, but without experimental proof of the validity of the concepts. I will give a brief overview of what these beginnings had led to in the design of mills for size and power, and show some of the results of this more detailed understanding of grinding processes. Concepts of Fracture Mills such as tumbling ball, rod, pebble and autogenous mills and vertical mills such as the Raymond, and E-type, apply compressive stress to lumps or particles relatively slowly. Compressive stress applied to a particle of an elastic brittle solid imparts overall strain energy to the solid and produces local regions of tensile stress, [(Fig. 1)] (Berenbaum and Brodie, 1959).

Jan 1, 1985

By M. King Hubbert, David G. Willis

A theoretical examination of the fracturing of rocks by means of pressure applied in boreholes leads to the conclusion that, regardless of whether the fracturing fluid be of the penetrating or non-penetrating type, the fractures produced should be approximately perpendicular to the axis of least stress. The general state of stress underground is that in which the three principal stresses are unequal. For tectonically relaxed areas characterized by normal faulting, the least stress should be horizontal; the fractures produced should be vertical with the injection pressure less than that of the overburden. In areas of active tectonic compression, the least stress should be vertical and equal to the pressure of the overburden; the fractures should be horizontal with injection pressures equal to or greater than the pressure of the overburden. Horizontal fractures cannot be produced by hydraulic pressures less than the total pressure of the overburden. These conclusions are compatible with field experience in fracturing and with the results of laboratory experimentation. INTRODUCTION The hydraulic-fracturing technique of well stimulation is one of the major developments in petroleum engineering of the last decade. The technique was introduced to the petroleum industry in a paper by J. B. Clark,' of the Stanolind Oil and Gas Co. in 1948, and since then its use has progressively expanded so that by the end of 1955 more than 100,000 individual treatments had been performed. The technique itself is mechanically related to three other phenomena concerning which ah extensive literature had previously developed. These are: (1) pressure parting in water injection wells in secondary-recovery operations, (2) lost circulation during drilling, and (3) the breakdown of formations during squeeze-cementing operations, all of which appear to involve the formation of open fractures by pressure applied in a wellbore. The most popular interpretation of this mechanism has been that the pressure had parted the formation along a bedding plane and lifted the overburden, notwithstanding the fact that in the great majority of cases where pressures were known they were significantly less than those due to the total weight of the overburden as determined from its density. Prior to 1948, this prevalent opinion had already been queried by Dickey and Andresen,2 in a study of pressure parting; and by Walker,"' who, in studies of squeeze cementing, pointed out that the pressures required were mostly less than those of the overburden, and inferred that the fractures should be vertical. J. B. Clark,' in his paper introducing hydraulic fracturing, and later Howard and Fast," and Scott, Bearden, and Howard," also of Stanolind, postulated that the entire weight of the overburden need not be lifted in producing horizontal fractures, but that it was only necessary to- lift an "effective overburden," requiring a correspondingly lower pressure. Hubbert,7 in discussing the paper by Scott and associates, pointed out that the normal state of stress underground is one of unequal principal stresses; and in tectonically relaxed areas, characterized by normal faults, the least stress should be horizontal. Therefore, in most cases, fracturing should be possible with pressure less than that of the overburden and, moreover, such fractures should be vertical. Harrison, Kieschnick, and McGuire,8 also on the expectation that the least principal stress should be horizontal, argued strongly in favor of vertical fracturing. Scott, Bearden, and Howard6 observed that, when

Jan 1, 1958

By A. S. Appleton, M. B. Bever

The stored and expended energies of an Au-Ag alloy deformed in torsion were investigated as functions of strain, strain rate, and temperature. At room temperature, the stored energy, measured by tin solution calorimetry, increased linearly up to an intermediate strain, above which it was nearly constant; at 78 °K the relation was linear up to fracture. Approximate values of the energy stored at 4 OK indicate a rapid increase below 78 OK. The energy stored by specimens deformed successively at 4o and 78oK changed slowly with increasing second strain from the curve for energy stored at 4 OK (measured after annealing at 78oK) to that for deformation at 78°K. The ratios of the stored to the expended energy at room temperature and 78 oK were unusually large. The additional amount of energy stored at the lower temperature is attributed mainly to an increase in the number and energy of dislocations. The stored energy as a function of strain rate at 4oK (measured after annealing at 78oK) increased to a sharp maximum over a narrow range of strain rates, whereas the expended energy changed only slightly. These effects are explained by the generation of point defects resulting from the enforced nonconservative movement of jogs in fast moving screw dislocations. The effect of straining successively at two strain rates at 4oK was also investigated. Stress-strain curves for 4o, 78oK, and room temperature and the results of micro examination of specimens deformed at these temperatures are also reported. THIS investigation is a continuation of a research program primarily concerned with the effects of variables on the energy relations involved in the deformation of metals. The variables investigated in earlier work have included strain, strain rate, temperature of deformation (room temperature and 78oK), and the type of deformation process. Most of this work has been carried out with Au-Ag alloys. In the investigation reported here, the deformation process was torsion, which had not been used previously in the program. The stored and expended energies involved in the deformation process were measured as functions of strain, strain rate and the temperature of deformation in the range 4" to 373" K. The microhardness and microstructure of deformed specimens were also investigated. The relation between stored energy and strain differs for different deformation processes.1,2 In general, the stored energy increases with strain and may attain a nearly constant value at large strains. This value and the rate at which it is reached, however, depend on the process of deformation. This dependence has been attributed, at least in part, to differences in the effective strain rate and rise in temperature of the specimen during deformation. The effects of the temperature of deformation on both the stored energy and the energy expended in the working process have been the subject of only a few investigations.3-6 Values of these quantities at low temperatures are of special interest as they yield information on point defects which, for many metals, is not obtainable at higher temperatures. The effect of the strain rate on the stored and expended energies has received little attention The strain rate has usually been measured only in relative terms, and hence the interpretation of its effect on the stored energy has been mainly comparative.

Jan 1, 1963

By Albert G. Anderson

I'm sure many of you have wondered what the term "Cross Flow Filtration" means. I shall explain it so that you more clearly understand it, and I hope you will be able to make it a useful tool for you in the future. In the simplest form of cross flow filtration contaminated fluid. or feed stream, is fed at a high velocity down the center of a porous tube, clear fluid weeps through the porous wall of the tube while the suspended solids concentrate in and are swept away by the feed stream. (See figure A). The high velocity flow of the feed stream tends to keep the inside of the tube free from a build up of contamination. By utilizing this cross flow approach to filtration, we have been able to concentrate such things as paint pigments, titanium dioxide, iron oxide, precipitated metals such as uranium yellowcake or silver, cadmium, some metal hydroxides, as well as gelatinous particles such as yeast cells. You will recognize that cross flow filtration lends itself to a recirculation loop system due to the necessity of maintaining a high velocity on the inside of the tube. This means that we generally have a storage tank from which we are recirculating, and the solids will concentrate in the storage tank. (See figure B). The storage tank design becomes an integral part of the total system and is very important to the total function of the system. The cross flow filter system requires only low pressure, and will work with centrifugal pumps, diaphram pumps or in some cases no pump at all. As an example, if we wanted to filter a 10 GPM (37.8 liters/min.) stream that contained 1% of submicronic particulate, we would flow the 10 GPM (37.8 liters/min.) into the storage tank and recirculate through the filter at a flow of anywhere from 30 GPM to 100 GPM (113.6 to 378.5 liters/min.) depending upon the filterability of the particles from the liquid. As the contaminated fluid recirculates, the solids will build up in the storage tank and loop until they reach a level, which is dependent upon the design of the tank. Then the settling rate in the tank becomes equal to the amount of the solids being fed to the tank, and as a result you reach an equilibrium in the amount of solids being presented to the inside of the porous tubing. We have in effect created a dynamic settling tank or thickener. A flow of 10 GPM (37.8 liters/min.) of clear filtrate comes off for either recycle or disposal with periodic removal of thick mud from the bottom of the tank. The thick mud can be either processed through a centrifuge, a rotary drum filter. a flat bed pressure type filter, or other means. Often the small volume of filtrate from the mud de-watering stream can be recycled back into the recirculation tank so that any fines in the feedback are again concentrated, thus closing the loop on the continuous process 100% removal system. Cross flow filtration will not do efficient molecular separation. We work primarily with discrete solids but we do separate them totally from the liquid carrier fluid. I know some of you will raise your eyebrows when I say that we remove 100% of the particulate, but in measurement tests, we have measured the filtrate cleanliness with instruments that record as low as .05 ppm and the instruments read zero, so that it is less than .05 ppm.

Jan 1, 1980

By Anton Sawatzky, Erich Vogt

By means of mathematical solutions to the appropriate diffusion equations, we describe the kinetics of the thermal diffusion of hydrogen in Zircaloy-2 for the various temperatures and concentrations encountered in a heavy water moderated reactor. When the hydrogen concentration is below terminal solid solubility only the a Phase is present. Redistributions are then described in terms of the characteristic functions of the difhsion equation. For higher concentrations both the a and 6 phases are present. We assume the two phases to be always in equilibrium. For moderately small hydrogen concentrations exact solutions of the two-phase equation approach the the approximate solutions derived by Sawatzky and for all concentrations the exact solutions exhibit the qualitative features of his result: the two-phase concentration increases with time, everywhere; in the absence of a hydrogen current at the hot end of the sample an a-phase region always exists there; the interface of the a + 6 , a-phase boundary moves toward the cold md of the sample and the hydrogen concentration is discontinuous at the interface. Simultaneous solulions of the a and a + 6 hydrogen distributims and of the concomitant interface motion are obtained and compared to the observations of Sawatzky and Markowitz. The kinetics of the hydrogen diffusion process are shown to lead to m apparent heat of transport of the a phase which is lower than the actual value (even for samples with long anneals) thus resolving at least partially the disparity between experimental measurements of this quantity. A number of recent papers1"4 have reported measurements on the diffusion and redistribution of hydrogen in Zircaloy-2 under temperature and concentration gradients. These studies were instigated by problems arising from increasing use of Zircaloy-2 as a fuel element cladding material in pressurized-water power-reactors. The Zircaloy-2 picks up hydrogen during the operation of the reactor: the consequent precipitation of zirconium hydride in the Zircaloy-2 has pronounced effects on its mechanical properties. The Purpose of the present Paper is to describe the kinetics of the thermal diffusion of hydrogen in Zircaloy-2 for the various hydrogen concentrations and temperatures likely to be encountered in reactors. When the hydrogen in the Zircaloy-2 is entirely in the solid solution phase (a phase), the differential equation for thermal diffusion is well known and the redistribution can be described by standard mathematical methods some of which are given in Section 11 below. The treatment of the a + 6 (hydride) region differs from the earlier treatment of sawatzky3 by taking into account several modifications of the two-phase diffusion equations as suggested to us by Markowitz4 and Kidson.5 Even with the modifications the results obtained in our more complete treatment are substantially the same as those found by Sawatzky. As we show, the difference between the results of sawatzky3 and kIarkowitz4 is largely due to the difference in the geometry of their experiments. In the next section we derive the modified two-phase equation for arbitrary geometry. The exact solution of this equation is given in Section III for linear and cylindrical geometry. It is shown that Sawatzky's approximate solution is quite accurate for almost all the temperatures and hydrogen concentrations which are actually encountered. In Sawatzky's approximate theory the hydrogen concentration everywhere in the two-phase region increases continuously. As his paper pointed out, this result, together with hydrogen conservation, implies that in a sample with no hydrogen current flowing into the hot end a two-phase region is always accompanied by a single-phase region at the hot end of the sample. The net hydrogen gain in the two-phase region is supplied by the decrease of hydrogen in the single-phase region and by the movement of the (a, a + 6) boundary toward the cold end of the sample. Hydrogen conservation at the (a, a + 6) boundary leads to a discontinuity in the concentration and its derivative there. In Section IV it is shown how these qualitative features of the thermal diffusion kinetics arise from the simultaneous solution of the a-phase differential equation and a + 6 phase equation. Methods derived by us previously for the solution of the redistribution in both phases and the accompanying motion of the boundary are used to obtain approximate solutions valid for large times. The solutions account well for the observed redistributions of Sawatzky and Markowitz. It is shown how the continuing motion of the (a, a + 6) boundary leads to measured values of the heat of transport lower than the actual value, even for specimens with relatively long anneals.

Jan 1, 1963

By J. U. Messenger

A technique is described whereby the resistance of an emudian to breaking can be quantitatively determined. Produced ailfield emulsions are usually the water-in-oil type and, accordingly, do not conduct an electrical current. However, there is a threshold of A-C voltage pressure above which an emulsion will break and current will flow. The more stable an emulsion, the higher the required voltage. A Fann Emulsion Tester, modified so that low voltages (0 to 10 v) can be accurately measured, is suitable. This technique has application in evaluating the effect of a demuksifier on the stability of an emulsion. Emulsions can, in essence, be titrated with demulsifiers by adding a quuntity of demulsifier, stirring, and measuring the voltage required to cause current to flow. Any synergistic effect of two or more materials added simultaneously can be followed accurately. A demulsifier that significantly lowers the threshold voltage (from 100 to 400 v to 0 to 10 v for the emulsions in this study) is effective and can cause the enlulsion to break. A demulsifier that will bring about this drop in the threshold voltage at low concentration ir very desirable. The technique is also well adapted for rapidly screening demulsifiers. INTRODUCTION Stable emulsions in produced reservoir fluids resulting from certain well stimulation and completion procedures are common problems. The use of suitable demulsifiers can often mitigate these difficulties. At the present time, a rapid and efficient method for selecting satisfactory demulsifiers is not available. It is badly needed. Reliance is now placed primarily on trial-and-error procedures. A new test method has been developed which permits a more rapid and precise selection of demulsifiers. It involves measuring the electrical stability potential of an emulsion before and after a demulsifier has been added. This paper describes this method and shows where it should have application in field emulsion problems. NATURE OF OILFIELD EMULSIONS Two immiscible components must be present for an emuhion to form; we are concerned here with crude oil and water. An emulsifier must be present for tin emulsion to be stable. J Emulsifiers can be substances which are soluble in oil and /or mter and which lower interfacial tension. They can be colloidal solids such as bentonite, carbon, graphite, or asphalt which collect at the interface and are preferentially wet by one of these phases. Unrefined crude oils can contain both types of emulsifiers, A popular theory is that, of the two phases in an emulsion, the dispersed phase will be the one contributing most to the interfacial tension.' Usually this phase contains the least amount of emulsifier. The stability of a water-in-oil emulsion is affected by the fol1owing: l) viscosity; (2) particle or droplet size; (3) interfacial tension between the phases; (4) phase-volume ratios; and (5) the difference in density between the phases. A stable emulsion is usually characterized by high-viscosity, small droplets, low interfacial tensions, small differences in density between its phases, and slow separatian of the phases. It also has low conductivity (high electrical stability potential). Water-in-oil and oil-in-water emulsions"' are both common; however, oil field emulsions are predominantly water-in-oil emulsions. The emulsions which commonly occur during oompletion and stimulation operations contain a combination of several of the following: acids, fracturing fluids (oil, water, acid), and formation water and oil. Produced emulsions usually contain formation water and oil. Emulsions form in oil wells because oil and water are mixed together at a high rate of shear in the presence of a naturally occurring or unavoidably produced emulsifier. During the completion and stimulation of productive zones, and while formation fluids are being produced, oil and water are very often commingled. These mixtures are formed into emulsions by agitation which occurs when the fluids are pumped from the surface into the matrix of the formation or produced through the formation to the surface. Restrictions to flow (such as perforations, pumps, and chokes)".'" increase the level of agitation; tight emulsions are more likely to form under these conditions. Often an emulsified droplet is an emulsion itself.'" Therefore, emulsion-breaking problems can be quite complex. The complexity can be even greater if a third phase (gas) is included. Demulsifiers operate by tending to reverse the form of the emulsion. During this process, droplets of water become bigger, viscosity is lowered, color becomes darker, separation of the phases faster and electrical stability potential approaches zero. Any of these effects could be followed as a means of determining emulsion stability. However, electrical stability potential is the most reproducible and most easily measured parameter for following the stability of a water-in-oil emulsion.

Jan 1, 1966

By C. R. Barrett, N. N. Singh Deo

The high-temperature steady-state creep behavior of a series of dilute copper-silicon alloys was studied to determine the effect of stacking fault energy on the creep-rate. The steady-state creep rate is, when taken at equivalent diffusivities decreases with decreasing stacking fault energy. The stress and temperature dependencies of is suggest that creep is a difusion controlled dislocation climb process. Electron microscopy studies of the creep substructure revealed: 1) the subgrain size is not a function of the stacking fault energy in these alloys, 2) the dislocation density not attributed to the subgrain walls seems to be higher during primary creep and decreases to a lower steady value during steady-state creep, and 3) the dislocation density during steady-state creep decreases with decreasing stacking fault energy. In the past few years numerous investigators have studied the influence of stacking fault energy on high-temperature creep strength. Most of these investigators have confined their attentions to studying the relationship between steady-state creep rate, is, and stacking fault energy, ?, when samples are tested under conditions of comparable stress and temperature. For the case of fcc metals, it was initially shown by Barrett and Sherbyl and since confirmed by many others2"4 that is decreases with decreasing ?, often following an empirical relation of the form i ?m where m is a constant about equal to 3. The application of theory to explain this observation has not been entirely successful. One of the main difficulties has been the almost complete lack of structural information (dislocation density, subgrain size, and so forth) for samples with different stacking fault energies, tested under high-temperature creep conditions. weertman5 has attempted to explain the stacking fault energy dependence of is on the basis of a dislocation climb mechanism. Assuming that both the rate of dislocation core diffusion and the ease of athermal jog formation decreases as ? decreases Weertman has argued that the rate of dislocation climb and hence the creep rate should also decrease as ? decreases. One questionable aspect of Weertman's analysis is the assumption that core diffusion down extended dislocations is slower than core diffusion down unextended dislocations. The only experimental work done in this area, by Birnbaum et al.6 on nickel and Ni-60 Co, has shown the core diffusivity to increase with decreasing ?. Theories of steady-state creep based on the diffusive motion of jogged screw dislocations often seem unable to predict even the qualitative nature of the es- relationship. Assuming that Weertman is correct in his assumption that the dislocation jog density decreases with decreasing ? then the jogged screw theories predict an increasing dislocation velocity with lower ?. It is usually assumed that the increase in dislocation velocity implies a corresponding increase in creep rate. However, two other factors must be considered before such a statement can be made. That is, we must know how both the mobile dislocation density and the effective stress (the difference between applied stress and internal stress) vary with ?. Significant changes in either one of these factors could outweigh any change in dislocation velocity accompanying a change in ?. And with the slower rates of recovery expected in low stacking fault energy materials it seems likely to expect both mobile dislocation density and effective stress to be dependent on ?. Sherby and Burke7 have suggested that stacking fault energy influences the creep rate in an indirect way. These authors cite evidence that the steady-state subgrain size generated during high-temperature creep is a function of ? decreasing with decreasing ?. Assuming the creep rate to be proportional to the area swept out by each expanding dislocation loop and that subgrain boundaries are good barriers to dislocations, then the creep rate should be proportional to subgrain area, hence increasing as ? increases. A critical evaluation of any of the above theories requires more quantitative information concerning the dislocation substructure generated during high-temperature creep. Accordingly this investigation was undertaken with an aim of studying the influence of stacking fault energy on tbe steady-state creep characteristics of a series of dilute copper-silicon alloys. Special emphasis was placed on studying the strain dependence of both the dislocation configuration and density. MATERIALS AND PROCEDURE Dilute copper-silicon alloys of the compositions shown in Table I were tested in tension at constant stress. The relative stacking fault energy of these alloys has been determined and is shown in Table 11. An Andrade-Chalmers lever arm was used to maintain constant stress and testing was carried out in a water

Jan 1, 1970

By W. A. Griffith

Fractional replication of factorial design, a general method for planning experimentation and for analysis of data obtained, is described as applied to a flotation investigation. This method leads to determination of optimum process conditions with minimum experimental effort. Its advantages over simple factorial design are demonstrated. A METHOD for planning experimentation and for analyzing the data secured will be demonstrated. This method, termed fractional replication of factorial design, employs a relatively small number of individual experiments to determine which of a large number of process variables are controlling, to determine which combination of levels of these variables is most likely to produce optimum results, and also to predict what results will be obtained with a particular combination of conditions not already tested. Although the general method is not new, having been developed by Finney in 1945,' the extent to which it can improve the effectiveness of applied research generally has not been recognized by metallurgists. The fractional replication procedure is particularly useful in flotation experimentation and an example from such an investigation will be used in the paper. However, it has equal value in any investigation in which similar experimental difficulties are encountered. In developing a flotation process for a particular mineral separation, the investigator is inevitably confronted with the following difficulties: 1—There are a large number of potentially important process variables. 2—Results of individual experiments are not highly reproducible, due in part to the difficulty in precisely controlling all the variables. 3—Considerable effort is expended in conducting individual experiments. 4—There are practical limits on the number of individual experiments which can be made. In situations of just such a type, modern statistical methods of planning experimentation and analyzing data have their greatest value. Applications of one such technique, called factorial design, to problems of this type have been described by Dorenfeld and others.'-' The simple factorial design is an efficient procedure when the investigator hopes to provide a comprehensive understanding of the effects and interrelationships of a small number of variables over a limited range. In applied research, this is seldom the main objective. Rather, the investigator usually wishes to determine which of the many potentially important variables are in fact controlling, which levels of the controlling variables will provide opti- mum metallurgical results, and what these results will be at optimum conditions. Interest in detailed trends is limited to the controlling variables and to levels in the region of optimum conditions. Simple factorial design has serious deficiencies for such objectives and is not the most efficient method of experimental design. Deficiencies of Factorial Design In a simple factorial design, an experiment must be made at every possible combination of each level of every variable, once these have been chosen and the levels of each to be included have been decided upon. As the number of variables or levels of each increases, the experimental program quickly reaches prohibitive size. For example, an investigation of only four variables, each at four levels, requires 256 individual experiments. Often upon completion of such an extensive program, it is found that one or more of the variables has metallurgically unimportant effects or that a poor estimate has been made as to the appropriate range of levels to be investigated. The result is that only a small proportion of the data obtained falls in the range of real metallurgical interest. Indeed, it frequently can be anticipated that certain combinations of levels of variables will not produce results of interest, but they still must be included so that the essential balance, or orthogonality, of the design will be retained. It may be true that factorial design will provide the greatest amount of information from a given number of experiments, but it does not necessarily follow that it will lead to the desired information with a minimum number of experiments. Much of the information provided may be of little value. Advantages of Fractional Replication The disadvantages of simple factorial design are overcome to a great extent by a modification known as fractional replication. This is a technique for sampling systematically the potential data of a full factorial experiment, that is, the data which would have been obtained if the complete factorially designed program had been completed. Only a fraction of the total array of experiments is made, but the experiments are chosen in such a way that the important advantages of factorial design and the accompanying analysis of variance are retained. The data obtained from the first group of experiments are used to determine which of several variables are controlling and which levels of these variables are most likely to produce the desired result. Unimportant variables and levels of variables then may be

Jan 1, 1956

By W. Hurst, R. E. Leeser, W. C. Goodson

Three aspects of gas deliverability are presented in this paper. The first treats with the gas deliverability or availability of a normal depletion-type dry gas field. Such encompasses not only the period of stabilized constant rate, but more so, the "tailings" when a fixed abandonment pressure is reached and the rate by necessity must decline. A comprehensive work plot is offered, developed from mathematics herein included, that removes the triai-and-errnr computations that attended such undertakings in the past. The second part treats with the discount factor of the open flow potential constant from what is observed initially in testing a gas well to what is evidenced when stabilization is reached. This prevails in tight formations, such as the Kansas Hugoton field which is offered as the example. The means of establishing this factor are pressure build-up curves which, as sustained by analytical deductions, reproduce this entire period of transient flow under conditions of a constant rate influx. Finally, what is offered in this paper is the deliverability performance of an exceedingly rich gas condensate field producing from a tight formation. The example shown is the Knox Bromide field in Oklahoma, producing from the Bromide formations. The results are ominous, showing early reduction in permeability to gas pow, due to the retrograde condensate forming in the pore space, with the attending early logging-up of these wells. The analytics of lowered permeability are incorporated in the gas deliverability formula along with the PVT data that gives the increased condensate liquid saturation as the gas flows to the well bore. This paper would not be complete without a critique oflered at the end. With the many gas wells now in production and those that have completed their life, there has been no factual information collected by any source as to what constitutes that permeability range where a gas well would be unimpaired in its gas deliverability by the presence of rich condensate content, and the lowered range where such would be harmful. This question confronts all producers. INTRODUCTION Various aspects of gas deliverability are presented in this paper that includes depletion-type reservoirs, deteriora- tion factor of the gas deliverability constant, and the performance of a rich gas condensate reservoir producing from a tight sand. With respect to the presentation of gas deliverability and its tailings for depletion-type gas reservoirs, one notes that this is essentially the information offered by every gas transmission company and producer appearing before the Federal Power Commission for Letters of Conveyance in the dedication of reserves. In the ordinary procedure, as many engage upon this study, trial-and-error calculations are included, particularly as apply to the tailings. For many years one of the writers has employed mathematical analyses to perform this step and avoid the complexities so associated. In the preparation of this paper these analyses have been amplified to include any slope n for the open flow potential relationship for which the tailings can be determined from Fig. 1. With reference to the deterioration or discount factor of the open flow potential constant as such occurs in the gas deliverability formula, this for the most part has been an unexplored subject. Although the issue first appeared in the Kansas Hugoton field, where such was surmised but only recently resolved, this situation of a deterioration of the gas deliverability constant can occur wherever dry gas production from a tight sand is encountered. The first concerted attacks upon this problem were the presentations of Hurst' and Goodson' before the Kansas Corporation Commission to show that transient fluid flow and unsteady-state flow formulas prevailed. This was amplified later before the Federal Power Commission3 to show that this deterioration factor could be identified from pressure build-up curves. This has been reported by McMahon.4 Its importance to the industry merits the review of these essential features in completing the program on the aspects of gas deliverability. Finally, as illustrated here, for a low permeability formation such as the Knox Bromide field where the gas is rich, representing some 165 bbl of condensate per MMcf of effluent gas, the gas deliverability can be of limited extent in the life of the field, leaving substantial amounts of condensate and gas unrecovered. In cases such as this, gas cycling is mandatory. This is particularly revealed by the fluid mechanics introduced here, employing factual field as well as laboratory data, to show this-restriction upon gas deliverability. PRESSURE DEPLETION What will now be offered is the study of gas deliver-

By D. Havlena, A. S. Odeh

The material balance equation used by reservoir engineers is arranged algebraically, resulting in an equation of a straight line. The straight line method of analysis imposes an additional necessary condition that a successful solution of the material balance equatiott should meet. In addition, this algebraic arrangement attaches a dynamic ineuning to the otherwise static material balance equation. The straight line method requires the plotting of one variable group vs mother variable group. The sequence of the plotted points as well as the general shape of the resulting plot is of utmost importance. Therefore, one cannot progrm the method entirely on a digital computer ar is usually done in the routine solution of the material balance equation. If this method is applied, then plotting and anaIysis are asential. Only the appropriate equations and the method of analysis and interpremtion with comments and discussion are presented in this paper. Illustrative field examples for the various cases treated are deferred to a subsequent writing. INTRODUCTION One of the fundamental principles utilized in engineering work is the law of conservation of matter. The application of this principle to hydrocarbon reservoirs for the purpose of quantitative deductions and prediction is termed "the material balance method of reservoir analysis". While the construction of the material balance equation (MBE) and the computations that go with its application are not difficult tasks, the criteria that a successful solution of the MBE should fulfill have always been a problem facing the reservoir engineer. True and complete criteria should embody necessary and dcient conditions. The criteria which the reservoir engineer uses possess a few necessary but no sufficient conditions. Because of this, the answers obtained from the MBB are always open to question. However, the degree of their acceptability should increase with the increase in the number of the necessary conditions that they should satisfy. Generally, the necessary conditions commonly used are (1) an unspecified consistency of the results and (2) the agreement between the MBE results and those determined volumetrically. This second criterion is usually overemphasized. Actually, the volumetrically determined results are based on geological and petrophysical data of unknown accuracy. In addition, the oil-in-place obtained by the MBE is that oil which contributes to the pressure-production history,' while the volumetrically calculated oil-in-place refers to the total oil, part of which may not contribute to said history. Because of this difference, the disagreement between the two answers might be of paramount importance, and the concordance between them should not be overemphasized as the measure of correctness of either one. In this paper, a third necessary condition of mathematical as well as physical significance is discussed. It is not subject to any geological or petrophysical interpretation, and as such, it is probably the most important necessary condition. It consists essentially of rearranging the MBE to result in an equation of a straight line. This straight line method of the MBE solution has invalidated a few long time accepted concepts. For instance, it has always been advocated that if a water drive exists, but one neglects to take it into account in the MBE, the calculated oil-in-place should increase with time. The straight line method shows that in some cases, depending on the size of the neglected aquifer, the calculated oil-in-place might decrease with time. The straight line method requires the plotting of a variable group vs another variable group, with the variable group selection depending on the mechanism of production under which the reservoir is producing. The most important aspect of this method of solution is that it attaches a significance to the sequence of the plotted points, the direction in which they plot, and to the shape of the resulting plot. Thus, a dynamic meaning has been introduced into the picture in arriving at the final answer. Since the emphasis of this method is placed on the interpretation of the sequence of the points and the shape of the plot, one cannot completely automate the whole sequence to obtain "the best value" as normally done in the routine application of the MBE. If one uses the straight line method, then plotting and analysis are musts. The straight line method was first recognized by van Everdingen, et al,2 but for some reason it was never fully exploited. The advantages and the elegance of this method can be more appreciated after a few cases are carefully treated and worked out by it. SOLUTION OF THE MATERIAL BALANCE EQUATION SATURATED RESERVOIRS The MBE for saturated reservoirs written in AIME symbols is

By Ralph W. M. Lai, D. W. Fuerstenau

This paper presents the results of a study of the principles which control the distribution of ultrafine particles between an oil phase and an aqueous phase. Alkyl sulfonates were used to control the wettability of fine alumina particles and to stabilize the water-isooctane emulsion system in the extraction process. The effects of sulfonate chain length, sulfonate adsorption, pH, and oil-water ratio on the recovery of solids in the oil phase were investigated. Higher extraction was found under conditions which yield higher hydrophobicity and large interfacial area of the oil-water emulsion. Today there is considerable interest in developing methods for the recovery of fine particles from disseminated ores. One such method might be the extraction of fine particulate solids from aqueous suspensions into an oil phase, and one example of a successful industrial process based on oil-water distribution phenomena is the emulsion flotation process used in concentrating manganese ore.' The oil agglomeration process of Farnand, Smith, and Puddington has been suggested for concentrating tin ore. Patents 3-4 for oil flotation processes extend back to that of William Haynes in 1860. The basis of all these processes is directly related to the stabilization of emulsions with solid particles. The first extensive study of oil-water emulsions stabilized by finely divided solids was made by Pickering. Later, Briggs,6 weston, 7 Chessman and King,8 and Leja and Schulman9 experimentally demonstrated how solid particles control water-oil and oil-water systems. Hildebrand10 suggested that fine particles can stabilize an emulsion because the fine particles are partially wettable by both oil and water, that is, that the contact angle must be finite. Takakuwa and Takamori 11 recently applied emulsion inversion to the study of collector adsorption in flotation systems. At the Royal School of Mines, Shergold and Mellgren12 are currently carrying out a comprehensive study of the extraction of fine particles into an organic liquid phase. In the present study of the liquid extraction or oil flotation of ultrafine particles, the recovery of hydro-philic alumina particles of 0.1-µ average diameter was determined after separating an isooctane-water mixture (emulsion) into its two phases. Alkyl sulfonates of different chain lengths were used to change the wettability of the alumina and to effect emulsifi-cation of the oil-water mixture. The effects of the oil-water volume ratio, sulfonate adsorption, and pH on the recovery of the fine particles were investigated. Since the surface charge on an oxide such as alumina is determined by the adsorption or interaction of H+ and OH- with the surface, pH control is extremely important. Furthermore, since the zero point of charge (zpc) of alumina13-14 occurs at pH 9 to 9.4, major adsorption effects with anionic collectors in these systems will occur below pH 9. The results of the foregoing investigation are discussed in terms of the recovery of the fine particles, the oil-water ratio, contact angles, surface tension, the adsorption of sulfonate onto the solid surface, and the emulsifica-tion of oil. MATERIALS AND EXPERIMENTAL METHODS Materials: The particulate solid used in all the experiments was Linde "A" alumina, a high-purity, readily available form of a-alumina. The surface area of the alumina, as evaluated from measurements of stearic acid adsorption from benzene, was 14.5 sq m per gm, based on a value of 20.5 sq Å for the cross-sectional area of stearic acid.13 All aqueous solutions were prepared from triply distilled conductivity water. Other materials included high-purity crystals of sodium alkyl sulfonates, isooctane, and, for pH adjustments, reagent grade hydrochloric acid and sodium hydroxide-Because sulfonic acids are strong acids, they are

Jan 1, 1969

By Henry A. Wriedt, John Chipman

Equilibrium in the reaction of hydrogen gas with oxygen in liquid nickel, iron, and their alloys has been studied at temperatures of 1500° to 1700°C. The equilibrium con^stant, 0/p, [% O], is greater in nickel then in iron by a factor of 45 at 1594 C. In the alloy steel range, the activity coefficient of oxygen is slightly increased by the presence of nickel. NICKEL is more resistant to oxidation than is iron, in both the solid and the liquid state. The solubility of oxygen,' however, is greater in nickel than in iron. The two metals and their liquid alloys constitute an interesting series of solvents in which the chemical behavior of oxygen can be studied by methods similar to those used for liquid iron.' Experimental Method The furnace assembly, including the tube for preheating the entrant gas, was the same as that described by Gokcen and Chipman. " controlled mixture of water vapor, hydrogen, and argon preheated to approximately the melt temperature impinged upon the surface of the liquid metal. The charge, weighing about 35 g, was held in an alumina crucible at constant temperature in the gas stream for a prolonged period, then quenched in a stream of cold helium, sampled, and analyzed for oxygen by the vacuum fusion method. The temperature of the liquid metal was measured by a Leeds and Northrup optical pyrometer which had been checked against a similar instrument certified by the National Bureau of Standards. The freezing point of electrolytic iron in hydrogen, taken as 153j°C, was used for frequent checks. This point and the emissivity determined by Dastur and Gokcen" established the transmissivity of the optical system. The emissivities of Fe-Ni alloys at 1535" are known from the experiments of Smith and Chipman,h nd it was assumed that the variation with temperature is parallel with that of iron. Recorded temperatures are probably accurate within k5"C. Preparation of Gas Mixtures—Since the ratio H,O:H, had to be much higher for nickel than for iron, it was found convenient to construct two separate systems for gas preparation. In both cases, the ratio of argon to the sum of the other two gases was at least 4, the proportion of hydrogen and water vapor being varied to produce H,O:H, ratios from 0.2 to 70. The primary purpose of the argon was to diminish the errors caused by thermal separation' in the gas phase. In addition, its presence lowered the saturation temperature required for a given H,O:H, ratio and minimized the difficulty of preventing condensation. Both systems were modifications of the one used by Dastur and Chipman.' For H,O:HI ratios greater than 2, welding-grade argon purified by passage over magnesium chips at 620°C was mixed with hydrogen produced in the electrolytic cell shown in Fig. 1. The design was based on a cell used by Bodenstein and Pohl," the important features of which area nickel wire cathode and sodium hydroxide solution as the electrolyte. The rate of hydrogen evolution was controlled by a calibrated ammeter and a rheostat in the circuit, the current efficiency being taken as 100 pct. Precise control of the argon flow rate at about 300 ml per min was achieved by means of a capillary flowmeter having a controlled pressure drop. The molar flow rate was measured volumetrically before each run, and a constant flowmeter setting was maintained thereafter. The A-H, mixture passed through a water saturator of the type previously described," the saturation temperature being maintained automatically with a fluctuation of k0.04"C. The tube between saturator and furnace was heated by a resistance winding. The system was of all Pyrex construction with the exception of the stainless steel tube containing magnesium chips. This was connected with the glass tubing at the upstream end by

Jan 1, 1957

By Robert E. Brand, Strathmore R. B. Cooke

For a study of the bonding that occurs in magnetite pellets during oxidation, cubes of magnetite, hematite, and quartz were prepared, each with one surface polished and nearly optically plane. These faces were fired together under pressure and at various temperatures to determine whether bonding would occur. Sections of surfaces that bonded were examined microscopically. An attempt was made to measure the strength effected by this treatment. EARLIER work1,2 has shown that the strength developed by fired magnetite pellets containing silica or silicates is due to three separate effects. Two of these are associated with the formation of slag in the interior of each pellet at temperatures above 1000°C, where the atmosphere is substantially though temporarily neutral. The third is independent of the silica content and is responsible for the development of a hard shell around the pellet during the initial stages of firing, a condition which permits the pellet to assume the increasing load of the overlying burden and prevents extrusion of the slagged and semi-liquid interior. Pellets made from hematite concentrates do not show these effects unless specific additives are incorporated, and it is these differences in behavior which in large part account for the difficulty in producing satisfactory pellets from such concentrates. Elsewhere it has been stated that "the most important factor responsible for the development of strength in fired magnetite pellets containing no additive is the oxidation of magnetite to hematite and the subsequent recrystallization and grain growth of this constituent."' The object of the investigation reported in the present paper was to study in detail the mechanism described above and to evaluate its contribution to the strength of the fired pellet compared with other possible mechanisms, excluding slag bonding. This was effected by placing together polished faces of two cubes prepared from specific minerals and then heating them under constant pressure at different temperatures both in air and in nitrogen. The materials used were magnetite, hematite, and quartz. The magnetite was coarsely crystalline material from western Montana, the hematite was Itabira ore from Brazil, and the quartz was of optical grade but of unknown origin. A high-grade Lebanon (Pennsylvania) magnetite concentrate was used to prepare briquets for the investigation of bonding strength. Analyses of the iron oxide materials are given in Table I. The specimens for furnace treatment were prepared with a carborundum cut-off wheel, which cut the minerals into cubes approximating % in. along each edge. Care was taken to make the sides as parallel as possible. The edges of the faces to be polished were beveled, and the cubes were mounted in brass rings with sealing wax, this being essential to prevent the specimens from rocking during polishing. To prevent excessive friction, the wax was recessed so that only the brass ring and the face of each specimen would be in contact with the polishing lap. Each specimen was then ground with No. 600 carborundum to produce a level surface before polishing. The Graton-Vanderwilt process was not used for polishing because of the long time required by this method, but several other methods were tried. Cloth laps using levigated alumina or chromium oxide did not produce the desired polish. Satisfactory results were obtained by means of a pitch lap: which was used for all the mineral cubes, with levigated alumina as the polishing medium. The lap was operated in a machine similar to those used in the preparation of single optical surfaces. It was run at low speeds and moderate pressures, and polishing was stopped occasionally to prevent overheating of the lap. A good polish with flat surfaces could be obtained in about 5 min. After polishing the sealing wax was removed either in hot water or in a solvent such as turpentine.

Jan 1, 1955

By J. C. Deppe

A method is presented for calculating approximate injection rates in secondary recovery operations. The method can be applied to cases of unequal fluid mobilities, irregular well patterns and boundary patterns. The steady-state pressure distributions for the four flood patterns reported by Muskat and for five additional patterns reveal that most of the difference in pressure between the injection and producing wells occurs in regions around the wells which can adequately be described as regions of radial flow. This leads to a method of calculating injectivity by approximating the flood pattern with radial flow elements (or a combination of radial- and linear-flow elements for some patterns such as the direct line drive). Irregular and boundary patterns can also be approximated by radial and linear elements. Each of these elements can be described by radial-and linear-flow equations and the results combined as series flow resistances to give an approximate equation for the initial injection rate. The mobility ratio does not affect the initial rate; therefore, if the well pattern is one of those regular patterns for which theoretical rate equations have been derived for unity mobility ratio, the approximate initial rate equation can be improved by adjusting it to match the theoretical equation. The available theoretical rate equations are listed, including five new cases. As the flood progresses, the injectivity changes because in general the flood front will divide the pattern into areas of different fluid mobilities. Simple shapes can be assumed for the flood front so that both areas can be divided into radial- and linear-flow elements. Radial- and linear-flow equations are applied to these elements to account for the change in flow resistance behind the front. To calculate injection rates after breakthrough, it is necessary to know the sweep efficiency at breakthrough. Areal sweep data are available in the literature for a number of patterns, and sufficiently accurate breakthrough sweep efficiency can be estimated from these data if it is not otherwise available. Again, simple shapes for the flood front can be assumed after breakthrough, and injection rates can be calculated for the remainder of the flood. When areal sweep data are available, these data can be used to cbeck. the calculated injection rates after breakthrough. Injection rates calculated by this approximate method compare favorably with available model data. INTRODUCTION When the fluid mobilities in the swept and un swept regions are equal, the injectivity will not change as the flood front advances; and for regular patterns it can be calculated by mathematical formu1as.l When the fluid mobilities in the swept and unswept regions are not equal, the injectivity will increase or decrease as the flood front advances. In this case, the injectivity has not been calculated by analytical means for any practical well pattern; and, furthermore, scale model and analog4,6 results have been published only for the five-spot pattern. The main object of this paper is to present an approximace method for calculating injectivity for the case of unequal mobilities. The method can be applied to regular, irregular and boundary patterns. Before the approximate method is discussed, the analytical formulas for equal mobilities will be summarized, and the analytical solution for radial flow with unequal mobilities will be used to show how the injectivity changes as the flood progresses. EQUAL MOBILITIES, REGULAR PATTERNS In a flooding operation the injectivity, the rate at which fluid can be injected per unit difference in pressure between injection and producing wells, depends on basic physical properties and, in addition, on three variables to be treated here—the area of the swept region, the fluid mobilities in the swept and unswept regions, and the well geometry (pattern,

By C. C. Martin, M. C. Fuerstenau, R. B. Bhappu

Experiments revealed that quartz could not be floated in conductivity water at any pH with a long-chained sulfonate as collector. Various cations, Fe+++, Al+++, Pb++, Mn++, Mg++, Ca++, are shown to function as activators when the pH is such that the cation hydrolyzes. As the solubility product of the various cation-sulfonates was exceeded at the concentrations involved, precipitated cation-hydroxy-sulfonates, e.g. Fe (RSO3)2OH must be functioning as the collector in these systems. A number of interesting and important phenomena were revealed during a study of the response of beryl to sulfonate flotation.' Perhaps the most important observation is the effect of various cations on flotation response. Experiments showed that in the vicinity of pH 3.0, about three ppm Fe+++ increased recovery of leached beryl from about 60 pct to 100 pct, whereas about four ppm Ca+ + reduced the recovery from 60 pct to 2 pct. With the exception of ferric iron, all of the other cations investigated reduced flotation recoveries drastically at this pH. As postulated in the paper Sulfonate Flotation of Beryl,' even with low concentrations of polyvalent cations and sulfonate (on the order of one ppm), cation-sulfonates are precipitated immediately in solution. That is, blueish-white clouds were noted to form after given additions of salts to known solutions of sulfonate. Furthermore, it was also shown that one-hundred-fold more ferric iron than necessary to form the precipitated cloud had to be added to effect complete flotation. In other words, some form of solid ferric sulfonate was found to function as the collector in these systems. The form of the compound was suggested to be Fe(RCO3)2OH and that bonding to the surface occurs through the hydrogen of the hy-droxyl. As the form of this compound was shown to be dependent on the extent of hydrolysis of ferric iron, it is to be expected that ferric iron will be unique in some pH region. Further, if the presence of hydroxyl in the neutral precipitated compound is of primary importance, then other cations should function similarly but at different values of pH, since they will hydrolyze at different concentrations of hydroxyl ion. To determine whether the premise of hydrolysis is correct, the study was extended to quartz systems in the presence of various cations at various values of pH. EXPERIMENTAL MATERIALS AND METHOD Sodium alkyl aryl sulfonate was chosen as the collector for the experimental work. The reagent, sup plied as a courtesy sample by the Shell Chemical Co., has the following physical properties: Physical form Solid — Finely ground Sulfonate content, pct wt 95 to 97 Molecular Weight 450 to 470 Number of carbon atoms in hydrocarbon chain 25 to 30 Melting Point, OF 250 to 260 Solubility Soluble in water with gel formation at concentrations higher than 25 pct by weight sulfonate This particular reagent was chosen because of its high active sulfonate content, long-chain length, and solid form. With the exception of the frother and sulfonate, all other chemicals were of reagent grade quality. Conductivity water was used in the experimental work. This water, made by passing distilled water through an ion exchange column, had an average measured conductivity of one micromho. Experiments were conducted at room temperature in a small glass flotation cell with the following procedure: 1) A redetermined amount of water and salt solution (e.g. CaC12 . 2H2O) were combined and the pH adjusted to a given value. 2) A given amount of sulfonate was added so that final solution volume was 130 cc. 3) One drop of Dow 250 frother was added. 4) Five g of quartz were added and the system conditioned for five minutes. 5) The pH of the system was measured (termed flotation pH).

Jan 1, 1963

By L. A. Wood, A. G. Winslow

Subsidence has occurred in several areas of the upper Gulf Coast region of Texas, although in most cases this is not evident without precise instrumental leveling. As referred to in this report, the region comprises about 7200 sq miles of the Gulf Coastal Plain and includes all or parts of the following counties: Harris, Fort Bend, Waller, Brazoria, Galveston, Liberty, Chambers, Orange, and Jefferson (Fig. 1). Topographic relief is low, ranging from sea level at the Gulf of Mexico to 300 ft or more in northwestern Harris County. In this upper Gulf Coast region, which embraces the most heavily populated and industrialized parts of Texas, industrial and municipal water supply is obtained chiefly from wells. Largest center of ground-water withdrawal is the Houston-Baytown area, where about 200 million gallons per day was pumped in 1956 for municipal and industrial use. Other large centers of ground-water withdrawal include the Alta Loma sector and the heavily industrialized Texas City area in Galveston County, where about 23 mgd was pumped in 1956. In the Freeport area in Brazoria County about 6 mgd was pumped in 1956 for city and industrial supply. In addition to withdrawals in metropolitan districts, great quantities of ground water are supplied by irrigation wells to thousands of acres of rice. These wells are widely dispersed over many miles, principally in Harris, Waller, Fort Bend, Brazoria, and Liberty counties, and effects of pumping are not so readily noticed as in the industrial and metropolitan sectors. Subsidence: Starting with the 1905-1906 line from Smithville, Tex., to Galveston, the U. S. Coast and Geodetic Survey has extended a network of first and second-order level lines throughout the upper Gulf Coast region. In 1918, when a line was run from Sinton, Tex., to New Orleans, La., a tie was made to the 1905 leveling at Houston. In 1932-1933 a line was run from Palestine, Tex., to Houston, and a good deal of leveling was carried out in the vicinity of Houston. At this time the entire network was adjusted to the Sea Level Datum of 1929. In 1935-1936 and 1942-1943 several new lines were run and considerable releveling was performed. Supplementary adjustments were made after both the 1936 and 1943 levelings in order to reduce the elevations to the Sea Level Datum of 1929. These adjustments probably masked some of the subsidence that had occurred up to that time, although a certain amount of subsidence was noticed in the 1943 leveling. When several lines in the vicinity of Houston and Galveston were releveled in 1951, the regional nature of the subsidence became apparent.' Direct comparison with the 1943 leveling was possible because the datum was not changed. Comparisons with prior leveling were unreliable because of differences in the early and subsequent datums. In 1953-1954 many lines in the region were re-leveled and in 1957 the entire network, including the elevations previously determined, was readjusted to determine the changes in elevations as indicated by various levelings between 1905 and 1953-1954. In this instance the adjusted elevations were held fixed near the boundaries of the region where the bench marks seem to have remained stable. The only sea-level connection held fixed was the connection of the 1905 leveling to the 1903-1906 tidal series at Galveston. Almost the entire region has been at least slightly affected by the subsidence, some areas to a great extent. The profiles of subsidence (Fig. 2) and the map showing the contours (Fig. 3) illustrate both regional and local changes. A large area of subsidence is centered in the Houston-Baytown district, where the maximum recorded is about 2.7 ft between 1905 and 1954. This probably falls short of the actual maximum, since land surface subsided as much as 2.5 ft in the industrial section alone between 1943 and 1954. Sub-

Jan 1, 1960

By W. W. Mein

IN March 1949 the Calaveras Cement Co. was sued by five landowners whose properties are located in the vicinity of the plant. These landowners—all of them cattle ranchers—sued for dust damages of $120,338 and for an injunction preventing the company from casting dust upon their properties in injurious quantities. The issue was split, the jury deciding the amount of damages and the court handling the injunction features. The action came to trial before a jury in April of that year and resulted in damages being awarded in the amount of $7,508. A month later the court issued its injunction requiring recovery of stack dust to an 87 pct minimum. The Calaveras plant, shown in Fig. 1, is located three miles from San Andreas in the heart of the Mother Lode country, about 45 miles east of Stockton, Calif. Since it is in a small town, the company's responsibilities to its employees and to the community are different from those of plants located in large centers of population. Construction of the plant was begun in 1925, and production was started the following year. Standard gray Portland cement as well as several specialty cements are produced. Production facilities have been continually enlarged so that the present capacity is more than double what it was 25 years, ago. The most important item in the company's expansion program was the purchase in 1945 of an 11 ft 3 in. x 360 ft Allis-Chalmers kiln which had been declared surplus by the Defense Plant Corp. The company is now able to produce about 2,500,000 bbl of cement a year, and employment at the plant is approximately 300. The production of cement has always involved a dust problem, solution of which has been attempted in a number of different ways. By far the most efficient method of catching dust discharge is electric precipitation. Soon after the company went into production in the 1920's, the landowners near the plant began complaining about dust discharge. The company settled the matter by paying them a total of $27,000 for dust damages incurred during the period from 1929 to 1938. In the latter year the two-unit Cot-trell Electric Precipitator shown in Fig. 2 was installed and collected over 85 pct of the dust resulting from operation of the two 11 ft 3 in. x 10 ft x 240 in. kilns. From 1938 until 1946 the landowners registered no objection and were apparently satisfied with that dust recovery. In 1946, however, the substantially higher production resulting from operation of the new kiln immediately overtaxed the facilities of the existing two-unit precipitator. In a short time claims for dust damages were again received from the landowners. The company informed them that because of postwar shortages it would take several years to install a new precipitator. The company stated that it intended to purchase an additional precipitator and offered in the meanwhile to pay the landowners a fair amount to cover their damages. This offer was rejected and the landowners filed suit in the Superior Court of Calaveras County asking for total damages of $120,-338 and for the injunction outlined above. Case for the Ranchers The attorneys for the plaintiffs produced three different types of witnesses: first, the landowners themselves and members of their families; second, other cattle raisers who were sympathetic to the plaintiffs; and third, several types of technical experts: veterinarians, soil chemists, and an aerial photographer. The testimony of these witnesses was designed to establish that the flue dust from the plant had not only damaged the land and the forage, but had also caused fluorine poisoning to their animals. By a combination of these causes, the plaintiffs attempted to prove substantial loss of profits.

Jan 1, 1952

By J. S. Levine, M. Prats

Several methods are available for calculating the performance of solution-gas-drive reservoirs from the PVT properties of the oil and from the relative permeability and other properties of the formation. These methods require a number of simplifying assumptions. The present method of computation has made use of a high-speed computer to solve simultaneously the nonlinear partial differential equations that describe two-phase flow by solution-gas drive in order to calculate the performance of a reservoir. Some of the results obtained by the nonlinear partial differential equation solution are compared with those obtained with an approximate method, which has been called the semisteady-state solution. The pressure and saturation profiles from the wellbore to outer boundary calculated by the two methods are compared for one constant-terminal-rate case and two constant-teminal-pressure cases. The agreement in these profiles, as well as in the values of average reservoir pressure and cumulative recovery, leads to the conclusion that, for most engineering calculations, the semisteady-state method will give a reasonable approximation to the numerical solution of the differential equations describing solution-gas drive. An unfavorable (as regards ultimate oil production) set of relative permeability curves was used in the calculations in the belief that the effect of the parameters which were studied would be emphasized to a greater degree. Furthermore, the reservoir was assumed to be completely homogeneous, and these results should not be considered applicable to any other type of reservoir. Gravity effects are not considered. The absolute permeability was varied from 25 to 0.5 md. At an economic limit of 2 B/D, the recovery for a 25-md reservoir is about 1.8 times as great as that for a 0.5-md reservoir. The effect of permeability on the producing gas-oil ratio is minor. Once PVT properties of the oil and the relative penneability properties of the reservoir are fixed, the pro- during gas-oil ratio is found to be a function of the fraction of oil recovered. Well spacings of 10, 40 and 80 acres were considered. For the assumed homogeneous-reservoir properties, the effect of spacing on recovery at an economic limit of 2 B/D was very slight. Certain dimensionless groups can be used to extend the results to other fields having different penneabilities, spacings, reservoir thicknesses, well radii and porosities, so long as the PVT and relative permeability properties are similar to those used in this paper. INTRODUCTION An important aspect of reservoir engineering is the prediction of the performance of the reservoir based on the limited information normally available early in the life of a field. Usually, soon after a field has been discovered, it is necessary to decide upon the spacing to be employed and the production program to be used for most efficient utilization of reservoir energy. Where these decisions have not been dictated by legal or political considerations, they have frequently been based on experience factors which indicate that two fields with similar characteristics will probably have similar performances. To the extent that the industry has been successful in economically exploiting oil reservoirs, this rule-of-thumb method has been found to have merit. A large amount of theory has been developed with which it is possible to predict the performance of a reservoir by using certain known properties of the oil and the formation. Because of the difficulties of the mathematics involved for a drainage area having square or rectangular boundaries, the problem usually has been simplified and idealized by assuming each drainage area to have radial symmetry and the field to be represented by a number of these drainage areas. Muskatl was one of the first to formulate the theory for two-phase flow and to solve the equations for a few cases by numerical integration, using the relatively slow means of computing available at that time. His equations describing

By Peter A. Allen

INTRODUCTION Lac Minerals Ltd., a significant North American producer of gold and other metals for over three decades, recently completed the financing for its third discovery in an eight year period. From 1975 to 1982, three important gold deposits in Quebec and Ontario: namely Bousquet, Doyon and Hemlo, were located and tested by company geologists. When current development plans are complete, Lac will have injected $400 million into these deposits. A discussion of these projects and their financing reveals some extraordinary challenges faced by top mining management. Responses to these challenges developed by Lac management may interest executives planning to proceed with new North American gold or base metal mines. MINE DE BOUSQUET In 1976, we explored a property containing a large low grade surface gold deposit. The geologial setting was exciting and we tested recent rethinking about the volcanigenic genesis of gold and base metal occurrences. These ideas helped pus to persist in drilling a nearby conductor and 600,000 ton, 0.22 ounce/ton orebody was outlined. With gold at $Cdn.150 per ounce the feasibility was approved and work started February 1978. Availability of a company custom mill 35 miles away reduced our immediate capital needs to $8.5 million. First of all, understanding the geological significance of results at Bousquet was critical. Mine and surface exploration geologists at Lac are continually and closely working together to better understand the nature of gold deposits in the Shield. This management emphasis on both pure and mine geology meant that a much larger ore reserve could confidently be predicted. (A proposal to quickly develop and ship the ore to a nearby company mine was discarded.) Instead of designing just for the initial plan of 500 tons per day we therefore included headframe, shaft and hoisting capacity for 1,700 tons and provided a compartment for future shaft deepening. These prebuilt additions were welcomed within two years. The type of collaboration whereby surface and mine geologists make study and field projects of each other's problems sets the foundation for this predictive ability. At Lac, the traditional responsibility for and supervision of the mine geologist by the mine manager is modified by this collaboration. This adds interest and scope to the geologist's work and lets him think beyond production requirements. Besides such management considerations, dedicated, professional geologists with integrity and lots of practical experience are needed to plan successful developments like Bousquet. Second, projecting future bullion revenue was critical and our experience selling production for future delivery gave us a way to ensure that loan repayment at term would not be jeopardized by falling metal prices. Output from a single mine is often reduced or curtailed for numerous reasons, so "selling forward" is done best and with least risk by companies like Lac with more than one or two producing locations for that metal. At first, a forward marketing program was considered unwise, speculative and unorthodox by analysts and other producers. Declining markets in the early 1980s prolonged gold project capital recoveries however, and critics were better able to appreciate the usefulness of such a tool. The notion that gold miners shouldn't sell forward is rooted in the belief that long-run gold prices will rise as modern social theory implementation in democracies inevitably erodes currency values. Selling against this trend goes against this industry belief or "ethic" and is at odds with the company shareholder who expects the gold price to rise not fall. Nonetheless, an efficient tool is available to the gold producer (not to most other metal producers) to quickly disengage from adverse price trends or to lock in spot prices instantaneously for two years or more of profits. Gold had recently retreated from U.S.$200 to $104 and now had reversed to $150. Lac had been forced to close several of its producers during the late 1960s because central bank-controlled

Jan 1, 1985