Search Documents

By H. G. Doll

A new electrical logging method called Laterolog is described which provides for better recording of formation resistivity. In this method a current, preferably of constant intensity, is forced into the formation perpendicularly to the wall of the hole as a sheet of predetermined thickness by means of a special electrode arrangement and of an autonlatic control system. The essential advantages of the Laterolog over the regular resistivity logs are explained. With the Laterolog, the mud column has much less influence than with conventional methods. The Laterolog is, therefore. particularly useful in wells drilled with highly conductive mud. Moreover. when the electrode system is located opposite a bed of moderate or small thickness, the disturbing effect of the adjacent formations becomes practically eliminated, provided the bed is thicker than the sheet of current. which itself is usually a few inches to a few feet thick. The sequences of beds are. therefore, much more sharply differentiated, and, in many cases, it is possible to read directly from the logs values close to the true resistivities of the formations without further corrections. Field examples are shown to illustrate the results of the method in various types of formations. INTRODUCTION As soon as it was introduced into the petroleum industry, electrical logging roved to be a powerful instrument for the delineation of the strata traversed by bore holes and for correlation of such strata." Research and engineering efforts were mostly devoted at this early time to the design of the techniqurs which could provide the best data for these purpoes. Later, attention was directed toward the quantitative analy5i.s of electrical logs, in order to obtain information on. reservoir characteristics, and, in particular, on oil and water saturation.3,4 Efforts have been made, accordingly, with a view to improving the already existing methods, or to finding new methods which could satisfy these new requirements." Quantitative analysis implies a determination as accurate as possible of the true resistivity of the formations. In fact, however; the "apparent resistivity" recorded in front of a given bed by a conventional electrical logging method is frequently quite different from the true resistivity of this bed. This is due to the combined influence of the mud column, of the adjacent formations both above and below the bed, and of the invaded zone, in which the original fluid has been more or less replaced by mud filtrate. When the beds are thick enough, and, as is frequently the case in sand and shale formations, when their resistivity is not much greater than that of the mud, it is generally possible to obtain their true resistivity directly on the electric log, or at least to obtain an approximate value which can be corrected without undue difficulty by means of recently published charts known as "departure curves."" When the beds are thin, and particularly when, in addition. their resistivity is substantially greater than that of the mud column, the apparent resistivities recorded on conventional logs are seriously affected, and it is generally very difficult, if at all possible, to estimate their true resistivity with a good enough accuracy. In the most severe cases, such as are encountered in holes drilled through hard formations with muds of high salinity, the conventional logs are so distorted and rounded that they do not even show clearly the houndaries of the different beds. In order to overcome the effect of the mud column and of the adjacent formations, it is possible to make use of appro. priate systems whereby the measurements involve a comp3ratjvelv thin slice of space at the level of the measuring device. A mention is due, in this respect, to a measuring system. called the "guarded electrode;" invented by C. Schlumberger over 20 years ago. The system, schematically represented in

Jan 1, 1951

By Fred H. Poettmann

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

Jan 1, 1951

By M. Muskat

Calculations are reported on the differential sensitivity of the updip invasion of oil strata of varving permeability to the driving pressure differential. It is assumed that the water-oil interfaces advance with sharp fronts, or that the microscopic displacement efficiency is independent of rate and capillary pressure effects. It is found that, under such conditions, limitations of withdrawal rates will not greatly inhibit the "fingering" and bypassing tendency through high permeability strata unless the differential fluid head across the original oil column is of the order of or exceeds about 80 per cent of the driving pressure differential. In practice, such restrictions will be feasible only for oil reservoirs and aquifers of high permeability and dip, but may develop automatically under combined gas injection and water drive operations. INTRODUCTION A commonly expressed view regarding the intrusion of edge-waters into oil producing formations is that the geometry of such advance is rate sensitive. Specifically, it has been suggested that high rates of advance, stimulated by rapid withdrawals, will lead to "fingering" and an irregular type of invasion, whereas at low withdrawal rates the water encroachment will be more uniform. Two factors tend to support this view. One is the resultant effect of capillary pressure and gravity gradients in the displacement processes' involved in water intrusion, which are magnified at low rates of flow so as to improve the recovery efficiency. The other is the direct retarding effect of the rising water head as it advances upslope into the producing area. Capillary pressure forces give rise to a rate sensitivity of the oil recovery by water displacement in uniform strata as well as differential effects in non-uniform formations. These are basically independent of the effect of gravity and would also occur in horizontal beds or when the oil and water have the same density. Unfortunately, however, the fundamental equation describing the displacement process under the action of capillary pressures can be solved only by numerical or graphical procedures, and moreover the solutions are sensitive to the character of the capillary pressure and relative permeability curves.* On the other hand, the direct effect of gravity can be given a formal analytical solution if the simplifying assumption is introduced that the water advances with sharp fronts in the individual uniform strata. It depends on the absolute permeabilities of the water invaded zones and hence results in a rate sensitivity of the differential invasion in stratified media. Although a tieatment of this factor alone will not give a complete evaluation of the effect of rate on water intrusion in non-uniform pays, it should provide at least a semiquantitative measure of the sensitivity of "fingering" to withdrawal rate. Such a treatment will be presented in the following sections. GENERAL THEORY The basic element of the theory is the history of water advance updip in a uniform stratum (Fig. 1). In contrast to the assumption of constant flow rate, generally made in the analysis of displacement processes, we shall assume here that

Jan 1, 1951

By H. E. Stamm, H. H. Spain, R. C. Rumble

This paper presents a reservoir analyzer study of the performance of the Woodbine formation in the East Texas basin. The study was made possible by the compilation of available information on the configuration. thickness and the pressure drawdown of the formation. The investigation was made of the pressure variations in the basin incident to production from the several Woodbine reservoirs. In addition, the apparent compressibility of Woodbine water was evaluated so that the potential water yield of the formation could be determined. The distribution of permeability of the Woodbine sand and the interference between producing areas were also investigated. In conjunction with this basic study of the Woodbine sand, an acceptable match of the production-pressure relationship in the East Texas Field was established on the analyzer. INTRODUCTION When production is taken from a reservoir contiguous to an aquifer, the resultant pressure gradient causes a water influx into the reservoir. For some years it has been the practice to predict the performance of water drive reservoirs with mathematical equations that relate the water influx into reservoirs to their pressure behavior. A rigorous solution of these equations is quite complex: hence, the conventional method requires the simplifying assumptions that the formation have constant thickness, permeability and poiosity, and that it have a known and regular shape. It is necessary to use performance records of the reservoir in question to determine the equation constants that make the results of the mathematical calculations duplicate the reservoir history. The values of the constants determined in this manner sometimes differ so greatly from theoretical or measured values that the constants appear unreasonable. The equation constants comprise such factors as the thickness, porosity and permeability of the formation, and the viscosity and compressibility of water. Usually the thickness, porosity and viscosity are known within a reasonable degree of certainty. whereas the formation permeability and the water compressibility are known with considerably less certainty. Hence. permeability and compressibility values are arbitrarily selected that will make the mathematical equations duplicate the reservoir behavior. In many reservoir studies the value of the compressibility of water thus selected has exceeded the reported value by a factor of from severalfold to manyfold. Whether this excessive value results from the inadequacy of the mathematical approach to take into account the shape irregularities and the heterogeneities of the aquifer, or is due to some type of formation compaction and decrease of porosity has not been established. Because of the significance of the potential water yield of formations with reference not only to oil reservoirs but to water resources as well, it is very important that a representative value for the effective water compressibility be determined. One of the best known examples of a field for which it has been possible in the past to predict the future reservoir pressure is the East Texas Field. Until about 1943 it was possible by means of mathematical equations to duplicate past pressure behavior and to predict with a high degree of accuracy the reservoir pressure of the field. Since 1943, however, the calculated pressures have deviated progressively from the observed pressures to such a degree that in recent years, it has

Jan 1, 1951

By B. J. Dotson, P. N. McCreery, R. L. Slobod, James W. Spurlock

A core sample porosity-check program is described. A number of laboratories participated in the investigation. which comprised measuring the porosities of ten selected natural and synthetic core samples. Each laboratory employed its own method (or methods) of measurement. These included Boyle's Law, water-saturation. and organic liquid-saturation methods. Brief descriptions of each method employed are given. In general, the measured values of porosity were in reasonably good agreement. The average deviation of porosity from the mean or average values for the group of samples was ±0.5 porosity per cent. The results of these measurements are compiled and presented in graphical form. INTRODUCTION The total hulk volume of a porous sample is derived from the sum of two volumes, the void or pore volume and the solid* or grain volume. The porosity fraction of the sample ia defined as the ratio of the void volume to the bulk volume. Multiplication of this ratio by 100 converts it to per cent porosity. Distinction is made here between total porosity and effective porosity. The former includes all void space, both interconnected and isolated, while the latter includes only those voids which are interconnected. This investigation deals only with the effective porosity.* Numerous cases had been noticed in which there appeared to be appreciable discrepancies in the results of porosity determinations on core samples. Differences of as much as 10 to 15 porosity per cent had been reported from different sources

Jan 1, 1951

By Walton J. Greer, System By John E. Owen

The guard electrode system measures the resistivity of formations by employing a thin disk of current which is caused to flow perpendicular to the bore hole. The control of this current disk is obtained in a brute force manner through the use of relatively long equipotential electrodes above and below the measuring electrode. The log obtained from the system is free of "lag," "plateau," "shadow," "reflection" and other distortions evident on conventional logs. The proportionate contributions of the hole size and invaded zone resistivity to the apparent resistivity reading are reduced when the mud and filtrate resistivities are lower. Thus, the guard system favors the use of conductive muds. The mathematical development of the theory is given in an appendix. Equations are derived to evaluate the effects of mud resistivity, hole size, invaded zone resistivity and depth, and the true formation resistivity. Sample logs from several provinces are reproduced and compared with conventional electric logs to illsutrate the much greater detailed lithology that can be presented by the guard electrode system. INTRODUCTION With increased ability of those in the petroleum industry to analyze formation characteristics from electrical resistivity data, there has become apparent the need for more precise measures of that parameter. The measurement of the resistivity of a homogeneous and isotropic material can be made simply when the body of the material is infinite in all dimensions and the measuring electrodes can be placed within the body. In such a case, the current leaving an electrode follows a spherical distribution and the mathematical expression for the potential level at any distance from that electrode can readily be derived.'.' If it is considered that the active electrodes are within a bore hole filled with mud of a resistivity different from that of the formation. the current does not flow in a spherical pattern near the current source or sink, and the calculations of potential levels become more difficult. The additional consideration of a section of formation invaded by mud filtrates, and thus having resistivity different from either that of the mud, the adjacent formations or of the undisturbed formation, makes the mathematical problem extremely complex. Should the vertical dimensions of the formations be considered as finite. the problem becomes so complex as to preclude reasonable solution. Yet, even these considerations have not described adequately the conditions encountered in a large proportion of the wells to be studied. In efforts to evaluate the effects of the various conditions. combinations of twn, three. and four electrode logging con- figurations of the conventional types have been used. The resulting logs have left much to be desired in some provinces. It has become apparent that the use of such configurations will not result in adequate measurements being made, particularly in cases where the formation is composed of thin sections. If all the current of a resistivity measuring circuit is made to flow through the formation opposite the exploring electrode, many of the shortcomings of the conventional measuring system are eliminated. In the guard electrode method of electrical well logging to be discussed in this paper, the current is so controlled by causing it to flow in a thin disk normal to the axis of the logging tool. EARLY HISTORY The guard electrode method of electrical well logging in its early years was a joint development by the Geophysical Research Corp. and the Amerada Petroleum Corp. The method was first experimentally investigated by means of model studies. The results of these studies indicated that current focusing systems properly applied to electrical well logging

Jan 1, 1951

By R. M. Leibrock, J. E. Huzarevich, R. G. Hiltz

The various factors considered in recommending the initiation of a gas injection project in the southern portion of the Cedar Lake Field are discussed. Performance history under gas injection operations is reviewed and these data are analyzed, utilizing both the material balance method and the fractional flow and frontal advance expressions. Results of the analysis of the performance data indicate that the injected gas has contacted and affected at least 60 per cent of the reservoir and a substantial increase in ultimate recovery can reasonably be expected. By holding the reservoir pressure appreciably above the bubble point, the well productive capacities have been maintained substantially above the level predicted for primary operations. The analysis of the Cedar Lake project suggests that in certain limestone reservoirs, at least, the probable success of gas injection cannot be predicted simply from ohservation of permeability distribution throughout the pay section, as indicated by core analysis data, on either one or a number of wells. Further, the performance of this particular project fails to indicate any basis for classifying carbonate reservoirs in general as being inherently unsuited to a dispersed type gas injection program, thus indicating that each reservoir should be considered on its own merits, regardless of the composition of the reservoir rock. INTRODUCTION Early in the life of the Cedar Lake Field, an extensive data gathering program was initiated to provide an accurate record of reservoir performance characteristics. From the study of these data it was apparent that there was a critical need for supplementing the natural reservoir energy in order to maintain well productivities and obtain the maximum ultimate oil recovery. Accordingly, detailed engineering studies were made of the various methods of secondary recovery which might be applicable. As a result of these investigations, the decision was made to initiate a gas injection program of sufficient intensity to maintain reservoir pressure at approximately 600 psia, or some 274 lb above the bubble point pressure of 326 psia. A full scale dispersed type gas injection program has been in operation on leases of the Stanolind Oil and Gas Co. in the southern portion of the field for nearly five years, and sufficient performance data are now available to evaluate the benefits which have been derived from this project. It is the primary purpose of this paper to analyze the performance data for the Cedar Lake gas injection project and to point out the significance of the ohserved behavior with respect to certain hypotheses which have been advanced in recent years concerning the probable success of gas injection projects in limestone reservoirs. This paper properly should be regarded more on the order of a progress report, inasmuch as some revision in interpretation will undoubtedly be required from time to time as additional performance data are obtained, although the satisfactory performance of the project to date leaves little doubt as to the ultimate success of gas injection in the Cedar Lake Field. As a result of the success of the project to date, a unit was formed in the southern part of the field, effective March 1, 1951, for the purpose of continuation of the gas injection program. Participants in this unit are the Mid-Continent Petroleum Co. and Stanolind Oil and Gas Co. GEOLOGY AND STRATIGRAPHY The Cedar Lake Field is located in the northern portion of the Midland Basin area as shown in Fig. 1. The southwest portion of the field lies within a playa, or dry salt lake, which covers an area of approximately eight square miles. As might be expected, it was this lake which furnished the inspiration for the name of the field. Except for its value as a salt water disposal pit, this lake has succeeded only in magnifying the difficulties in developing this portion of the field. Typical of this section of West Texas, the area in general is relatively flat and has a semi-arid climate. The localized structure which favored the accumulation of oil is an anticline with approximately 100 ft of closure. The major axis of the structure extends in a general southeast-northwest direction. Originally this structure was defined by seismograph data, which have been subsequently confirmed by development. In general, the geologic column is typical of that found throughout the basin. From the surface to depth of approximately 1,800 ft, surface sands and undifferentiated red beds. probably Triassic. are encountered. Below this point to the producing horizon, all formations are of the Permian age.

Jan 1, 1951

By M. R. J. Wyllie, F. Morgan, P. F. Fulton

The importance of formation factor, F, not only in electric logging but as a fundamental rock parameter has recently been stressed.',: The desirability of investigating the range of variation of the resistivity index exponent, n, in the relationship I = S ;", where I is the resistivity index and Sw the water saturation as a fraction of the void volume of a porous medium, has also been urged.3 The magnitude and variation of n with saturation and rock texture is a subject not only of theoretical interest but also one of prime importance in the interpretation of electric logs. A simple technique has recently been developed which enables both F and u to he measured with high accuracy and which may also find acceptance as a convenient method for the determination of irreducible saturation attainment in the restored state method of core analysis. Experience has taught that reproducible measurements of F are possible only if the resistance measuring electrodes are so arranged with respect to a plane face on a porous medium that they are able to make electrical contact with substantially all entry pores in that plane. In practice this may be achieved by using a platinized-platinum gauze electrode backed by some absorbent material (such as felt) which has been saturated with a fluid identical with that used to saturate the porous medium. Applicatiorl of pressure to the electrode and absorbent material then forces the gauze against the plane face of the porous medium and simultaneously squeezes saline solution through the meshes of the gauze. By this means the electrode is in continuous aqueous contact with all pores and satisfactory and reproducible low resistance contact with the porous medium is achieved. Clearly this method, although satisfactory for measurements of F, cannot be applied to the making of continuous resistance measurements on a porous medium while capillary pressure desaturation is being carried out. However, accepting the principle that for satisfactory results electrical contact must be made between a measuring electrode and all pores adja- cent to that electrude, methods of bringing electrodes into intimate contact with the surfaces of porous media were investigated. Two methods were ultimately found to be satisfactory: in the one, the metal electrode is formed on the required portion of the porous medium by the use of a metal spray gun, while in the second the electrode is painted on with an ordinary camel's hair brush. The first method has the advantage of permitting the use of any metal which can be sprayed, but has the disadvantage of requiring rather elaborate and expensive equipment. The second method is presently limited to silver electrodes although in principle other metals, e.g. platinum or gold, could be used. Moreover, the method is so simple and cheap, and has been found to be so satisfactory that it will be described in some detail. The core being investigated is cut into a right circular cylinder and is extracted and dried in the usual manner. The ends are then lightly painted with silver conducting paint* of the type used in printed electrical circuits. The quantity of paint used is not critical but the recommended, minimum compatible with entirely coating the core ends is recommended, particularly on the end that contacts the barrier. The core is then dried at atmospheric temperature for one hour or for shorter periods at any suitable elevated temperature up to about 110°C. It will be found that silver coatings so prepared are not only strongly adherent but also permeable and the core may be the core may be desaturated by the ordinary capillary pressure technique through one of the coated faces. The same permeability is characteristic also of thin metal coatings formed using the spray-gun technique. An ordinary Lucite capillary pressure desaturation cell has been adapted to a form suitable for measuring the resistivity of the saturated silver faced cores both at 100 per cent saturation (i.e., F) and at intermediate saturations down to the irreducible minimum. This has been achieved as follows: A Coors porcelain barrier, having a displacement pressure of c 30 psi was grooved across a diameter. Dimensions of this groove were c 1 mm deep and 1 mm wide at the top. The groove was then painted thickly with silver conducting paint, the paint in the groove being extended lightly over the edges of the groove for a distance of c 1 mm on each side. A 30 gauge silver wire was then arranged in the groove in a form of a spring bow, each end of the silver being held at diamet~ically opposite ends of the groove by means of plastic cement. The arc of the bow at its highest point was arranged to be a millimeter or so above the face of the barrier, while one end of the bow wire was led by means of a pressure-tight connection through the wall of the capillary pressure cell. The groove in the barrier was then Surrounded by suitably cut portions of Kleenex in the conventional manner so as to ensure capillary continuity from core to barrier, and the core placed on the barrier. The weight of the core distorted the silver spring bow and good electrical contact was thereby made between the outside of the cell and the lower painted silver electrode. Electrical connection to tile top painted silver

Jan 1, 1951

By Walter Rose

It is the purpose of this note to call attention to the circumstance complicating the attainment of uniform gas-liquid distributions in multi-phase flow systems, and especially in those of the so-called Hassler type.' Although the complication arises directly as a consequence of gas compressibility and the dependence of fluid distributions on interfacial curvature phenomena, little reference to the problem has yet appeared in the literature.' In fact, it is only the paper of Geffen et al.3 which gives any experimental evidence that a problem exists. On the other hand, previous authorsl,4,5 have made the tacit, albeit fallaceous assumption that uniform gas-liquid distributions automatically are established in linear flow systems by the expedient of maintaining the pressure gradients equal in the flowing gas and liquid phases. Such a device, it can be shown, will succeed only if the immiscible fluids are both equally compressible or both essentially non-compressible. Another necessary condition required before uniform fluid-fluid distributions can be achieved is that the porous matrix is isotropic. Uniform fluid saturation distribution conditions obtain in isotropic media only when the curvature of each interface of contact between wetting and non-wetting fluids is everywhere the same throughout the interspaces. This is a necessary condition which follows from the consideration that variation in interfacial curvature gives rise to finite capillary pressure gradients and therefore to variation in the saturation distribution. That this is not a sufficient condition follows from the consideration that hysteric possibilities allow for different saturations even though the capillary pressure gradient is zero. In any event, it will- be recalled that it is the intent in the Hassler scheme of relative permeability determination' to obtain initially (by the capillary pressure drainage or imhibition process) uniform conditions of fluid distribution in the core sample, and then to maintain this uniformity during mixture flow so that the resultant fluid mobilities which are calculated will refer to steady-state transfer under fixed conditions of uniform saturation. If the flowing fluids are incompressible, it is recognized that employment of the same value of the pressure gradient in each fluid will result in maintenance of the initially obtained condition of zero capillary pressure gradient. On the other hand, the consequence of gas compressibility is that the gradient in the gas pressure varies from point to point in the flowing stream. and since the pressure gradient in an incompressible liquid is constant it is thus impossible to maintain zero capillary pressure gradient during the simultaneous flow of gas-liquid mixtures. Therefore, one would expect to observe a variation in the gas-liquid saturation in the directions of the capillary pressure gradients, depending in magnitude on the particular way saturation changed with capillary pressure. In order to formulate the manner gas compressibility effects prevent the attainment of uniform gas-liquid distributions during mixture flow, consider a linear flow system of unit length. As noted, it is necessary to assume an isotropic medium (viz. one where permeability is independent of position and direction). One then can arbitrarily set the pressure

Jan 1, 1951

By H. F. Dunlap, R. R. Hawthorne

A method of calculating formation water resistivities from chemical analyses is presented which is somewhat faster and more accurate than previously described methods. For 26 formation water samples taken from wells in the Texas Gulf Coast, the average deviation of calculated from measured water resistivities was 3.3 per cent, with an average bias of + l.l per cent in the calculated values. This compares with average deviations of 4.6 per cent and 7.8 per cent, with corresponding average biases of +3.2 per cent and -6.3 per cent, for the two methods previously described in the litera-

Jan 1, 1951

By Leonard B. Lipson

Boyer, Morgan and Muskat,' with improvements by Morgan, McDowell, and Doty,l and Laird and Putnam,' have described a scheme for the determination of fluid saturatic in porous media by measuring the x-ray linear absorbing power of particular fluids in the media. This method is one of several which depend on the linear absorbing power of the fluids in the porous medium. Other examples of linear absorption methods involve measurement of dielectric constant, gamma ray absorption, and ultrasonic attenuation. It is the purpose of this note to demonstrate by theortical considerations that the absorption versus saturation curve for a linear absorption method is not necessarily a single valued function of the saturation. Such non-unique behavior must be considered and evaluated before suitable saturation determinations can be made by a linear absorption method. To illustrate the discussion, examples will be given only for the simplest cases of combinations of two different dielectrics. The dielectric constant method was chosen for illustration because it was thought to demonstrate most graphically the point in question. According to Argue and Maass,4 the resultant dielectric constant of a combination of two different dielectrics depends on the orientation of the dielectrics in the electric field of the test condenser used to measure the dielectric constant. One ease involves alternate slabs of dielectrics E1 and E2 perpendicular to the electric field of the test condenser, as shown in Fig. 1A. For this case: l/'E, = A/E, + B/E,.......(1) where E,,, is the dielectric constant of the combination, A is the ratio of the path length through dielectric E, to the distance between the plates of the condenser (assuming the dielectrics have the same cross sectional area as the condenser), and B is the ratio of the path length through dielectric E2 to the distance between the plates of the condenser. It will be noted that the constants A and B are numerically equal to the fractional "saturation" of the region between the test plates of the condenser by the two dielectrics E1 and E2, respectively. Another case involves alternate columns of the two dielectrics parallel to the electric field and forming a continuous path between the condenser plates, as shown in Fig. 1B. For this case: E., = AE1 + BE2.........(2) where E., is the dielectric constant of the combination, A is the ratio of the cross sectional area of dielectric E, to the cross sectional area of the condenser plates (assuming the lengths of the dielectrics equal the distance between the condenser plates) and B is the ratio of the cross sectional area of dielectric E2 to the cross sectional area of the condenser plates. As before the constants A and B are the effective "saturation" of the condenser by the two dielectrics. Still another case is described by Lichtenecker5 where dielectric E2 is shown in the form of small particles distributed at random in dielectric E1. This is shown in Fig. 1C. For this case: log Em = A log E1 + B log E2.....(3) where Em is the dielectric constant of the combination. A is the ratio of the volume of dielectric E1 to the volume of the condenser, i.e., the effective "saturation" of the condenser by dielectric E,. B is the ratio of the volume of dielectric E2 to the volume of the condenser, i.e., "saturation" of dielectric El. Equations (l), (2) and (3) are plotted in Fig. 2 for E1 = 2.5 (oil) and E2 = 80 (water). Equations (1) and (2) may he thought of as representing the limiting cases of the configuration of two fluids which might be expected to exist in a porous medium. Equation (3),

Jan 1, 1951

By Henry A. Stiff

The classification and correlation of water analysis data presents many problems which can be solved by graphic methocis. The pattern system, a new type of graphic procedure described in this communication, is believed to have several advantages over older methods. Examples of the application, of the pattern system to the solution of problems encountered in petroleum production are given. INTRODUCTION Several graphic methods for presenting analytical water data have been developed and are now in use.1,2,3,4 It is believed, however, that a recently developed type of graph called the "pattern" offers several advantages over other methods. This system presents a better picture of the total salt concentration than is usual in such graphs. The effect of dilution or concentration has been reduced to a minimum, and at the same time distinction between various types of water has been improved. The system is extremely versatile, yet so simple it can be plotted on ordinary graph paper and adapted to almost any type of filing system. DISCUSSION The essential feature of the pattern system is the graph shown in Fig. 1. Horizontal lines extending right and left from a vertical line at zero form the graph. Positive ions are plotted to the left while negative ions are plotted to the right. The figure immediately beneath each ion gives the scale. Most oil field waters can be plotted on a scale where 100 milli-equivalents of sodium and chloride and 10 milliequivalents of each of the other ions are represented by one scale unit. For highly concentrated brines, 1,000 milliequivalents of sodium and chloride and 100 milliequivalents of each of the

Jan 1, 1951

By John R. Spencer

A capillary tube variable volume cell is described, which has operated satisfactorily over a range of 100°F and 3,500 psi to 550°F and 1,500 psi. The cell contents are entirely visible over the length of the capillary with the exception of a space at the top amounting to three per cent of the total volume. Observations on the critical behavior of simple systems have been greatly facilitated by the apparatus, it being possible to establish the critical conditions for a particular mixture in the course of a day's run. INTRODUCTION Many types of apparatus have been designed for obtaining data on the phase relationships of hydrocarbon systems. Such apparatus can in general be subdivided into two classifications; (l), indirect measurement of phase volume, and (2), direct visual observation of phase volume, each of which has outstanding advantages and disadvantages. Indirect measurement of phase volume requires the use of rather elaborate apparatus including internal probes in which are located the device for determining the locus of the interface between phases. Such equipment is by its very nature costly, and is complex from the standpoint of operating technique. Probing techniques can be used, however, over wide ranges of temperature and pressure with a high degree of reproducibility. Apparatus for visual observation of phase behavior obviously has the advantage of enabling the experimenter to see exactly what takes place, but because the material to be tested is contained in a cell constructed of both metal and glass, is limited severely as to the temperature range of operation. Furthermore. the maximum pressure usually is limited by the choice of material used for sealing the glass to the metal, and by the strength of the glass itself. In the interest of increased efficiency of observation, as much of the cell as possible should be visible to the eye. Glass capillary tubes have been used successfully but have always given difficulty at the metal-to-glass seal. An apparatus has been devised to overcome the problems of sealing and resistance to temperature effects that have heretofore been the limiting factor with glass apparatus in the observation of phase volumes at elevated temperature and pressure. APPARATUS The objective which dictated the general arrangement of the apparatus (Figs. 1 and 2) was simplicity of operation. The air bath, A, was constructed of an iron framework to which was attached approximately one inch of glass wool bats which were covered on the outside with "Masonite." A roughly octagonal center piece was installed to provide forced circulation of the air by the fan, B, past the heaters, D, cell, F, and thermoregulator, C, respectively. The entire top of the box is removable to facilitate charging of the cell and to permit easy maintenance. Two windows are provided before and behind the cell to allow the operator to observe the contents against the background fluorescent light, E. The heaters,

Jan 1, 1951

By M. R. J. Wyllie

In a recent publication1 wetting phase relative permeability was exvressed as: and it was stated that a similar expression applied, mutatis mutandis, to non-wetting phase relative permeability; i.e., In Equations (1) and (2), Sw and S0 are respectively the wetting and non-wetting phase saturations as a fraction of the pore volume; PC is the capillary pressure; I; the wetting phase electrical resistivity index; IN, the analogous non-wetting

Jan 1, 1951

By Kefton H. Teague

Attempts to locate domestic supplies of sillimanite have been unsuccessful until recently. This paper describes recent discoveries of sil-limonite-bearing schists in the Southeastern States, with emphasis on geology of the deposits, mineral composition, origin, and reserves. Statements concerning possible methods of mining, milling, and utilization of sillimanite are given. Results from refractory tests on a laboratory scale are summarized. SINCE 1917, when synthetic mullite was first used in spark plug porcelains, there has been a continued search for natural sillimanite. During the recent war this search was intensified, not only for sillimanite but also for any materials which could be used satisfactorily as a substitute for the imported kyanite from India. The Indian massive kya-nite does not decrepitate during calcination as does our domestic bladed kyanite, thus a coarse grog can be obtained from it. Some of the most promising domestic materials suited for high temperatures are limited in quantity. Massive topaz from Chesterfield County, S. C., showed promise of possessing desirable refractory properties; but the action of combined fluorine, liberated during calcination, upon equipment and surroundings, as well as a limited supply of ore and the presence of impurities discouraged its use. Residual lumps and boulders of massive kyanite, similar to the material imported from India, occur in Georgia. This material, like the massive topaz, is of limited quantity. Andalusite from California and dumortierite from Nevada are satisfactory high-temperature refractory materials. Their uses have been restricted to specialized products since the location and character of the deposits are such that large production is impractical. Sillimanite, like kyanite and andalusite, has the theoretical chemical composition of Al2SiO5. It has a hardness of 6 to 7; vitreous, silky to subadamantine luster; specific gravity of 3.2 to 3.3; and is gray to bluish gray in color. Under the binocular microscope, the crystals are transparent. The mineral may occur as dense, fibrous mats (fibrolite) composed of fibrous, sometimes radiating hair-like crystals in schist which are associated with igneous intrusions. Deposits of this type have widespread distribution, but none discovered to date appears to have commercial possibilities. Sillimanite also occurs as bundles of crystals disseminated in biotite schist, as noted in the Hart-Elbert-Madison County, Ga., area and in some of the South Carolina and North Carolina deposits. This type of material has been called prismatic sillimanite by the U. S. Bureau of Mines. Deposits of this variety appear to offer commercial possibilities. When heated above 1650°C, sillimanite expands about 6.5 pct and is converted into a mixture of mullite (3Al2O3 . SiO2) and vitreous silica. This is a stable mixture which has a specific gravity of about 3.15. General Geology: Sillimanite (fibrolite) is a mineral common to the pre-Cambrian metasedi-ments of the Southeast. It occurs throughout the crystalline schists which lie between the Fall Line on the southeast and the established Paleozoic rocks on the northwest (fig. 1). In addition to the meta-sediments, this area contains pre-Cambrian igneous

Jan 1, 1951

By S. A. Falconer

This paper deals with the practical application of the Dutch State Mines cyclone separator for fine-coal cleaning. The more important operating variables are discussed, and results of a number of continuous-scale tests on various fine coals with a 6-in. cyclone separator are given. Descriptions of the Cyanamid Pilot Plant and the 50 tph Cyclone Separator Plant at the Emma Washery in Holland are included. TODAY, as never before, the coal producers of this country are faced with the necessity for obtaining maximum yield of marketable coal of all sizes down to and including the finest. Particular interest is being displayed in a more effective means of recovering the finer sizes which in the past so often have been discarded to waste or treated by inadequate equipment. The reasons for this interest are easily understood. Rising mining costs have resulted in the increased use of mechanized equipment underground for full seam mining. This, in turn, means more refuse in the run-of-mine and also more fines in the feed going to the cleaning plant. Stringent antistream pollution laws now make it necessary at least to impound the fines, and usually it is a paying proposition to recover the coal from these fines before sending them to the storage area. Finally, the market and price of the finer sizes has increased, but at the same time specifications have tightened and the demand is for higher grade products. During the past four years a great deal of publicity has been given to this new process and apparatus for treating fine sizes of ores as well as coal. Excellent discussions of the apparatus itself and the theoretical aspects of the basic operating principles involved in its operation have been presented by M. G. Driessen,1-3 M. R. Geer and H. F. Yancey,4,5 J. W. Hyer6 and more recently by D. A. Dahlstrom8 as well as others.7,9-11 In view of this, it has been decided to include in the present paper only a brief description of the cyclone separator and its operating principles. This paper deals primarily with some of the more practical commercial aspects of the process and equipment required for its operation; the variables in the operation of the cyclone, their effect and means for control; and the results obtained with a 6-in. semicommercial size separator on a variety of fine coals, bituminous and anthracite, using finely ground magnetite for the separating medium solids. A description is included of the commercial-size cyclone separator pilot plant installed by the Dutch State Mines at its Emma mine in Holland, where finely ground slate is utilized as medium solids. Field of Application: The Dutch State Mines cyclone separator operates most efficiently on ores and coal in the size range of % in. to about 48 mesh. On some types of feed, effective separation can be made on sizes as small as 100 mesh, but up to the present, this may be considered exceptional. Depending upon the type of medium used, the specific gravities of separation can be varied within wide limits, up to as high as 2.70, if magnetite is used as medium. Early Development: The development of the Dutch State Mines cyclone as a separator, or cleaner, was the outgrowth of early experimental work by Dries-sen' and others in Holland, just prior to World War 11, in the course of which was discovered the effectiveness of the cyclone in reclaiming and thickening the dense medium used in the "Loess" process of

Jan 1, 1951

By William B. Stephenson

The paper describes efficient and effective methods of automatically controlling sand pump installations. Particular reference is made to liquid-level controls actuating variable speed pump-driving units. Included is a discussion of various combinations of equipment, with schematic drawings. HERETOFORE many sand pump installations in the mining and metallurgical field involved guesswork as to the size of the pump and then simply direct-connecting it to a constant speed motor which was often large enough to do twice the work required. Considerable thought was given to mill location in an effort to avoid as much pumping as possible. Hillside design of mills with gravity flow of tailing was much to be desired. Power costs were low, construction was relatively cheap, and labor in many places meant nothing in comparison to the values derived from the ores being treated. Today most of these conditions are reversed. Ores are becoming more complex with lower values per ton of ore treated. Power is more expensive, construction is higher, and labor definitely takes a very important place in the profit and loss statement. Each piece of equipment incorporated in the design of a new mill is selected to provide highest efficiency based on careful metallurgical analysis of the ore to be treated both in good and poor markets. Each is designed with its flexibility in mind and with provision for adjustments in its performance to meet actual conditions in the operation of the mill. Sand pumps have taken an important position in mill design. Operators the world over are relying more and more on the performance of a sand pump in the efficient operating of their plant. Many ordinary jobs are being done today with pumps driven through adjustable pitch diameter V belt drives in order that changes in pump speed can be made easily to meet actual operating conditions. Operating costs have long since proved the folly of operating a pump at too high a speed for the volume and head to be overcome. Metallurgical results have shown the benefits derived from a smooth flow as contrasted with a surging flow from an uncontrolled pump. Construction costs have risen to such an extent that surge bins and other means formerly employed must be omitted from today's flowsheet. Several installations of considerable size and involving the use of automatic pump speed control have proved the adaptability of such controls since they were placed in service in recent years. In the following presentation each system discussed and illustrated schematically has been installed in one or more instances and has proved itself reliable and definitely economical from the standpoint of capital investment, operation and maintenance. Variable Volume-Head Speed Control: Probably the most popular, effective and efficient automatic control setup for sand pumping installations is the one shown schematically in fig. 1. Commonly used

Jan 1, 1951

By C. E. Lesher

THIS paper reports a study of the reactivity of 950°F and 1650°F cokes as measured by relative rates of reduction of iron oxides at temperatures up to 2200°F. Previous work cited shows general acceptance of the theory that reduction by carbon is a gaseous reaction, and that kind and character of carbon as well as particle size have measurable effect on the velocity of reaction. As will be shown, the data obtained in this study confirm those conclusions. The work was not designed to examine iron oxide reduction equilibrium, but if reaction velocity be defined as the speed with which "a reaction tends to approach conditions of equilibrium," the data here presented may be considered as a study of reaction rates, and the relative degree of reduction to metallic iron as the measure of reactivity. Three standardized combinations of Lake Superior brown iron ore with carbon were tested by similar procedures. One combination was a mechanical mixture of carefully sized high-temperature coke (1650°F) with the ore. The second was a mechanical mixture of the ore with Disco* obtained by carbonizing the identical coal at 950 °F. The third was an agglomerate prepared by carbonizing the coal and ore at 950°F, premixed in proportions to give as nearly as possible the same relative amounts of carbon and ore as the mechanical mixtures. This agglomerate, obtained by heating the finely divided ore (through 30 mesh) with coking coal through the plastic temperature range so as to form solid aggregates, gives a product in which the oxide particles are impregnated with, and intimately bound together with low-temperature coke. The agglomerate-—ore-Disco—was most active in oxide reduction; the mechanical mixtures of Disco and ore next in order, with coke the least reactive. General Discussion: Carbon exists in many forms and it is well known that the form or nature of the carbon used in reduction of oxides is related to the critical temperature of reduction. Sugar carbon, charcoal, and lampblack are forms of carbon that will reduce oxides at lower temperatures than high-temperature coke, and coke will, in turn, give a lower critical reduction temperature than graphite. There have been many investigations of this characteristic of carbons. Johnson' reported a difference of 130°F (70°C) in the critical reduction temperature of zinc oxide as between charcoal 1891 °F (1033°C) and Acheson graphite turnings 2021°F (1105°C) with zinc oxide. Bodenstein2 using charcoal and coke, found a difference of 138°F (77°C) comparing an experimental figure of 2066°F (1130°C) for coke and 1928°F (1053°C) for charcoal, in the reduction of zinc oxide. He concluded that this is very marked and observed that the "type of carbon merely raises or lowers the temperature at which rapid reaction takes place." Comparing the effectiveness of types of carbon in reduction of zinc oxide, it was found that a "brown coal coke" gave 97 pct zinc elimination at 1832°F (1000°C), as compared with 48 pct with "hard coal coke."' A wide range of metallic oxides was studied by Tammann and Sworykin,4 who found that the temperature at which decomposition of oxides begins depends on the nature of the carbon used. Carbon in the form of graphite, lampblack, and sugar carbon was investigated. Sugar charcoal will reduce Fe2O3 to Fe3O4 at 842°F (450°C) as compared with 1112°F (600°C) for coke, according to Meyer." Direct reduction of iron oxides by charcoal begins at 1382°F (750°C), but "first becomes intense" at 1652°F (900°C), whereas with coke, direct reduction begins at 1742°F (950°C), and "first becomes appreciable" at 2012°F (1100°C).6 he total reduction of the sample under certain conditions when heated in a current of CO with charcoal was about 100 pct for limonite and about 77 pct for magnetite. Using coke under the same conditions, the respective percentages were 75 and 47. In a study of processes for sponge iron7 by the Bureau of Mines, the conclusion was reached that a low-temperature char from noncoking subbituminous coal is the most satisfactory solid reducing agent. In a critical study of zinc smelting from a theoretical viewpoint Maier8 concluded that the reduction is by CO, that the reaction between ZnO and CO is intrinsically more rapid than the subsequent reduction of CO2 by C, which is limited by diffusion rates, which in part effectively limits the smelting process. Maier said that the operation is improved with the activity of the reducing carbon. An active carbon, he said, is one maintaining a low CO, content in the retort. Reactivity of Carbon: One form of carbon is more potent in reducing oxides than another. A carbon that reacts faster than another at a given temperature is said to be more reactive. Reactivity is measured by several methods, using carbon dioxide, air, or steam as reactants.9 ebastian and Mayers" have developed a method for the determination of absolute reaction rates between coke and oxygen by a study of ignition points under certain conditions. These and other investigators have established the relative reactivity of types of carbon. Lignite, charcoal, bituminous coal, cokes in the ascending order

Jan 1, 1951

By Fred C. Bond, Jen-Tung Wang

Comminution energy is principally energy of deformation before breakage, which appears as heat An empirical equation is presented which covers the entire comminution range. The new strain-energy theory considers comminution from the known principles of mechanics and the reduction ratio. Energy requirements according to the different theories the romnnred. THE present status of the theory of comminution is extremely unsatisfactory. The amount of various ores and rock products crushed and ground annually approaches the staggering sum of perhaps one half billion tons. Yet the basic concepts underlying an operation of such magnitude are unknown, and actual knowledge of what takes place in comminution is almost entirely empirical. The responsibility for this situation cannot be charged to apathy or lack of appreciation of the importance of the problem. Many attempts have been made to evolve a workable hypothesis, but none has been completely successful. General Considerations In order to break rock, it must be subjected to a stress which strains the rock beyond its critical breaking point. The stress imparts energy to the rock, most of which is released in the form of heat when the stress and resulting strain are removed. If the strain induced exceeds the critical strain, the energy is released by breaking; if the induced strain does not reach the critical point, the energy is released with removal of the stress as energy of resilience. Rock is a brittle material and it is assumed that it breaks approximately at its yield point, so that no permanent deformation results from strains below the critical. When rock is broken the total energy input is accounted for by the heat liberated, and by the surface energy of the new surface produced. A small and probably negligible amount is released as noise. Rock is commonly broken under compression and the applied stress is ordinarily compressive. However, for breakage to occur, it is only necessary that the induced strain exceed the critical value, and this strain may represent the resultant of compressive, shearing, and tensile forces. The stress-strain diagram of rocks under compression, and the modulus of elasticity as determined therefrom, should be fundamental considerations in the development of any theories of rock breakage. Crushing in jaw or gyratory crushers results primarily from a squeezing action, and reduction in a hammer mill is the result of impact. The action in rod and ball mills is a combination of impacting, squeezing, and wearing away by attrition, or rubbing. In ordinary crushing and grinding, the forces are applied at protruding points and are not distributed evenly throughout the rock. This is one of the features of comminution which theoretical considerations have not covered adequately. Another feature is the effect of impact velocities, which may greatly influence the total energy required for breakage. Impact velocities may vary from 1 or 2 fps in crushers to 10 or 20 fps for ball mills, and as high as 100 or more for hammer mills. However, according to the findings of geophysicists, the velocity of sound and of compression waves in stone is many times greater; the primary longitudinal compression waves travel at perhaps 15,000 fps, and the secondary transverse waves, which may reach the surface and cause cracks to form thereon, travel at perhaps 5000 fps. When a crack tip forms, the total surrounding stresses are concentrated in this tip, which rapidly extends throughout the rock particle.1 It would seem that the energy required to deform the rock beyond its critical strain, resulting in the formation of a crack tip, represents practically all of the energy required. Most of the primary crack tips presumably form on the surface. A simple method of designating feed and product size is necessary for the evaluation of crushing and grinding results. Taggart's suggestion2 that the 80 pet passing size be used has been found very practical. In this paper, the square opening screen size which 80 pct of the feed passes is designated as the feed size F, and the size which 80 pct of the product passes, or the product size, is designated as P. The reduction ratio at the 80 pet passing size, or F/P is designated as n. The feed and prod-

Jan 1, 1951

By B. H. Caudle, R. L. Slobod, E. R. Brownscombe

An earlier publication has discussed three methods lor obtaining relative permeability data on small core sample-: and the apparatus and technique for the capillary pressure displacement method. This paper describes the apparatus and technique for the solution gas and the dynamic displacement methods and presents a routine procedure for obtaining oil-gas and water-oil relative permeability data. Theoretical and experimental considerations are presented to show that the end effect commonly associated with the dynamic flow mechanism is extremely small where constant rates are employed in the flowing phase. An effect of flow rate on relative permeabilities obtained by the dynamic system is found only when gas is one of the flowing phases and this effect is ascribed to a form of channeling in the capillary system. The apparatus and procedure* used to obtain relative permeability data with three phases flowing are described and some preliminary results of the use of this method are shown. INTRODUCTION In a previous publication from this laboratory there appeared a brief discussion of the concepts behind three basic methods for obtaining relative permeability data.' These three methods were called the capillary pressure displacement method. the solution gas displacement method, and the dy. namic displacement method — the names being suggestive of the type of process used for obtaining the desired saturation prior to making the permeability measurements. In the same publication, the apparatus and technique for the capillary pressure displacement method were described and some typical results obtained by this method were presented. This paper will present: 1. The routine procedure used to obtain complete relative permeabilitv data (In small core samples with two phasec flowing; 2. The experimental techniques for the solution gas and dynamic displacement methods for obtaining relative permeability data; 3. The results of some studies on the mechanism of fluid flow through consolidated porous media; and, 4. The preliminary results on the determination of relative permeability for a system in which three phases are flowing. TECHNIQUE AND APPARATUS FOR THE DETERMINATION OF RELATIVE PERMEABILITY WITH TWO PHASES FLOWING In order to obtain relative permeability data on a small core sample which will describe the various phases of reservoir depletion. it is necessary that the methods used duplicate, as closely as possible, the various mechanisms involved in rescbrvoir depletion. It is important that. in any type of relative permedbility determination. the core be originally at virgin reservoir conditions with respect to the fluids within the pore *paces. (for example. oil used to displace water from core until connate water saturation estahlished). and that the subsequent saturation changes in the core be analagous to tile reservoir depletion mechanism for which the data are to he used. PROCEDURE OF DETERMINING RELATIVE PERMEABILITY DATA WITH TWO PHASES FLOWING A standardized procedure has been set up for obtaining relative permeability data and other information on small core samples. This procedure essentially consists of 11 steps as shown below: 1. Determination of tile capillary pressure curve.'! This information is used to help select samples for relative permeability determinations and to charactrrize. in general. the pore size distribution of the sample. 2. Determination of the specific permeabilities to air, oil. and water. These values are those on which the relative pern~eahility calculations are based. 3. Reduction of the core to connate water by an oil drive using the high speed centrifuge. 1. Determination of the critical gas saturation value and the associated relative permeability to oil by the solution gas drive method. 5. Determination of the relative permeabilities to oil and gab in the region where both phases are flowing by the capillary presqure method. 6. Reduction of the sample to residual oil in the centrifuge by an air drive. and the determination of the relative permeability to gas at this point.

Jan 1, 1951