Search Documents

By W. H. Justice

Cycling operations have been successfully conducted in the La Gloria Zone for eight years. The cycling of this zone is now virtually completed. It is proposed that this report shall review the cycling operations of this reservoir so that the industry can more accurately evaluate the cycling process as a production method. The La Gloria project is one of the largest cycling operations in the State of Texas. The gas rights in the field are unitized in three operating units. By mutual agreement Magnolia Petroleum is the designated operator of the project for the 12 operators owning the gas rights. The La Gloria Corp, is the owner and operator of the cycling plant. The development and operation of the field are under the direction of an Operators Committee comprised of a representative of each operator. The Operators Committee is assisted by an Engineering and Geological Sub-Committee. The La Gloria gas field is situated approximately three miles north of Fal-furrias, Tex., and lies principally in the extreme south end of Jim Wells County, with its outer limits extending into the north end of Brookes County, Tex. It is a simple elliptically-shaped dome which has a closure of 350 ft. The field encompasses a surface area of 5,500 acres. The gas reserves are found in 20

Jan 1, 1950

By F. M. Stewart, F. H. Callaway, R. E. Gladfelter

In this paper oil initially in place is calculated by the various methods commonly used for analysis of water drive fields using data available as time progresses. Rate and pressure are predicted by means of the same methods using data available at the end of two years' history. These predictions are then compared with subsequent performance permitting an evaluation of the methods, which are Modified Schilthuis, Simplified Hurst, van Everdingen et al (1952), and the Electric Analyzer. INTRODUCTION Several variations of the material balance equation, differing in the term used to evaluate water influx, are in common use for estimating oil in place, predicting rate and volume of water influx, and forecasting producing rates and reservoir pressure in water drive reservoirs. The electric analyzer with an electrical network set up to simulate the oil reservoir and aquifer and their contained fluids may be used for the same purpose. Examples of the application of one or another of these methods have appeared in the literature. The accuracy and reliability of the material balance for calculating oil in place have been covered.' However, there have been few articles in which more than one method has been applied to an actual field allowing a comparison. It is felt that a field example, simple in nature, where the Modified Schilthuis, Simplified Hurst, van Everdingen et al (1952), and the Electric Analyzer have been applied would be of benefit in providing an evaluation of these methods as to their comparative usefulness and reliability in calculating initial oil in place and predicting performance. The Torchlight Tensleep reservoir is a small field in a large aquifer. The reservoir oil is greatly under-saturated, having a bubble point of only 55 psi and 2 cu ft/bbl of gas in solution. Changes in the reservoir volume factor in the range of pressures encountered are small, meaning little error in converting production to reservoir conditions. The pressure has declined from 1,618 psi initially to 137 psi during the first four and one-half years. Water production has been negligible. For Torchlight, the material balance equation as discussed later is considerably simplified by the omission of terms having to do with free gas. Thus with simpler equations, less variables to consider, the large pressure drop, the absence of a gas phase, and rapid pressure-rate adjustments, which mean an unusual amount of history in a short time, the reservoir appears suitable for comparing the methods listed above. GENERAL INFORMATION ON THE FIELD The Torchlight field is located on the eastern side of the Big Horn Basin, Big Horn County, Wyo., approximately 3 miles east of the town of Basin. The Tensleep reservoir was discovered in September, 1947. It lies at a depth of about 3,000 ft and produces oil of 35" API gravity. A total of five wells have been completed in the Tensleep, but only four are now producing. The reservoir is well defined in areal extent in that a total of 13 wells have penetrated the Ten-

Jan 1, 1955

By C. R. Criss, R. J. McCormick

ThiS paper .summarizes the development and perform. ance of the Coldwater Oil field, Isabella County, Mich. Production is obtained at a depth of 3750 ft from a dolomite reservoir, which is probably both Rogers City and Dundee. There appears to be no persistent barrier to vertical fluid movement between the Rogers City and Dundee, and both formations are believed to have a common original oil-water contact. No we1ls are open to production at the level of the Dundee. The discovery of the field in 1944 was followed by the drilling of 81 producing wells and nine dry holes. By the end of 1952, the field had produced 12,763,000 bbl of oil, at which time 69 active wells were producing 3600 B/D oil and 21,000 B/D water. The reservoir has an effective natural water drive, and bottom-hole pressure has declined 75 psi from an original pressure of 1453 psi. Oil is undersaturated at reservoir conditions with a solution GOR of 512 cu ft/bbl and a saturation pressure of 1190 psi. INTRODUCTION Although the Coldwater Field is now in its ninth year of production, it ranks seventh in cumulative recovery, and second in current daily oil production in the State. Following complete development of the field in 1946, it has consistently contributed approximately 11 per cent of Michigan's daily production. Coldwater appears to have a greater pay thickness than most fields producing from the Rogers City forma- tion in the west central part of the Michigan basin. Except as influenced by the extent of closure and the thickness of the oil column, reservoir performance is believed to be representative of Rogers City fields in this part of the basin. For this reason, it is felt that the following presentation of factual data for the field will be of general interest. HISTORY AND DEVELOPMENT Discovery of the Coldwater field during Aug., 1944, followed favorable subsurface geological correlations and core tests drilled to a shallow stratigraphic marker horizon. The field is located in the west central part of the Michigan Basin in Coldwater Township (T16N-R6W), Isabella County. It is 60 miles west-northwest of Bay City, and near the geographic center of the southern peninsula of Michigan. The discovery well flowed 200 bbl of oil the first 24 hours from the Rogers City dolomite at a depth of approximately 3750 ft. Field production is principally from the Rogers City with a minor amount assumed to come from the underlying Dundee dolomite, although none of the wells are open to production at the level of the Dundee. Both the Rogers City and Dundee formations belong to the Devonian system (See Geologic Column-Fig. I.) In addition to Rogers City-Dundee oil production, sweet gas is produced from the Michigan Stray sand of the Mississippian system at a depth of 1400 ft. The gas reservoir has approximately the same areal extent as the underlying oil reservoir. Although gas production is an important feature of the Coldwater field, this report is concerned only with the oil reservoir. The productive area is approximately 3200 acres. The pool is developed on a 40 acre spacing pattern with wells drilled in the center of the north half of 40 acre tracts, except for a few off-pattern wells drilled along

Jan 1, 1955

By D. R. McCord

Data of a known statistical quality has been successfully used in a system of conventional fluid mechanics equations essentially free of empirical "conformance factors to prediet early detailed field 03 in a reservoir where gravity drainage is important and oil and rock properties vary con. siderably. The system is useful where comparative recovery predictions for one or more cases of natural depletion and gas cap pressure maintenance are needed and the reservoir is not too heterogeneous. INTRODUCTION Many authors, have examined various facet* of the gravity drainage-gas cap pressure maintenance problem over the last several years. Two of them. Burtchaell and Elkins, have published practical studies. In the case described in the present paper the best development and operating plan for one of the important reservoirs of the Bolivar Coastal Field in Western Venezuela was urgently needed. It was requested that sound comparative economics be presented along with an engineering evaluation of no additional drilling, maximum drilling. several cases of gas pressure maintenance, of water flooding and combination cases. A rapid analysis of the effect on recovery of the range of crude and rock properties indicated that none of the published methods using average properties would yield adequate comparative predictions. This paper covers the predictive system developed to deal with natural depletion and gas cap pressure maintenance. The large amount of data available for this field permitted statistically reliable correlations to be made of most rock and oil properties, and the evident smooth areal variations of these properties encouraged the construction of a lineal mathematical model of material balance and fluid flow equations, which could be manipulated by a modified iterative procedure. The manipulation is much more adaptable to computing machines than manually operated calculators and subsequent studies have been prepared for the IBM Card Program Calculator. A matrix as opposed to linear model has been tentatively worked out by other members of our group to deal with more heterogeneous reservoirs. It is more adaptable to computers such as IBM-701 or the Remington-Rand Univac. In 1950, when the system reported below was developed, no rapid computing machines were available to the author. Their existence and promise encouraged an initial approach by hand. Some details of oil and rock Property correlation are included because there are a surprisingly large number of medium-sized reservoirs in the world that show systematic property variations which are significant for a rigorous analycon. All large reservoirs (over half a billion barrels recovery) with which the author is familiar contain at minimum a significant temperature variation. Using the rock and oil property correlations, past pressure and production, the gas-oil ratio history of the reservoir was computed and found to be in close agreement with field measurements. The future performance under cases of natural depletion, gas injection of 100 MMcf/day, and of 130 MMcf/day, was then computed; a 70 MMcf/day case was also carried far enough to provide data for comparative economics. The 130 MMcf/day case proved to be the most attractive. A 50,000 horsepower gas turbine driven centrifugal compressor plant to take this quantity of adjacent field gas from 30 to 1,950 psia and inject it into the reservoir is under construction. The total cost of the project will exceed 20 million dollars.

Jan 1, 1953

By D. L. Garthwaite, F. M. Stewart, F. K. Krebill

Pressure has been maintained in the Elk Basin Ten-deep reservoir since the initiation of inert gas injection in Sept., 1949. Oil is being produced under conditions favorable for gravity drainage, including high angle of dip, appreciable structural closure, and fairly good per~neability. Results being obtained are high per cent recovery of oil-in-place, sustained field productivity, reduced operating problems, and recovery of plant products. The high recovery has been calculated by alpplication of the gravity drainage theory described in an earlier AIME paper.' It is also confirmed from field performance by comparing nil recovery to date with gas cap spuce voided. GENERAL INFORMATION The Elk Basin field lies in Park County, Wyo., and Carbon County, Mont., approximately 50 miles east of Yellowstone Park. Located at the northern end of the Big Horn Basin. the field is situated on an elongated asymmetrical anticline (Fig. 1). The pressure maintenance project covered by this paper involves the Embar-Tensleep reservoir in this multi-pay field. Discovered in Nov., 1942, it has a proved productive area of over 6,300 acres. It consists of the 210-ft thick Tensleep sandstone of Pennsylvanian Age overlain by 40 ft of Embar dolomite of Permian Age. The two formations apparently are in communication with each other and are produced as a common source of supply. As 98 per cent of the reserves are contained in the Tensleep sandFtone, in this paper the Embar-Tensleep reservoir will he called simply the Tensleep. The reservoir is found at an average depth of 4,900 it, is approximately 7 miles long and 2 miles wide, and has a maximum oil productive closure of 2,330 ft. The strata dip an average of 21" on the west flank and 45" on the east. There are presently 132 Tensleep wells in the field drilled on 40-acre spacing; of these, 17 are shut in because of high gas-oil ratio (GOR) and eight are gas injection wells. Rock characteristics and related information for the pay section follow: Porosity, avg 10.7 per cent Permeability (air) from cores, avg 118 mds Permeability from PI'S 91 mds Connate water, avg 8 per cent Initial pressure (-400 ft) 2,234 psia Pressure April 1954 (-400 ft) 1.327 psia At high structural positions much of the better sandstone has connate water content of only 2 to 5 per cent as indicated by the coring of wells with oil in nearby Tensleep fields. In special imbibition tests, samples readily imbibed oil but repelled water. In tests on crushed samples, the sand grains rapidly settled in oil and carbon tetrachloride but formed a stable unsinking floc in water. This information indicates that the Tensleep sand is preferentially oil wet. The oil produced averages 29' API gravity and has a high sulfur content. Extensive sampling and testing have established the fact that characteristics of reservoir oil in this field vary greatly with structure, see Fig. 2. Additional detail on the variation of oil characteristics with structure at Elk Basin may be found in Reference 2. At original reservoir pressure the oil was undersatu-rated with gas, the bubble point being 1,250 psi at the crest but only 500 psi at the lowest elevation sampled. As noted on Fig. 2, other fluid properties—reservoir volume factor. solution ratio, gravity, viscosity, and hydrogen sulfide content—likewise vary greatly. These variations complicate mathematical analysis, for it is necessary to average these characteristics in order to enter them into reservoir engineering expressions. To further complicate analysis, sand-face pressure also varies greatly because of the increase in fluid head down structure. As an approximation, average properties under initial conditions may be read from Fig. 2 at the -450-ft elevation.

Jan 1, 1956

By M. A. Selim, W. H. Somerton

Thermal expansions and heats of reaction of three typical sandstones were measured in the temperature range of 25°to 1,000°C. The significance of these data in subsurface heat-transfer calculations is discussed. Linear thermal expansions were measured both parallel and perpendicular to the bedding planes. Volume expansions are reported as the perpendicular expansion plus two times the parallel expansion. Expansion behavior of the sandstones was found to be controlled by the expansion characteristics of the quartz content. Differential expansion of the quartz grains and other minerals included in the sandstones caused permanent deformation of the heated samples after they were cooled to room temperature. Structural damage resulting from heating is probably an important cause of the reduction of thermal conductivity of heated samples. Measurements of heats of reaction were based on the known heat required for a- quartz inversion. Thermal reactions, which probably include deby-droxylation of clay minerals and decomposition of carbonate minerals contained within the samples, were found to require more than one-fourth of the amount of heat necessary to raise the temperature of the rock alone. In shales and limestones, the reaction heat could be substantially greater than that required from specific heat considerations alone. INTRODUCTION In earlier work, l-3 methods of measuring thermal properties were developed and thermal data for several sedimentary-rock types were reported. Data included thermal conductivity, specific heat and thermal diffusivity. In addition, calculations of heat transfer in subsurface formations may require data on the thermal-expansion characteristics of rock and the heats of reaction of mineral constituents. Thermal expansion of rocks is relatively small in magnitucle and, from the standpoint of change in bulk density, has only minor effects on heat-transfer cliaracteristics. Thermal-expansion behavior, however, may have significant effects on the structure of rocks. Differential expansion of different minerals and along different crystal-lographic axes may result in structural damage which could effectively reduce the thermal conductivity of the rock. Many mineral constituents of sedimentary rocks undergo phase or related changes when heated to sufficiently high temperatures. Although some of these reactions (such as the a-ß quartz inversion) are reversible and the heat absorbed is returned to the system upon cooling, many of the major reactions are irreversible. Through certain temperature ranges, the additional heat energy required to complete the reactions may be nearly as great as that necessary to raise the temperature of the rock in the absence of thermal reactions. If the heat requirements for these reactions are not considered, serious errors in heat-transfer calculations may result. Thermal-tzxpansion and thermal-reaction data have been obtained for three typical sandstones. Methods of measurement and results obtained are reported and discussed in the following. THERMAL EXPANSION OF SANDSTONES The linear expansions of three outcrop sandstones (Bandera, Berea and Boise) have been measured in the temptbratue range of 25° to 1,000° C. Expansion measurements were made in the directions parallel and perpendicular to the bedding. The differential expansion apparatus used is that described by Mitoff and Pask.4 The test sample is heated in an electric furnace with a temperature rise of 6°C/minute. The lengthening of the samples upon heating is compared with the small and known expansion of a fused silica rod. The change in length is transmitted to an X-Y recorder by means of a Statham transducer. A maximum error of ± 1.5 per cent has been obtained with this apparatus

By R. E. Jenkins

The Bell Canyon member of the Delaware Mountain group has yielded quite a large number of fields in which completion and production problems have been numerous and complex. Reserves are difficult to estimate due to the problem of evaluating the formation water saturation, and the feasibility of water flooding is questionable. Extensive laboratory investigations were undertaken to determine if normally measured rock and contained fluids properties could account for many of the peculiar characteristics of the formation. A statistical study or routine core-analysis data was made. Capillary-pressure data were obtained by several techniques. Fresh-water and brine permeabilities were measured and relative permeabilities of oil, water and gas were determined. Some qualitative wetta-bility tests were performed, and the petrography of a few thin sections was studied. The results of this work are presented, with typical or average data shown for each test. The permeability -to- water characteristic of the Delaware was found to be abnormally poor. Correlations of laboratory and field performance data indicate high water saturations in much of the formation. Pore geometry is highly uniform. INTRODUCTION The Delaware formation has undergone extensive exploration and development in the past few years. This has resulted in the discovery of quite a number of fields—producing oil for the most part. A large percentage of the productive wells have produced water along with the oil. Often there is no apparent pattern for the percentage of water cut experienced. In several instances vertical displacments of over 200 ft have been established between communicating wells in which the formation appears to be about the same, yet the structurally high wells yield a much higher water cut than the lower wells. Drill-stem tests have contributed little to well completion because they nearly always show small amounts of water (usually mud filtrate), regardless of what is produced later in the life of the well. Other formation evaluation methods have had no more than limited success in indicating ultimate productivity. Consequently, predicting the type of fluid productivity of Delaware wells has been very difficult. Hydraulic fracture treatments are almost universal, and unstimulated Delaware wells have not exhibited the productivity that normally would be expected of wells with the specific permeabilities encountered. Some operators have indicated that their Delaware wells need frequent re-stimulation. The uncertainty of reservoir water saturations reported for the Delaware has made difficult the task of evaluating oil in place in the reservoirs and estimating primary recovery. The unreliability of this data also will cause feasibility studies for water flooding or other secondary-recovery methods to be problematical. Several opinions have been voiced to explain why the Delaware performs as it does. Chronically poor well completions, dynamic water and tilted water tables, capillary inequi-libria, peculiar relative permeability characteristics and oil-wetness are among the suggested reasons. The dual purpose of this paper is (1) to describe the principal rock characteristics, and (2) to present an evaluation of those characteristics so that approaches can be made to explain the irrational behavior of wells completed in the Delaware. STUDY OF ROUTINE CORE-ANALYSIS DATA A statistical study of routine core-analysis data has been made to relate fundamental rock properties and characteristics of the formation to well productivity. For the most part, well production data were obtained from operators and augmented by production data from commercial information services. Only data which were available in detail and which were considered to be accurate were used. The core data used were obtained on samples from the perforated or open interval only and, thus, may not always represent the total productive interval in a well. All the core-analysis data included in this study are from the Bell Canyon member of the Delaware and were obtained by "plug-type" techniques, as opposed to whole-core or full-diameter techniques. Porosity values were obtained by the summation-of-fluids technique, permeability values are air permeabilities corrected for slippage to an equivalent liquid permeability, and the oil and water contents were determined by means of a high-temperature downdraft retort. All the procedures used are discussed at length in Ref. I. The core-analysis data were sorted

By R. N. Upadhyay, M. Maleki, I. Fatt

Conventional laboratory core analysis tests on samples of two limestone reservoir rocks indicate that about 20 per cent of PV is in dead-end pores. These tests (electric logging formation factor, mercury injection capillary pressure and miscible displacement) were carried out on 1/4-in. diameter test plugs. Test results show a clear dilference between these samples and sandstone or homogeneous limestone reservoir rock. Although the amount of dead-end pore space can be only roughly estimated, the presence of such pore space seems clearly indicated. Pressure transient studies also show presence of deadend PV. Although they do not give quantitative results, pressure transient data yield a reasonable estimate of the size of the neck connecting dead-end pores to the main flow channels. INTRODUCTION Equations conventionally used to describe reservoir flow behavior contain the implicit assumption that all connected pore spaces contributed to both porosity and permeability. Several authors3-7 have pointed out the changes in pressure transient behavior and in electric log interpretation that may result if this assumption is incorrect and, instead, dead-end or cul-de-sac pores are present. There is a need for laboratory tests that can detect presence of dead-end pores in core samples. With such information on hand the petroleum engineer can make more rational use of the mathematical tools now available for analysis of reservoir flow behavior. This paper describes laboratory studies designed to detect and, if possible, give a quantitative measure of dead-end PV in laboratory-size core plugs. Three reservoir rocks were used, two of which were limestones suspected of having deadend pore spaces and a well-known sandstone, used as a comparison standard, in which there is believed to be little or no dead-end pore space. All the studies were designed to measure the natural dead-end PV; i.e., the pore space which is dead-ended because of rock structure. During multiphase flow in a rock without dead-end pores, some parts of one of the phases can become surrounded by the other, thereby giving (for certain flow behavior) an effective dead-end pV8,9. Such behavior will not be described here.

Jan 1, 1967

By L. L. Handy

Two equations are compared for calculating imbibi-tion rates of water into porous media. A linear relation was found to apply when the square of the volume of water imbibed was plotted against time for water imbibing in natural rock materials. The effect of adsorbed organic materials on surface wettability was determined from measurements of rate of imbibition, water saturation and water permeability. INTRODUCTION An understanding of the mechanisms by which water displaces oil from porous media requires an understanding of the role of capillary forces in the displacement process. The capillary forces are determined by properties of the water-oil-solid surfaces. Consequently, some measure of the surface properties of reservoir rocks is necessary if displacement behavior in oil reservoirs is to be predicted from laboratory measurements. Several experimental procedures have been proposed to assign quantitative wettability indexes to reservoir rock surfaces. The most recent of these proposals are those of Bobek, Mattox and Denekas' and Amott.' These indexes are designed to show a continuous variation from the preferential oil-wet to the preferential water-wet systems. They require measuring some property of the rock which is a function of surface wettability. The quantities are measured on unaltered core material and compared with values obtained for known oil-wet and water-wet extremes on the same material. These methods are useful but are semi-empirical in nature. They have the disadvantage that the measured quantities may be functions of other variables in addition to surface wettability. Some understanding of the behavior of water-oil-rock systems may be obtained from studies of the simpler water-air-rock systems. Surfaces of a given wettability are more easily maintained for the latter. Also, certain approximations can be made for the flow of water and air which are not justified for water-oil flow. The theoretical and experimental work discussed in this paper is part of a study of the role of capillarity in the dynamic displacement of air by water. The expec- tation is that observations for the water-air system will lead to generalizations which are applicable to water-oil systems. THEORY Imbibition is defined as that process by which one fluid displaces another from a porous medium as a result of capillary forces only. Differential equations have been derived to describe the displacement of oil by water during imbibition, but these equations have not been solved.".' With certain assumptions, however, solutions are possible for water-air imbibition.

By C. Gatlin, N. E. Garner

Results of impulsive wedge penetration tests on two synthetic, plastically deforming rocks are presented. Basic data obtained were force-time, displacement-time, and force-displacement curves for the impacts, plus the crater geometry. Wedge geometry and blow frequency were varied over a considerable range. The synthetic rocks consisted of wax-sand mixtures; two waxes of diflerent ductilities were used to provide variable "rock" characteristics. Conventional triaxial tests showed that these synthetic rocks exhibited force-deformation curves and Mohr envelopes quite similar to real rocks, except that strengths were much lower. Measured forces from static penetration tests agreed closely with theoretical values; however, dynamic force values were much higher than the static. These latter disparities are attributed to the viscous nature of the waxes. Thus the utility of these or similar rock models must depend on the scaling of rock viscosity, which is as yet unknown for impulsive loadings at elevated stress states. It appears, however, that some macroscopic, static phenomena may be studied with wax-sand rock models. INTRODUCTION The resistance of solid materials to indentation or perforation by projectiles or other penetrators has been studied by workers in many areas. Despite these efforts no universally accepted laws or formulas are available for describing experimental observations. In the metals field the force-deformation behavior of impacting bodies is often analyzed by the Hertz law for elastic collisions, the Meyer law if plastic deformations occur, or some combination of both.' The similarities of these expressions to empirical drilling formulas of the oil industry are apparent. Beginning with the basic contributions of Simon and co-workers at Battelle,' a number of experimental papers concerning the reaction of rocks to vertical impact have appeared in the U. S. mining and petroleum literature.'-' Most published data have, to date, been obtained at atmospheric pressure, although some early high pressure information was reported by Payne and Chippendale.8 Maurer" has recently utilized available brittle impact data to develop a drilling rate equation based on the experimentally observed proportionality between crater volume and blow energy. His result agreed with earlier efforts by both Somerton, who used dimensional analysis, and Outmans, who used plasticity theory. It has long been known that rocks exhibit different modes of failure depending on the state of stress. The literature in this area is considerable; however, papers by Bredthauer, Handin and Hager,13 nd Robinson", are adequate to illustrate the point. Since rocks flow plastically at certain triaxial stress conditions, the mathematical theory of plasticity has been used to analyze the rock drilling problem. Cheatham'" has altered the wedge identation solution of Randtl to rocks, and has developed useful equations for penetrator forces under a variety of conditions. Outmans" has utilized Hill's solution in a similar manner to develop a drilling rate equation. Both Cheatham and Outmans used the linear Mohr-Coulomb rule to relate rock strength and confining pressure. The actual stress at the hole bottom is not easily ascertained, although photoelastic studies by Galle and Wil-hoit," plus the analytical treatment of Cheatham and wilhoiti8 provide some insight. Consequently it is not clear to what extent the highly idealized rheological model of a perfectly plastic solid can be realistically applied to the rock drilling problem. This paper is the first report on a long range experimental study of crater formation in rocks at elevated stress states. The data presented here are from the first phase of the project. Data obtained from impulsive wedge impacts on two synthetic, plastically deforming rocks are presented. MODEL ROCKS Geologists have long been faced with modelling the behavior of the earth and, as a consequence, have studied scaling problems in some detail.' In general, their main problem is handling the wide disparity between laboratory and geologic time. In our studies the time effects (blow velocity or rate of loading, blow duration, etc.) were essentiafly the same for both. model and prototype, as were were geometry and tooth penetration. Thus application of available scaling laws suggests that Similarity is obtained if the stress-strain curves of model and prototype are similar." For this reason Hubbert and Willis''

By P. D. White, J. T. Moss, J. S. McNiel

Results of experimental work on the in situ combustion process were first published in this country in 1953' when Kuhn and Koch described results of a three-well test in Jefferson County, Okla. Shortly there-after Grant and Szasz2 described field studies in Nowata County, Okla., in which several modifications of the basic combustion process were investigated by injection of various oxidizing gas mixtures other than air. In the past five years speculation concerning this method of oil recovery has continued, with a number of interesting papers on various aspects of in situ combustion appearing in the literature3,4,5,6,7. However, only one of these' provided actual engineering information on pattern-type field experiments. The primary objectives of the Oklahoma three-well test were to check operation of specially designed electrical ignition equipment and evaluate problems associated with initiation of combustion in a natural oil sand reservoir. Immediately following successful completion of this test in the Pontotoc sand, plans were activated for a second field experiment to be conducted in the same sand a short distance away. The objective of this second field test1 as to establish a relationship between field observations and information derived from laboratory experimentation and computations. Secondary objectives of the test were to develop operating tech- niques and define operating problems, and to obtain information that would assist in the economic and engineering appraisals of the process. Careful planning was required since this test was designed with the hope of providing basic information on the characteristics and operation of a combustion drive. It was necessary to include test facilities which would not only permit precise control and measurement of the process variables but would also provide sufficient flexibility of control to insure that variables in the combustion process could be studied during the course of the test under a wide range of conditions. In order to accomplish these objectives elaborate instrumentation and several control wells were included in the design of the test installation. This type installation was felt to be well suited for an experimental study, but was in no respect intended as a counterpart of a commercial-scale operation. PRELIMINARY LABORATORY WORK An important objective of the field experiment, as stated previously, has been to obtain reliable data on operating variables which might be compared with similar data obtained in laboratory combustion tests. The value of successful correlations between quantities that can be measured in the laboratory and those obtained in a full-scale commercial venture has been well established for other industrial processes. The difficulty of following a process underground should make such correlations even more essential here. Also. approaching the field test phase of development of in situ combustion, it was considered desirable to precede design of test equipment with laboratory combustion studies under a variety of conditions to establish likely ranges of variables. These included temperature, air injection rate, frontal advance rate and produced gas composition.

By P. F. Gnirk, J. B. Cheatham

An experimental investigation has been made to study the effects of tooth angle and confining pressure on the force required for indentation of dry rock samples under confining pressure. In these experiments static single-tooth indentations were made on two standstones, a limestone, a schist and a slate under confining pressures of 5,000 to 15,000 psi. Sharp wedge-shaped teeth with included angles ranging from 30°to 105°were used. In all cases the force-displarenzent curves are approxinlately linear for rocks tested under confining pressure, as contrmted with the discontinuous curves obtained for silnilar rocks at atmospheric conditions. Photonzicrographs of thin sections through the craters indicate very little evidence of chip formation during indentation. In general the experimental data lie between calculated values for perfectly rough and perfectly smooth tooth-rock interfaces. INTRODUCTION The variables associated with ctrilling at depth are numerous and, in many instances, are not easity isolated and controlled. Also, many of these variables are interdependent; however, if the variables which are sufficiently independent of each other are held constant and are appropriately introduced, the dependent variables can be more easily studied and evaluated. It is the purpose of the present investigation to study experimentally the simplest case of bit tooth-rock interaction. In essence, the experiments have involved the static loading of sharp wedge-shaped tools on the surfaces of effectively semi-infinite dry rock samples subjected to confining pressures. During this investigation the following quantities have been maintained constant: 1. Pore pressure: atmospheric (i.e., the rock is dry). 2. permeability: although the rocks may be permeable, they are "jacketed" in such a manner that the fluid is not permitted to flow through them. 3. Temperature: 75F. 4. Rate of loading: essentially static (approximately 0.002 in./second). 5. Bit tooth: a sharp wedge-shaped tooth loaded normal to the rock surface. 6. Rock surface: smooth and flat. 7. Drilling fluid: hydraulic oil. In addition to the above conditions, the rock surface was free of a layer of cuttings. There was no indexing to previously indented surfaces, and the influence of corners of the borehole was neglected because the rock samples were effectively semi-infinite. The following quantities have been varied, and the resulting effects have been studied: 1. Rock a) Type (e.g., limestone, sandstone, etc.). b) Mechanical properties ) Necessarily dependent c) Porosity on type of rock. 2. Bit tooth (included angle of wedge-shaped tooth). 3. Depth of penetration. 4. cobfining pressure. On the basis of laboratory experiments which have shown that many rocks exhibit ductile behavior at confining pressures corresponding to depths usually reached during drilling,'.' a theory has previously been advanced which predicts the force required for penetration by a single sharp wedge tooth into a ductile* material.' The material is assumed to be isotropic, homogeneous and rigidly plastic (no elastic deformation) and to obey the Coubomb yield criterion (Figs. 1 and 2). Equations governing the force required for penelration are given for two extreme cases:** 1. A smooth, perfectly lubricated tooth-rock interface (fridionless boundary), in which the solution gives a minimum value for the penetrating force: in which F = force of the tooth, b = length of tooth edge,

By W. F. Hower, W. Brown

Large-scale sand consolidation tests were conducted in an effort to determine the reasons for the successes and failures of this method of sand control. Several different consolidating materials were used in treating both clean and bentonitic sands that were packed in a chamber having a capacity of 3.3 cu ft. The results were essentially the same for all of the different consolidating materials, The data show that low-viscosity consolidating materials pumped at a relatively slow rate gave the best results. Where the formation has produced sand, the treating fluids can compress the formation, thus permitting the channeling of fluids to another horizon. Pressure-packing these zones before attempting to consolidate is recommended. Sands containing more than 4 per cent of water-swelling clays are not good candidates for consolidation. It is indicated that loose sand, particularly when it is bentonitic, can be fractured during the placement of the treating fluids. INTRODUCTION Sand production in oil and gas wells has plagued the industry for many years, and numerous cures for this problem have been suggested. Most methods have been successful to a certain degree, but the great variety of well conditions that exist in the different areas has magni- fied the problem and limited the successful use of the various systems. Four review papers1-4 present a wealth of information concerning the degrees of success that have been obtained by the different sand-control methods. The bridging of sand grains by the use of gravel packs and screens has been quite successful. However, these methods do not leave the casing clear for all types of multiple completions, and the cure does not last for the production life of the well in some instance:;. The control of loose sands by sand consolidation with resins has never been as successful as desired. It has always been hoped that such a treatment would eliminate all sand problems for the life of the well, but. initial applications, starting in the middle 1940's, were only moderately successful. Lott, et a1,3 reported a success ratio of approximately 50 per cent and made the following conclusions. The highest percentage of successes were obtained where: a. Consolidation of a zone was made at the time of initial completion or prior to the production of sand. b. The interval treated was less than 12 ft in length. c. Between 30 and 50 gill plastic/ft of producing interval was displaced through the perforations. REASONS FOR SAND CONSOLIDATION FAILURES Our own experiences in the field of sand consolidation point toward the following conditions as the major reasons for the failure of sand consolidation attempts. 1. Mud-plugged perforations and mud invasion of the formation. 2. Sand in the casing covering all or part of the perforations. This sand could be either formation sand or one of the coarser sands used as propping agents in hydraulic fracturing. 3. Holes in the casing. 4. Channels behind the casing. 5. Attempting to treat too long a perforated section. 6. Too high a percentage of water-swelling clays in the formation. 7. Formations that have produced sand. Recent attempts were made to treat perforated sections ranging from 10 to 30 ft, in wells that have produced sand, by using a straddle packer that was raised and lowered through the perforations as the consolidating material was being pumped. In most instances, the pressure required to pump fluid into the formation varied considerably as the tool was raised and lowered. This suggested the possibility that significant differences in permeability were present or that only part of the formation had produced sand. There were times when a sudden break in pressure indicated that a fracture was being formed. Research conducted several years ago concerning the problem of the control of water in air and gas drilling indicated that shale sections could be fractured quite easily. In addition, it was determined that it was easier to pump fluids into shale bodies by fracturing the shale itself, or the interface between the shale and sand, than to pump into a fluid-saturated formation. Formations that produce sand are usually adjacent to shale bodies and frequently have shale streaks of various thicknesses inter-bedded in the sand. Therefore, where shale is exposed to fluid pressure it

By C. Gatlin, M. E. Chenevert

The effects of bedding plane orientation on the elastic constants and the yield strengths of three laminated rocks (one sandstone and two shales) and one isotropic rock (a limestone) were studied. The directional dependence of the elastic properties of these rocks was experimentally evaluated using a triaxial compression cell and auxiliary stress-strain measuring equipment. Symmetry of Poisson's ratios within the bedding plane suggested that horizontal isotropy exists, but the bedding planes do give rise to an appreciable difference between properties in the horizontal and vertical directions. For the three bedded rocks studied, Young's modulus was lower normal to bedding than along bedding. Yield strengths were determined at confining pressures from 0 to 12,000 psi in a triaxial compression cell. The rocks studied showed strength reductions as high as 40 per cent when the test specimen was oriented at 20°-30° to the bedding planes. The mechanical behavior of these rocks suggested that the rock properties of shear strength and/or coefficient of internal friction can vary with direction, depending on the particular rock tested. Tensile strengths were also measured and found to be lowest when failure occurred along bedding. This work shows that bedded formations exhibit sizable directional variations in both their elastic constants and yield strengths. It is suggested that these variations may be accounted for by using the "elastic laminate" model and the "variable coefficients" failure model. INTRODUCTION The nature of rock deformation at elevated pressures has been studied by many workers; papers by Handin,l,2 and Robinson3 illustrate the present state of knowledge. Most investigators have either chosen rocks which were as isotropic as possible (in order to avoid complications of data interpretation and analysis) or they have oriented their samples so that the effects of anisotropies (such as bedding planes) have been avoided. Of the studies performed, few were concerned specifically with mechanical anisotrop-ies. Griggs4 has presented limited data for specimens cut parallel, normal, and 45° to the bedding plane; he relates the strength anisotropy observed to the fabric (bedding) anisotropy. His tests were primarily concerned with large deformations (20 per cent strain); thus, no directional values of the elastic constants were reported. Handin has reported the results of similar experiments. Bott5has discussed rock strength anisotropies due to faults, cleavage, or bedding. He was concerned primari ly with determining the shear stress on such planes and did not mention the effect of friction. Jaeger6 later generalized Bott's work by taking friction into account and presented a limited theory for the failure of rocks having a "single plane of weakness", and also for rocks having a constant coefficient of friction, and a shear strength which varies with bedding plane orientation. Donath and Cohen,7 and Donathg have evaluated rock strengths from shale and slate specimens cut normal and parallel to bedding. A dependence of cohesive strength (7,) on the specimen orientation was also shown. Adler9 has also studied this problem and lists similar results. He assumes that all bedded rocks behave according to Jaeger's single-plane-of-weakness theory. Kalinin and Belorussov10 list results for strength tests parallel and perpendicular to bedding and use this information as a basis for hole deviation analysis. From the literature it is apparent that sedimentary rocks have been tested under widely varying conditions of stress; however, the assumption of isotropy is generally, but not always, made. Since geologic sedimentation often deposits sediments in very definite lavers.,. it seems that more svstema tic attention should be given to the possible effects of this natural bedding. Bedding, as used here, refers to visible regularities of grain size or orientation resulting from depositional processes. The mechanical strength anisotropies, which bedding may present, is not the only anisotronic property to be considered, since the elastic

Jan 1, 1966

By V. J. Sikora, T. S. Hutchinson

This paper presents a method for making a water-rlrive ana1gsis without prior knowledge of aquifer geometry and uniformity using a standard desk calculator. Although it is necessary to know the initial oil in place to use this method, this is a minor limitation to the scope of (he method since in most reservoirs it is possible to obtain reasonably accurate volumetric oil-in-place estimates. An equation is developed which relates pressure at the water-oil contact to the water influx rate as a function of time by a factor called the resistance function. The resistance function introduces the composite effect of the aquifer geometry and flow resistance disrribution. It is pointed out that the characteristic shape of this function makes it possible to start with on approximation of the function and successively improve the approximation until the correct resistance function curve is obtained. In this fashion water-drive estimates can be made without the limitation of assuming simplifed aquifer shapes and flow distributions. This is the novel feature of this development. Methods are given for extrapolating the final curve to calculate future aquifer behavior. Equations are developed for adjusting the pressures and water influx rates where it appears that possible errors in these quantities make it difficult or impossible to obtain a useuble resistance function curve without this adjustment. Application of the pressure build-up analysis techiiique to estimate. some of the aquifer properties is also presented INTRODUCTION A great many oil pools are the result of oil accumulation in some type of trap in an otherwise large and continuous porous stratum. The void space of this stratum outside of the oil pool itself is filled with water or brine. In analyzing performance of the oil pool surrounded by this water aquifer, it is quite necessary. in most cases, to include behavior of the water. When the pressure at any point in a fluid system is lowered, such as by opening a well in an oil sand, fluids in the immediate neighborhood of this point will begin flowing towards this lower pressure sink. As pressure in this area drops due to flow towards the sink, fluids from farther out will start to flow towards the lower pressure. As more and more fluid is removed from the system, the distance from the sink or well within which flow is occurring will continually increase; that is, thc region of disturbance will grow. If some rigid boundary such as a fault is reached by the disturbance, this area will cease to grow; but if some movable boundary is reached, such as a water-oil contact, the area will grow on out into the water, although rate of growth may and almost always does change. The relation between amount of pressure drop and amount of fluid flowing at any point and at any time in an aquifer depends on such factors as compressibility and viscosity of the fluids, the porosity and permeability of the rock, geometry of the whole system and withdrawal rate or pressure drop. With these factors known, it is theoretically possible to calculate the pressure-flow behavior of the system, However, in practice true solutions to the problem are next to impossible due to complexity of reservoir systems. A number of approximate methods of solution have been developed based on various simplifying assumptions. One frequent assumption is that the water-oil contact can be located and equations defining oil reservoir and aquifer solved separately by assigning values to various parameters so as to match past pressure and production history. Usually the properties of the aquifer are not known, since few, if any, wells are drilled through the aquifer; but the water influx and the pressure at the reservoir boundary are known over some time. If it can now be assumed that the aquifer is circular, pie-shaped, or linear, and that it has uniform properties, it is possible to fit a theoretical dimensionless curve to the past aquifer performance history and therefore calculate future water drive. These theoretical curves are available in the literature,1,2,3 Of course, the assumption of uniform aquifer properties is almost always somewhat in error. Fortunately, however, moderate variations from the average have little effect on behavior of the system; and, hence, the best fit of a theoretical curve is frequently satisfactory. In other cases variations in aquifer properties and geometry are large enough that none of the available theoretical curves will give an acceptable fit of the data. In these cases, methods of fitting the data with various electrical analyzers have been developed.1,2' It is the purpose of this paper to present an approach

By T. K. Perkins, L. E. Bartlett

Values of rock surface energy (i. e. energy required to form a unit area of new surface) are useful for interpretation of drilling and fracturing phenonema This paper describes the adaptation of a cleavage technique (splitting of large rocks under controlled conditions) for measuring surface energies. Values of surface energy obtained at room temperature and at low crack extension rates are reported for a variety of sandstones and limestones. Elastic moduli, tensile strengths and compressive strengths of the samples are also reported. The values of Young's moduli (measured in bending) are relatively low, probably because of the low stress levels imposed. Measured values of surface energy me generally higher than anticipated and probably reflect inelastic yielding near the extending tip of the crack. INTRODUCTION The technology of petroleum production includes many operations such as drilling and fracturing, which involve rock breakage. One of the fundamental properties of a rock needed to fully explain breaking and fracturing phenomena is the surface energy (i. e. the energy required to form a unit area of new surface). In this paper we describe an adaptation of the cleavage method (splitting of the rock under controlled conditions) to utilize large rock samples. Application of the theory of elastic cantilever beams and the Griffith concept of mechanical stability permits the calculation of surface energies from simple measurements. Values of surface energies obtained at room temperature on dry samples and at low crack extension rates are reported for a variety of sandstones and limestones. Measured values of Young's Moduli, compressive strengths and tensile strengths are also reported for the rocks. THE EXPERIMENTAL METHOD Equipment designed to cleave large rock samples is shown on Fig. 1. This equipment consists essentially of two steel blocks cemented to a rock sample and forced apart with a steel rod. A ball joint and a pivot prevent bending stresses in the rod. A differential thread arrangement permits the rod to advance only 0.012 in./revolution of the driving nut. The axial force in the rod is measured with type FA-50-12 Baldwin Lima Hamilton strain gauges and a BLH type 20 strain gauge reader. Separation of the blocks is measured with a Starrett 656617 dial gauge. The massive nature and design of the equipment minimize the amount of elastic energy stored in the steel parts as the stress in the rod is increased. This will permit the rock to split under controlled conditions as will become clear later. Fig. 2 shows a sketch of the rock sample during cleavage. Guide slots are cut longitudinally along the rock sample as suggested by Berry.1 An initiating slot is also cut across the top of the block to a depth of 1 to 2 in. During cleavage, elastic energy is stored in the flexed beams of rock of length C. If the length of steel rod L is held constant and the crack extends in length, then the elastic potential energy stored in the rock decreases. However, surface energy must be supplied as the crack extends in length. Griffith6 has proposed the following criteria for mechanical stability. The crack will extend as long as the potential energy available is greater than the surface energy required. However, the crack will stop extending when the change in potential energy with change in crack length is just equal to the change in surface energy with change in crack length. See Eq. 1.

By L. L. Handy

The development of a theory for miscible liquid displacement requires evaluation of the variables which affect growth of the mixing zone between solvent and displaced oil. Factors which appear to be important are individual fluid viscosities, viscosity ratios, flood rate, fluid densities, flow characteristics of the porous medium and molecular diffusion coefficients of the fluid components. The primary purpose of this paper is to evaluate diffusion effects. Theoretical treatments to date have been limited to floods for which the viscosity ratio is one. Two principal theories have been proposed. Von Rosenberg adapted for porous media a theory derived for capillary tubes by G. Taylor.' ,' In this theory molecular diffusion perpendicular to the direction of flow is a primary factor in sharpening the flood front. Slow floods give sharper fronts for a given distance traveled than fast floods. An alternative theory considers miscible liquid displacement as a statistical problem.3,4,5 Diffusion is not an important factor in this theory, but it leads to the same general type of equation as von Rosenberg's. Both theories predict S-shaped concentration profiles with the same dependence on distance traveled. The statistical or "dispersion" theory predicts rate independence, however. To supplement rate studies a direct measurement of a diffusion effect would be helpful in evaluating which of the two proposed mechanisms best describes miscible liquid displacement for one-to-one viscosity ratio systems. No quantitative theory has been proposed for floods in which a low-viscosity fluid displaces a high-viscosity fluid. It might be anticipated, however, that the extensive fingering observed in floods with adverse viscosity ratios would increase opportunities for an exchange of components between displaced and displacing liquids by a diffusional process. Even if molecular diffusion were not an important mixing mechanism for one-to-one viscosity ratio systems, it could be significant in systems with adverse viscosity ratios. METHOD FOR EVALUATING DIFFUSIONAL MIXING In solvent flooding the displacing liquid can become mixed with the displaced liquid by a number of processes; in general, either mechanical or diffusional in character. To evaluate the contribution of diffusion, a method was sought that would distinguish diffusional mixing from mechanical mixing. Suppose two soluble tracers are added to the displacing liquid and that one tracer has a diffusion coefficient much greater than the other. Then, if diffusion is important during miscible liquid displacement and if diffusional transport is primarily transverse to the direction of flow, substances with high diffusion coefficients will have sharper concentration profiles than those with comparatively low molecular diffusion coefficients. In order to have material balance the profiles for the two tracers must intersect. For one-to-one viscosity ratios the intersection would occur at the 50 per cent concentration point. It is generally assumed that diffusion in the direction of flow is negligible. This assumption appears reasonable because flow velocities are ordinarily orders of magnitude greater than diffusional velocities. If this is not a valid assumption, however, the shapes of the fronts will be further degraded by longitudinal diffusion. The effects of rate and molecular diffusion coefficients are the opposite for longitudinal diffusion to those for lateral diffusion. Higher rates or lower diffusion coefficients would tend to give sharper fronts. The double tracer method offers several unique advantages in evaluating effect of molecular diffusion. First, the method gives a direct measure of the effect of diffusion on mixing zone size. Second, the effect of a difference in diffusion coefficients is determined from a single experiment. Comparison of several experiments at different rates is not required. Third, the method is applicable for any viscosity ratio. Fourth, the method is unaffected by density differences between the two fluids which might result in rate dependent gravity effects. EXPERIMENTAL PROCEDURE Suitable tracers are substances whose diffusion coefficients differ as widely as possible, but whose adsorptions on rock surfaces is negligible. The two tracers selected were methanol and sucrose in water solutions. The approximate diffusion coefficients for methanol and sucrose, as given in the International Critical Tables, are 1.3 X 10-5 and 0.3 X 10-5 m2/sec, respectively. Concentrations in the water solutions were 200 gm/liter of solution for both methanol and sucrose. These high

By R. R. Harvey

Recent evidence has indicated that the films forming at crude oil-water interfaces are the result of naturally occurring su,rface-active components in the crude oil rather than oxidation products resulting from exposure of the oil to air. Since the presence of such films may strongly affect the displacement of oil from porous media, a knowledge of their composition and formatiorz is of fundamental importance. In this study, without the assumption of a specific mechanism, film formation is considered to involve the transfer of a weakly surface-active species in solution in the oil to the oil-water interface. The film gas law and the Gibbs adsorption equation may thus he applied to the system. The rigid films which form upon decrease of the oil-water interfacial area are considered a result of the very slow rate of re-solution of the film material into the oil. When the shrinkage of the interface is very rapid, the amount of film material leaving the interface is negligible, and the system may then be treated in the same manner as an insoluble film. The combination of these treatments leads to a new concept of the film factor in which it may he quantitatively related to the concentration of film-forming material at the interface. Suitable experiments can then yield the following information. 1. Original concentration of surface-active species in the oil. 2. The change of concentration with time. 3. Concentration of surface active species at the interface. 4. The change of this concentration with time. 5. Interfacial tension of the system at any time and its variation with time. 6. Equilibriunz concentrations for a given set of conditions. 7. Concentration requirements for solid film formation. 8. Molecular diameter of the adrorbed species. These calculations may be made solely from data taken from measurements on the original oil. While using the pendant drop method of interfacial tension measurement of crude oil-water systems, Bartell and Niederhauser' remarked on the thin, transparent, pliable, solid films or membranes which were often formed at the interface between crude oil and water. As a result of their observations of these films, Bar-tell and Niederhauser concluded: 1. The films are formed by material from the oil adsorbed at the interface rather than by a reaction at the interface. 2. The adsorption of the film-forming material is a rather slow process. 3. There are good indications that the process is reversible; i.e., the film constituents may re-dissolve into the oil. 4. The film will spontaneously spread over any new surface formed. 5. The process of adsorption at the interface is accompanied by a decrease in the interfacial tension of the system. Although they had proved that film formation was not due to a reaction at the interface, Bartell and Niederhauser were unable to demonstrate conclusively whether these film-forming constituents were present in the original reservoir oil or whether they were formed as oxidation products after exposure of the oil to air. Recent evidence seems to point to their presence in the original oil. Although the true nature of the film-forming constituents is still uncertain, these conditions would be fulfilled by the assumption of the presence of a slightly soluble surface-active component in the original oil. This component would migrate to a freshly formed interface because it would thereby lower the free energy of the system. The component at the surface constitutes an expanded monolayer which would become compressed, upon retraction of drop, due to the slow rate of re-solution of the film components into the bulk of the oil. It is possible, by involved extraction processes, to separate or at least concentrate this film-forming fraction. However, using the assumptions made above, the technique of monolayer film study may be applied to reveal much of the nature of the adsorbing material without resorting to an actual separation of the film-forming constituents.

By G. D. Boozer, W. H. Somerton

RESEARCH CORP., LA HABRA, CALIF. Thermal diffusivities of some typical sedimentary rocks have been measured by a rapid unsteady-state technique. Thermal data including diffusivity and conductivity are presented for the temperature range of 200" to 1,800°F. Diffusivities calculated from steady-state conductivity rneasurernents compared favorably with unsteady-state data at the lower temperatures. Unsteady-state data are less reliable at higher temperatures; reported values may be as much as 20 per cent high at temperatures in excess of 1,500°F. This is caused partly by the simplified method used in calculating diffusivity from experimental data and partly by the many thermal reactions which occur at the higher temperatures. The marked decrease in thermal diffusivity and conductivity at elevated temperatures indicates the importance of considering the temperature dependence of thermal properties in subsurface heat-transfer calculntions. Knowledge of the effects of liquid saturation, overburden and pore pressures, and heats of reaction of reservoir materials will be required before complete subsurface thermal analyses may be made. INTRODUCTION Heat-transfer calculations are becoming increasingly important in the analysis of in situ combustion and underground nuclear explosions. Thermal data for natural rock materials at elevated temperatures are generally unavailable. Calculations at high temperatures may be grossly inaccurate if temperature-dependent thermal properties of rocks are not considered. It is known that the heat capacities of rocks increase with temperature. Heat capacities may be calculated from chemical or mineralogica1 analyses of the rock system. Thermal bulk expansion of rocks is believed to be small, although few data are available. Thermal conductivities of most natural rock materials decrease with increased temperature. Thermal diffusivity, defined as the thermal conductivity divided by the product of the heat capacity and the bulk density of the material, is used in unsteadv-state heat-transfer calculations. It is apparent that thermal diffusivity decreases markedly with increased temperature. A new method, which gives a nearly continuous record of thermal diffusivity within the temperature range of 200° to 1,800°F and requires less than 1.5 hours for the run,' was used to measure the diffusivity of several sedimentary rocks. Thermal conductivities were calculated from measured bulk densities, heat capacities calculated from chemical analysis and the measured diffusivity values. Thermal conductivities were measured at lower temperatures by a steady-state technique1 to test the agreement of data obtained by the two methods.

By J. M. Smith, D. Kunii

Effective thermal conductivities of sandstones filled with stagnant fluids were measured using a steady-state technique. Data were obtained for seven sandstone samples, taken from four different locations and ranging in permeability from 18 to 590 md. The measurements with gases (helium, nitrogen, air and carbon dioxide) covered a pressure range from 0.039 psia to 400 psig. Data were taken for four liquids — n-heptane, methyl alcohol, 79.8 weight per cent glycerol-water solution and pure water at atmospheric pressure. The experimental results were used to evaluate the theoretical equations for predicting stagnant conductivities developed earlier 4 The low-pressure measurements permitted evaluation of the consolidation parameter bpDp/ks (necessary to utilize the theory) for the various types of sandstones. Using these characteristic values, the tbeoretical equations correlated well with the experimental conductivity data for the several fluids and rock samples. INTRODUCTION An aspect of heat transfer in solid-fluid systems of considerable current interest is the effective thermal conductivity of porous media. The stimulus for study of the subject arises from the need for sound procedures for designing thermal methods of petroleum production. The general system occurs when there exists a flow of fluid through the pores of the solid material. However, a logical starting point in developing a theory for predicting the effective thermal conductivity in the general system is to attack the special case when the porous solid is filled with stagnant fluid. Since the flow rates anticipated in thermal production processes are very low, such stagnant conductivities k,O are also of practical significance. Kunii and Smith4 recently reviewed the existing data for stagnant conductivities and proposed a theory for heat transfer in porous rocks. This leads to equations for predicting the stagnant conductivity as a function of the properties of the fluid and solid phases. The existing experimental information1-3,5,7,8 covered a range of porosities and agreed well with the theory. However, the available data were insufficient to examine critically the effects of two important properties on koe, — the thermal conductivity of the fluid and the pressure. The objectives of the present study are twofold: (1) to measure the effect on koe of (a) the pressure for gases, and (L) the thermal conductivity of the fluid kg for bloth gases and liquids; and (2) to determine for rocks of different porosities and permeabilities the apparent solid conductivity ks and the consolidation parameter bpDp/ks required in applying the theoletical equations. Measurements extrapolated to zero pressure are helpful in evaluating ks and the consolidation parameter. However, the extrapolations must be based upon data obtained at low pressures because in small pores the molecular conductivity of the fluid may be unusually depressed. This occurs when the mean free path of the molecules is of the same order of magnitude as the pore diameter, that is, in small-diameter pores at low pressures. The phenomenon has been studied recently by Schotte6in beds of unconsolidated fine particles. With particles in the size range of 45 to 200 microns, a significant reduction of kg with pressure was found even at absolute pressures greater than 1 atm. In the sandstones employed in the present investigation, pore diameters lower than those experienced by Schotte are likely, so that the effect would presumably be greater. SCOPE OF MEASUREMENTS Stagnant conductivities were measured for four types of sandstones using a total of seven samples. The characteristics of the samples, including porosity and permeability values, are given in Table