Search Documents

By F. Brons, W. C. Miller

Pressure information for use in material-balance calculations is obtained, where possible, from pressure build-up surveys in shut-in wells. Using proper extrapolation methods, static pressures are obtained which, by averaging, give the field static pressure. Often, particularly in fields with a large number of wells, only spot pressure readings are available. A graphical method is presented here which enables the determination of the corrections which must be applied to the spot readings to give an estinzate of the static pressure. Although the static pressures thus obtained are not so reliable as those from build-up surveys, they are more nearly correct than the spot readings themselves for use in material-balance calculations. INTRODUCTION The average static pressure in a bounded reservoir is determined from a number of individual-well pressure measurements. Where pressure build-up curves are available, extrapolation methods such as described in Ref. 3 may be used to determine the static pressure in the drainage area of each individual well. The possibility of using pressure build-up curves presupposes that each well being surveyed is closed in during a certain period and that, during this time, a pressure recorder is left suspended in the hole. For surveys in fields with a large number of wells, this is usually not feasible for economic and other reasons. Either a large number of pressure recorders would have to be used simultaneously or, using a limited number of recorders, the time interval needed to complete a field-wide survey would be too long. The data obtained for this incomplete coverage are insufficient to permit reliable averaging of the few calculated static pressures over the entire reservoir. In many such cases another type of pressure information is available, i.e., spot pressure readings. During this type of pressure survey all wells are closed-in simultaneously and, after a reasonable shut-in time (e.g., three days), pressure readings are taken by running the recorder in and out of each successive well. These surveys, which are less expensive, are often taken at regular intervals, say once or twice a year. In fact, such surveys are required by the state regulatory bodies in many portions of the united States. These spot readings are normally lower than the static (fully built-up) pressure due to the relatively short shut-in time. No direct means of extrapolating such a reading to the static pressure is available because no pressure build-up record is obtainable. Consequently, these data are inferior to pressure build-up information as a basis for estimating the stabilized reservoir pressure in the drainage area of a specific well. For the common case where abundant spot readings but only a few pressure build-up surveys are available, our studies suggest that it may be possible to design a program which will combine the two types of information and provide an estimate of the average reservoir pressure that is more reliable than that obtainable solely from either set of data. in time) gives a point marked p. The actual static pressure at the time the well was shut-in is marked by p.. The illustration shows that, even if a build-up curve is available and can be extrapolated, a correction

By J. D. Testerman

A statistical technique to identify and describe naturally occurring zones in a reservoir and to correlate these zones from well to well is described. The technique is particularly useful in describing a reservoir where cross flow between adjacent strata is important in determining reservoir behavior. Although it has been used primarily for permeability zonation, the technique is general and can be used to correlate any reservoir property or related data, such as the information contained in well logs. INTRODUCTION One of the first problems encountered by the reservoir engineer in predicting or interpreting fluid displacement behavior during secondary recovery processes is that of organizing and using the large amount of data available from core analysis. Permeabilities pose particular problems in organization because they usually vary by more than an order of magnitude between different strata. Due to the sheer volume, it is almost always necessary to group data and to use an average value to represent a number of measurements. Perhaps the most common method now used to group permeability data is the capacity-fraction technique, which ranks permeabilities in order of magnitude, regardless of the physical location of the permeabilities within the reservoir. The cumulative per cent capacity is plotted against cumulative per cent thickness. This plot is divided into an arbitrary number of zones, generally of equal thickness. Five zones (or averaged groups of data) usually are obtained, each of which is treated as homogeneous in subsequent calculations. The division so obtained has no physical meaning; strata in the same zone, calculation-wise, are usually not adjacent in the reservoir. Reservoir engineering techniques being developed will handle crossflow that occurs between adjacent communicating reservoir strata because of imbibition and gravity segregation. Since crossflow occurs between physically adjacent layers within the reservoir, a new zonation technique recognizing the actual location of strata within the reservoir is necessary. Similarly, the recognition of natural zones is important for predictions of oil recovery by processes involving diffusion. One such process is miscible displacement, where predictions of lateral diffusion within the reservoir must recognize the actual location of the invaded zones in relation to the rest of the formation. Natural zones must also be adequately recognized to account for heat transfer within the reservoir during thermal exploitation. Because of the complexity of the problem, statistics appear to offer the only practical hope of dividing a reservoir into physically-meaningful natural zones. This paper presents a statistical technique for identifying these natural zones and for ascertaining which ones are likely to be continuous between adjacent wells. The zones defined have minimum variation of permeability internally and a maximum variation between zones. The technique is general and can thus be applied to reservoir properties other than permeability. The method will guide the reservoir engineer in estimating which zones are likely to be continuous between wells. However, a statistical correlation based on permeabilities in two different wells is no guarantee that the zones so defined are, in fact, continuous. Rather, the assumption of continuity must be consistent with geological data concerning the depositional environment, as well as justified on the basis of engineering judgment in combination with statistics, just as judgment is required with conventional zonation methods. CALCULATION PROCEDURE The reservoir zonation technique is a two-step operation. The steps are individually described, and a sample calculation is presented in the Appendix. ZONATION OF INDIVIDUAL WELLS First, the set of permeability data at a single well is zone, into Zones. These zones are selected so that variation is minimized within the zones and maximized between the zones. The equations4,6 used to zone the data are where B = the variance between zones, , = the number of zones, i = the summation index for the number of zones, j = the Sumation index for the number of data within the zone, mi = the number of data in the ith zone, k,. = the mean of the permeability data in the ith zone, k . = the over-all mean of the data in the well, W = the pooled variance within zones, N = the total number of

By N. T. Cotman, C. H. Pickering, P. B. Crawford

A beat-conduction model has been developed to study the flow of fluids in a stratified oil reservoir which is being subjected to unsteady-state depletion. To simulate stratification, plates of different metals were joined together so that cross flow could occur between strata. A description of the model, plate construction and necessary instruments is included. The data indicate that stratified oil reservoirs experiencing cross flow may be treated as a single, homogeneous, producing sand having properties intermediate between those of the layers making up the system. A technique for averaging core-analysis data to arrive at the proper average value of rock permeability and porosity is presented INTRODUCTION Although oil has been produced from reservoirs for 100 years and reservoir injection processes have been undertaken for 40 years, the industry is still inadequately informed on the problem of reservoir inhomogeneities. Stratification is perhaps the best-known type of reservoir heterogeneity. Urenl discussed the possible significance of this nonuniformity in his 1927 paper on the theoretical aspects of water flooding. Since this time the problem of vertical permeability variations has been the subject of much discussion and numerous papers. Levorsen2 simply describes the origin of layered beds in his Geology of Petroleum. Some beds are field-wide and may be traced from well to well, while others apparently pinch-out and become discontinuous. In addition, the physical characteristics of a layer may vary from place to place. For these reasons, a stratified system is not easily defined. In calculations involving oil reservoirs in which permeabi1ities vary, methods presently used contain many simplifying assumptions. The reservoir frequently is considered as consisting of permeable layers, each of uniform vertical and lateral permeability and each behaving as if separated from adjacent layers by thin impermeable layers. Other assumptions frequently made in reservoir calculations are steady-state flow conditions, unit mobility ratio, linear geometry, average or homogeneous porosity, and negligible capillary-pressure and gravity effects. Muskat3 has. developed solutions to the stratification problem for mobility ratios other than unity for both linear and exponential permeability distributions. He assumes lateral uniformity and continuity of all productive strata and does not consider porosity. This study will take into account transient effects, cross flow and varying porosity, as well as the variations in permeability, by means of a heat analogy. Muskat5 has shown the similarity of the equations which describe the flow of fluids in a porous medium to the flow of heat through a solid. Landrum, et al,6 made an analogy between the flow of fluid in an unsymmetrical reservoir and the conduction of heat in a similarly shaped metal plate. This analogy is summarized briefly as follows. If fluid is produced from an oil reservoir above the bubble point, the resulting pressure distribution will be identical to the temperature which would result in a metal or heat-conducting plate if heat were similarly removed. This means that a study of heat flow in stratified thermal models will behave like fluid flow in stratified oil reservoirs. For mathematical correspondencies, set- Ref. 6. The fluid mobility k/p is proportional to the thermal conductivity, and the porosity-times-compressibility product is proportional to the density times specific heat of the thermal plate. EXPERIMENTAL EQUIPMENT AND PROCEDURE The stratified reservoir model used in this study is illustrated schematically in Fig. 1. A Leeds and Northrup Speedomax Type G temperature recorder

By M. R. Tek, M. L. Katz

The flow of fluids in stratified porous systems, in which adjacent layers of the system are in communication over part or all of their common interface, in many instances involves crossflow between layers. A theoretical study has been made of the unsteady-state flow of slightly compressible fluids during depletion of bounded stratified porous systems. Consideration has been restricted to single-phase flow, with negligible capillary and gravity effects. Analytical expressions have been derived for pressure distribution and for (luxes at the producing face in two-layer systems with complete fluid communication between layers. These expressions occur in the form of double Fourier and Fourier-Bessel series. The method of extending the solution to any number of layers is demonstrated by deriving the equation for pressure distribution in a three-layer system. On the basis of the exact results for two-layer systems, certain generalizations regarding the behavior of stratified systems in general are presented. The performance of stratified systems, in terms of cumulative flux as a function of time, is shown to lie at all times between bounds established from single-layer theory. The upper bound is given by treating the system as a single layer with arithmetic average physical properties. The corresponding lower bound is the summation of the fluxes from each layer as treated individually. Variation in performance of stratified systems between the upper and lower bounds is determined by the extent to which crossflow occurs in the system. The principal parameters affecting the amount of crossflow are the vertical permeability of the system, the ratio of system thickness to "radius of drainage" and, for radial systems, the ratio of system thickness to interior radius. At very early times, the system performance is near the lower bound. As time and the radius of drainage increase, the per- formance tends toward the upper bound. INTRODUCTION One of the principal problems confronting the reservoir engineer today is the incorporation of the effect of reservoir heterogeneity into predictions of reservoir performance. The first part of this problem is the characterization of the actual reservoir heterogeneity in a mathematically useful manner.1, 2 It is then necessary to calculate the effect of this heterogeneity on the system performance. It is frequently possible to estimate this effect by comparison with the known behavior of an idealized system. One such idealized system is the stratified or "layer-cake" model. Many studies have been made of the behavior of stratified systems, from both a theoretical and experimental standpoint. The great majority of the theoretical studies have been concerned with systems of separated layers, with no fluid crossflow. Among the more noteworthy of such theoretical papers is that recently published by Lefkovits, et al,3 in which a comprehensive analysis was presented for the behavior of bounded stratified systems which are unconnected except at the wellbore. Three recent references have dealt with stratified systems including the effect of crossflow. Jacquard4 has presented mathematical derivations pertaining to depletion of two layer systems producing at constant rate. A doctoral dissertation by Katz5 contains analytical solutions and tables of results for pressure distribution and fluid flux during constant-terminal-pressure depletion of two-layer and three-layer systems with crossflow, together with tables of eigenvalues necessary in the computations. The dissertation also presents results from experimental studies of stratified systems using a heat-transfer analog model. This paper is largely based on the material in this dissertation. A third paper has been presented by Russell and Prats6 since completion of the work resented herein, containing mathematical derivations for and discussion of the behavior

By H. N. Hall

Various factors must be considered in an engineering evaluation of gravity-drainage reservoirs. Among these are: (1) the effect of producing rate on total oil recovery; (2) the effect upon well productivity and ultimate recovery of the pressure level maintained during the producing life of the reservoir; (3) the economic advantage of full or partial pressure maintenance; and (4) estimate of the rate of gas production and injection and the possible purchase of gas under conditions of full pressure maintenance to ascertain compressor facilities needed. All of these factors can be evaluated only when a reliable method is employed for determining reservoir performance in gravity-drainage reservoirs. The purpose of this paper is to present a general method for calculating the performance of a gravity-drainage reservoir. This method is applicable for conditions of complete pressure maintenance, partial pressure maintenance and normal pressure depletion. Provisions are made to take into account variations throughout the reservoir of reservoir configuration, changes in permeability and fluid composition. Based on the method presented in this paper, an IBM 650 computer program has been developed. The past performance of an actual gravity-drainage reservoir producing under conditions of declining pressure and no gas injection was duplicated using this program. INTRODUCTION In tilted reservoirs the production of oil is influenced by drainage of oil from upstructure to downstructure locations. When this downstructure drainage of oil is sufficient to cause effective segregation of the gas and oil in a reservoir, the reservoir is usually classified as a segregation drive or gravity-drainage reservoir. (Discussion will be restricted to gravity-drainage reservoirs which have no encroachment of edge water.) The important feature in gravity-drainage reservoirs is the density difference between reservoir oil and gas. These phases tend to segregate in the reservoir with the result that in the gas cap the oil saturation is maintained at a higher level by drainage of oil from the gas-cap area. Oil can be produced from the oil zone at a low gas-oil ratio and reservoir energy is thereby conserved. The standard material balance in not adequate for predicting gravity-dramage reservoir performance because it does not take into account the difference in saturation above and below the gas-oil contact. Several authors'.' have presented methods for calculating the performance of gravity-drainage reservoirs in which reservoir pressure is maintained constant by gas injection into the gas cap. Using some simplifying assumptions, these methods can be employed with a desk calculator to give acceptable results. The problem of predicting the performance of gravity-drainage reservoirs under the conditions of declining reservoir pressure is many time more complex than that of constant pressure. fierefore, attempis to develop a method suitable for desk calculation have required excessively simplified assumptions. In the past several years, highspeed digital computers have become more widely available for reservoir engineering problems. These corn puters are well suited to problems such as the prediction of the performance of gravity-drainage reservoirs with pressure decline. Many of the simplifying assumptions necessary for hand computation can be eliminated so that a realistic approach to the gravity-drainage process can be made. CONCEPTUAL PICTURE OF OIL MOVEMENT IN GRAVITY-DRAINAGE RESERVOIRS Before attempting to develop an analytical treatment for conditions occuring in a gravity-drainage reservoir, a concept should be formed concerning the movement of fluids in the reservoir as oil is produced. A review of the literature'.' shows that it is customary to classify gravity-drainage operations into two categories—(1) with complete pressure maintenance, and (2) with declining pressure. The same line of reasoning will be followed in presenting the concept of the movement of fluids in the reservoir because it is easier to visualize the movement of fluids under conditions of complete pressure maintenance. After discussing complete pressure maintenance, an analogy will be made between that and the case of declining pressure. It should be kept in mind throughout that the final aim for the problem of solving gravity-drainage performance with digital computers will be to develop a general program for any kind of gravity-drainage reservoir. COMPLETE PRESSURE MAINTENANCE One feature which is generally common in gravity-drainage reservoirs is a gas cap located at the top of the structure. This is shown in Fig. 1(a). Fig. l(b) shows oil saturations that might occur through the reservoir. In the gas cap, oil saturation is lower than

By D. D. Smith, D. P. Squier, E. L. Dougherty

A system of differential equations describing the temperature behavior of fluid injected at constant surface temperature in a well is derived and .solved analytically. A formula for the fluid temperature at any time and depth is given, as well us a special formula valid for very large times. These formulas are used to calculate temperatures for several typical cases. The results indicate that, initially, the temperature of the water entering the formation is considerably lower than the injection temperature. This condition lasts for only a short period— less than three days for most cases of practical interest. Following this highly transient period, during which the temperature of the fluid entering the formation builds up to about 50 to 75 per cent of the injection temperature. the system enters a quasi-steady state in which the temperature changes are very slow. After severl years, the bottom-hole temperature will still be 50" to 100°F lower than the injection temperature, hilt the heat losses may he tolerable. INTRODUCTION Predicting the behavior of a hot-water flood requires that the temperature of the water entering the injection interval be estimated. This report describes the development and solution of a system of equations which describes the temperature behavior of the injected water in the wellbore with certain simplifying assumptions. The only previous means known to the authors for describing such a process is that of Moss and White.' Their results appear to be close to those obtained by our method in the practical cases which were compared; this agreement is largely due to the fact that in our method temperature soon approaches a quasi-steady state, as was assumed in their method throughout. However, our model covers all times, is continuous (whereas the Moss-White model depends on breaking the depth into discrete intervals) and. we feel. more closely describes the physical problem. FORMULATION OF THE PROBLEM PHYSICAL SYSTEM AND ASSUMPTIONS The injection procedure consists of pumping water at a fixed surface temperature T., down an infinitely long cylindrical well or tubing of inner radius Any material exterior to the water column such as mud, casing, or cement is regarded as part of the formation. The general behavior of the system may be described qualitatively as follows. When the hot water is first introduced into the system, the temperature difference between the formation and the water is large, resulting in a high rate of heat transfer. As a result, the temperature adjacent to the wellbore rises very quickly. Because the segment of the formation adjacent to the wellbore largely controls the heat transfer rate, the heat transfer rate will become relatively constant when this portion has reached a temperature close to that of the water opposite it. The temperature of the water and formation then increase very slowly with time. The length of the initial highly transient period and the temperature of the water at its conclusion will be functions of depth, injection rate, injection-string radius, surface injection temperature and the physical properties associated with the water-formation system. The following additional assumptions were made. 1. There is no heat transfer by radiation in the system. 2. There is no heat transfer by conduction in the vertical direction in either the injection stream or the formation. 3. The linear volumetric and mass flow rate of the water is constant throughout the injection stream. 4. No horizontal temperature gradient exists in the injection stream. 5. The product of density and heat capacity is constant for both the water and the formation, and the formation thermal conductivity is constant. 6. Initially, both the water in the wellbore and the reservoir are at a temperature given by the (constant) ambient surface temperature plus the product of depth and geothermal gradient (assumed constant). At large distances for the wellbore (r m), the formation will remain at this temperature. 7. The water temperature and the formation temperature at r — r,, are equal for all depths D. DERIVATION OF EQUATIONS The differential equation satisfied by the fluid temperature T,(D, t), which is obtained by writing a heat balance on a cylindrical differential of volume dV of the injection string between the depths D and D i dD, is

By W. D. Kruger

Methods for analyzing flooding or cycling projects by means of two-dimensional flow calculations are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The accuracy of such analyses depends on the basic reservoir data employed, especially the absolute permeability distribution in the reservoir; however, absolute permeability data normally are determined from well flow tests and/or core measurements and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. A calculation procedure is presented for determining the areal permeability distribution in the reservoir. Use of the procedure allows verification of the basic reservoir data through the matching of past reservoir conditions, thereby providing proper data for reliable predictions of future operations. The calculation procedure is a numerical technique based on a mathematical model of the reservoir and is designed for solution by means of an electronic digital computer. Results from field application of this procedure are presented. Permeability distributions from measured data are compared to those calculated to indicate the need for the analysis. Comparisons of field-measured and calculated pressure distributions are given to establish the validity of the calculation procedure. The conclusions of this paper are (I) use of the calculation provides data that allow a match of past reservoir conditions and thus reliable prediction of future operations; and (2) the procedure should, where applicable, be used to provide adequate reservoir data for use in two-dimensional flow calculations and related reservoir analyses. INTRODUCTION One means of providing more precise engineering control over reservoir operations during all stages of depletion is to develop some type of model to simulate the reservoir, using information developed from the behavior of the prototype to predict the performance of the actual reservoir. The use of mathematical rather than physical models for this purpose has been gaining in importance recently because of the application and availability of electronic computers. Mathematical models avoid instrumentation difficulties and allow easier handling of large-scale problems. At present, analyses employing two-dimensional models are common, and thought is being given to the use of models in three dimensions. The use of multi-dimensional models is necessary to show not only what is occurring in the reservoir, but also where in the reservoir the phenomena of interest are taking place. Methods for analyzing flooding or cycling projects by means of two-dimensional numerical models are presented in the literature. The use of these methods allows the determination of optimum operating conditions for such projects and provides information on pattern efficiencies, relative flow rates, pressure distributions and floodout histories. The applicability of numerical models and the results of using such analyses have been presented through field examples. The references indicated here discuss fully the data necessary to build the mathematical reservoir model and the way the model is used to study or predict the future performance of the reservoir. The accuracy of the analysis of such projects depends on the basic reservoir data employed in the construction of. the model. In order to be sure the model represents the reservoir, it must be used to match some known set of reservoir conditions, i.e., to match past reservoir performance. The accuracy of the model is directly dependent upon the absolute permeability distribution in the reservoir; however, absolute permeability data are normally determined from well flow tests and/or core measurements, and are so limited in quantity and quality that they do not adequately describe the reservoir flow properties. There is always the need for better data and for verification of existing data. The purpose of this paper is to show the use of the mathematical model to match past performance in flooding or cycling projects, the use of numerical techniques to adjust the basic absolute permeability data if necessary so that a match is obtained, and in this way to generate the effective areal permeability distribution of the reservoir. A review of the basic concepts of the numerical model will be given, followed by a discussion of the data-adjusting technique for matching past performance. TWO field examples will be cited to show

By R. A. Greenkorn

This project originated in a practical problem—we needed five tracers that could be used together to locate flow paths in a pilot flood. While tracers for subsurface liquids have been used since the turn of the century,1-15 none of those reported in the literature appeared to be either consistent or quantitative enough for our purposes. Most were used in field systems without controlled experiments to determine the accuracy and precision of analysis, and many were tracers that could not be used collectively. The ideal tracer, of course, would follow the fiuid of interest exactly, traveling at the same velocity as the fluid front. But the ideal is impractical to attain because adsorption-desorption effects cause the tracer to lag behind the front; these effects, plus diffusion-dispersion effects, cause the tracer front to spread more than the fluid front. Thus, our objective was not to locate a tracer that would be ideal for all circumstances but rather, to find one that would approximately follow the fluid, or one that under controlled conditions could be corrected to calculate the movement of the fluid front. We considered tracers satisfactory—(1) if they were easy to analyze; (2) if their breakthrough-elution curves were not too different from those for the chloride ion, a tracer believed to follow the fluid front closely; and (3) if we could calculate from the curves a material balance, at 1.25-PV (pore volume) injected, within 5 per cent of that calculated from chloride curves. Of a possible 35 materials, we selected 13 tracers that could be quickly and easily identified and whose analysis was claimed to be accurate within 5 per cent. Only one of these, tritiated water, was a radioactive tracer. Radioactive tracers are easy to detect even in small quantities, but they require special handling and special equipment. Also, those that can be used together are limited because special equipment is required to separate emissions from the various tracers. Two of the original 13 tracers were eliminated in static tests to determine how accurately and precisely they could be analyzed, and to check on gross adsorption. The remaining 11 were flowed through a 9-ft linear sandstone model, and breakthrough-elution curves were obtained. Finally, three tracers were field-tested as breakthrough tracers. These tests are described in the following sections. The 13 tracers considered in these experiments were EDTA (ethylene diamine tetra acetic acid), fluo-rescein, picric acid, salicylic acid and ammonium, boron (as borate), bromide, dichromate, iodide, nitrate and thiocyanate ions, plus chloride ion and tritiated water. All but the chloride ion and tritiated water were subjected to static sand tests to eliminate the tracers that could not be analyzed quickly and accurately and to eliminate those that showed excessive adsorption. All but two of the tracers, EDTA and salicylic acid, qualified for flow tests on this basis. The procedures used in the static tests (see Appendix A for analytical details) were as follows. Each tracer (in the amount shown in Table 1) was dissolved in 1,000 ml of water. Then, 800 ml of this solution was added to 500 gm of sand, and the remaining 200 ml was reserved as a control standard. The sand-tracer mixture was agitated once a day over a 10-week period, and during this period we analyzed five duplicate samples from this mixture, along with five duplicate samples from the control solution. From the control-sample analyses, we constructed a control chart on which the sand-tracer analysis results were plotted. The control chart, and its use to determine accuracy and precision of analytical results, as well as the amount of adsorption, is described in Appendix B. The results of the static tests are summarized in Table 1, which lists the concentration c (concentration at time t) divided by c, (initial concentration). If there was no adsorption, c/c, = 1. The acceptable deviation from this value varied from tracer to tracer, depending upon the limits of analytical precision established for the different tracers. The best results were obtained with boron, bromide and dichromate, with all values falling within limits of accuracy and precision. Ammonium, iodide, nitrate and picric acid also were satisfactory, although the ammonium results were erratic. Flourescein and thiocyanate were unsatisfactory ini-

By J. M. Smith, D. Kunii, P. Adivarahan

Effective thermal conductivities were measured for seven samples of porous rocks through which gases or aqueous salt solution were flowing, parallel and countercurrent to the flow of heat. The results are analyzed in terms of the contributions of the solid and fluid phases to the total energy transfer. The solid contribution for the sandstone samples studied was found to be a Junction of the porosity and an effective thermal conductivity for the solid phase. The fluid contribution was a function of the Peclet number. Based upon the equations presented, the effective conductivity of the system can be estimated from the properties of the fluid, the flow conditions, and the porosity and pore-size distribution of the rock. INTRODUCTION Heat-transfer characteristics of porous media containing fluids are needed for the design and evaluation of thermal methods of petroleum production. The study reported here considers the problem of countercurrent flow of fluid and heat in sandstone rocks, and is a continuation of previous studies concerned with parallel and radial flow of energy and fluid in both rocks517 and beds of fine particles.6,10 In addition to these earlier works, Asaad,1 Kimura,4 Preston,8 Zierfuss and vlietll and waddams9 have reported effective conductivities for systems filled with stagnant fluids. There appears to be no published data on rocks through which fluids are flowing. The purpose of this study was to present such information and, particularly, to investigate the effects of properties of the rock and the flow rate on the results. The rock samples were naturally occurring sandstones from different locations and were similar to petroleum-bearing formations. The properties of the seven samples showing the range in porosities and permeabilities are given in Table 1. The mean pore sizes were calculated from mercury penetration data by averaging the radius with respect to pore volume. The details of this method and comparison with other procedures for estimating 7 are given elsewhere. Four fluids, nitrogen, carbon dioxide, helium and 10 weight per cent sodium-chloride solution were used with each sample. For Bartlesville sandstone, BA-5, runs were made at mean pressures from 15 to 250 psig using nitrogen. The other measurements were carried out at a mean pressure of 100 psig. The maximum temperature range eacountered was 70° to 300°F. The flow rate of fluid ranged from about 1.0 to 130 Ib/(hr)(sq ft), and this corresponded to a modified Reynolds-number range of 5 x 10-4 to 1.0. The variation in Prandtl number was 0.70 to 4.1. EXPERIMENTAL WORK The equipment consisted of apparatus to supply a constant flow of gas, or salt solution, at constant pressure to the bottom of a vertical test section, and metering and flow-regulation apparatus after the test section. The flow of salt solution was

By S. W. Churchill

A digital computer was used to obtain an exact numerical solution for the transient behavior of the insulation and earth adjacent to an isothermal, submerged flat surface for a single set of parametric values. Comparison of the computed results with analytical solutions for limiting conditions revealed that a complete and general solution for all parametric values could be constructed from these limiting solutions. Complete and general solutions for insulated spheres and cylinders can also be constructed from limiting solutions. The procedure is illustrated in part for an insulated spherical tank. The heat flow was found to fall rapidly to a pseudo-steady state; after some time, it then decreased slowly to zero for a flat plate and cylinder, and to a low steady-state rate for a sphere. The accumulative heat flow during the initial falling-rate period may be a significant fraction of the heat flow during the entire first year. INTRODUCTION Underground storage of liquefied natural gas and liquefied petroleum gas has received considerable recent attention. The rate of heat flow from the earth to the storage cavity and the resulting temperature field in the earth are important factors in the technical and economic evaluation of such storage facilities. The objective of this paper is to indicate how complete solutions can be developed for the transient flux and temperature field for various geometrical configurations. The representative properties and dimensions, and the resulting parameters utilized in the illustrative results, are indicated in Table 1. The results presented herein are for dry earth. The effect of the latent heat of solidification of water in the soiI will be described in a subsequent paper. The temperature field in the insulation and earth is determined by energy balances, boundary conditions and initial conditions. The physical problem can be described mathematically with virtually no idealizations insofar as physical properties are known. It appears possible to solve the equations by numerical integration on a high-speed digital computer for any geometrical configuration and conditions. For complex situations, however, the computations are expensive and the results are highly specific. Analytical solutions have been developed for a few simple but important geometrical configurations and conditions, including one-dimensional heat transfer from (or to) earth at an initially uniform temperature to (or from) flat plates, spheres and circular cylinders. The solution for an insulated flat plate has also been derived but is in the form of a slowly converging series involving taklated functions. It was planned to use a computing machine to evaluate this series for a number of representative conditions. However, upon examination of the results of preliminary computations it was discovered that a complete, general and accurate solution could be developed by interpolation between several much simpler solutions for limiting conditions. This technique then was used to develop a a complete solution for an insulated sphere. The equations presented herein are derived in Carslaw and Jaegerl and other books on applied mathematics, or they are simple extensions of these previous results. Hence, all derivations are omitted. Emphasis is instead placed on the results and their

By N. Marusov, D. K. Lowe, J. C. Karp

This paper considers, from an engineering viewpoint, several factors involved in creating, designing and locating horizontal barriers for controlling water coning. This is an effort to consolidate new concepts with previous information so that a reasonable selection of barrier materials, dimensions and vertical position can be made. Coning theories previously developed are briefly reviewed and an effort i,s made to reduce the results of coning-theory calculations to a point where routine calculations can be made with a desk calculator. It is expected that these simplified calculations will give usable predictions of the amounts of improvement attainable with barriers of various dimensions. Apparatus and procedures used for testing the suitability of commercial cements are described and test results are presented. INTRODUCTION Just about as soon as the phenomenon of water coning was recognized as an oilfield problem, there were suggestions that the production of water coning could be controlled or completely suppressed by means of horizontal barriers. It was also suggested that natural barriers such as shale streaks were helpful in restricting bottom-water production. The implication was that wells not having the benefit of continuous shale streaks to suppress bottom waters should, in some way, be supplied with an artificial barrier. With regard to the selection of barrier materials, there are many references which give the merits of liquid barriers.',' These liquids include surface-active agents," recipitates' and emu1sions. There is no real need to review the merits of these materials again. This paper deals with the use of solid barriers which could be made from various cements. The two major restrictions placed on these barrier cements are that (1) they must be commercially available in large quantities, and (2) they must be applicable by means of standard tools and equipment with techniques that have been fully developed and are in general use at the present time. A schematic drawing of a horizontal barrier in a reservoir is shown in Fig. 1. Two maximum stable water cones are illustrated-—-one for the wellbore diameter and the other for a barrier. A maximum stable cone is one that has attained its maximum volume at the critical production rate and is on the threshold of producing water. The volume displaced by each water cone and the rise in water table are indications of the amount of water-free oil produced. Corrections for porosity, connate water and residual oil, of course, must be applied. The water-oil interface which defines a maximum stable cone is called a free surface. It is a surface over which the pressure is constant. This free surface has the shape of a limiting streamline at and below which there is no movement or transport. The shape of the water cone describes the pressure drop in the reservoir. The effect of a completely impermeable barrier on the cone shape is essentially the same as extending the wellbore out to the barrier radius. The placement of a horizontal barrier requires that a horizontal fracture be created using a single-point entry technique. Two methods may be employed for filling this fracture. One method consists of propping the fracture before filling with cement. The other consists of fracturing with fracture fluid, then immediately following with cement

By L. D. Mullins, W. R. Bartlett, R. H. Barham, W. L. Wahl

This paper describes an electrical model and its applicution to the analysis of four reservoirs in Saudi Arabia. The model has 2,501 mesh points and represents 35,000 sq miles of the Arab-D member. Details of modeling such as mesh size, control problems and standards of performance in matching reservoir history are discussed. The particular performance match achieved for the Arab-D rnember is presented. Details such as permeability barriers, aquifer depletion and iwerference between oil fields are given. The performance match realized in the Abqaiq pool is presented in detail. INTRODUCTION The resistor-capacitor network and associated control equipment described in this paper comprise an electrical analog of a reservoir system. Similar equipment has been used to study the transient response of reservoirs for many years. The unique feature of the model and application to be described is the extremely large size of the model and reservoir system, and the detail observed in simulating the reservoir with the model. The Arabian American Oil Co. first became interested in analog computers for simulation of oil reservoirs in 1949. Since that time, several models have been developed, each more elaborate and refined so that the reservoir system might be more closely simulated. The current model is the latest in a series designed, built and operated by the Field Research Laboratory of Socony Mobil Oil Co. in collaboration with Aramco. It has been and continues to be used to study the regional performance of the Arab-D member limestone reservoir. The Arab-D member is one of the Middle East's most prolific producing horizons.. THE MODEL The theory of simulating a reservoir system with an electrical system has been presented in the literature.1-5 Therefore, this paper will not discuss the theoretical aspect of the problem except to point out the correspondence between the fluid system and electrical system, as shown in Table 1. In general, the complete model is made up of input devices, output devices, central control and a resistance-capacitance (RC) network. At times, the RC network alone is referred to as the "model". However, it should be evident from the text which meaning is attached to the word "model". A discussion of the equipment follows. THE RESISTANCE-CAPACITANCE NETWORK The RC network consists of 2,501 capacitance decades interconnected through 4,900 resistance decades. The components are arranged to form a rectangular network of 2,501 mesh points in a 41- X 61-mesh array. Imposing the mesh grid system on the continuous reservoir system divides the reservoir into discrete areal segments. These discrete segments may be of various sizes. More precisely, the mesh size need not be uniform throughout the model. The RC network is fabricated in two sections which are connected at the top. An inside view of the "tunnel" formed by the two sections is shown in Fig. 1. The height and width of the tunnel are shown in the figure. Numerals appear along the bottom and along the back opening of the tunnel. These numbers denote the x and y co-ordinate positions of mesh points. Fig. 2 presents a rear view of one-half the model. The length dimensions of the model, as well as a rear view of the capacitor decades, are shown in this figure. The control dials used in adjusting the resistance and capacitance values on the model can be seen in the enlarged portion of the model shown in Fig. 3. The electrical capacity at any mesh point can range from 0 to 1.0 microfarads set to the nearest tenth of a microfarad. The electric resistance connecting any two mesh points can range from 0 to 9,990,000 ohms set to the nearest 1,000 ohms. External capacitors may be added to any or all mesh points if the need arises. The values of electrical resistance and capacitance are adjusted manually by manipulating the two types of decade units. INPUT EQUIPMENT A considerable quantity of equipment is used to control the input to the RC network. These input devices are

By M. R. Tek

A new method for correlating the data on multiphase flow through vertical pipe is presented. The correlation is based on a "two-phase f factor" concept which was developed and successfully applied to horizontal multiphase flow by Bertuzzi, Tek and Poettmann. Field data previously published on several flowing and gas-lift wells have been used as the basis for the developed correlation. The application of the method to actual two-phase flow problems indicate that this method is capable of predicting the pressure distribution in vertical multiphase flow strings well within the accuracy range usually desired for common engineering and design calculations. A new working chart developed for calculation of two-phase pressure gradients and a graphical step-by-step procedure for the computation of pressure distribution are presented along with an example problem. INTRODUCTION Multiphase flow through vertical pipe is encountered in many engineering installations. In petroleum, chemical process, nuclear engineering and many other industries, problems associated with simultaneous flow of two or more phases through vertical pipe have been of interest for a long time. This interest has increased considerably during recent years due to applications to new processes in petroleum production and refining and to problems of steam generation and heat removal from nuclear reactors. One prominent example of vertical two-phase flow is provided by the gas-lift process where oil, water and gas flow simultaneously. If the pressure profile in a gas-lift well can be predicted within reasonable accuracy, it would be possible to get good estimates of the power required to lift the oil, the optimum depth, the pressure and the rate at which to inject gas. Furthermore, the effect of production rate and tubing sizes on these quantities can be evaluated before any design decision is made on the installation and operation of the flow string. The majority of experimental work available in the literature deals with two-component systems where individual phase flow rates in and out of the pipe remain constant. The general problem of prediction of pressure drops in multiphase flow systems is very complicated. The co-existence of numerous flow patterns of widely different geometry and mechanism, conditions of surface instability and the nature of force fields acting upon the system are among the major difficulties commonly encountered. The classical approach of fluid dynamics which would be based upon the formulation and solution of Navier-Stokes equations has been found by many investigators completely devoid of any hope of success-—not only because of inherent nonlinearities but also because of insurmountable analytical difficulties standirng in the way of setting up the boundary conditions. The presence and effect of interfacial forces on mu1tiphase flow systems further complicate the theoretical approach. For these reasons, many investigators choose to adopt semi-empirical if not purely empirical approaches in order to obtain solutions of engineering utility. A prominent practical solution has been given by Poettmann and Carpenter' in the form of an empirical correlation. In their paper, total flowing densities of fluids and solubility effects of gas in oil have been taken into account for the correlation of field data covering a wide range of operating conditions for oil wells. They treated the gas. oil and water as a single phase of combined propcrties and correlated the multiphase friction factor as a function of the product (pvd), i.e., density X velocity X diameter of the flow string. However, because the product (pvd) is dimensional while the frictior factor is dimensionless, the generality of their result's is somewhat restricted. It seems that the omission 01 viscosity effect may be one of the reasons for the scatter of data as shown in their correlation chart. The concepts of the "two-phase f factor" and the "two-phase Reynolds' number function" were recently developed and successfully applied to correlate horizontal multphase flow by Bertuzzi, Tek and Poettmann.' Recently, two new methods of correlation by Tek and Chan' have been presented on simultaneous flow of liquid and gas through vertical pipe. These two methods of correlation and the working equations in field units necessary for their application are included in this paper. The extension of the conccpt of two-phasc Reynolds' number function successfully developed for horizontal flow into vertical multiphase flow systems, development and evaluation of working charts permitting calculation of two-phase pressure gradients and

By H. D. Attra

The effects of phase behavior on reservoir economics is a function of both fluid composition and flow properties of the reservoir rock. In some operations such as gas injection into volatile crude-oil resertoirs, the stock-tank recovery derived from vaporization of the reservoir oil by the displacing gas can approach the recovery derived from the displacement mechanism. In other cases of gas-cycling reservoirs, the dry injected gas can effectively vaporize a retrograde liquid phase, allowing a substantial reduction in cycling pressure with no appreciable loss in recoverable liquids. Because of these effects, it is essential that phase behavior be considered in developing the economics of any secondary recovery program in which gas is the injected fluid. The inclusion of phase behavior in conventional calculation procedures has been retarded because of the complexity of the calculation and the large volume of laboratory data required. However, by applying high-speed computers and new laboratory techniques, the phase-behavior displacement calculation has been developed to economically include this basic fundamental of gas displacement. The calculation is performed by reducing the reservoir to a series of one-dimensional segments. Phase changes, saturation distributions and changing fluid properties are evaluated for each segment as a function of time. The relative volumes of equilibrium liquid and gas displaced from each segment are determined by a trial-and-error procedure to meet the requirements of phase equilibrium and two-phase flow. The produced fluids are flashed through conditions of surface separation approximating those in the field to obtain actual values for stock-tank liquid recovery and producing gas-oil ratios. Based on the results of successful field studies, the technique appears to have application to reservoirs producing under a nonmiscible gas-displacement mechanism and should substantially increase the reliability of this type of analysis when phase behavior is an important factor. Predicting reservoir performance under gas-injection operations involves three basic principles of reservoir engineering. First, a material balance exists between the displacing gas and the displaced reservoir fluid. Second, the oil and gas flow in relationship with the two-phase flow property of the reservoir rock; and third, phase equilibrium exists between the displacing gas and the native reservoir fluid. To illustrate these three principles, consider the one-dimensional reservoir shown in Fig. 1. The reservoir is divided into a series of segments representing a one-dimensional cross section between an injection area and a producing area. So represents the oil saturation, Cw represents the connate-water saturation and Sg represents the displacing gas saturation. The gas is moving into Segment 1 and the reservoir oil is being produced out of Segment N. Consider the case when gas has just moved into Segment 1 and the displacing gas-oil contact is at position a-a The volume of oil pushed out of Segment 1 is a function only of mobility and the two-phase flow properties of the reservoir rock, that is, Kg/Ko and µo/µg. Behind the front if the injected gas is foreign to the reservoir oil (a nonequilibrium gas), there will either be vaporization of the oil into the gas or condensation of the gas into the oil. In either case the net result is a mass transfer and a re-distribution of composition and liquid and gas saturation. If the displacement is at constant pressure and is continued, the gas moves into new segments, proceeding succesively into Segment 2, 3, etc. The total volume of oil produced from the reservoir

By L. W. Emery

Undergrorlnd combustion operations were initiated in a 60-acre Bartlesville sand "shoe-string" reservoir in Allen Connty, Kans., in 1956. Tests in separate patterns were conducted using various co~nbinations of air and recycle gas to propagate combustion fronts from the injection toward the producing wells. These patterns were made up of 6 injection and 20 prodrrcing wells Gas and liquid prorluctiorz from each pattern was measured on an individual-well basis, and comparisons were made between the three patterns to ascertain the relative effects of injected gas composition on production behavior. Breakthrough of the combustion front at a well was characterized by an increase in water production from the well followed by an increase in bottomhole temperatrrre to approximately 250" F. After burning fronts had broken through at five producing wells, operations were terminated in 1960. From the total project approximately 79,000 bbl of oil were produced during thermal operations at a cumulative produced GOR of 23 Mcf/bbl. No appreciable change in the character of the produced crude was observed. Combustion in the reservoir was maintained with injected gas compositions ranging fronz 6 per cent oxygen in recycle gas to 100 per cent air. lnjectiotz of large quantities of recycle gas resulted in higher producing GOR's from offset wells than were measured from a pattern into ~vhich straight air ~vas injected. The air required to move the combustion front through I acre-ft of reservoir was computed to be 20 MMscf. This valrre was found to be relatively independent of the quantities of recycle gas injected. The recovery efficiency from the swept area was esti~izated to be about 59 per cent. Areas swept were similar in shape to tlzose obtained with a laboratory potentiometric model. Samples of sund taken from behind the burning front by coring indicated almost total oil removal from the sand. Petrographic analysis of the core samples indicated that the sand had been heated to peuk temperature of rlbout 1,200" F. No rignificant difference in peak temperature was forrnd in two areas where compositions of injected gas were quite different. Compression costs for thermal recovery were estimated to be $1,20/bhl of produced oil. INTRODUCTION The use of the "forward combustion" process as an oil recovery method has received a great deal of attention. This method involves ignition of the formation in an injection well, followed by propagation of a combustion front through the reservoir. Combustion is maintained by the injection of an oxygen-containing gas to react with reservoir hydrocarbons. As the flame front progresses through the reservoir, oil and formation water are vaporized, driven forward in the gaseous phase and recondensed in the cooler part of the formation. In turn, the condensed fluids push oil into the producing wellbores. Completed field tests of the process were first reported by Kuhn and Koch,' and by Grant and Szasz.' Results from other tests have since been reported by Walter,3 by Moss, White and McNeil,' and by Gates and Ramey." ach of these tests essentially utilized a single injection well surrounded by four or more producing wells. Sinclair Research, Inc., elected to do field experimental work using a number of test patterns in a single field in order that comparisons between various operating schemes could be made. The site selected and purchased in 1955 for this experimental work was a 60-acre Bartlesville sand reservoir located in Allen County, Kans. Combustion operations were initiated in mid-1956. Between that time and termination of the project in mid-1960, combustion fronts were propagated from injection wells to producers in three separate well patterns, using different mixtures of air and recycle gas. The test was terminated before sweep of the three patterns was complete so that information about the effect of combustion on the swept areas could be obtained by coring. Results from the test in the form of injection and producing well performance have been carefully recorded, and these form the general basis for this paper. DESCRIPTION OF RESERVOIR The reservoir in which the combustion tests were conducted is a Bartlesville sand "shoe-string", typical of a number of small reservoirs in Southeastern Kansas. Average reservoir characteristics are shown in Table 1. Fig. 1 is an isopachous map of the producing sand showing the reservoir to be approximately 400-ft wide and 2,500-ft long. Maximum net productive sand thickness is 21 ft. Fig. 2 shows a typical core analysis obtained by coring with water-base mud. The reservoir has no appreciable dip and is closed on the sides by degradation of sand into shale. The main body of sand is heavily laminated with shale stringers, which are not continuous between wells. The main reservoir is overlain by 30 to 40 ft of laminated low-permeability sand and shale streaks. No information is available on the original properties of

By P. E. Baker

Approximate solutions are presented for the heat-flow equations in a loss-free linear system with a moving source and with heat transfer by convection and conduction, representing in situ combustion in an oil reservoir. The rate of heat generation and the velocity vf of the source (flame front) are assumed constant. A cold front or convection front caused by the cooling effect of in-jected air develops behind the flame front and moves with velocity vc, which is determined by the thermal properties of the air and rock. Because of thermal conduction the fronts are not sharp; i.e., they do not represent temperature discontinuities. The following solutions apply after a sufficiently long time. where Q is heat generated per unit volume of formation, k is thermal conductivity, CF is volumetric heat capacity, a is thermal diffusivity, and subscripts 1 and 2 refer to values behind and ahead of the flame front, respectively. The parameter CF2 is contained in v . Te is initial reservoir temper- . C2 e ature and air-injection temperature. Other mathematically possible solutions are presented. The different solutions result from the use of different values for various physical parameters. The solution just presented is thought to be most commonly applicable. INTRODUCTION In this paper, approximate analytical expressions are derived for the temperature profile developed in the underground combustion method of oil recovery. A qualitative description of the process may be found in a paper by Tadema.1 The flame front is regarded as a moving source with a constant rate of heat generation. Heat transfer is by conduction and convection. Previous publications, of which a few are cited here, 2-5 have also attacked this problem and have presented solutions pertaining to different methods of approximating or idealizing the process studied. The present paper also treats an idealized case: linear flow is assumed for heat and fluids, and vaporization and condensation effects are neglected. The results obtained are similar to those derived by Bailey and Larkin,3 but are obtained by an entirely different procedure. The method used here should give considerably greater insight into the physics of the in situ combustion process. Four types of cases are considered, distinguished by thermal and fluid-flow properties of the system. A constant rate of heat generation and constant rate of flame-front advance are assumed. Flame-front temperature is determined for each case using a heat balance which involves the entire temperature profile. The applicability of the different solutions to real systems is discussed after they are all presented. Again, this is done for the purpose of understanding better the physics of the process. Consideration of the nonapplicable solutions and of the reasons for their not applying is significant. In a paper by Cooperman,4 the only case solved is one considered here as not applicable. Fluid flow, as such (pressures, saturations, oil recovery, etc.), is not treated. This subject has been covered in a paper by Wilson, et al.6 Also, we are not concerned with requirements for continuance of the process; we only attempt to determine the temperature history for cases where the process has continued for some time. PHYSICAL SYSTEM AND APPROACH TO PROBLEM A linear flow (heat and fluid) system is assumed and the x-axis is taken parallel to the flow direction. The combustion zone (flame front) is vertical,

By O. G. Kiel, G. W. Swift

The concept of "a continuous succession of steady-states", which has been applied successfully by Aromfsky and Jenkins to obtain a solution for the nonlinear partial differential equation describing the transient Darcy flow of gas through porous media, is demonstrated to be equally valid for transient non-Darcy flow. A mathematical model, which numerically solves the partial differential equation, is used to check the validity of the succession of steady-states solution. Comparison of sand-face pressure histories compelted by the two methods shows excellent agreement. The utility of the succession of steady-states solution in predicting performance of gas wells rests in the fact that no special computation equipment is required. The development of the succession of steady-states solution leads also to a practical method for determining and analyzing field test data. A method for taking gas-well test data under constant-rate conditions is presented. Experimental data obtained in the field by employing the constant-rate method are presented and analyzed in accordance with the succession of steady-states solution. Analysis of data in this fashion is demonstrated to give direct "in situ" information for reservoir permeability, porosity and turbulence or non-Darcy coefficient. INTRODUCTION The economics of gas production are dependent upon the transient behavior of flow within the reservoir. For production from a finite reservoir, the transient flow behavior can be subdivided into two parts. At first, the transient caused by the movement of the pressure "wave" into the reservoir is of importance. Later in the production history, the pressure-wave movement ceases and the second transient stage of material depletion becomes controlling. For reservoirs of relatively high permeability, it can be shown that the pressure wave moves into the reservoir and stabilizes quite rapidly. In the case of relatively impermeable reservoirs, quite the opposite is true. Although it is theoretically possible to compute the production capability of a well from the properties of the reservoir as determined by static tests and core analyses, more reliable information is obtained by conducting flow tests on the well and thereby obtaining some measure of "in situ" formation properties. For gas wells, there are two basic types of tests in existence: the flow-after-flow method,' and the isochronal method.' Both of these techniques are tailored to obtain data that can be analyzed in accordance with the empirical performance equation: In addition, the isochronal method makes provision for the sluggish nature of pressure-wave movement in "tight" formations by requiring pressure build-up between flows and by stipulating that data obtained on successive flows be analyzed at equal elapsed flow times. It can be demonstrated that either test is valid for reservoirs of high permeability. Further, since it has been pointed out that the pressure wave stabilizes rapidly for reservoirs of this type, tests of relatively short duration will give stabilized information on the performance of a well. Further decline of sand-face pressure and/or production rate may be determined by employing material-balance techniques. Cullender' points out that for relatively impermeable reservoirs the flow-after-flow method gives invalid results. (See Appendix C.) If the isochronal method of testing is used, there are two alternatives: (1) the tests must be conducted for a sufficient length of time to obtain stabilized information (which may require months to accomplish); or (2) some method for extrapolating the results of short-term isochronal tests must be employed. The first alternative is impracticable because of manpower, conservation and economic considerations. Recourse to the second alternative requires some assurance regarding the reliability of the extrapolation technique. Poettmann and SchilsonJ present an empirical method for predicting stabilized performance. The present investigation was originally initiated to determine the reliability of this technique. To do this, a mathematical model was developed to simulate the Darcy and non-Darcy flow of gas through porous media. The model consisted of a finite-difference approximation of the nonlinear partial differential equation which was solved on an IBM 7090 computer. Long-term production histories were simulated by the model and compared against predictions obtained from the Poettmann-Schilson method. As the work progressed, it became apparent that a straightforward predictive equation could be developed by utilizing the concept of a succession of steady-states. As a result, the emphasis of the work was redirected to exploit the advantages of the new method.

By G. Thodos, E. M. Phillips

A reduced density correlation for the gaseous and liquid regions of ethane has been developed from all the experimental PVT data available in the literature for this substance. Saturated vapor and liquid densities reported for methane, propane and n-butane were utilized to establish the rectilinear-diameter line for ethane. This line was then used with the saturated-vapor data for ethane to produce the saturated envelope for this substance. Most of the saturated liquid data reported for ethane were found to be inconsistent with the saturated-liquid curve resulting from this study. The correlation covers the regions 0.65 < TR <6.0 and 0 < PR < 15. For regions for which no experimental data have been reported, the isobars of the correlation were extended by the use of the theorem of corresponding states. This correlation is presented both on rectilinear and log-log coordinates. INTRODUCTION At the present time, there is considerable interest in the exact PVT behavior of pure hydrocarbons. Matschke and Thodos12 recently have developed a comprehensive compilation of experimental vapor and liquid density data for methane. Similar studies have been conducted for other pure substances.4, 7 In the present study the density data available in the literature for gaseous and liquid ethane have been critically reviewed and utilized to develop a reduced density correlation for this substance. SATURATED LIQUID AND VAPOR DENSITIES The rectilinear-diameter rule proposed by Cailletet and Mathias5 states that the sum of the saturated liquid and vapor densities of a pure substance is a linear function of temperature which passes through the critical point. This law, although empirical, has been well established by continued verification since its conception in 1886. This rule has been applied to the saturated density data of methane,6 propane 17 and n-butane8 to produce the straight lines presented in Fig. 1 which pass through the critical point TR = 1.00, Pr = 2.00. The rectilinear-diameter relationship for ethane was established from the lines for the other substances in Fig. 1 by plotting the quantity p~, + PR against zc for different reduced temperatures. A critical analysis of the saturated vapor and liquid densities reported for ethane by Kuenen,10 Maass and McIntosh,11 and Sage, Webster and Lacey18 revealed that the values of PR1 + pr resulting from these studies were inconsistent with the rectilinear-diameter line obtained from the reliable density data of the other substances, as shown in Fig. 1. In the analysis of these data, reduced quantities were calculated using the critical con-

By H. J. Ramey

As fluids move through a wellbore, there is transfer of heat between fluids and the earth due to the diflerence between fluid and geothermal temperatures. This type of heat transmission is involved in drilling and in all producing operations. In certain cases, quantitative knowledge of wellbore heat transmission is very important. This paper presents an approximate solution to the wellbore heat-transmission problem involved in injection of hot or cold fluids. The solution permits estimation of the temperature of fluids, tubing and casing as a function of depth and time. The result is expressed in simple algebraic form suitable for slide-rule calculation. The solution assumes that heat transfer in the wellbore is steady-state, while heat transfer to the earth will be unsteady radial conduction. Allowance is made for heat resistances in the wellbore. The method used may be applied to derivation of other heat problems such as flow through multiple strings in a wellbore. Comparisons of computed and field results are presented to establish the usefulness of the solution. INTRODUCTION During the past few years, considerable interest has been generated in hot-fluid-injection oil-recovery methods. These methods depend upon application of heat to a reservoir by means of a heat-transfer medium heated at the surface. Clearly, heat losses between the surface and the injection interval could be extremely important to this process. Not quite so obvious is the fact that every injection and production operation is accompanied by transmission of heal between wellbore fluids and the earth. Previously, the interpretation of temperature logs&apos;,&apos; has been the main purpose of wellbore heat studies. The only papers dealing specifically with long-time injection operations are those of Moss and White3 and Lesem, et al.&apos; The purpose of the present study is to investigate wellbore heat transmission to provide engineering methods useful in both production and injection operations, and basic techniques useful in all wellbore heat-transmission problems. The approach is similar to that of Moss and White:&apos; DEVELOPMENT The transient heat-transmission problem under consideration is as follows. Let us consider the injection of a fluid down the tubing in a well which is cased to the top of the injection interval. Assuming fluid is injected at known rates and surface temperatures, determine the temperature of the injected fluid as a function of depth anti time. Consideration of the heat transferred from the injected fluid to the formation leads to the following equations. For liquid, Eqs. 1, 1A and 2 are developed in the Appendix. These equations were developed under the assumption that physical and thermal properties of the earth and wellbore fluids do not vary with temperature, that heat will transfer radially in the earth and that heat transmission in the wellbore is rapid compared to heat flow in the formation and. thus, can be represented by steady-state solutions. Special cases of this development have been presented by Nowakl and Moss and White.3 Both references are recommended for excellent background material. Nowak&apos; presents very useful information concerning the effect of a shut-in period on subsequent temperatures. Since one purpose of this paper is to present methods which may be used to derive approximate solutions for heat-transmission problems associated to those specifically considered here, a brief discussion of associated heat problems is also presented in the Appendix. Analysis of the derivation presented in the Appendix will indicate that many terms can be re-defined to modify the solution for application to other problems. Before Eqs. 1, 1A and 2 can be used, it is necessary to consider the significance of the over-all heat-transfer coefficient U and the time function f(t). Thorough discussions of the concept of the over-all heat-transfer coefficient may be found in many references on heat transmission. See McAdams5 or Jakob," for example. Briefly, the over-all coefficient U considers the net- resistance to heat flow offered by fluid inside the tubing, the tubing wall, fluids or solids in the annulus, and the casing wall. The effect of radiant heat transfer from the tubing to the casing and resistance to heat flow caused by scale or wax on the tubing or casing may also be included in the over-all coefficient. According to McAdams, on page 136 of Ref. 5>

By W. E. Brigham, P. W. Reed, J. N. Dew