Search Documents

By E. W. Hough, B. B. Wood, M. J. Rzasa

An apparatus for the determination of fluid-fluid interfacial tension by the pendant drop method has been constructed. The apparatus is refined beyond those previously described in that the samples are introduced into an evacuated windowed cell without contacting air or mercury. Interfacial tension data are reported for the water-methane system between 74 and 280°F and 15 and 15,000 psi. The variation of wettability of the stainless-steel dropper tip with pressure and temperature was observed incidentally to the measurements. INTRODUCTION The most refined previous apparatus --for the determination of fluid-fluid interfacial tension at reservoir conditions was built by Hauser and Michaels.' These authors used the pendant drop method, previously investigated by Andreas, Hauser, and Tucker.' Recently, Fordham3 has made the pendant drop method an absolute one by numerical integration of the differential equation of drop shape. It is necessary to know the difference in density of the two fluids and the acceleration due to gravity in order to evaluate the interfacial tension from measurements of drop dimensions. Contact angle is not involved in the measurements by this method. Some precise tests of the method have been made at atmospheric pressure and 77°F by Douglas," who found agreement with established values to within about 0.1 per cent for water and air, cyclohexane and air, benzene and air, and benzene and water. No values of interfacial tension from the recent pendant drop apparatusJ have been published by Hauser and Michaels. However, values for the water-benzene and water-decane systems are given in Michaels' thesis5 or Dressures from 15 .to 10,000 psi and temperatures from 73 to 268ºF. Some determinations on gas-oil and gas-water systems in the range of 15 to 4,000 psi and 78 to 178°F have been made by Hozott,' by the drop-volume method. The capillary rise method has been used to investigate the interfacial tension in gas-oil systems at 88°F to several thousand psi by Swartz.7 Other determinations in gas-oil systems at pressures under 1,000 psia have been made by Jones: and Beecher and Parkhurst.9 The scarcity of data for pressures above 5.000 psi is apparent. APPARATUS An assembly drawing of the working section of the interfacial tenqion apparatus is given in Fig. 2. 'The endanant drop (A) is housed in a pressure chamber having glass windows (B). An agitator (6) is provided. The metal parts of the apparatus coming into contact with the sample are made of AISI: Type 316 steel. Drain grooves are provided around the whdows so that the sample in contact with the O-rings is not returned to the working section when the pressure is reduced. A schematic arrangement of the apparatus is given in Fig. 3. The working section containing the pendant drop (A) is filled by samples from bombs (B) with the aid of the vacuum system (C). The O-ring sealed pistons (D) used to separate the sample from the hydraulic fluid are shown together with the pressure production and measurement system (E). A front view of the apparatus is shown in Fig. 1. Most of the hydraulic system is at the top of the panel. The pump and displacement piston are at the lower right. The Variacs for the air bath heaters. the camera. and the vacuum gauge are a1so visible on the panel. The working section together with a window closure and a drop-forming tip closure are shown in Fig. 9. The carbon arc light source with filter is visible in the upper part of the picture. Two sides of the air bath with the impeller and a part of the convection shield can also be seen. The assembled air bath is shown in Fig. 10.

Jan 1, 1951

By H. C. H. Darley

Particles of montmorillonite-type clays consist of plate-like clay crystals stacked face to face. Under certain conditions in a suspension, these particles disperse to individual platelets with the formation of diffuse double layers round each platelet. This action has a profound influence on the rheo-logical and other properties of the suspension. This paper describes a test for degree of dispersion and uses this test to correlate changes in dispersion with changes in viscosity and gel strength. The test is of practical value in determining and controlling dispersion in drilling muds. INTRODUCTION It has long been recognized that the hydration and dispersion of clays, particularly those of the montmoril-loilite type, strongly influences the rheological and other properties of drilling muds, but little has been done to measure either of these phenomena under conditions that exist in drilling muds. The hydration of clays as such has been extensively studied by those interested in clay mineralogy, but in drilling muds it is to be expected that dispersion is at least as important a factor as hydration because of the much higher water-to-clay ratios and the very high degree of mechanical agitation. The term dispersion as generally applied to colloidal systems means the subdivision of large particles to small. Particles of montmorillonitic clays, however, consist of plate-like particles stacked face to face. Dispersion of these aggregates to individual platelets may be expected to have an unusually large influence on the properties of the suspension because of the shape factor, the enormous increase in specific surface, and the pronounced surface phenomena exhibited by the crystal surfaces. It is this form of dispersion that is particularly considered in this paper, -and the term degree of dispersion is used to imply the number of platelets per particle. The term aggregation is used to denote the opposite of dispersion, the association of individual platelets face to face to form stacks. It is important, as Garrison' has pointed out, to distinguish between the deflocculation-flocculation phenomenon and the dispersion-aggregation phenomenon. Flocculation of clay suspensions is characterized by the formation of voluminous sediments that are stable under quite high centrifugal fields. This indicates an edge to plate association of particles to form an open network type of structure, which obviously will have quite a different influence on the properties of the suspension than aggregation of the platelets face to face. The work described in this paper was primarily concerned with developing a test for degree of dispersion in clay suspensions and by means of this test to follow changes of dispersion brought about by addition of electrolytes and to show the influence of dispersion on the rheological properties. The Clay-Volume Test This test is based on the assumption that the greater the number of particles per unit weight of clay in a suspension the greater will be the volume of sediment when it is flocculated under standard conditions and centrifuged. After considerable experimentation the following procedure was established: 1 cc of the test mud vigorously shaken with 100 cc of 1-N NaCl in an ASTM oil centrifuge tube and centrifuged for five minutes in an International No. 2 centrifuge at 2,000 rpm. When drilling fluids of different solids content are being compared, the clay volume is expressed in cc per gram of clay solids, but generally changes in dispersion in a given mud are being compared and it is more convenient to report as cc of sediment. The test is strictly a comparative test, so that it is not necessary to adhere to the conditions specified above, e.g., with field muds a stronger flocculant such as Al2(SO4)3 may be preferable to NaCl.

Jan 1, 1958

By U. N. Bhrany, N. Arbiter

F. C. Bond (Consulting Engineer, Processing Machin-Dept., AllisChalmers Manufacturing Co., Milwaukee) This is a very comprehensive paper. It deals with 1) size distribution functions, 2) energy-particle size relationships, and 3) criteria for grinding efficiency. The method of analysis used consists of tracing the increase and decrease of successive Tyler scale screen fractions as grinding progresses. The treatment is based upon the assumption that the work done in grinding is concentrated entirely upon the coarse size fractions in the mill. It is postulated that when a feed consisting entirely of particles of one size fraction is batch ground that finer fractions increase at a constant rate as grinding progresses, until finally each in turn begins to be ground and to decrease in amount. It is shown that under these conditions the empirical Gates-Gaudin-Schuhmann log-log size distribution plot can be logically explained. It is difficult to imagine any condition of ball mill operation, or any probability analysis, which would result in no grinding contacts on the finer particles present. However, the authors' postulate must be accepted if their conclusions are to be sustained. They conclude that the equation expressing the for mation rate of each size fraction can be resolved into the G.-G.-S. size distribution equation, and that the curvilinear portion of the log-logsize distribution line represents the only particle sizes which are being ground as well as produced. The author's statement that the G.-G.-S. Eq. 4 applies only to the finer portion of the size distribution while the Rosin-Rammler Eq. 6 applies over the whole size distribution range certainly supports the view that an exponential equation of the Rosin-Rammler type corresponds more closely to actual size distribution than does the power law G.-G.-S. equation. However, the G.-G.-S. equation is used to derive Eq. 8 and Eq. 9, which state that the energy required to grind a ton through a particular mesh varies inversely as the particle size to the exponent -a , whic h is the slope of the extended G.-G.-S straight line. This is stated to apply to both batch and continuous grinding. The G.-G.-S. slope exponent thus becomes very important in all crushing and grinding analyses by this method. The authors include some data which can be used to test their conclusions, as described below: We will consider grinding a feed composed of 14 mesh (1190u) to -200 mesh (74u). The average of 17 Allis-Chalmers grindability tests on quartz gives a work index (Wi) of 12.77. The -200 mesh product should have 80 pct passing 41.6u or P = 41.6. The equivalent 80 pct passing size of the 14 mesh feed particles is considered to be the particle diameter, as is true for spherical particles. Thus F is 1190u. In the work index equation 10 wi 10 wi 127.7 127.7 P f 41.6 1190 = 16.1 kw-hrperton. In Table I the ball mill exponent for quartz is 1.0 and from Eqs. 9 and 10: w=kw-hr/ton =Constan/(do/di) = 16.1 ton so that (1190) Constant/ Constant and the Constant equals 259. The tons per horsepower-hour value is 0.0464, in agreement with Table IV. The exponents for other minerals are listed in Tabl e 11. In Table IIA, the kilowatt hours per ton required to grind these minerals between the same sizes as the quartz, both by the authors' method and the work index equation, are listed. According to the authors kw-hrperton = 259/(16.1)".

Jan 1, 1961

By J. N. Ong, J. P. Ratledge, J. H. Boyce

SINCE German operators opened the Tsumeb mine in the early 1900's, continuous operation has been interrupted only by enforced shutdowns during two world wars and the depression of the 1930's. Original metallurgical practice was blast furnace smelting. This was applicable only to low zinc ore, leaving large tonnages of high zinc ore to be stockpiled on the surface. A gravity concentrator using tables, buddles, and jigs operated for several years, but differential flotation of the ore did not get beyond the experimental stage. When the property was put up for sale by the South African Custodian of Enemy Property there was an estimated 434,000 tons of ore on surface dumps. The dump material, assaying 4.48 pct Cu, 15.8 pct Pb, and 10.9 pct Zn, had become highly oxidized during a stay of many years on the surface. But these dump ores had to be used in evaluating possible metallurgy because the shaft became flooded after the 1940 shutdown and underground ore was unavailable. Many companies examined Tsumeb but quickly lost enthusiasm because of the formidable metallurgical problems presented by the dump ore. However, on the basis of extensive test work, the present company decided that differential flotation would be successful and purchased the property after competitive bidding. Geology and Mineral Association is Complex The Tsumeb mine lies near the northern boundary of South West Africa at an altitude of 4200 ft. The village of Tsumeb is 335 air miles northeast of Wal-vis Bay, territorial port of entry. The orebody is a hydrothermal replacement deposit associated with an irregular intrusion of aplite cutting folded pre-Cambrian dolomite beds. In plan the deposit is a lens of 250x600 ft maximum dimensions. In section the orebody dips 55" south from the surface to 2000 ft, below which there is a reversal in dip to 75" north. Mineralization is localized in the aplite and in the adjoining brecciated dolomite. In the upper levels the ore is largely confined to a single lens, but below the 1890 level massive sulphide ore forms north and south lenses that partially enclose a pipelike body of medium to low grade disseminated ore, part of which is of too poor a grade to mine. Drilling below the 24th, or bottom level, has disclosed that a high percentage of the mineralization from the 3000 to 3300-ft horizons is comprised of copper and lead oxides. The deep oxidation may be local, since it appears due to circulating waters in a major fault cutting the orebody. Tsumeb ranks as one of the highest grade base metal deposits in the world. It is geologically and mineralogically famous for 41 copper, lead, zinc, vanadium, and germanium ore minerals, but its complex mineral assemblage is not attractive to metallurgists. In addition, there is a gradual change in oxidation intensity from the upper levels to the lowest level now operating, which is the 24th. This is well illustrated by the copper minerals. In the upper levels malachite and azurite are predominant, in the intermediate levels chalcocite, and in the lower levels tennantite. Throughout the mine there are sufficient heavy metal soluble salts to activate zinc, and large amounts of depressants are required even for so-called clean sulphide ores. When the same sulphide ore is exposed to prolonged oxidation aboveground, as was the case with the dump ores, reagent requirements become fantastic. Plant Construction First staff members arrived at Tsumeb in March 1947, and by May essential services had been restored and mine unwatering began. In 2 1/2 months company employees built a temporary jig and sorting plant for treating existing dumps to provide a saleable product during high metal prices. Meanwhile, concentrator design and construction was contracted to a South African firm, with subcon-

Jan 1, 1956

By Douglas R. Piteau

Douglas R. Piteau (Engineering Geologist, De Beers Corp., Kimberley, South Africa; Presently Consulting Specialist at Nchanga Mine Open Pit, Chingola, Zambia)-Mr. Lutton is to be commeneded on his analyses and presentation of such practical data which can be used for predicting a preliminary optimum slope angle in rock. It appears that he assumes all slopes analyzed to be straight. Although there are overriding factors which affect the general slope configuration, such as structural discontinuities, ground water, and other environmental factors such as weathering conditions, regional stress, etc.,8 the affects of plan geometry under certain conditions also can be of outstanding importance. In investigation of natural slopes, for example, the effects of plan geometry are not entirely appreciated and it is in this respect that discussion is given. Experience in southern Africa has shown that the overall geometry of surface openings in rock, particularly those with dimensions of the order of 1500 ft or less, i.e., where radii of curvature of the slopes are less than, say, 700± ft, must be considered in rational slope-stability analyses. That is, the slope angle is not only a function of slope height, strength of the mass, and geometry of structural discontinuities, but also is dependent upon the plan geometry of the local slopes and the plan shape of opening. Technically speaking, adequate analysis of the slope must be three-dimensional, the third dimension being that of overall plan configuration and the analysis must incorporate its subsequent effects. However, current theories of slope stability all deal with the slope as a two-dimensional problem, i.e., a slice of unit length of an infinitely long slope is considered for a condition of plane strain. In this the tacit assumption is made that the plan radii of the crest and toe of the slope are infinite and no attention is given to the plan geometry of the slope. Clearly this is not the condition encountered in practice, particularly in open pit mining where excavations, relatively speaking, can be small and where these radii of curvature accordingly can have important affects on the slope stability. This is illustrated by the analyses of Piteau and Jennings 9 of the natural slopes which had developed around pipes of four diamond mines in the semiarid Kimberley area of South Africa. Although these slopes consist of dolerite overlying shale, the values of the slope angles of all the pits were corrected for a height of slope of 320 ft assuming that no dolerite occurred. This brought the data onto a comparable basis. These data were analyzed statistically and a clear empirical relationship, as is shown [in Fig. 4], is seen to have formed between angle of slope and curvature of slope, i.e., the angles of slope are greatest where the radii of curvature are least. Also, it can be seen that the slope of each line is different, indicating that the plan shape of the open pits also are significant, the curve being steepest for the open pits which are more circular in shape in contrast to those which are more elongated. The effect of slope curvature on the average angle of slope for the two extreme radii of 200 ft and 1000 ft for all four open pits is 39.5° and 27.3°, respectively. Hence, the concave curvature effect represents an average increase in slope angle of 12.2° from the essentially straight slope. In quantitative terms in Mr. Lutton's Fig. 2, this represents displacing the curve for a 320-ft pit from the second curve from the bottom (i.e., shale in subtropical climate) to that of the third from the bottom (i.e., quartzite in temperate environment).

Jan 1, 1972

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 Richards R. Vandervoot

T HE creep behavior of many polycrystalline metals at moderate stresses can be described by the empirical relationship: where Em is the minimum creep rate, A is a constant, a is the applied stress, E is the average elastic modulus at the test temperature, and D is the self-diffusion coefficient.&apos; This equation is generally applicable for strain rates of less than around 10-5 per sec and temperatures greater than about one-half the melting point in absolute units. Under these conditions atomic mobility plays an important role in the deformation process, and steady-state creep occurs when a dynamic equilibrium is reached between the generation and recovery of dislocation substructure. Although several theoretical mechanisms have been proposed to explain the creep of polycrystalline solids,2"4 the dislocation climb model2 is in accord with most experimental observations. The majority of creep data used to develop Eq. [1] came from metals with either the fcc or hcp crystal structures. However, on the basis of results for a and 7 iron as well as those for a and a thallium, Sherby concluded that the high-temperature deformation behavior of close-packed metals and bcc metals is probably the same.&apos; The greater creep resistance of fcc and hcp metals and alloys on a relative temperature scale was attributed to lower diffusivities in the close-packed structures. To further compare the creep properties of fcc, hcp, and bcc metals, creep data for all six bcc refractory metals were analyzed according to Eq. [I]. The results are shown in Fig. 1. The sources of information for the quantities used to prepare this figure are listed in Table I, except for high-temperature elastic moduli. The latter data were obtained from Refs. 28 and 29. The dashed lines in the figure represent the range and breadth of data analyzed by Sherby for fcc and hcp metals. In general, the deformation behavior of all of the refractory metals, ex- cept possibly that of tungsten, can be represented by Eq. [1]. This correlation is strong evidence of the similarity that exists between creep of pure close-packed metals and creep of pure bcc metals. At (Em/D) < 108, the results of Flagella and Tarr on tungsten (W-1 in Fig. 1) agree with the five-power stress law, while those of Green (W-5 in Fig. 1) show a seven-power stress dependence. Where (Em/D) > 108, the stress exponents for other data on tungsten were in the range 10 to 15. It is possible that the creep rate dependence on stress in this region is not a power law but exponential (i.e., Em - eßo). However, it was not feasible to differentiate between these two cases because of the steep slope and scatter in the data. Note the inflection point in the curve for the data of Flagella and Tarr (W-1 in Fig. 1). The slope of the line changes from 5 to around 10 at Em/D = 108. The breakdown of the fifty-power stress law at high stresses is not uncommon in other metals, and this phenomenon may result from nonequilibrium conditions in the material such as excess lattice vacancies produced by rapid straining.30 The mechanical strength of a polycrystalline metal is related to its microstructure. Attempts to establish the influence of grain size on the creep rate of metals and alloys have led to a wide variety of results. For example, the following effects have been observed: creep rate decreases as grain size increases;31&apos;32 creep rate decreases as grain size decreases;&apos;,33 creep rate is a minimum at some critical grain size and increases with either an increase or decrease in

Jan 1, 1969

By Leonard B. Lipson

Boyer, Morgan and Muskat,&apos; with improvements by Morgan, McDowell, and Doty,l and Laird and Putnam,&apos; 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/&apos;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 R. P. Alger, C. A. Doh, M. P. Tixier

The principle, the equipment and field operation of sonic logging are described. The tfio-receiver system produces logs independent of hole size and mud. Field experience is given and forms the basis for the interpretation of the log. The derivaiion of porosity values from measured velocities is discussed, according to the type of formations (hard formations, compacted sands, and unconsolidated formations). The time-average equation proposed by M. R. J. Wyllie, A. R. Gregory and L. W. Gardner is uscd as a basis for the computation of porosity in limestones, cemented sandstones and compacted sands. Variations in the lithologic character of limestones do not seem to change the porosity calibration markedly. The compaction of sands is related to the compaction of shale adjacent to them; and, thus, the shale velocity is used for the establishtnent of empirical relations for the computation of porosity in un-consolidated formations. Various forrnu.las are tentatively presented to account for shale and fluid content. Field experience demonstrates that considerable attenuation of sonic energy takes place in unconsolidated formations, particularly when gas bearing, and in fractrured formations. Unusually large attenuation produces skipped cycles, a feature easily recognized. The application of sonic logging to structura1 studies is featured, showing its possible integration with the dipmeter. Correlation and its application to the irrterpretation of seisrmic surveys are reviewed. The paper is illustrated with field examples. Sonic logging is the recording of the time required for a sound wave to traverse a definite length of formation. Sonic travel-times are inversely proportional to the speed of sound in the various formations. The speed of sound in subsurface formations depends upon the elastic properties of the rock matrix, the porosity of the formations and their fluid content and pressure. Below the "weathered" or low-velocity layer extending from 50 to 100 ft or so below the surface, sound velocities may range from about 6,000 ft/sec in shallow shales to as much as 34,000 ft/sec in dolomites. In hard formations (well cemented and/or compacted), the sonic log reflects the amount of fluid in the formations; hence, it correlates well with their porosity. In unconsolidated formations, which are usually of fairly high porosity, the sonic log gives an approach to porosity determination, when its readings are corrected for lack of compaction, shaliness, and fluid content. In such formations the sonic log may also indicate the presence of gas and may distinguish between oil- and water-bearing beds. Field experience with the sonic log in soft formations is more limited than that in hard formations. As a result, interpretation in soft formations is not as well developed at the present time. In all types of formations the sonic log is of considerable value for correlation and for the more detailed and accurate interpretation of seismic data. EQUIPMENT AND PRINCIPLE FIELD OPERATION The basic features of the sonic log" are a two-receiver system and two available short spacings. The log is recorded on film using a standard field unit. An SP curve is run along with the sonic curve to give a more interpretable log and to establish absclute depth control for comparison with other logs run in the same well. If desired, a gamma-ray curve may be recorded simultaneously with the sonic curve. Fig. 1 shows the sonic-gamma ray combination sonde. The sonic portion consists of a sound transmitter and three receivers mounted 3, 4 and 6 ft from the transmitter. When the first and second receivers are used, a sonic log of I -ft spacing (distance between receivers) is obtainea. Use of the first and third receivers provides a 3-ft spacing. The transmitter emits pulses at the rate of l0/sec. After the transmission of a pulse the first arrival of sound energy at each receiver triggers its response system. The differ-

By H. I. Kaplan, C. Laird

Pulsating contpresslon experiments have been carried out on coppev single crystals in order to test the adequacy of mechanisms which have been suggested for stage I cvack grouth when tension-compression loading is the mode of- applying stress. Although dejorrrzation continues slowly for many millions 01 cycles under Pulsating compression stresses both within slip bands formed in the first cycle and by the creation of new bands, the morphology of these bands is not typical of fatigue slip bands in general and neither initiation of cracks nor propagation of artificially induced cracks was observed. It is therefore concluded that 1) initiation and growth of cvacks in fatigue is a consequence of repeatedly reviersed, microplastic deformation rather than intermittent, unidirectional deforma-tion, 2) gvowth of- cracks in polycrystalline specimens FAILURE of a plain specimen subjected to cyclic tension-compression stresses is known to take place by two stages of cracking.1,2 The first stage is characterized by propagation along slip planes roughly at 45 deg to the stress axis, in the outermost layers of the specimen, and at a rate of the order of angstroms per cycle. The average macro orientation of a stage I crack is normally at right angles to the stress axis,3 however, except in a few cases, e.g., single crystals in suitable orientation, when both the micro and macro directions of propagation are at 45 deg to the axis. After a stage I crack has penetrated to such a depth that the stress on the remaining section of metal is significantly increased, stage II crack propagation succeeds stage I, at right angles to the stress axis both micro and macroscopically. Direct evidence of the processes occurring during a subjected to pulsating compression stresses is caused by the setting-up of local tensile stresses due to interactions between grains, and 3) both unslipping mechanisms and operation of the "plastic blunting process " within a single slip band are adequate descriptions 01- stage I growth, being fundamentally similar. The slip mode of a material, however, is considered to control whether unslipping or the plastic blunting process takes place. Thus materials of planar slip mode fauor unslipping while those hacing a wavy mode, and capable of cross slip, can fail in stage I by the plastic blunting process. Similarly, and perhaps more importantly, reversed torsion and tension-compression modes of applying stress may be expected to cause the unslipping mechanism and the plastic blunting process, respectively. single stress cycle4,5 shows that stage II crack propagation takes place by a repetitive plastic rounding and closing of the crack tip. This mechanism has come to be called the "plastic blunting process".6 On the other hand, the extreme smallness of the phenomena occurring in stage I growth prevents elucidation of its mechanism by direct observation. Many suggestions have been proposed for this mechanism,&apos;"-" most of them involving some kind of random or systematic unslipping process. It is possible, however, that the mechanism of stage I growth is essentially the same as that of stage 11, but occurring on a much smaller scale.6 In this situation, it has been observed that cracks propagate along slip bands comprising dislocation cell structures and bounded on either side by regions of different dislocation density and distribution." Since there is evidence that such bands are softer than their matrix, * it may be expected that the plastic blunting process occurs within a single slip band during stage I growth, rather than on two gross bands as in stage 11, when the acting stress is greater.= The blunting of the crack tip allowed by coordinated

Jan 1, 1968

By J. L. Huitt, S. R. Darin

The use of a partial monolayer of propping agent to obtain a high flow capacity for a hydraulically induced fracture is discussed. From the results of laboratory work it was shown that a modified form of the Kozeny-Carman relation could be used to describe the flow in the partial monolayer propped fracture. With equations presented in the paper, the density pattern of the propping agent (number of particles per unit of fracture surface) that results in the maximum flow capacity for the fracture can be determined. The maxitnum flow capacity obtained with a partial monolayer is often an order of magnitude greater than the flow capacity obtained in greater width fractures containing multilayers of the propping agent. INTRODUCTION One of the predominant factors controlling the success of a hydraulic fracturing operation is the propping of the fracture. The trend in hydraulic fracturing recently has been an increase in the ratio of propping agent-to-fluid. The primary purpose of this increased ratio is to sustain a propped fracture of greater width. It was shown in a recent study&apos; that for some formations low concentrations of propping agent would result in "closing" or "healing" of the fracture. In some cases, in an effort to insure sufficient propping agent concentration, the fracturing operation is designed to obtain a pack of the propping agent in the fracture. The placing of a pack of propping agent in the fracture provides a fracture of maximum width; however, the width alone does not control the flow capacity of the fracture. The flow capacity is dependent on the permeability of the fracture, as well as the fracture width. Thus, increasing the permeability to obtain larger flow capacities is as important as obtaining a greater width fracture. By increasing the concentration and particle diameter of the propping agent, fractures of sufficient flow capacities are obtainable for effective well stimulation in a number of formations. However, in some of the recovery techniques, the fracture capacity obtained by commonly used practices in fracturing is inadequate. One such recovery technique is gravity drainage via a horizontal fracture placed near the base of the productive zone. For this recovery technique, the fracture capacity desired is usually one to two orders of magnitude greater than the fracture capacity normally obtained in fractures propped with the size of sand commonly used (20-40 mesh). By the use of particles of larger diameter, the permeability of a pack2 is increased; however, even a multilayer pack results in a fracture flow capacity lower than that desired for some formations. For example, a fracture of 1/4-in. width, packed with 8-12 mesh sand of about 2,000 darcies permeability, would result in a fracture flow capacity (permeability X fracture width) of only 40,000 md-ft. If the fracture flow capacity needed is an order of magnitude greater than that cited in the example, an approach different from the formation of a pack in the fracture is needed. One approach is a partial monolayer of large-diameter propping agent. In working with this type of fracture system, it was found that the flow problems were unique to the system; therefore, the study discussed in this paper was made. THEORY AND RELATED STUDIES Previous studies, although related to that described in this paper, have not pertained directly to the flow in a fracture propped with a partial monolayer of propping agent. However, studies such as flow between parallel plates,&apos; flow in packed columns,&apos; fracture flow capacity&apos; and flow in simulated fractures reduce the scope of this flow problem considerably. A fracture propped with a partial monolayer of dropping agent may be considered to be between two extremes — between an open fracture and a bed of particles packed In a conduit of fracture geometry. Thus, it seems appropriate to review the extreme cases.

By O. D. Gonzalez, R. A. Oriani

The migration of hydrogen and of deuterium dissolved in a iron and in nickel induced by an applied electrical potential has been measured over a range of temperature. In all cases the intevstitial solute migrates to the cathode. The deduced charges of transport are positive, independent of or slightly dependent on temperature, and markedly larger for the beavier mass. These results are interpreted as showing predominating momenlum transfer between electron holes and the acticated complex of the elementary diffusion act. THE migration of components of a metallic system under the action of an applied electrical potential is called electromigration or electratransport. This phenomenon has been known since 1861 and is becoming of increasing interest to metal scientists because of the opportunity it offers to explore the connections between mass transport and electron transport. Until recently, reliable data have been very scarce because of the difficulties inherent in the measurements. In addition, theoretical understanding has been either nonexistent or rudimentary and at present this situation is only moderately improved. Whereas the prevalent notion had been (and is still widely held) that electromigration directly reflects the state of ionization of a component in its equilibrium state in the metal, it is now becoming clear that the major, if not the sole, factor is the momentum transfer between the charge carriers, themselves in motion because of the applied field, and the metallic components in their activated state. It becomes important therefore to obtain reliable data over a range of temperature on electrornigration in systems where only one component moves, that is, either in pure metals (vacancy motion) or in interstitial solid solutions. Our choice of systems was dictated by our measurements&apos; of the Soret effect in solid solutions of hydrogen and deuterium in iron and nickel. Indeed, the main reason for making measurements on electromigration at all is that we suspected a mechanistic connection between thermal diffusion and electromigration. In our former paper&apos; we proposed that what is common to both phenomena is the momentum transfer between the charge carriers in the metal and the atom as it jumps from one equilibrium position to another. In Soret diffusion, however, there is in addition the dissipation of the energy of the jumping atom via pho-nons. We have since found that similar ideas for this interconnection have been developed by Fiks.2 In the present paper we will discuss in more detail our understanding of the mechanism of electromigration itself, but we turn first to a presentation of the experiments and results. EXPERIMENTS AND RESULTS Let us consider a metallic system in which only one component, m , can diffuse at the temperature in question, and to which are applied simultaneously an electrical potential gradient, grad cp, and a chemical potential gradient, grad p,. One may write3 for the material flux Jm and the electronic charge flux Je by applying the condition, grad = 0, to Eq. [I] and evaluating grad(um/T) for the isothermal case and for activity coefficient independent of composition over the relevant composition range. ere, c is the concentration and D the diffusivity of the mobile component. The experimental technique that we have used to measure Z* for the isotopes of hydrogen individually dissolved in each of two transition metals has been developed from a method used by Wagner and Heller.4 (A somewhat similar method was also used by Herold.5) A wire or thin rod of the metal under study connects two volumes of gaseous hydrogen at identical pressures; hence grad um = 0. A direct current is made to pass through the specimen, serving to maintain the specimen at the desired temperature and also to produce the electromigration of the hydrogen dissolved in the metal in equilibrium with the gaseous hydrogen. Under these conditions, the flux of hydrogen from one volume to the other is given, from Eq. [I] for grad( um/T) =0by

Jan 1, 1968

By C. G. Dunn

This note describes a method for revealing micro-structures in tantalum through an {hkl} effect which was found during a study involving amorphous and crystalline tantalum oxides. Three steps are involved. First, an amorphous film of Ta2O5 is formed anodically1 on the surfaces of the tantalum. The amount needed depents on the thickness of the substrate tantalum; 1600Å of oxide, however, is generally satisfactory in the 1/2 -mil (12 µ) to 3-mil thickness range. The film at this point is uniform in color. Second, the anodized specimen is heated to 800°C in. vacuo during a period of a few minutes and held at temperature for 5 min. During such a treatment some oxygen leaves the oxide through dissolution into the substrate tantalum (e.g., at 10-8 Torr oxide films thinner than about 800 completely disappear on 5-mil strip). Considerable crystallization of the amorphous oxide also occurs adjacent to the metal/oxide interface but not throughout the oxide thickness. The thickness of oxide needed as a residue is not too critical; values between 500 and 1000Å (optical estimate&apos;) may be most useful. (The presence of a residual oxide layer after the heat treatment usually is not apparent; in fact, in most instances the surfaces have a metallic appearance. This can be the case also for thicker residual films needed in the detection of a crystalline oxide phase by surface X-ray diffraction—e.g., if the initial oxide film is 3000Å thick, enough crystalline oxide is formed for X-ray detection.) The third step eliminates the metallic appearance at the surface and reveals (or produces) a varied colored film of tantalum oxide. Two methods can be used and both will be illustrated. If the specimen is dipped in 48 pct HF for a few seconds, the tantalum is found to be covered with a multicolored oxide. Moreover the colors correlate with the grain structure of the substrate tantalum through an orientation dependence of the type illustrated in Fig. 1, where results of an orientation analysis by the Laue X-ray method have been plotted. Other film conditions would have provided a different set of colors for the same orientations. Any such correlation, however, provides a basis for a qualitative assessment of the distribution of orientations in a large sample. The cold-rolled structure of tantalum foil produced by a recrystallization anneal and a cold-rolling reduction of 95 pct is not easily obtained by usual metal-lographic methods.&apos; The structure, however, is easily revealed by the present technique, Fig. 2. The alternative third step is to reanodize the specimen. The optical nature of the final anodic film depends on the presence of the crystalline oxide phase; color variations indicate optical thickness variations in the film. As the reanodization voltage is increased from zero, straw and light brown colors appear (as usual), followed now by variations in color from grain to grain depending on {hkl} orientation. Suitable end points, which are determined by voltage and time, are easily established. Fig. 3, for example, shows the grain structure developed in 5-mil tantalum strip during a 30-min re-crystallization and grap-growth anneal at 1800°C. An anodic film of 1600Å thickness was added prior to the special 800°C anneal. Reanodization was done at constant current density (near 0.2 ma per sq cm) in an aqueous solution of H3PO4 (0.01 pct) and was terminated after the initial brown colors were replaced by a deep blue at approximately 30 v. Although good grain contrast can be obtained with bright field illumination, even greater contrast for black and white micrographs may be obtained with polarized light, as in Fig. 3.

Jan 1, 1967

By C. Wagner

When an alloy solidifies and the composition of the solid differs from the composition of the liquid, atoms of the alloying elements rejected at the solid-liquid interface have to diffuse toward the bulk liquid. Diffusion may be supported by convection. Theoretical calculations have been made for a liquid without convection, for a liquid involving natural convection, and for solidification of a liquid alloy at the surface of a rotating disk as an example of forced convection. WHEN an alloy solidifies, the equilibrium composition of the originating solid phase differs in general from the composition of the liquid phase. Thus, segregation may take place, i.e., the composition of the solid phase formed at different times differs from the bulk composition. The degree of segregation depends on several factors such as the partition ratio of the alloying elements between the liquid and the solid phase, the rate of solidification, and convection. To clarify, a theoretical analysis has been worked out. The following presuppositions are made: 1—Calculations are confined to the solidification of binary alloys in which the concentration of component A (solvent) is much greater than that of component B (solute). 2—The equilibrium concentration of the alloying element B is assumed to be lower in the solid phase than in the coexisting liquid phase. 3—At the solid-liquid interface, partition equilibrium is assumed. 4—The dependence of the partition ratio of the solute between the liquid and the solid phase on concentration and temperature is disregarded since only small concentrations of the solute are considered. 5—The concentration of the alloying element B in the solid phase is supposed to be lower than the saturation concentration for the formation of a second solid phase. 6—diffusion in the solid phase is disregarded, since diffusion rates in the solid state are comparatively low. When the partition ratio of the solute between the liquid and the solid phase is greater than unity, a part of the solute is rejected at the solid-liquid inter- face and has to diffuse toward the bulk liquid. Thus there is necessarily a concentration difference in the liquid between the solid-liquid interface and the bulk liquid as has been pointed out by Hayes and Chipman,&apos; Rutter and Chalmers,&apos; Tiller, Jackson, Rutter, and Chalmers, and others (see Fig. 1). In view of the concentration gradient in the liquid phase, the ratio of the concentration of the solute in the solid phase to that in the bulk liquid differs from the equilibrium ratio, although at the solid-liquid interface, partition equilibrium is assumed. Under certain limiting conditions, the solid may even have virtually the composition of the bulk liquid as is shown below. General Equations Let u be the linear crystallization velocity in cm per sec, D the diffusion coefficient of component B in the liquid phase in cm2 sec-&apos;, c&apos;,, and c", the concentrations of component B in the liquid and the solid phase, respectively, in mol per cm", and s the distance from the solid-liquid interface toward the bulk liquid at a given time t. The partition ratio is assumed to be equal to a partition constant K. Thus. The difference in the specific volume between the liquid and the solid phase may be disregarded. In view of the discontinuity of the concentration of the solute at the interface according to Eq. 1, u(c&apos;, — c",), , , = tic&apos;, (S = 0)(1 — K) gram-atoms B per unit area per unit time are rejected at the interface and have to diffuse toward the bulk liquid. From Fick&apos;s first law it follows that uc&apos;,,(s = 0) (1 ~K) == - D(dc&apos;,,/?>s). _-„. [2] The concentration profile in the liquid near the interface depends on the interplay of diffusion and convection as is discussed in the following sections. A schematic concentration profile is shown in Fig. 1. The "effective thickness of the diffusion boundary

Jan 1, 1955

By P. J. Root, J. E. Warren

An idealized model has been developed for the purpose of studying the characteristic behavior of a permeable medium which contains regions which contribute significantly to the pore volume of the system but contribute negligibly to the flow capacity; e.g., a naturally fractured or vugular reservoir. Unsteady-state flow in this model reservoir has been investigated analytically. The pressure build-up performance has been examined in some detail; and, a technique for analyzing the build-up data to evaluate the desired parameters has been suggested. The use of this approach in the interpretation of field data has been discussed. As a result of this study, the following general conclusions can be drawn: 1. Two parameters are sufficient to characterize the deviation of the behavior of a medium with "double porositym from that of a homogeneously porous medium. 2. These parameters can be evaluated by the proper analysis of pressure build-up data obtained from adequately designed tests. 3. Since the build-up curve associated with this type of porous system is similar to that obtained from a stratified reservoir, an unambiguous interpretation is not possible without additional information. 4. Differencing methods which utilize pressure data from the final stages of a build-up test should be used with extreme caution. INTRODUCTION In order to plan a sound exploitation program or a successful secondary-recovery project, sufficient reliable information concerning the nature of the reservoir-fluid system must be available. Since it is evident that an adequate description of the reservoir rock is necessary if this condition is to be fulfilled, the present investigation was undertaken for the purpose of improving the fluid-flow characterization, based on normally available data, of a particular porous medium. DISCUSSION OF THE PROBLEM For many years it was widely assumed that, for the purpose of making engineering studies, two param- eters were sufficient to describe the single-phase flow properties of a producing formation, i.e., the absolute permeability and the effective porosity. It later became evident that the concept of directional permeability was of more than academic interest; consequently, the degree of permeability anisotropy and the orientation of the principal axes of permeability were accepted as basic parameters governing reservoir performance. 1,2 More recently, 3-6 it was recognized that at least one additional parameter was required to depict the behavior of a porous system containing regions which contributed significantly to the pore volume but contributed negligibly to the flow capacity. Microscopically, these regions could be "dead-end" or "storage" pores or, macroscopi-cally, they could be discrete volumes of low-permeability matrix rock combined with natural fissures in a reservoir. It is obvious that some provision for the inclusion of all the indicated parameters, as well as their spatial variations, must be made if a truly useful, conceptual model of a reservoir is to be developed. A dichotomy of the internal voids of reservoir rocks has been suggested.7,8 These two classes of porosity can be described as follows: a. Primary porosity is intergranular and controlled by deposition and lithification. It is highly interconnected and usually can be correlated with permeability since it is largely dependent on the geometry, size distribution and spatial distribution of the grains. The void systems of sands, sandstones and oolitic limestones are typical of this type. b. Secondary porosity is foramenular and is controlled by fracturing, jointing and/or solution in circulating water although it may be modified by infilling as a result of precipitation. It is not highly interconnected and usually cannot be correlated with permeability. Solution channels or vugular voids developed during weathering or burial in sedimentary basins are indigenous to carbonate rocks such as limestones or dolomites. Joints or fissures which occur in massive, extensive formations composed of shale, siltstone, schist, limestone or dolomite are generally vertical, and they are ascribed to tensional failure during mechanical deformation (the permeability associated with this type of void system is often anisotropic). Shrinkage cracks are the result

By G. A. Vissac

The only fine coal washer with proved automatic controls, the feldspar jig is capable of good efficiencies even at low separating gravities, handles a variety of products, and treats 150 tph and over. IN continental Europe the feldspar jig is used almost exclusively for cleaning the fine sizes of coal; it operates with an artificial bed, made up generally of pieces of feldspar or of any other hard and heavy rock with good cleavages. A primitive form of this machine requiring close supervision was used in this country years ago, but output was small. The modern machines are generally air-pulsated Baum-type jigs, with capacities up to 150 tph, and treat sizes from % in. to 0, usually previously dedusted or deslimed. These jigs are equipped with very accurate and sensitive automatic controls insuring satisfactory operation even at low separating gravity and with feeds irregular in quality and quantity. This paper will deal more particularly with coal cleaning separation, but the technical studies apply equally to all minerals. Theory of the Feldspar Jig To understand fully the operation of the feldspar jig and to appreciate its potential value, it is necessary to measure and analyze the motion of small particles of various specific gravities, falling in various media, under various conditions. As no practical formulas were available for that purpose, the writer has developed new theories and established new sets of formulas, which cannot be demonstrated here. They are applicable to all minerals and to many problems of coal or mineral concentration. These formulas do not pretend to absolute accuracy, but results, tested by comparison with experimentally recorded data, are sufficient for practical purposes, as they give definite indications of trends, relative values, and order of magnitude of the many factors involved. For pulsating motions no earlier formulas are available. In the case of continuous currents or motions, such formulas as Rittinger&apos;s, Allen&apos;s, or Stokes&apos;, over a hundred years old, apply only to special shapes and sizes of particles, involve factors and coefficients difficult to measure or appreciate, and deal only with free motions. They are of little use in solving practical problems. Units used with the formulas presented here are in the CGS system. Lateral friction and suction have been taken into account according- to methods currently used in aerodynamics, and only hindered motion is dealt with. As the feldspar jig is primarily intended to create a fluid bed condition, the mixture of fluids and solids obtained is homogeneous, and the apparent specific gravity of the separating medium is calculated accordingly. For instance, if 20 pct of the volume is coal and shale, with an average of 1.60 sp gr, in a liquid of 1.10 sp gr, the actual specific gravity of the separating fluid is 0.80 x 1.1 + 0.20 X 1.6 = 1.20. The fundamental equations are as follows: The general motion of a particle falling in a medium of specific gravity equal to u, hindered motion. is dv/dt = g/Ls [ L (s-u) - u/g v2 [1] g = 981 cm per sec, acceleration of gravity, and L = thickness of the particle of specific gravity (s) and velocity (v). The initial acceleration dv/dt = g (1-u/s ) [2] is a function of the specific gravity of the particle, only for a given medium, and independent of its dimensions—an ideal separating condition. The terminal velocity Vt = gxL- (s-u)/u, [3] for a given value of u, is a function of s and L. Two particles, L, s, and L2 s2, are said to be equivalent if they have the same terminal velocity, and accordingly L1/L2=s2-u/s1-u. [4] Two such particles cannot be separated by continuous currents, in any concentrating device. Because the initial accelerated motion, which provides for efficient separation, is followed by a constant velocity, which has a limiting effect, it is essen-

Jan 1, 1956

By J. Geertsma

The relationship between porosity and the speed of propagation of acoustic waves in fluid-saturated porous rocks as measured by the Sonic log and by ultrasonic tecbniques is analyzed. Biot&apos;s continuum theory 1 is used to explain the difference in acoustic wave propagation between a dry and a liquid-saturated porous material. The porosity is a variable in this theory. However, the acoustic wave propagation in the dry rock depends too on porosity, and this dependence is not predicted by the theory. Frequently in dry sandstones, a nearly linear relationship between reciprocal acoustic wave velocity and porosity is observed in the low-porosity range. The physics behind this behavior is outlined. An empirical relationship of the form, 1/v = A + F, applies accordingly for many porous dry rocks, provided the porosity is the only variable. The presence of a liquid in the pores changes the value of B, and this change is found to be in agreement with the Biot theory. The time-average relation introduced some years ago2 results in an equation of the same type — 1/V = F/Vf + (1 - F)/Vr — but is not based on a sound physical picture. Still, this relation sometimes predicts approximately correct A and B values. Carbonate rocks with their complicated pore structures sometimes show a different relationship between wave velocity and porosity, unfavorable for log interpretation. Examples are presented. The simultaneous presence of calcite, dolomite and anhydrite, with their different grain densities and matrix compressibilities, complicates acoustic-log interpretation in carbonate rocks still further. Other complicating effects of acoustic-log interpretation are discussed. They are related to the influence of shale streaks and natural fractures on the average wave velocity observed by the logging tool and to the effect of adsorption phenomena on wave propagation in unstressed rocks particularly in sandstones. INTRODUCTION The velocity of propagation of sound waves in porous sediments as a function of depth is a quantity frequently measured. Velocity loggers of various d-sign may be used, for which velocity is defined as the. distance between a wave generator and a wave detector (one-receiver system) or the distance between two wave detectors (two-receiver system), divided by the shortest time required for a vibrational pulse to cross this distance. Velocity loggers ale used to assist in seismic prospecting, to differentiate between the various types of sedimentary rock layers, and also to determine rock porosity and pore fluid content. This paper is especially concerned with the relation between velocity and porosity and the effect of the pore fluid. Most velocity-log interpretation in terms of porosity is based on a formula advocated by Wyllie, et al,2 suggested earlier by Hughes and Jones,3 and known as the "time-average relation". This formula app1ies for a model consisting of a layered system of parallel alternating slices of a solid and a liquid, crossed by the wave path perpendicular to the interfacrbs. This must be an unsuitable model for the derivation of wave propagation properties of liquid-saturated porous media. It suggests that only rock matrix and fluid properties influence the wave velocity. The surprising fact, nevertheless, is that applications of this formula to clean sandstones under specific conditions (sufficient effective stress) results in porosity predictions which are close to reality. Other authors, such as Gassmann,4 White and Sengbush,5 Brandt,6 and Hicks and Berry,7 make use of more realistic models to arrive at a wave velocity formula for porous sediments. Their models consist of various liquid-saturated packings of frictionless spherical grains and are, therefore, also highly idealized. A number of experimental results obtained in the laboratory can be better explained with these model theories than by the time-average relation. A study of a general character, not depending a priori on a special model of the porous structure, was still lacking.

By Donald D. Carr, Lawrence F. Rooney

Perhaps no other mineral commodity in this volume has as many uses as limestone and dolomite. These carbonate rocks are the basic building blocks of the construction industry, the material from which aggregate, cement, lime, and building stone are made. Carbonate rocks, and their derived products, are used as fluxes, glass raw material, refractories, fillers, abrasives, soil conditioners, ingredients in a host of chemical processes, and much more. Although the use of limestone and dolomite goes back before recorded history, new uses and products are continually being found. Carbonate rocks form about 15% of the earth&apos;s sedimentary crust and are widely available for exploitation. Found extensively on all continents, they are quarried and mined from formations that range in age from Precambrian to Recent. Reserves of carbonate rock are large and will last indefinitely, although high-purity deposits may be absent or have limited availability in certain areas. Limestone and dolomite are so useful and so abundant that in 1981 about 644 million tons were produced in the United States. In fact, about 75% of all stone quarried or mined in the United States was carbonate rock. Sand and gravel was the only mineral commodity produced in greater quantity. In this chapter we have sought to provide the general reader an overview of carbonate rocks, their composition, distribution, production, and uses, with emphasis on developments since 1960. In doing so we borrowed extensively from our report on the applied geology of limestone and dolomite (Rooney and Carr, 1971). The literature on the limestone and dolomite resources of the United States and Canada is so vast and diverse that one can easily become mired down in print. We tried to list those references that are easily accessible and that provide the economic geologist with a starting point for more intensive study. The table listing selected references on limestone and dolomite resources in Canada and the United States reflects our attempt to be selective rather than exhaustive and to direct the reader to the best sources of detailed information. Geology Mineralogy The principal carbonate rocks used by industry are limestone and dolomite. Limestones are sedimentary rocks composed mostly of the mineral calcite (CaCO3) and dolomites are sedimentary rocks composed mostly of the mineral dolomite (CaCO3 • MgCO3). Aragonite (CaCO3), which has the same chemical composition as calcite but a different crystal structure, is economically important only in modern deposits, such as oyster shells and oolites. Aragonite is metastable and alters to calcite in time. Some other carbonate minerals, notably siderite (FeCO3), ankerite (Ca2MgFe (CO3)4), and magnesite (MgCO3), are commonly found associated with limestones and dolomites but generally in minor amounts. Because of their similar physical properties, the carbonate minerals are not easily distinguished one from the other. Specific gravity,

Jan 1, 1983

By M. C. Leverett, M. E. True, W. B. Lewis

This paper states the theory underlying the design of two kinds of dimensionally scaled models of parts of idealized oil fields. One of these simulates an oil well and its surrounding sand for a distance of 16 ft. radially from the well. The other model simulates linear flow through layered sands. Construction and operation of the models are described and typical data are given. The unique features in the design of the models are: (I) the treatment of the permeability-viscosity quotient as a single variable and (2) the use of previously reported experimentally developed relations between relative permeability and saturation and capillary pressure, saturation, porosity, permeability and interfacial tension. The models at present simulate flow of oil and water only, through unconsolidated sands. They are superior to unscaled models, which may give rise to erroneous conclusions. The models are designed to study the desirability of various methods of well completion and the effects of various factors on recovery of oil from layered sands. Introduction Recent work on the fundamental behavior of flowing mixtures of fluids in sands has resulted in the development of some new theory in this field, much of which has been presented in previous papers.l-4 However, the application of this theory to the behavior of oil wells or oil fields has proved too complicated for ordinary mathematical techniques, except in very highly simplified cases. A practical solution to this problem appears therefore, to lie in the use of dimensionally scaled models of elements of oil fields. The object of this paper is, first, to present the theoretical and experimental basis of the design of such models, and, second, to describe their construction and operation. Experiments now under way on the models already built will yield material for other reports, which will be devoted to interpretation of the results in terms of oil-field behavior. So far, simple dimensional models of two kinds have been built and are being operated. These models represent different parts of an idealized oil field but obey the same underlying theory. The only fluids represented are oil and water, since certain questionable approximations are necessary for gas-liquid systems. General Theory Dimensional models have been used in other fields for a long time, and their theory is adequately given in the literature. 5 Here, it suffices to state that each physical element entering into a dynamic relation may be given a definite description in terms of three fundamental, but arbitrarily chosen, dimensions. These usually are mass, length, and time, symbolized here by the small letters m, I, and t. Two systems are said to be dimensionally similar or dimensionally scaled if for each element in one system there is a corresponding element in the second system and the ratios of the numerical magnitudes of corresponding elements are constant for all elements of the same dimensional type.

Jan 1, 1942

By W. B. Lewis, M. E. True, M. C. Leverett

This paper states the theory underlying the design of two kinds of dimensionally scaled models of parts of idealized oil fields. One of these simulates an oil well and its surrounding sand for a distance of 16 ft. radially from the well. The other model simulates linear flow through layered sands. Construction and operation of the models are described and typical data are given. The unique features in the design of the models are: (I) the treatment of the permeability-viscosity quotient as a single variable and (2) the use of previously reported experimentally developed relations between relative permeability and saturation and capillary pressure, saturation, porosity, permeability and interfacial tension. The models at present simulate flow of oil and water only, through unconsolidated sands. They are superior to unscaled models, which may give rise to erroneous conclusions. The models are designed to study the desirability of various methods of well completion and the effects of various factors on recovery of oil from layered sands. Introduction Recent work on the fundamental behavior of flowing mixtures of fluids in sands has resulted in the development of some new theory in this field, much of which has been presented in previous papers.l-4 However, the application of this theory to the behavior of oil wells or oil fields has proved too complicated for ordinary mathematical techniques, except in very highly simplified cases. A practical solution to this problem appears therefore, to lie in the use of dimensionally scaled models of elements of oil fields. The object of this paper is, first, to present the theoretical and experimental basis of the design of such models, and, second, to describe their construction and operation. Experiments now under way on the models already built will yield material for other reports, which will be devoted to interpretation of the results in terms of oil-field behavior. So far, simple dimensional models of two kinds have been built and are being operated. These models represent different parts of an idealized oil field but obey the same underlying theory. The only fluids represented are oil and water, since certain questionable approximations are necessary for gas-liquid systems. General Theory Dimensional models have been used in other fields for a long time, and their theory is adequately given in the literature. 5 Here, it suffices to state that each physical element entering into a dynamic relation may be given a definite description in terms of three fundamental, but arbitrarily chosen, dimensions. These usually are mass, length, and time, symbolized here by the small letters m, I, and t. Two systems are said to be dimensionally similar or dimensionally scaled if for each element in one system there is a corresponding element in the second system and the ratios of the numerical magnitudes of corresponding elements are constant for all elements of the same dimensional type.

Jan 1, 1942